Pfadanalyse mit R lavaan
10. Mediation im Längsschnitt (2 Messzeitpunkte)
Arndt Regorz, Dipl. Kfm. & MSc. Psychologie, 02.01.2023
Dieses ist eine Begleitseite zum Video-Tutorial über die Schätzung eines Pfadmodells mit Längsschnitt (2 MZP).
(Hinweis: Mit Anklicken des Videos wird ein Angebot des Anbieters YouTube genutzt.)
R-Code aus dem Video-Tutorial
Hier ist der gesamte Code aus dem Video-Tutorial:
library(lavaan)
# Mediationsmodell (2 MZP)
med_modell <- '
# Regressionspfade
IV2 ~ IV1
MED2 ~ MED1 + a* IV1
DV2 ~ DV1 + b * MED1
# Korrelation Disturbances endogene Variablen (t2)
IV2~~ MED2 + DV2
MED2 ~~ DV2
# Benutzerdefinierter Effekt
ind := a * b '
med_fit <- sem(med_modell, data= mediation.längsschnitt.2.mzp.lavaan)
summary(med_fit, fit.measures = TRUE)
modindices(med_fit, minimum.value = 10, sort=TRUE)
# Schätzung mit Bootstrapping
med_fit_b <- sem(med_modell, data= mediation.längsschnitt.2.mzp.lavaan,
se="bootstrap", bootstrap = 5000)
summary(med_fit_b, fit.measures = TRUE)
parameterEstimates(med_fit_b)