SEM/CFA: z-Test oder Likelihood-Ratio-Test (Likelihood-Quotienten-Test)?
Arndt Regorz, Dipl. Kfm. & M.Sc. Psychologie, 14.11.2022
(Hinweis: Mit Anklicken des Videos wird ein Angebot des Anbieters YouTube genutzt.)
Hier finden Sie den R-Code, mit dem die Werte aus dem o.g. Video berechnet worden sind:
library(lavaan)
# 1 Testung Strukturparamter mit z-Test
# 1a Identif. über ersten Indikator jeweils
strukturmodell_1a <- '
# Faktorladungen
Industrialisierung1960 =~ x1 + x2 + x3
Demokratie1960 =~ y1 + y2 + y3 + y4
Demokratie1965 =~ y5 + y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_1a <- sem(data = PoliticalDemocracy, model = strukturmodell_1a)
summary(model_fit_1a, fit.measures = TRUE, standardized = TRUE)
# 1b Identif. über zweiten Indikator jeweils
strukturmodell_1b <- '
# Faktorladungen
Industrialisierung1960 =~ NA * x1 + 1 * x2 + x3
Demokratie1960 =~ NA * y1 + 1 * y2 + y3 + y4
Demokratie1965 =~ NA * y5 + 1 * y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_1b <- sem(data = PoliticalDemocracy, model = strukturmodell_1b)
summary(model_fit_1b, fit.measures = TRUE, standardized = TRUE)
# 1c Faktorvarianzen gleich 1
strukturmodell_1c <- '
# Faktorladungen
Industrialisierung1960 =~ x1 + x2 + x3
Demokratie1960 =~ y1 + y2 + y3 + y4
Demokratie1965 =~ y5 + y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_1c <- sem(data = PoliticalDemocracy, model = strukturmodell_1c,
std.lv = TRUE)
summary(model_fit_1c, fit.measures = TRUE, standardized = TRUE)
# 2 Testung Strukturparamter mit LR-Test
# 2a Identif. über ersten Indikator jeweils
# Vergleichsmodell
strukturmodell_2a <- '
# Faktorladungen
Industrialisierung1960 =~ x1 + x2 + x3
Demokratie1960 =~ y1 + y2 + y3 + y4
Demokratie1965 =~ y5 + y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + 0 * Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_2a <- sem(data = PoliticalDemocracy, model = strukturmodell_2a)
summary(model_fit_2a, fit.measures = TRUE, standardized = TRUE)
lavTestLRT(model_fit_1a, model_fit_2a)
# 2b Identif. über zweiten Indikator jeweils
strukturmodell_2b <- '
# Faktorladungen
Industrialisierung1960 =~ NA * x1 + 1 * x2 + x3
Demokratie1960 =~ NA * y1 + 1 * y2 + y3 + y4
Demokratie1965 =~ NA * y5 + 1 * y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + 0 *Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_2b <- sem(data = PoliticalDemocracy, model = strukturmodell_2b)
summary(model_fit_2b, fit.measures = TRUE, standardized = TRUE)
lavTestLRT(model_fit_1b, model_fit_2b)
strukturmodell_2c <- '
# Faktorladungen
Industrialisierung1960 =~ x1 + x2 + x3
Demokratie1960 =~ y1 + y2 + y3 + y4
Demokratie1965 =~ y5 + y6 + y7 + y8
# Strukturmodell/Regression
Demokratie1960 ~ Industrialisierung1960
Demokratie1965 ~ Demokratie1960 + 0 * Industrialisierung1960
#Kovarianzen
y1 ~~ y5
y2 ~~ y6
y3 ~~ y7
y4 ~~ y8
'
model_fit_2c <- sem(data = PoliticalDemocracy, model = strukturmodell_2c,
std.lv = TRUE)
summary(model_fit_2c, fit.measures = TRUE, standardized = TRUE)
lavTestLRT(model_fit_1c, model_fit_2c)