我想知道是否有人知道在 R 中运行多个中介模型的方法。我知道中介包允许多个简单的中介模型，但我想运行一个同时评估多个中介模型的模型。

我假设我可以在 SEM 框架（路径分析）中做到这一点，但想知道是否有任何新的软件包可以计算多个中介的中介分析的典型统计数据（间接效应、通过中介的总效应比例等），并且可以利用自举。我知道这是一个长远的目标，但我认为在投入时间从头开始开发之前我应该​​问一下。

更新：（2013 年 11 月 11 日）

自从几年前提出这个问题以来，我已经学会了使用美妙的 R 包 lavaan 来做多重中介。

这是示例代码：

model <- '
# outcome model 
outcomeVar ~ c*xVar + b1*medVar1 + b2*medVar2

# mediator models
medVar1 ~ a1*xVar 
medVar2 ~ a2*xVar

# indirect effects (IDE)
medVar1IDE  := a1*b1
medVar2IDE  := a2*b2
sumIDE := (a1*b1) + (a2*b2)

# total effect
total := c + (a1*b1) + (a2*b2)
medVar1 ~~ medVar2 # model correlation between mediators
'

请注意，a1、a2、b1、b2 和 c 是标签。然后运行模型：

fit <- sem(model, data=dataframe)

并查看输出：

summary(fit, fit.measures=TRUE, standardize=TRUE, rsquare=TRUE)

最后，生成引导置信区间：

boot.fit <- parameterEstimates(fit, boot.ci.type="bca.simple")

有关详细信息，请参阅 lavaan 网站：http: //lavaan.ugent.be/