# \donttest{
# note: iter = 250 for demonstrative purposes
# data
Y <- bfi[complete.cases(bfi),]
# males and females
Ymale <- subset(Y, gender == 1,
select = -c(gender,
education))[,1:10]
Yfemale <- subset(Y, gender == 2,
select = -c(gender,
education))[,1:10]
# fit model
fit <- ggm_compare_estimate(Ymale, Yfemale,
type = "ordinal",
iter = 250,
prior_sd = 0.25,
progress = FALSE)
plot(summary(fit))
# }
Run the code above in your browser using DataLab