Obtain raw ordinal predictions

How can I obtain the posterior probabilities for all ordinal outcomes after fitting an ordinal model? I am looking for the same output as fitted(mod,summary=F), but including the observational variance. When I do predict(mod,summary=F) I get the predicted ordinal categories, not the probabilities for each category.

Interesting question! What are you looking to estimate/quantify/achieve with the draws? Are you examining population- or group-level effects?

Population effects—