compare_fit: compare model fit for 'modsem' models

In modsem: Latent Interaction (and Moderation) Analysis in Structural Equation Models (SEM)

compare model fit for modsem models

Description

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

Usage

compare_fit(est_h1, est_h0, ...)

Arguments

est_h1	object of class modsem_da or modsem_pi representing the alternative hypothesis model (with interaction terms).
est_h0	object of class modsem_da or modsem_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 to lavaan::lavTestLRT. currently only used for modsem_pi models.

Examples

## Not run: 
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)

## End(Not run)

modsem documentation built on June 13, 2025, 9:08 a.m.