如何(正确)使用nlme 4/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。我将使用read (如果有的话)将回归方程转换为语法。到目前为止我还没有找到任何线索。

Answer 1

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

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

如果这被认为是最佳实践，我可以保留这个Q/A或者删除。