ordinal example

ordinalexample.R

library(datasets)
library(sjmisc)
library(sjPlot)
library(ordinal)

set.seed(111)
ab <-datasets::airquality
ab$orddep<- as.factor(ab$Month)

ord.1 <- clm(orddep ~ Temp  , data=ab)

# # in logits so not so easy to explain
summary(ord.1)

# # helpful since its in odds ratio format and not logits

exp(coef(ord.1))

# # helpful since its in odds ratio format and not logits but also looks better
tab_model(ord.1)

# # helpful to visualize ** this is how i understood it better.
plot_model(ord.1,show.intercept = TRUE,show.values = TRUE,
           show.p = TRUE)