如何（正确地）使用 lme4/nlme 指定一个跨关卡效果答案

【问题标题】：How to (properly) specify a cross level effect with lme4/nlme如何（正确地）使用 lme4/nlme 指定一个跨关卡效果
【发布时间】：2018-01-26 17:20:16
【问题描述】：

我正在尝试使用 lmer() 和 lme() 函数正确指定跨级别效果。

我正在使用修改后的 sleepstudy 数据版本，我在其中添加了一个 2 级变量：score 给每个参与者：

library(lme4)
data(sleepstudy)
score <- gl(n = 5, k = 10, length = 180)
sleepstudy$score <- as.numeric(score)

我的目标是找到可以复制以下公式的正确 lme4::lmer()（和 nlme::lme() 类似物）语法：

1 级：
$Reaction = \beta_0 + \beta_1(天)+ r$

2 级：
$ \beta_0= \gamma_{00} + \gamma_{01}（分数）+ u_0$
$ \beta_1= \gamma_{10} + \gamma_{11}（分数）+ u_1$

编辑：我不确定为什么数学方程没有正确渲染，所以这里有一个 img 的样子：

由于某些未知原因，我很难将正规方程转换为 lme4/nlme 语法。

PS。我会使用读取（如果有的话）将回归方程转换为语法。到现在都没找到。

【问题讨论】：

【解决方案1】：

通过反复试验，我找到了给定示例的正确解决方案：

lmer(Reaction ~ score*Days  + (Days|Subject), data = sleepstudy)

如果认为这是最佳做法，我可以保留此 Q/A 或删除。

【讨论】：