BPSpost2mcmc: Function to Convert BPS Posterior Samples into MCMC Objects
In BayHaz: R Functions for Bayesian Hazard Rate Estimation

Description Usage Arguments Value Note See Also Examples

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.

1	BPSpost2mcmc(sampost)

sampost

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

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

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

BayHaz-package, BPSpostSample

# 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)