glmGetSamples: Get posterior samples for a specific model configuration for...
In hypergsplines: Bayesian model selection with penalised splines and hyper-g prior

Description Usage Arguments Value Author(s) Examples

Get posterior samples for a specific model configuration for generalised response

1 2	glmGetSamples(config, modelData, mcmc = getMcmc(), computation = getComputation())

`config`	the model configuration vector
`modelData`	the data necessary for model estimation, which is the result from `glmModelData`
`mcmc`	MCMC options, result from `getMcmc`
`computation`	computation options produced by `getComputation`

A list with samples from the shrinkage hyperparameter t = g / (g + 1) and the (linear and spline) coefficients, and further information.

Daniel Sabanes Bove daniel.sabanesbove@ifspm.uzh.ch

## get some data
attach(longley)

## get the model data
md <- glmModelData(y=as.numeric(Employed > 64),
                   X=cbind(GNP, Armed.Forces),
                   family=binomial)

## and get samples for one specific model
samples <- glmGetSamples(config=c(2L, 3L),
                         modelData=md,
                         computation=
                         getComputation(higherOrderCorrection=FALSE))

str(samples)