I am a beginner for stan
, and brms
, and I'm having problems with removing the posterior predictive distributions. Say I have a simple logistic regression
fit = brm(y ~ x, family="bernoulli", data=df.training)
where y
is binary and x
continuous. For test data (or even training data), I thought that now I can get a predictive distribution for the Bernoulli probability by p
changing probs
to
predict(fit, df.test, probs=seq(0, 1, 0.1))
However, although the output of this command gives me estimates that are continuous in the range [0,1]
(this makes sense), the values of the confidence interval seem binary (this does not make sense to me) ... How do I get all the subsequent predictive distribution for p
?
source
share