# lme model
fit1 <- lm_model(
data = popular,
response_variable = popular,
predictor_var = c(sex, extrav)
)
fit2 <- lm_model(
data = popular,
response_variable = popular,
predictor_var = c(sex, extrav),
two_way_interaction_factor = c(sex, extrav)
)
compare_fit(fit1, fit2)
# see ?measurement_invariance for measurement invariance example
Run the code above in your browser using DataLab