Question

这是关于lme()语法的问题。我的回复变量是“响应”。我的固定变量是＆＃39; year＆＃39;。我有两个随机变量：＆＃39;学生＆＃39;它嵌套在＆＃39; school＆＃39;。

中

我希望包含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)