posterior_predict {bayesGAM}

Posterior predictive samples from models fit by bayesGAM

Description

Draw from the posterior predictive distribution

Usage

posterior_predict(object, ...)

## S4 method for signature 'bayesGAMfit'
posterior_predict(object, draws = NULL, ...)

## S4 method for signature 'glmModel'
posterior_predict(object, draws = NULL, results = NULL, ...)

Arguments

object

Object of type bayesGAMfit generated from bayesGAM.

...

Additional arguments for postrior_predict

draws

An integer indicating the number of draws to return. The default and maximum number of draws is the size of the posterior sample.

results

Matrix of HMC posterior samples

Value

a list of D by N matrices, where D is the number of draws from the posterior predictive distribution and N is the number of data points being predicted per draw.

References

Goodrich B, Gabry J, Ali I & Brilleman S. (2020). rstanarm: Bayesian applied regression modeling via Stan. R package version 2.19.3 https://mc-stan.org/rstanarm.

Jonah Gabry, Ben Goodrich and Martin Lysy (2020). rstantools: Tools for Developing R Packages Interfacing with 'Stan'. https://mc-stan.org/rstantools/, https://discourse.mc-stan.org/.

Examples

f <- bayesGAM(weight ~ np(height), data = women, 
              family = gaussian, iter=1000, chains = 1)
res <- posterior_predict(f, draws=100)