跟着Nature学绘图——分组直方图！

 

点击上方

“科研后花园”

关注我们

复现图片如下：

主要复现的图片是上图中的图b-i类似的样式，具体代码如下：

1、设置工作环境并加载所需R包：

#设置工作环境
rm(list=ls())
setwd("D:/分组直方图")

#加载包
library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics

2、加载绘图数据（根据原文图片随机编写数据，无实际意义）：

#加载数据
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)

PS: 以上内容是小编个人学习代码笔记分享，仅供参考学习，欢迎大家一起交流学习。

代码及源数据获取：https://github.com/wzsBio/code

推荐阅读

1

2

3

4

5

6

7

温馨提示

如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。