BPSpost2mcmc {BayHaz}R Documentation

Function to Convert BPS Posterior Samples into MCMC Objects

Description

A function to convert a first order autoregressive BPS posterior sample into an MCMC object, so that package 'coda' can be used for output diagnostics.

Usage

BPSpost2mcmc(sampost)

Arguments

sampost posterior sample of BPS hazard rates (as generated by BPSpostSample)

Value

An MCMC object, complete with burn-in and thinning information.

Note

If package 'coda' is not available, a matrix with meaningful column names is returned.

See Also

BayHaz-package, BPSpostSample

Examples

# set RNG seed (for example reproducibility only)
set.seed(1234)

# select a BPS prior distribution
hypars<-BPSpriorElicit(r0 = 0.1, H = 1, T00 = 50, ord = 4, G = 3, c = 0.9)
# load a data set
data(earthquakes)
# generate a posterior sample
post<-BPSpostSample(hypars, times = earthquakes$ti, obs = earthquakes$ob)

# convert the posterior sample into an MCMC object
MCMCpost<-BPSpost2mcmc(post)

[Package BayHaz version 0.1-3 Index]