How to get conditional effects and posterior predictions for the zero inflated (zi) part of a zero_inflated_poisson model

Is it possible to get a posterior_prediction of the zeroinflated part of the zero_inflated_poisson model? I have a model that defines the poisson part and the zero inflated part as follows:

mdl <- brms::bf(y ~ year + s(age), zi ~ s(age))

I am able to get the predictions from the model just fine, but I would also like to explore just the zero inflated part of the model. I tired:

brms::posterior_predict(fit, resp = "zi")

but that did not seem to do anything.

A related issue is that I can’t work out how to plot the conditional_effects of the zi part of the model, but I am probably missing something here?

Try posterior_linpred(fit, dpar = "zi") or posterior_epred(fit, dpar = "zi")

1 Like