r - 嵌套 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)