如何解决stat_compare_means 比较不在 x 轴上的组
我一直在研究如何绘制分组箱线图的平均 p 值。
这是我的一部分数据:
ID<-c('E5b','R6',"S22","E5b","R6","S22")
value<-c(1.02048033657553e-05,7.03779542465882e-07,3.51889771232941e-07,5.69459095210849e-06,5.42341995438904e-07,1.08468399087781e-06,1.15124329576991e-05,2.34947611381614e-07,4.69895222763228e-07,1.02807349661977e-05,2.12704861369607e-06,1.90550741185218e-06,1.52440592948174e-06,1.23540828390671e-06,4.11802761302236e-07,5.22781921260155e-06,1.04556384252031e-06,1.71521997010029e-06,2.05826396412034e-06,4.18012063828162e-06,7.60021934233022e-07,2.93951950197596e-05,2.31458228502044e-07)
condition<-c("E","E","R","S","S")
family<-c("Unassigned","Unassigned","Siphoviridae","Siphoviridae")
df <- data.frame(ID,value,condition,family)
my_comparisons <- list(c("E","R"),c("E","S"),c("R","S"))
p <- ggboxplot(df,x = "ID",y = "value",color = "condition",palette = "jco",add = "jitter")+
facet_wrap(~family,scales='free_x')
p + stat_compare_means(comparisons = my_comparisons,bracket.size = .4,size = 8)+
stat_compare_means(method = 't.test',size = 3)
这仅显示一组的 p 值?
Warning message:
"Computation failed in `stat_signif()`:
missing value where TRUE/FALSE needed"
Apparently stat_compare_means 只能比较 x 轴上的组。有没有办法在所有条形之间绘制条形及其 p 值?
理想情况下,除了 x 轴上的 ID 外,我想创建如下内容。非常感谢任何帮助!
解决方法
需要在 stat_compare_means 中添加 (aes(group = condition)。
p <- ggboxplot(df,x = "ID",y = "value",color = "condition",palette = "jco",add = "jitter")+
facet_wrap(~family,scales='free_x')+
scale_y_continuous(labels = comma)+
stat_compare_means(aes(group = condition),size=3)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。