如何（正确）使用lme4 / nlme指定交叉级别效果

Question

我正在尝试使用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。我会使用读数（如果有的话）将回归方程式转换为语法。到目前为止，我没有找到任何。

Answer 1

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

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

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