我尝试使用“中介”包运行HLM中介:

med.fit <- glmer(M ~ treat + control  + (1|subject_id) ,family = binomial(link = "logit"), data = R1_data)

out.fit <- glmer(Y ~ M+ treat + control+ (1 + M|subject_id),family = binomial(link = "logit"), data = R1_data)

med.out <- mediate(med.fit, out.fit, treat = "treat", mediator = "M", sims = 1000)

我收到此错误消息:

[.data.frame(y.data,int.term.name [p])中的错误:选择了未定义的列

如何解决这个问题呢?这是原始数据和代码:

名字(R1_data)

[1] "subject_id"
[3] "Presented_is_solvable"
[5] "JOS"
[17] "Answer_JOS"
[23] "Matrix_Z_score"

库(调解)

med.fit < - glmer(JOS~Matrix_Z_score Presented_is_solvable(1 | subject_id),family =二项式(link =“logit”),data = R1_data)

out.fit < - glmer(Answer_JOS~JOS Matrix_Z_score Presented_is_solvable(1 JOS | subject_id),family =二项式(link =“logit”),data = R1_data)

med.out < - mediate(med.fit,out.fit,treat =“Matrix_Z_score”,mediator =“JOS”,sims = 1000)