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?