# From the help page of lavaan::cfa().
library(lavaan)
HS.model <- '
visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9
'
fit <- cfa(HS.model, data = HolzingerSwineford1939)
vec_rsquare(fit)
vec_sample_vcov(fit)
vec_sample_var(fit)
vec_est_var(fit)
vec_est_se(fit)
HS.model.sem1 <- '
visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9
textual ~ a * visual
speed ~ b * textual
ab := a * b
'
fit_sem1 <- sem(HS.model.sem1, data = HolzingerSwineford1939)
HS.model.sem2 <- '
visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9
textual ~ a * visual
speed ~ b * textual + cp * visual
ab := a * b
'
fit_sem2 <- sem(HS.model.sem2, data = HolzingerSwineford1939)
vec_def_var(fit_sem1)
vec_def_se(fit_sem1)
vec_lavTestLRT(fit_sem1, fit_sem2,
model.names = c("No Direct", "Direct"))
vec_lavTestScore(fit_sem1,
add = "speed ~ visual")
vec_lavTestWald(fit_sem2,
constraints = "cp == 0")
if (requireNamespace("semTools")) {
vec_compRelSEM(fit)
}
Run the code above in your browser using DataLab