Rでの複数のメディエーション分析

Rで複数のメディエーションモデルを実行する方法を誰かが知っているのではないかと思っています。メディエーションパッケージでは複数の単純なメディエーションモデルを使用できますが、複数のメディエーションモデルを同時に評価する1つのモデルを実行したいと思います。

私はこれをSEMフレームワーク（パス分析）で行うことができると仮定していますが、複数のメディエーターのメディエーション分析に典型的な統計を計算するパッケージの新しい人（間接効果、メディエーションによるトータル効果の割合など）、ブートストラップを利用できます。私はこれが長いショットであることを知っていますが、ゼロから開発する時間を投資する前に尋ねるべきだと思いました。

更新日：（11/11/2013）

数年前にこの質問をして以来、私は素晴らしい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のWebサイトを参照してください：http ://lavaan.ugent.be/

lavaanパッケージは、SEM用のRパッケージです。これを使用して、複数のメディエーション仮説をテストできます。ブースターがあります。