See this comment in paul-buerkner/brms#1230: Do we want the behavior of predict.refmodel(), proj_linpred(), and proj_predict() to be the same for "stanreg" and "brmsfit" fits? If yes, how should that common behavior look like: The current "stanreg" or the current "brmsfit" behavior?