请问各位小伙伴,R语言如何实现自变量为多分类,中介和因变量都是2分类的中介分析。非常感谢
set.seed(1)
X<-rnorm(100) #随机生成X,Y,M
M<-0.5*X+rnorm(100)
Y<-0.7*M+rnorm(100)
Data<-data.frame(X=X,Y=Y,M=M)
## 模型定义
model <- ' # direct effect
Y ~ c*X
# mediator
M ~ a*X
Y ~ b*M
# indirect effect (a*b)
ab := a*b
# total effect
total := c + (a*b)
'
fit <- sem(model, data = Data)
summary(fit)
我们可以得到直接效应(c),间接效应(ab),总效应(total),X-M(a)以及M-Y(b)的效应。