与 geom_signif 函数的多重比较，R

Question

该包ggsignif对于快速、轻松地指示ggplot图表中的显着比较非常有用。但是，该comparisons调用需要手动键入要比较的每对值。

例如。

library(ggplot2)
library(ggsignif)

data(iris)

ggplot(iris, aes(x=Species, y=Sepal.Length)) + 
  geom_boxplot() +
  geom_signif(comparisons = list(c("versicolor", "virginica"),c('versicolor','setosa')), 
              map_signif_level=TRUE)

我想知道如何通过立即引用所有可能的组合来避免这种情况？例如，expand.grid(x = levels(iris$Species), y = levels(iris$Species))，给出所有组合

           x          y
1     setosa     setosa
2 versicolor     setosa
3  virginica     setosa
4     setosa versicolor
5 versicolor versicolor
6  virginica versicolor
7     setosa  virginica
8 versicolor  virginica
9  virginica  virginica

但如何让这个被接受geom_signif(comparisons=...呢？

软件包信息可在此处获取https://cran.r-project.org/web/packages/ggsignif/index.html

Answer 1

根据 Adam Quek 的评论，您只需转置创建的矩阵并将每一行转换为列表：

split(t(combn(levels(iris$Species), 2)), seq(nrow(t(combn(levels(iris$Species), 2)))))

$`1`
[1] "setosa"     "versicolor"

$`2`
[1] "setosa"    "virginica"

$`3`
[1] "versicolor" "virginica" 

ggplot(iris, aes(x = Species, y = Sepal.Length)) + 
  geom_boxplot() +
  geom_signif(comparisons = split(t(combn(levels(iris$Species), 2)), seq(nrow(t(combn(levels(iris$Species), 2))))), 
              map_signif_level = TRUE)