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—