# NOT RUN {
# note: iter = 250 for demonstrative purposes
# data
Y <- bfi
Y <- subset(Y, select = c("E5", "N5",
"gender", "education"))
fit_mv_ordinal <- estimate(Y, formula = ~ gender + as.factor(education),
type = "ordinal",
iter = 250,
progress = FALSE)
regression_summary(fit_mv_ordinal)
# }
Run the code above in your browser using DataLab