如何用R中的mediation包做中介效应，需要添加协变量

11171

收藏 2019-06-15

目前的代码没有添加协变量，求问如何往里添加协变量
library(mediation)
model.m <- lm(A2M_p99NA_Cubic_spline~sta_pm2.5_day_lag00, data=data_mediation_na.omit)
model.y <- lm(SBP_average~A2M_p99NA_Cubic_spline+sta_pm2.5_day_lag00, data=data_mediation_na.omit)
model<- mediate(model.m, model.y,treat="sta_pm2.5_day_lag00", mediator="A2M_p99NA_Cubic_spline",sims=50)
summary(model)
plot(model)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

admin_kefu

2019-6-19 16:49:39

您好，如果您的求助没有解决，请到项目交易发布需求，会有更快更专业的用户帮助您 https://bbs.pinggu.org/prj/

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

njysta

2019-8-22 23:50:26

直接往模型里面添加协变量就可以啊
model.m <- lm(A2M_p99NA_Cubic_spline~sta_pm2.5_day_lag00+covariates, data=data_mediation_na.omit)
model.y <- lm(SBP_average~A2M_p99NA_Cubic_spline+sta_pm2.5_day_lag00+covariates, data=data_mediation_na.omit)