Мені цікаво, чи хтось знає про спосіб запустити модель множинного посередництва в Р. Я знаю, що пакет медіації дозволяє створювати кілька простих моделей медіації, але я хочу запустити одну модель, яка одночасно оцінює кілька моделей медіації.
Я припускаю, що я можу це зробити в рамках 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: http://lavaan.ugent.be/