使用 stat_compare_means() 显示统计检验 + 个体数量 (n = )

Question

我目前使用

stat_compare_means()

在箱线图上显示与

ggplot2

的比较测试。但是，我希望该函数显示测试名称和给定 p 值之间样本的个体数 (n = )。有没有办法将该信息添加到

stat_compare_means()

中？我可以轻松地用

geom_text()

或

annotate()

显示它，但这不太方便，因为我希望它写在 p 值之前。

data(iris)
ggplot(iris, aes(x=Species, y=Sepal.Length)) +
  geom_boxplot() +               
  geom_jitter(width = 0.1, shape=21, colour="black", fill="grey95", stroke=0.5, size=2) +
  stat_compare_means(label.y.npc = 0.93, method = "kruskal.test", size = 5, family = "serif") +
  theme_classic()

Answer 1

我们可以通过调整

n

的

label

美感来获得总的

stat_compare_means

:

ggplot(iris, aes(x=Species, y=Sepal.Length)) +
  geom_boxplot() +               
  geom_jitter(width = 0.1, shape = 21, colour = "black", fill = "grey95", 
              stroke = 0.5, size = 2) +
  stat_compare_means(label.y.npc = 0.93, method = "kruskal.test", 
                     size = 5, family = "serif",
                     aes(label = paste0(after_stat(method), 
                                        ', n =', nrow(iris),
                                        ', p = ', after_stat(p.adj)))) +
  theme_classic()