在R中的多个boxplot中添加pvalue

Question

我正在尝试将p值添加到我的boxplot并显示在我的情节中。我在互联网上找到了一些东西：

到目前为止，我可以获得pvalue并将其设置在标题中，但它看起来非常糟糕。

是否可以在我的绘图的右上角创建一个框并显示每个数据的p值？

这是我正在做的一个例子：

par(mfrow=c(2,2))

boxplot(x = df, main = round(t.test(df)$p.value, digits=3))
boxplot(x = df1, main = round(t.test(df1)$p.value, digits=3))
boxplot(x = df2, main = round(t.test(df2)$p.value, digits=3))
boxplot(x = df3, main = round(t.test(df3)$p.value, digits=3))

Answer 1

这是一个简单的例子：

# Fake data
set.seed(448)
dat = data.frame(x = c(rnorm(20), rnorm(20,1,1)), y=rep(c("A","B"), each=20))

boxplot(dat$x ~ dat$y, xlim=c(0.5,3))
rect(2.85 - 0.2, max(dat$x) - 0.12, 2.85 + 0.2, max(dat$x) + 0.12)
text(2.85, max(dat$x), paste0("p < ", round(t.test(dat$x ~ dat$y)$p.value, 3)))