r - R中的箱线图比较-如何在箱线图中以垂直线显示p值？

Question

我生成了一个箱线图并测试了向量 x1 和 x2 与wilcox.test. 如何在箱线图中实现测试结果？

> x1 <- rnorm(1000)
> x2 <- rnorm(1000) +10
> wilcox.test(x1,x2, paired=TRUE)

    Wilcoxon signed rank test with continuity correction

data:  x1 and x2
V = 0, p-value < 2.2e-16
alternative hypothesis: true location shift is not equal to 0

> boxplot(x1, x2)

感谢你们！

R版3.1.0（2014-04-10）——《春之舞》

平台：x86_64-w64-mingw32/x64（64位），Windows 7

Answer 1

尝试做：

text(x=1, y=10, labels="p-value<0.0001")

我是 R 的新手，但如果您输入x和y坐标并输入标签，您可以在图表上手动输入测试结果。

我相信您已经知道这一点，但\n将用于在您的标签中制止。

Answer 2

您可以按照此处所述执行此操作：

library(tidyverse)
library(ggpubr)

x1 <- rnorm(1000)
x2 <- rnorm(1000) +10
dt =  tibble(x1, x2) %>% pivot_longer(1:2)

my_comparisons <- list( c("x1", "x2") )
ggboxplot(dt, x = "name", y = "value",
          color = "name", palette = "jco")+ 
    stat_compare_means(comparisons = my_comparisons, method.args = list(paired=TRUE))