compare model fit for modsem
models
compare_fit.Rd
Compare the fit of two models using the likelihood ratio test (LRT).
hypothesis model, and est_h1
the alternative hypothesis model. Importantly,
the function assumes that est_h0
does not have more free parameters
(i.e., degrees of freedom) than est_h1
.
alternative hypothesis model
Arguments
- est_h1
object of class
modsem_da
ormodsem_pi
representing the alternative hypothesis model (with interaction terms).- est_h0
object of class
modsem_da
ormodsem_pi
representing the null hypothesis model (without interaction terms).- ...
additional arguments passed to the underlying comparison function. E.g., for
modsem_pi
models, this can be used to pass arguments tolavaan::lavTestLRT
. currently only used formodsem_pi
models.
Examples
if (FALSE) { # \dontrun{
m1 <- "
# Outer Model
X =~ x1 + x2 + x3
Y =~ y1 + y2 + y3
Z =~ z1 + z2 + z3
# Inner model
Y ~ X + Z + X:Z
"
# LMS approach
est_h1 <- modsem(m1, oneInt, "lms")
est_h0 <- estimate_h0(est_h1, calc.se=FALSE) # std.errors are not needed
compare_fit(est_h1 = est_h1, est_h0 = est_h0)
# Double centering approach
est_h1 <- modsem(m1, oneInt, method = "dblcent")
est_h0 <- estimate_h0(est_h1, oneInt)
compare_fit(est_h1 = est_h1, est_h0 = est_h0)
# Constrained approach
est_h1 <- modsem(m1, oneInt, method = "ca")
est_h0 <- estimate_h0(est_h1, oneInt)
compare_fit(est_h1 = est_h1, est_h0 = est_h0)
} # }