ggplot2构建误差线组图并添加组间显著性标记

R语言数据分析指南

发布于 2023-08-18 13:57:34

2520

发布于 2023-08-18 13:57:34

文章被收录于专栏：R语言数据分析指南R语言数据分析指南

欢迎关注R语言数据分析指南

❝本节来复现nature上的一张小图介绍如何使用「ggplot2构建误差线组图并添加组间显著性标记」，下面小编就通过一个案例来进行结果展示，添加了详细的注释希望各位观众老爷能够喜欢 ❞

论文

❝Lactate limits CNS autoimmunity by stabilizing HIF-1α in dendritic cells ❞加载R包

library(tidyverse)
library(ggsci)
library(rstatix)
library(ggpubr)
library(cowplot)

数据清洗

df <- read_tsv("data.txt") %>%
  pivot_longer(-Id) %>% 
  mutate(group = sub("\\.\\.\\..*", "", name),
         Id=as.factor(Id)) %>% select(-name)

构建绘图数据

plot <- df %>% group_by(Id,group) %>% 
  summarise(value_mean=mean(value),sd=sd(value),se=sd(value)/sqrt(n()))

绘制主图

p <- plot %>% ggplot(aes(Id,value_mean,fill=group,group=group,
                  ymin=value_mean-se,ymax=value_mean+se))+
  geom_errorbar(width=0.2)+
  geom_line(color="black")+
  geom_point(aes(pch=group,color=group,fill=group),size=3)+
  scale_shape_manual(values=c(25,21,22,24))+
  scale_fill_npg()+
  scale_color_npg()+
  labs(x=NULL,y=NULL)+
  theme_classic()+
  theme(plot.margin=unit(c(0.1,1,0.1,0.1),units=,"cm"),
        axis.line = element_line(color = "black",size = 0.4),
        panel.grid.minor =element_blank(),
        panel.grid.major =element_blank(),
        axis.text.y = element_text(color="black",size=9),
        axis.text.x = element_text(margin = margin(t=3),color="black",size=9,angle=0,vjust = 0.5),
        legend.key=element_blank(),
        legend.title = element_blank(),
        legend.text = element_text(color="black",size=8),
        legend.spacing.x=unit(0.1,'cm'),
        legend.spacing.y=unit(1,'cm'),
        legend.key.width=unit(0.7,'cm'),
        legend.key.height=unit(0.4,'cm'),
        legend.background=element_blank(),
        legend.position = c(0.2,0.95),
        legend.box.margin = margin(0,0,0,0))

统计分析

df %>%
  t_test(value ~ group) %>%
  adjust_pvalue() %>% add_significance("p.adj") %>% 
  select(2,3,p.adj.signif)

添加注释

ggdraw()+
  draw_plot(p)+
  draw_line(x=c(0.95,0.95),y=c(0.35,0.63))+
  draw_line(x=c(0.95,0.95),y=c(0.65,0.8))+
  draw_text(x=0.975,y=0.5,text="****",angle=90,size=20,color="red")+
  draw_text(x=0.975,y=0.72,text="****",angle=90,size=20,color="red")