Hello,
I have a mixed-effect models I fit with brms and was wondering if there is any way to extract the full epred from conditional_effects.
I wanted to get the full distribution as I wanted to make some bell curve plots for each value of one of my random effects.
For example, where region is a random_effect, and I want to have a plot for each region
conditional_effects(fit_model, effects='Region',re_formula=NULLE)
Cheers!