Как включить "значимость" HSD Тьюки непосредственно в графы ggplot2 в R?
У меня есть следующие данные (дата)
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.