我正在尝试使用中介程序包对多级数据(1-1-1模型)进行中介分析。我已经省略了丢失的数据,将所有变量集中在组内,并为我的预测变量和中介变量指定了随机截距和斜率。
我的lmer模型运行,但是当我尝试使用中介命令时出现以下错误:
Error in array(NA, c(J, K)) : 'dims' cannot be of length 0
这是我的代码:
model.m <- lmer(M ~ X + (1 + X|group), data = mydata)
model.y <- lmer(Y ~ X + M + (1 + X|group) + (1 + M|group), data = mydata)
med.out <- mediate(model.m, model.y, treat = "X",
mediator = "M", sims=50, group.out="group")
有人对为什么我遇到阵列错误以及如何解决它有想法吗?
谢谢!
劳拉