其他
跟着Nature学绘图——分组直方图！
点击上方
“科研后花园”
关注我们
复现图片如下：
主要复现的图片是上图中的图b-i类似的样式，具体代码如下：
1、设置工作环境并加载所需R包：
#设置工作环境rm(list=ls())setwd("D:/分组直方图")
#加载包library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics2、加载绘图数据（根据原文图片随机编写数据，无实际意义）：
#加载数据df <- read.table("data.txt",header = T,sep='\t')df$group <- factor(df$group,levels = c("Sal","Coca","Ket","LSD","MDMA"))df$sample <- factor(df$sample,levels = c("48 h","2 wk"))head(df)3、颜色：
col <- c("#000000","#575757","#c53a8e","#e79600","#a42422")col2 <- c("#767475","#000000","#872860","#f0b75b","#c97c7b")4、绘图：
ggplot(df,aes(sample,value))+ geom_bar(aes(fill=group),color="black",stat="summary",fun=mean,position="dodge",size=1)+ stat_summary(fun.data = 'mean_sd', geom = "errorbar", width = 0,linewidth=1)+ geom_point(aes(color=group),shape=21,size=4,stroke=1.5)+ geom_hline(yintercept = 0, linetype = 1, color = "black", size = 1)+ facet_grid(~group,scales = 'free_x',space = "free")+ scale_color_manual(values = col2)+ scale_fill_manual(values = col)+ theme_classic()+ theme(axis.line = element_line(size = 1), axis.text.x = element_text(color = "black", angle = 90,vjust = 0.5,hjust = 1,size = 15), axis.text.y = element_text(color = "black",size = 15), axis.ticks = element_line(color = "black",size = 1), legend.position = "none", strip.background = element_blank(), strip.text = element_text(color = "black",size = 18))+ labs(x=NULL,y=NULL)推荐阅读
1
2
3
4
扩增子测序数据分析还不会？小编整理的全套R语言代码助您轻松解决问题！
5
6
7
温馨提示
如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。