sample function

Obtain draws from a Bayesian model sampler

sample(x, n_save = 1000L, n_burn = 0L, mh = set_mh(), verbose = TRUE)

Arguments

• x: Bayesian model
• n_save, n_burn: Integer scalar with number of draws to save / burn.
• mh: Settings to tune the Metropolis-Hastings step. See set_mh.
• verbose: Logical scalar. Whether to print status updates.

Returns

Returns a numeric matrix with stored draws. The Bayesian model is modified in place.