if (FALSE) {
H0 <- "
 # Outer Model
 X =~ x1 + x2 + x3
 Y =~ y1 + y2 + y3
 Z =~ z1 + z2 + z3
 # Inner model
 Y ~ X + Z
"
estH0 <- modsem(m1, oneInt, "lms")
H1 <- "
 # Outer Model
 X =~ x1 + x2 + x3
 Y =~ y1 + y2 + y3
 Z =~ z1 + z2 + z3
 # Inner model
 Y ~ X + Z + X:Z
"
estH1 <- modsem(m1, oneInt, "lms")
compare_fit(estH0, estH1)
}
Run the code above in your browser using DataLab