線形混合モデルの事後テストを実行したいと考えています。私は最初に を使用しましたがlsmeans∷lsmeans、これは正常に動作します。問題は次lsmeansのとおりです。少なくとも私が知っていることではありません。そのため、グループを手動で並べ替えました。私の実験計画には多くの処理があり、多くの事後テストを実行するため、これはなんとなく苦痛です (そしてエラーが発生しやすくなります)。
私は を見てきました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)