如何(正确地)使用 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 原样,如果认为这是最佳做法,我可以将其删除。
我正在尝试使用 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 原样,如果认为这是最佳做法,我可以将其删除。