我目前正在将lavaan 包用于R
结构方程模型。我想计算每个路径系数的效果大小(即部分 eta 平方)。是否已经有一个包可以做到这一点?
c
例如,如何计算a
和b
回归系数的影响大小?
set.seed(1234)
X <- rnorm(100)
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)
理想情况下,该方法在基于潜在变量构建模型时也应该有效。