data(Boston, package = "MASS")
df <- within(Boston,
medv <- cut(medv,
breaks = c(0, 10, 15, 20, 25, 50),
ordered = TRUE))
fit(medv ~ ., data = df, model = POLRModel)
Run the code above in your browser using DataLab