Brms: How to keep it consistent for mu and shape in posterior_linpred (Weibull time-to-event)?

[paul.buerkner]

One option could be to call prepare_predictions directly, which will then internally create the samples for the new random effects levels. Then, the resulting object can be passed to posterior_epred(., dpar = "<dpar>", scale = "linear").