Специальный тест на смешанной модели - функция отображения групп
Я хочу запустить специальные тесты для линейной смешанной модели. Я изначально использовал lsmeans∷lsmeans, который работает отлично. Проблема: с lsmeansрезультат не "отсортирован" по группам (скажем, ab b c и т. д.), по крайней мере, я не знаю об этом. Итак, я отсортировал группы вручную. Поскольку в моем экспериментальном проекте есть много методов лечения и много специальных тестов для запуска, это несколько болезненно (и подвержено ошибкам).
Я смотрел на agricolae∷HSD.test, который дает группы для специальных тестов. Недостатком является то, что HSD.test не работает со смешанными моделями multcomp::cld, multcomp::hsdни multcomp::multcompLetters в этом отношении).
РЕДАКТИРОВАТЬ: ОТВЕТИТЬ
# here, the groups are sorted so that the highest mean is "a"
library(lsmeans)
comp_2_by_2<- lsmeans(model, pairwise ~ FactorA:FactorB)
groups_HSD<- cld(comp_2_by_2, alpha=0.05, Letters=c("a","b","c","d"), reverse=TRUE)