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上有很多很好的信息。