如何向多面箱线图添加统计信息？

Question

我正在尝试创建一个分面箱线图，它由集群分面。 x 轴是时间，y 轴是计数。 我正在尝试生成将第 0 天与第 1 天、第 0 天与第 2 天以及第 1 天与第 2 天进行比较的统计数据。

该数据是整个数据框：

structure(list(time = c("Day 0", "Day 0", "Day 0", "Day 0", "Day 0", 
"Day 0", "Day 0", "Day 0", "Day 0", "Day 0", "Day 1", "Day 1", 
"Day 1", "Day 1", "Day 1", "Day 1", "Day 1", "Day 1", "Day 1", 
"Day 1", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2"), count = c(1L, 4L, 
1L, 1L, 1L, 1L, 2L, 3L, 1L, 1L, 59L, 1L, 1L, 1L, 3L, 1L, 2L, 
5L, 3L, 1L, 1L, 2L, 1L, 1L, 1L, 54L, 3L, 6L, 1L, 1L, 1L, 2L, 
1L, 1L, 7L, 1L, 3L, 1L, 207L, 5L, 3L, 3L, 11L, 2L, 1L), cluster = c("C", 
"C", "C", "C", "D", "D", "D", "D", "D", "D", "A", "A", "D", "D", 
"D", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "B", "C", 
"B", "B", "C", "C", "C", "C", "C", "B", "B", "B", "A", "A", "A", 
"A", "A", "A", "A", "A")), row.names = c(NA, -45L), class = "data.frame")

对于统计数据，我使用的是 wilcox.test。 这是我用于测试的代码：

wilcox.test(CDR3_D0$count, CDR3_D1$count, alternative = "two.sided", exact = FALSE)
wilcox.test(CDR3_D0$count, CDR3_D2$count, alternative = "two.sided", exact = FALSE)
wilcox.test(CDR3_D1$count, CDR3_D2$count, alternative = "two.sided", exact = FALSE)

为了统计，我把完整的数据帧按照时间分成了3个数据帧。 以下是这些数据帧的代码：

Day 0:
structure(list(time = c("Day 0", "Day 0", "Day 0", "Day 0", "Day 0", 
"Day 0", "Day 0", "Day 0", "Day 0", "Day 0"), count = c(1L, 4L, 
1L, 1L, 1L, 1L, 2L, 3L, 1L, 1L), cluster = c("C", "C", "C", "C", 
"D", "D", "D", "D", "D", "D")), row.names = c(NA, -10L), class = "data.frame")

Day 1:
structure(list(time = c("Day 1", "Day 1", "Day 1", "Day 1", "Day 1", 
"Day 1", "Day 1", "Day 1", "Day 1", "Day 1"), count = c(59L, 
1L, 1L, 1L, 3L, 1L, 2L, 5L, 3L, 1L), cluster = c("A", "A", "D", 
"D", "D", "C", "C", "C", "C", "C")), row.names = c(NA, -10L), class = "data.frame")

Day 2:
structure(list(time = c("Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2", 
"Day 2", "Day 2", "Day 2", "Day 2", "Day 2", "Day 2"), count = c(1L, 
2L, 1L, 1L, 1L, 54L, 3L, 6L, 1L, 1L, 1L, 2L, 1L, 1L, 7L, 1L, 
3L, 1L, 207L, 5L, 3L, 3L, 11L, 2L, 1L), cluster = c("D", "D", 
"D", "D", "D", "B", "C", "B", "B", "C", "C", "C", "C", "C", "B", 
"B", "B", "A", "A", "A", "A", "A", "A", "A", "A")), row.names = c(NA, 
-25L), class = "data.frame")

我用来生成带有统计信息的分面箱线图的脚本是：

CDR3_clust_plot <- ggplot(CDR3_clust_plot_df, aes(y = count, x = Time, group = Time)) + 
  theme_bw() +
  geom_boxplot()+
  theme(legend.position = "none")+
  scale_y_log10(limits = c(1, 250)) +
  facet_wrap(vars(cluster), scales = "fixed")+
  stat_compare_means(method = "wilcox.test", aes(label = paste0(..p.format..)), label.y = 900, label.x = 1.45)+
  geom_signif(y_position = 875, xmin = 1, xmax = 2, annotations = "")+
  labs(y = "Count", x = "Time")

我最终得到了多面但没有统计信息的图，并收到以下警告：

警告信息：

1： stat_compare_means()中的计算失败：缺少参数“x”，没有默认值

2：删除了包含缺失值的 7 行 (geom_text)。

3：删除了包含缺失值的 6 行 (geom_signif)。