(Using brms v2.3.0)

Hi,

If I fit a model using `sample_prior='only'`

, i.e.,

fit <- brm(formula = tp ~ 1 + technique + category + subject,

data = d,

prior = c(set_prior(“normal(0,1)”, class=“Intercept”),

set_prior(“normal(0,1)”, class=“b”)),

family = poisson(link=log),

sample_prior = “only” #To ignore the likelihood and check if the model is sane

)

I would expect that I later could do:

`ppc_dens_overlay(y = d$tp, yrep = posterior_predict(fit, draws=25))`

but unfortunately, I get:

`Error in validate_yrep(yrep, y) : NAs not allowed in 'yrep'.`

Fake data to reproduce this can be found here:

I have a nagging feeling I’ve completely missed something here…