r - 在 R 中指定 3 级随机截距

Question

我正在使用 R 中的 lmer() 函数(lme4 包)分析一项纵向研究，我在其中测量了 120 名受试者，共测量了 6 次。首先，我指定了一个这样的模型:

library(lme4)
model1 = lmer(DV ~ 1 + X1*X2  + (1+X1|SubjectID), REML="false")

X1 是一个时变变量(level-1)，X2 是一个主题级变量(level-2)。

因为这些主题嵌套在几个团队中，所以我被建议在团队级别(级别 3)包括一个随机截取。但是，我只找到了如何同时包括随机截距和斜率:

model2 = lmer(DV ~ 1 + X1*X2  + (1+X1|TeamID/SubjectID), REML="false")

有谁知道如何只向模型 1 添加 3 级随机截距？

Answer 1

通过使用术语 (1|SubjectID)，您告诉模型仅针对 SubjectID 的不同实例期望不同的基线。为了告诉模型期望对固定效应 X1 的不同响应，我们使用 (1+X1|SubjectID)。因此，您只需要条款

(1|TeamID) + (1+X1|SubjectID)

在你的模型中。

顺便说一下，在 Cross Validated 上有很多关于此的有用信息.