用R进行中介效应分析

中介效应分析用mediation包来完成。分析之前需要先安装mediation包，安装方法，请在R控制台输入：install.packages("mediation")。

脚本与注释

library(mediation)   # 调用mediation包
b <- lm(job_seek ~ sex + age+ econ_hard+ treat, data=jobs)             # 自变量到中介变量的回归方程
c <- lm(depress2 ~ sex + age+ econ_hard+ treat + job_seek, data=jobs)  # 自变量+中介变量到因变量的回归返程
contcont <- mediate(b, c, sims=50, treat="treat", mediator="job_seek") # 采用准贝叶斯近似估计法
summary(b)                                                             # 显示b的结果
summary(c)                                                             # 显示c的结果
summary(contcont)                                                      # 显示中介效应分析的结果