跟着Nature Communications学作图:R语言ggplot2绘制带有条纹的分组柱形图
语言 分组 绘制 Nature ggplot2 跟着 作图 带有
2023-06-13 09:16:32 时间
论文
Pan-African genome demonstrates how population-specific genome graphs improve high-throughput sequencing data analysis
https://www.nature.com/articles/s41467-022-31724-3
本地pdf s41467-022-31724-3.pdf
论文中公布了大部分图的数据,但是没有公布对应的作图代码,没有关系,我们可以自己写代码试着模仿,今天的推文重复一下论文中的Figure 2A 带有条纹的分组柱形图
示例数据截图
image.png
这里实现条纹柱形图用到的是 ggpattern
这个R包
参考链接
https://coolbutuseless.github.io/package/ggpattern/index.html
https://github.com/coolbutuseless/ggpattern
安装
remotes::install_github("coolbutuseless/ggpattern")
因为是ggplot2的扩展包,除了把作图函数替换,其余的细节都可以用ggplot2的语法来调节
读取数据
library(readxl)
dffig2a<-read_excel("data/20220806/41467_2022_31724_MOESM4_ESM.xlsx",
sheet = "figure 2a")
dffig2a
library(tidyverse)
dffig2a %>%
pivot_longer(-'Super-population') -> new.dffig2a
作图代码
library(ggplot2)
cols<-c("#ffa657","#fd8011","#6cbe6c","#349734",
"#eba0d5","#da7dbd","#63a0cb","#1f7ab4",
"#d0d166","#bbbe21")
ggplot(data = new.dffig2a,aes(x=`Super-population`,y=value))+
geom_bar_pattern(stat="identity",
position = "dodge",
aes(pattern=name,
fill=name),
pattern_density=0.01,
fill=cols,
color="black",
show.legend = FALSE)+
scale_pattern_manual(values = c('Divergence'='stripe',
'Diversity'="none"))+
scale_y_continuous(expand = expansion(mult = c(0,0.1)),
labels = scales::percent,
limits = c(0,0.25/100),
breaks = seq(0,0.25/100,by=0.05/100))+
labs(x=NULL,y=NULL)+
theme_classic()+
theme(axis.line.y = element_blank(),
axis.ticks.y = element_blank(),
panel.grid = element_line(linetype = "dashed"),
panel.grid.major = element_line(),
panel.grid.minor = element_blank()) -> p1
p1
ggplot()+
geom_rect_pattern(data=data.frame(x=1,xend=2,y=1,yend=2),
aes(xmin=x,ymin=y,xmax=xend,ymax=yend),
pattern_density=1,
fill="white",
color="black")+
geom_rect_pattern(data=data.frame(x=1,xend=2,y=2.5,yend=3.5),
aes(xmin=x,ymin=y,xmax=xend,ymax=yend),
pattern="none",
pattern_density=1,
fill="grey",
color="black")+
theme_void()+
geom_text(data=data.frame(x=2,y=1.5),
aes(x=x,y=y),label="Divergence",
hjust=-0.1)+
geom_text(data=data.frame(x=2,y=3),
aes(x=x,y=y),label="Diversity",
hjust=-0.1)+
xlim(1,4) -> p2
p1+
annotation_custom(grob = ggplotGrob(p2),
xmin = 4,xmax = Inf,
ymin = 0.2/100,ymax=0.25/100) -> p3
p3
library(patchwork)
p3+p3
image.png
示例数据可以在论文中去下载,代码直接在推文中复制
相关文章
- 【R语言】常规除法,取整,取余,向上取整,向下取整,四舍五入
- 【R语言】dplyr对数据分组取各组前几行
- 【R语言】因子在临床分组中的应用
- 【说站】c语言中如何计算数组长度
- 2022-10-31:以下go语言代码输出什么?A:map[];B:nil;C:Panic;D:编译错误。 package main import “fmt“
- R语言简单随机分组/区组随机/分层随机
- R语言Apriori关联规则、kmeans聚类、决策树挖掘研究京东商城网络购物用户行为数据可视化|附代码数据
- Go语言中常见100问题-#1 Unintended variable shadowing
- 跟着Nature Genetics学作图:R语言ggplot2普通箱线图/分组箱线图/分面箱线图
- 跟着NaturePlants学作图:R语言ggplot2分组折线图完整示例
- 跟着NatureCommunications学作图:R语言ggtree根据分组给进化树上色
- R语言ggplot2做柱形图如何让起始位置不是0?
- 【C 语言】文件操作 ( 读取文件中的结构体数组 | feof 函数使用注意事项 )
- Linux和C语言:联合开发精彩软件(linux和c语言)
- 在 Fedora 上开启 Go 语言之旅
- 为什么使用 Go 语言?
- 用c语言实现MySQL中自增的操作(c mysql 自增)
- JavaScript语言对Unicode字符集的支持详解
- go语言工程结构