我有一个正在执行Kruskal-Wallis测试的 Dataframe 。我想将结果的紧凑字母显示添加到：
1.平均值、sd、Z统计、P值+字母表

1. A显示紧凑字母的结果图
   关于如何在TukeyHSD后进行此操作，有丰富的知识，但我无法使用其他模型（即Kruskal-Wallis Dunn检验、welch ANOVA配对t检验等）在Hocs后找出答案。我曾想过手动进行此操作，但它不起作用。
   这是我目前掌握的情况：

data <- iris
kruskal <- kruskal.test(Petal.Width ~ Species, data = data)
kruskal_ph <- dunnTest(Petal.Width ~ Species,
                  data=data,
                  method="bonferroni")

#Make CLD here ....


#OR

#First make comparisons table and mannually insert letters. E.g.

pvalue <- kruskal_ph[[2]][, "P.adj"]
    hi <- kruskal_ph[[2]][,"Comparison"] %>%
              as.data.frame()
      hi <- hi %>%
          rename(hi, comps = .) %>%
          separate(comps, c("group1","group2"))
      comparisons <- cbind(hi, pvalue)
    comparisons$Sig <- comparisons$Sig
    comparisons$Sig[comparisons$pvalue > 0.05] <- "ns"
    comparisons$Sig[comparisons$pvalue <= 0.05] <- "*"
    comparisons$Sig[comparisons$pvalue <= 0.01] <- "**"
comparisons

# Make CLD from comparisons here

# Add to a summary table (how?) and graph?

library(ggplot2)
ggplot(data, aes(x = Species, y = Petal.Width, colour = Species, fill = Species)) +
  geom_bar(stat = "summary", fun = "mean")