# NOT RUN {
# }
# NOT RUN {
# generate some data
income_options <- c("below_20", "20_to_40", "40_to_100", "greater_100")
income <- factor(sample(income_options, 100, TRUE),
levels = income_options, ordered = TRUE)
mean_ls <- c(30, 60, 70, 75)
ls <- mean_ls[income] + rnorm(100, sd = 7)
dat <- data.frame(income, ls)
# fit a simple monotonic model
fit <- brm(ls ~ mo(income), data = dat)
# summarise the model
summary(fit)
plot(fit, N = 6)
plot(marginal_effects(fit), points = TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab