在 R 中的多个箱线图中添加 pvalue

Question

我正在尝试将 p 值添加到我的箱线图中并显示在我的图中。我在网上找到了一些接近的东西：

到目前为止，我可以获取 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)))

在 R 中的多个箱线图中添加 pvalue

Adding pvalue inside multiple boxplot in R

r

boxplot