如何(正确地)使用 lme4/nlme 指定跨层效应

How to (properly) specify a cross level effect with lme4/nlme

我正在尝试使用 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 级:
$反应 = \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。我会使用阅读(如果有的话),其中回归方程被翻译成语法。到目前为止我还没有找到。

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

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

我可以保留此 Q/A 原样,如果认为这是最佳做法,我可以将其删除。