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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

目前的代码没有添加协变量，求问如何往里添加协变量
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)

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：mediation ATION Media 中介效应 IAT

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

直接往模型里面添加协变量就可以啊
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)

您好，想请教一下语句里的sims=50，这个50是怎么定义的呢？我看随着这个数值不同跑出来结果也不一样，是和样本量有关吗？

换个昵称 发表于 2019-9-16 19:48
您好，想请教一下语句里的sims=50，这个50是怎么定义的呢？我看随着这个数值不同跑出来结果也不一样，是和样 ...

就是做bootstrap的次数，一般推荐1000或5000

橙橙啊 发表于 2019-6-15 19:46
目前的代码没有添加协变量，求问如何往里添加协变量
library(mediation)
model.m

好的好的好的好的