r - 混合模型的事后测试 - 组显示功能

Question

我想为线性混合模型运行事后测试。我最初使用过lsmeans∷lsmeans，效果很好。问题是：使用lsmeans，结果没有按组“排序”（例如：a ab bc 等），至少不是我所知道的。所以我手动对组进行了排序。由于我的实验设计中有许多治疗方法，并且要运行许多事后测试，这在某种程度上是痛苦的（并且容易出错）。

我看过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)