Значения баров и звездочек для коробок

Question

Кто-нибудь создал цикл или функцию для рисования значимых баров (со звездочками *,, *) для нескольких блокпостов сравнения?

Я использую фрейм данных с несколькими переменными (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.

Я ценю любую помощь. Некоторое время я буду продолжать искать варианты.

Bests