Zastanawiam się, czy ktoś zna sposób na uruchomienie wielu modeli mediacji w R. Wiem, że pakiet mediacji pozwala na wiele prostych modeli mediacji, ale chcę uruchomić jeden model, który ocenia wiele modeli mediacji jednocześnie.
Zakładam, że mogę to zrobić w ramach SEM (analiza ścieżki), ale zastanawiałem się, czy ktoś nowy z pakietu, który obliczył statystyki typowe dla analizy mediacji dla wielu mediatorów (efekty pośrednie, proporcje całkowitego efektu poprzez mediację itp.), i może użyć ładowania początkowego. Wiem, że to długa szansa, ale pomyślałem, że powinienem zapytać, zanim zainwestuję czas w tworzenie od zera.
AKTUALIZACJA: (11.11.2013)
Od czasu, gdy zadałem to pytanie kilka lat temu, nauczyłem się korzystać z cudownego pakietu Lavaan z pakietu R do wykonywania wielu mediacji.
oto przykładowy kod:
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
'
Zauważ, że a1, a2, b1, b2 i c są etykietami. Następnie uruchom model:
fit <- sem(model, data=dataframe)
I spójrz na wynik:
summary(fit, fit.measures=TRUE, standardize=TRUE, rsquare=TRUE)
Na koniec wygeneruj przedziały ufności bootstrap:
boot.fit <- parameterEstimates(fit, boot.ci.type="bca.simple")
Zobacz stronę internetową Lavaan, aby uzyskać więcej informacji: http://lavaan.ugent.be/
źródło
Odpowiedzi:
Pakiet lavaan to pakiet R dla SEM. Można go użyć do przetestowania hipotezy o wielu mediacjach, a istnieje też boostrap.
źródło