当 2 个随机效应嵌套时，lme() 单个随机交互项的正确语法

Question

这是一个关于 lme() 句法。我的响应变量是“响应”。我的固定变量是“年”。我有 2 个随机变量：“学生”嵌套在“学校”中。

我想包括一个 year*school 互动，但我不想包括 year*student 一。

这是我到目前为止所使用的语法，但这似乎包括 2 个随机交互，而我只想要 1 个。

lme1 = lme(response ~ 1 + year, random = ~ year | school/student, 
            method = "REML", data = data)

Answer 1

这会更容易做到 lmer 来自 lme4 包裹。假设每个学生都有一个唯一的标识符，它看起来像这样（未经测试）。

lmer(response ~ year + (year | school) + (1 | student), data = data)