我正在尝试引导置信区间,以得出我对中介模型的间接影响的点估计。由于它实际上是(第一阶段,a路径)主持的中介,因此我实际上需要做三遍,以找出主持人(以下为b
)在平均值(财政)和±1标准偏差。
我用两种方法估算了模型:
1。使用undefined
程序包。
fiscal
2。使用mediation
程序包。
summary(model_m <- lm(unaccept_i ~ condition*fiscal, data=d1))
summary(model_y <- lm(redis_i ~ condition*fiscal + unaccept_i, data=d1))
mod_med <- mediation::mediate(model_m, model_y,
covariates = list(fiscal = 0),
treat="condition",
mediator="unaccept_i",
boot=T,
sims = 5000)
我可以尝试通过上述两种方法之一访问间接效果(是否可能?尚不清楚),也可以直接自己计算:
psych
唯一的问题是这种方法不会给我带来置信区间。
有人能弄清楚如何自举/蒙特卡洛模拟置信区间,以便我可以绘制条件间接效应(a * b)的平均值和主持人的+1 SD吗?
非常感谢!
PS:这是数据的一部分:
psy_med <- psych::mediate(redis_i ~ condition*fiscal + (unaccept_i),
data=d1, plot=F, n.iter=5000)