Sample an MCMC posterior — sample.nsbi_mcmc

Runs the sampler chosen by posterior.nsbi_nle() or posterior.nsbi_nre(). Draws are cached on the posterior object, so asking for the same or fewer draws from the same observation returns immediately instead of re-running a chain – which is what makes summary() and repeated calls tolerable.

Usage

# S3 method for class 'nsbi_mcmc_posterior'
sample(
  x,
  size = 1000,
  n = size,
  obs = NULL,
  refresh = FALSE,
  verbose = FALSE,
  seed = NULL,
  ...
)

Arguments

x: An nsbi_mcmc_posterior object (named x to satisfy the sample() generic).
size, n: Number of posterior draws (n is an alias for size).
obs: Observation to condition on (defaults to the posterior's x_obs).
refresh: Force a new run even when a cached one would do.
verbose: Report sampling progress.
seed: Optional integer seed for this run. It overrides the seed given to posterior(). A cached run is reused unless refresh = TRUE, so pass both to get a fresh, reproducible chain.
...: Unused.

Value

An n x dim matrix of posterior draws (class nsbi_samples), with convergence diagnostics attached as attribute diagnostics.