BayHaz (version 0.1-3)

BPSpost2mcmc: 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

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)