Function to Plot BPS Hazard Rates

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 |

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

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

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

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.