Analiza wielu mediacji w R.

19

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/

wmmurrah
źródło
Nie jestem pewien, czy semPLS czy plspm pozwolą na analizę mediacyjną, ale warto to sprawdzić.
chl
@wmmurrah będzie na barana na twoje pytanie, czy jedyną zaletą bootstrapowania jest uzyskanie przedziałów ufności?
lf_araujo
1
@If_araujo Jeśli jesteś w trakcie testowania hipotez, zamiast wartości p należy zastosować przedziały ufności ładowania początkowego, ponieważ te ostatnie wymagają często naruszonych założeń normalności. Efekt pośredni, który jest iloczynem dwóch współczynników ścieżki, jest zwykle wypaczony, co sprawia, że ​​założenia wartości p są wątpliwe, chyba że są duże. Więc nawet jeśli nie chcesz używać interwałów, są one lepsze od wartości p.
wmmurrah
@ If_araujo patrz: Patrz: Preacher, KJ i Hayes, AF (2008). Strategie asymptotyczne i resamplingowe do oceny i porównania efektów pośrednich w wielu modelach mediatorów. Metody badań behawioralnych, 40 (3), 879–891.
wmmurrah

Odpowiedzi:

5

Pakiet lavaan to pakiet R dla SEM. Można go użyć do przetestowania hipotezy o wielu mediacjach, a istnieje też boostrap.

użytkownik34609
źródło
Niedawno stałem się wielkim fanem lawy. Zobacz moją aktualizację powyższego pytania. Naprawdę mam nadzieję, że rozwój lawy będzie kontynuowany!
wmmurrah