require(ordinal)
require(emmeans)
model <- clm(rating~contact*temp,data=wine)
EMM <- emmeans(model,~contact:temp|cut,mode="linear.predictor")
# Probabilities
rating.emmeans(EMM)
# Cumulative probabilities
rating.emmeans(EMM,type="cumprob")
# Most probable rating
rating.emmeans(EMM,type="class1")
Run the code above in your browser using DataLab