次のデータがあります（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に呼び出すことなく、テストの結果（「重要性」）をグラフに直接組み込む方法（できるかどうか）と方法がわかりません。