0

我正在使用包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工作？

4

0 回答 0