Я пытаюсь запустить HLM-посредничество с пакетом "mediation":
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 = binomial (link = "logit"), data = R1_data) </p>
out.fit <- glmer (Answer_JOS ~ JOS + Matrix_Z_score + Presented_is_solvable + (1 + JOS | subject_id), семейство = binomial (ссылка= "logit"), data = R1_data) </p>
med.out <- посредник (med.fit, out.fit, Treat = "Matrix_Z_score", mediator = "JOS", sims = 1000) </p>