Please also provide the following information in addition to your question:

- Operating System: Ubuntu 19.04
- brms Version: 2.9.0 (Github)

I’m working with a categorical model of the form:

My model is:

mod <- brm(Reason ~ V1 + V2 +V3 + (1|ID),

data = Data ,

family = categorical(link = “logit”),

prior = priors)

I am using predict.brmsfit to return predicted probabilities of each ‘Reason’ with:

predict(mod, newdata = NULL, re_formula = ~ (1|ID),

summary = T)

and using cbind to attach the returned predicted probabilities to the original data frame for visualization.

Is it possible to also return 95% credible intervals surrounding the predicted probabilities?