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

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
# 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)