混合模型的事后检验 - 群体显示功能

Question

我想为线性混合模型运行事后测试。我最初使用lsmeans∷lsmeans，这很好用。问题是：使用lsmeans，结果不会在组中“排序”（例如：a 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)