# Fit structural equation model
HS.model <- '
visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9
'
fit1 <- lavaan::cfa(HS.model, data = lavaan::HolzingerSwineford1939,
estimator = "MLR")
fit0 <- lavaan::cfa(HS.model, data = lavaan::HolzingerSwineford1939,
orthogonal = TRUE, estimator = "MLR")
# Chi-square difference test
# lavaan::anova(fit1, fit0)
satorraBentlerScaledChiSquareDifferenceTestStatistic(
T0 = lavaan::fitMeasures(fit0)["chisq.scaled"],
c0 = lavaan::fitMeasures(fit0)["chisq.scaling.factor"],
d0 = lavaan::fitMeasures(fit0)["df.scaled"],
T1 = lavaan::fitMeasures(fit1)["chisq.scaled"],
c1 = lavaan::fitMeasures(fit1)["chisq.scaling.factor"],
d1 = lavaan::fitMeasures(fit1)["df.scaled"])
Run the code above in your browser using DataLab