sommer 包相当于 lmer(反应 ~ 天数 + (天数 | 主题), 数据 = sleepstudy)

Question

在lme4中，

lmer(Reaction ~ Days + (Days | Subject), data = sleepstudy)

指定了一个线性混合模型，具有按主题随机截距和天数随机斜率，以及相关随机截距和随机斜率。

在 sommer 包中，

mmer(Reaction ~ Days, random= ~ Subject + vsr(Days, Subject), data = sleepstudy), tolParInv = 1e-6, verbose = FALSE)

仅返回不相关的随机截距和随机斜率项，对应于

lmer(Reaction ~ Days + (0 | Subject) + (0 + Days | Subject), data = sleepstudy)

。

对于

lmer(Reaction ~ Days + (Days | Subject), data = sleepstudy)

，是否可以在 sommer 包中指定精确的等价物？

Answer 1

这段代码能给你你想要的吗？

fm2 <- mmer(Reaction ~ Days,
random= ~ Subject + vsr(Days, Subject),
data=DT, tolParInv = 1e-6, verbose = FALSE)

我直接从这里的文档中获取它：