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

【问题标题】：Proper syntax for lme() single random interaction term when 2 random effects are nested嵌套 2 个随机效应时 lme() 单个随机交互项的正确语法
【发布时间】：2014-06-03 14:33:02
【问题描述】：

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

我想包含一个year*school 交互，但我不想包含一个year*student。

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

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

【问题讨论】：

我看不出学生是如何与学校嵌套的：您能否展示一下您的数据是什么样的？
为了更准确一点，我会说你有两个 分组变量（学校和学生 [校内]），两个随机效应或三个取决于您如何设置问题（按学校、按年、按学生截取，其中前两个可以组合成向量值随机效应）

标签： r interaction mixed-models

【解决方案1】：

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

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

【讨论】：