Compute intervals from the posterior predictive distribution.

# S3 method for brmsfit
predictive_interval(object, prob = 0.9, ...)

Arguments

object

An R object of class brmsfit.

prob

A number p (0 < p < 1) indicating the desired probability mass to include in the intervals. Defaults to 0.9.

...

Further arguments passed to posterior_predict.

Value

A matrix with 2 columns for the lower and upper bounds of the intervals, respectively, and as many rows as observations being predicted.

Examples

if (FALSE) {
fit <- brm(count ~ zBase, data = epilepsy, family = poisson())
predictive_interval(fit)
}