我使用包 lme4 :
运行以下混合模型moms.rest.lmer <- lmer(rest ~ habitat + (1|year) + (1|ID))
其中rest是一个数值,而habitat是一个3级的因子变量。我试图使用包 lsmeans 中的lsmeans对3个栖息地水平进行成对比较:
moms.rest.lsm <- lsmeans(moms.rest.lmer, pairwise ~ habitat, adjust="Tukey")
我使用这种方法超过一年没有问题,但我需要重新运行模型,现在我突然遇到这个错误:
Loading required namespace: lmerTest
Failed with error: ‘there is no package called ‘lmerTest’’
Install package 'lmerTest' to obtain Satterthwaite degrees of freedom
我读到包 lmerTest 有自己的lsmeans版本(lsmeansLT),但我不想使用 lmerTest 包,因为它似乎没有做成对比较。如何使用 lsmeans 来使旧代码工作?