r - 如何使用 tidyverse 和 rstatix 显示 Tukey 组？

Question

在使用 Tukey 测试后rstatix，我想使用 Tukey 组向数据框中添加另一列。

library(tidyverse)
library(rstatix)

ggplot(data = iris,
       aes(x = Species,
           y = Petal.Length)) +
  geom_boxplot()

Tukey <- iris %>%
  tukey_hsd(Petal.Length ~ Species) %>%
  add_significance() %>%
  # ????? <------- Help here

我期待这样的事情：

> Tukey
# A tibble: 3 x 9
term    group1     group2     null.value estimate conf.low conf.high p.adj p.adj.signif group
* <chr>   <chr>      <chr>           <dbl>    <dbl>    <dbl>     <dbl> <dbl> <chr>      <chr>       
  1 Species setosa     versicolor          0     2.80     2.59      3.00 3e-15 ****     A       
  2 Species setosa     virginica           0     4.09     3.89      4.29 3e-15 ****     B   
  3 Species versicolor virginica           0     1.29     1.09      1.50 3e-15 ****     C

额外：如果可能的话，请帮我把组的字母放在箱线图的每个盒子上。

Answer 1

Tukey <- iris %>%
  tukey_hsd(Petal.Length ~ Species) %>%
  add_significance() %>% add_column(Group=c("A","B","C"))