我有以下数据（dat）
V W X Y Z
1 8 89 3 900
1 8 100 2 800
0 9 333 4 980
0 9 560 1 999
我希望对上述数据集执行 TukeysHSD 成对测试。从测试结果来看，我想将显着性比较合并到图中（显示显着组之间的“*”或“**”符号）。
这是尝试的代码：
library(ggplot2)
library(reshape2)
dat1 <- gather(dat)
ggplot(data = dat1, aes(x = key, y = value)) + stat_summary(fun.data = "mean_cl_normal", colour = "red", size = 1)
pairwise.t.test(dat1$value, dat1$key, p.adj = "holm")
我不知道是否（我可以）以及如何将测试的结果（“重要性”）直接合并到图中，而不将每个结果保存为外部数组，然后将其调用到 ggplot2 中。