ggplot2添加p值和显著性—科研工具箱-叨客学习资料网

请关注公众号【叨客学习资料】在使用网站的过程中有疑问，请来公众号进行反馈哦

数据可视化分析中我们经常需要进行数据间的统计分析，并进行显著性标记，虽然ggpur包被大佬吐槽制造混乱，但在进行显著性标记标记方面也是有其可取之处。

下面主要展示ggpubr与ggsignif包对数据进行统计分析，喜欢的小伙伴可以关注我的公众号R语言数据分析指南将分享更多实用文档，先行拜谢了

加载R包

library(pacman)
pacman::p_load(tidyverse,ggpubr,rstatix,ggsci,ggsignif,reshape2)

wilcox(秩合检验）

ggplot(iris,aes(Species,Sepal.Length,fill=Species)) + 
geom_boxplot()+geom_jitter(shape=16,size=2,position=position_jitter(0.2))+
  geom_signif(comparisons = list(c(\"versicolor\", \"virginica\"),
                                 c(\"versicolor\",\"setosa\"),
                                 c(\"virginica\",\"setosa\")),
              map_signif_level=T,
              textsize=6,test=wilcox.test,step_increase=0.2)+
  guides(fill=FALSE)+xlab(NULL)+theme_classic()

t_test分析

ggplot(iris,aes(Species,Sepal.Length,fill=Species)) + 
  geom_boxplot()+
  scale_fill_jco()+
  geom_jitter(shape=16,size=2,position=position_jitter(0.2))+
  geom_signif(comparisons = list(c(\"versicolor\", \"virginica\"),
                                 c(\"versicolor\",\"setosa\"),
                                 c(\"virginica\",\"setosa\")),
              map_signif_level=T,
              textsize=6,test=t.test,step_increase=0.2)+
  guides(fill=F)+xlab(NULL)+theme_classic()

anova分析

ggplot(iris,aes(Species,Sepal.Length,fill=Species)) + 
  geom_boxplot()+
  scale_fill_jco()+
  geom_jitter(shape=16,size=2,position=position_jitter(0.2))+
  stat_compare_means(method = \"anova\")+theme_bw()

kruskal分析

ggplot(iris,aes(Species,Sepal.Length,fill=Species)) + 
  geom_boxplot()+
  scale_fill_jco()+
  geom_jitter(shape=16,size=2,position=position_jitter(0.2))+
  stat_compare_means()+theme_bw()

对多组数据进行分面

iris %>% melt() %>%
  ggplot(aes(variable,value,fill=Species)) + 
  geom_boxplot()+
  scale_fill_jco()+
  stat_compare_means(label = \"p.format\")+
  facet_grid(.~variable,scales = \"free\",space=\"free_x\")+
  theme_bw()+theme(axis.text.x = element_blank())

iris %>% melt() %>%
  ggplot(aes(variable,value,fill=Species)) + 
  geom_boxplot()+
  scale_fill_jco()+
  stat_compare_means(label =  \"p.signif\", label.x = 1.5)+
  facet_grid(.~variable,scales = \"free\",space=\"free_x\")+
  theme_bw()+theme(axis.text.x = element_blank())

版权声明 1、本网站名称：叨客学习资料网
2、本站永久网址：https://www.leobba.cn
3、本网站的文章部分内容可能来源于网络，仅供大家学习与参考，如有侵权，请联系站长进行删除处理。
4、本站一切资源不代表本站立场，并不代表本站赞同其观点和对其真实性负责。
5、本站一律禁止以任何方式发布或转载任何违法的相关信息，访客发现请向站长举报
6、本站资源大多存储在云盘，如发现链接失效，请联系我们我们会第一时间更新。

THE END