Generic function to extract samples from objects.

## Generic.
samples(object, ...)

## Method for "bamlss" objects.
# S3 method for bamlss
samples(object, model = NULL, term = NULL,
  combine = TRUE, drop = TRUE, burnin = NULL,
  thin = NULL, coef.only = FALSE, ...)

Arguments

object

An object for which samples should be extracted.

model

Character or integer, specifies the model for which samples should be extracted.

term

Character or integer, specifies the term for which samples should be extracted.

combine

Samples stored as a mcmc.list, e.g., when a model is estimated on multiple cores, can be combined into one large sample matrix.

drop

If there is only one model for which samples should be extracted, should the list structure be dropped?

burnin

Integer, specifies the number of samples that should be withdrawn as a burn-in phase.

thin

Integer, specifies the step length of samples that should be extracted, e.g., thin = 10 mean that only every 10th sample is returned.

coef.only

Logical, should only samples of model coefficients be returned?

Other arguments.

See also

Examples

# NOT RUN {
## Generate data.
d <- GAMart()

## Estimate model.
b <- bamlss(num ~ s(x1) + s(x2) + s(x3), data = d)

## Extract samples for "s(x2)".
sa <- samples(b, term = "s(x2)")
head(sa)

## Trace plot.
plot(sa)
# }