Predict with brms fit object

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

  • Operating System: MacOS High Sierra
  • brms Version: 2.7.0

I’m running a multilevel model in Stan with one binary predictor and 1 factor-type random intercept and a truncated outcome variable - i.e y | trunc(lb=10, ub = 50) ~ x1 + (1|x2). The model runs fine.

However, I wanted to do some posterior checks. The data looks like this:

At first I ran predict() on the brms fit object and this was the distribution of the data

I don’t think this is an issue with model specification because the pp_check() on the same object looks like this

I thought predict worked on a brms object in same way posterior_predict works on the rstanarm fit objects - have I misunderstood?

Thanks! :)

You may use posterior_predict also with brms or predict(..., summary = FALSE).

Great thanks @paul.buerkner - I knew I was doing something silly!