我尝试使用“中介”包运行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)