Значения баров и звездочек для коробок
Кто-нибудь создал цикл или функцию для рисования значимых баров (со звездочками *,, *) для нескольких блокпостов сравнения?
Я использую фрейм данных с несколькими переменными (19) и следующий код:
: присущ рефлективные, вербальные
boxplot$STATUS - это переменные с 4 уровнями (A, B, C, D)
par(mfrow=c(5,4))
for (i in c(1:19)){
boxplot(log(cito[,i]+1)~cito$STATUS)
}
TukeyHDT пост-тест
TukeyHSD(aovp(log(cito$CRP+1)~cito$STATUS,center=T))
Итак, я хотел бы представить в цикле коробчатого графика результаты значимых различий между группами (cito$STATUS)
Получено на TukeyHSD.
Я ценю любую помощь. Некоторое время я буду продолжать искать варианты.