While I have used frequentist mixed model for the last several years, I have only delved into Bayesian approaches over the last few. Bayesian analysis and reporting are still somewhat rare in the fields I work in, so I’m unsure of the validity of some of my reporting techniques.

I typically provide a regression summary table with mean/median coefficient/parameter estimates and an 89% CI.

In cases of logistic models (binomial, bernoulli, ordinal), I have been using the fantastic `tidybayes`

package to generate fitted draws (response probabilities) and then with guidance from the research team designating ROPE boundaries as a way of presenting modelling results in an intuitive manner.

Is this approach valid?

For example:

This is from a bernoulli family model:

To produce this plot, I use `add_fitted_draws`

to add predicted response probabilities and then use `compare_levels`

to create the probability differences. Then I manually add the ROPE boundaries to the figure and accompany the figure with a table of percentages of the probability difference distributions inside the ROPE computed from in `bayestestR`

.