r - 如何将 Tukey'sHSD 的“重要性”直接合并到 R 中 ggplot2 的图中？

Question

我有以下数据（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 中。