我被困在如何绘制分组箱形图的平均p值上。
这是我数据的一部分:
ID<-c('E5b','R6',"S22","E5b","R6","S22","E5b","R6","S22","E5b","R6","S22","E5b","R6","S22","E5b","R6",
"S22","E5b","R6","S22","E5b","R6","S22","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, 0, 1.90550741185218e-06, 1.52440592948174e-06, 0, 1.23540828390671e-06, 4.11802761302236e-07,
0, 5.22781921260155e-06, 1.04556384252031e-06, 0, 1.71521997010029e-06, 0, 2.05826396412034e-06, 4.18012063828162e-06,
0, 7.60021934233022e-07, 2.93951950197596e-05, 0, 2.31458228502044e-07)
condition<-c("E","E","E","E","E","E","E","E","E","E","R","R","R","R","R","R","R","R","R",
"R","R","S","S","S","S","S","S","S","S","S")
family<-c("Unassigned","Unassigned","Siphoviridae","Unassigned","Unassigned","Siphoviridae","Unassigned","Unassigned",
"Siphoviridae","Unassigned","Unassigned","Siphoviridae","Unassigned","Unassigned","Siphoviridae","Unassigned",
"Unassigned","Siphoviridae","Unassigned","Unassigned","Siphoviridae","Unassigned","Unassigned","Siphoviridae",
"Unassigned","Unassigned","Siphoviridae","Unassigned","Unassigned","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外,我想创建如下内容。非常感谢您的帮助!
发布于 2021-02-15 18:55:17
需要在stat_compare_means中添加(aes(组=条件)。
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)
https://stackoverflow.com/questions/66196884
复制相似问题