次のように、lsmeans を使用して線形混合効果モデルを実行するとします。
model <- lmer (variable ~ A * B + (1|C), data = data)
3つのレベルを持つ「A」の重要な主効果を見つけたところ。次に、「lsmeans」パッケージを使用して 3 つのレベルを比較する事後テストを実行しました。
lsmeans (model, pairwise~A, adjust = "Tukey")
Aのレベルを比較するときに、この事後比較でもランダム係数「C」が考慮されるのでしょうか?