使用mediation包
首先构建X->M的模型,根据数据类型进行选择模型
a<-lm(M~X)
其次构建X->M->Y的模型
b<-lm(Y~X+M)
计算直接效应(ADE)、间接效应(ACME)和总效应,中介效应占比为Prop.mediated
set.seed(123)
result<-mediate(a,b,treat="X",
mediator="M",boot = T) #这里treat和mediator的名称,为字符型
summary(result)
首先构建X->M的模型,根据数据类型进行选择模型
a<-lm(M~X)
其次构建X->M->Y的模型
b<-lm(Y~X+M)
计算直接效应(ADE)、间接效应(ACME)和总效应,中介效应占比为Prop.mediated
set.seed(123)
result<-mediate(a,b,treat="X",
mediator="M",boot = T) #这里treat和mediator的名称,为字符型
summary(result)