我想用 ggplot2 绘制数据并使用 ggpubr 添加 p 值。compare_means()
但是绘制的 p 值与我用or计算的 p 值不匹配t.test()
。
这是我的数据:
Group <- rep(c('A'), each=12)
Drug <- rep(c('x','y'), each= 6)
RC <- c(13076, 10814, 14297, 13119, 13616, 5, 19671, 18318, 12058, 17624, 9565, 13689)
mydf <- as.data.frame(cbind(Group,Drug,RC))
mydf$RC <- as.numeric(as.character(mydf$RC))
这是情节的代码:
library(ggplot2)
library(ggpubr)
ggplot(mydf, aes(x=factor(Group),y=log10(RC)))+
geom_boxplot(aes(color=Drug))+stat_compare_means(aes(group = Drug), method = "t.test", label = "p.format")
显示的 p.value 为 p=0.32。但是,当我使用其他方法计算 p.value 时,我得到了 ~0.149
compare_means(RC~Drug,data=mydf, method = "t.test")
t.test(RC~Drug,data=mydf, exact= FALSE)
这里已经解决了一个类似的问题,但在这种情况下,与和compare_means()
相比,给出了另一个结果。我确保拥有最新版本的 ggpubr (ggpubr_0.1.6.999)。stat_compare_means()
t.test()
我错过了什么?