如果我使用lsmeans执行线性混合效果模型,如下所示:
model <- lmer (variable ~ A * B + (1|C), data = data)
我在哪里发现了&#34; A&#34;的重要影响,它有三个等级。然后我进行了一项事后测试,比较三个级别,使用&#34; lsmeans&#34;包:
lsmeans (model, pairwise~A, adjust = "Tukey")
我想知道这种事后比较是否也采用随机因素&#34; C&#34;在比较A?
的水平时要考虑到