на основе этого замечательного сообщения, я пытаюсь запустить собственное модерируемое посредничество, но моя модель имеет дихотомический модератор только на пути b.
Код выглядит следующим образом:
library(mediation)
Model.m<-lm(mpg ~ drat, data = mtcars)
Model.y<-lm(disp ~ mpg*vs + drat, data = mtcars)
Mod.Med.vs1 <- mediate(model.m = Model.m, model.y = Model.y,
covariates = list(vs = 1),
boot = TRUE,
boot.ci.type = "bca", sims = 2000, mediator="mpg", treat="drat")
# Here I get the error
Mod.Med.vs0 <- mediate(model.m = Model.m, model.y = Model.y,
covariates = list(vs = 0),
boot = TRUE,
boot.ci.type = "bca", sims = 2000, mediator="mpg", treat="drat")
Mod.Med <- mediate(model.m = Model.m, model.y = Model.y,
boot = TRUE,
boot.ci.type = "bca", sims = 2000, mediator="mpg", treat="drat")
Я получаю сообщение об ошибке сразу после запуска Mod.Med.vs1
:
Error in `[.data.frame`(pred.data.t, , vl) : undefined columns selected
Важно, что если я умножу drat
на vs
также на Model.m
, он работает нормально, но опять же, я не заинтересован в исследовании влияния модератора до посредника.
Есть предложения?