0

我想用 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()

我错过了什么?

4

1 回答 1

1

您正在对日志值进行测试:

t.test(log10(RC) ~ Drug, data = mydf, exact = FALSE)
# 0.3237
于 2018-02-22T19:18:57.513 回答