BayHaz (version 0.1-3)

CPPevalHR: Function to Evaluate CPP Hazard Rates

Description

A function to evaluate a (prior or posterior) sample of CPP hazard rates on a grid of time points.

Usage

CPPevalHR(time, sample)

Arguments

time

vector of time points where the hazard rates in the sample should be evaluated

sample

sample of CPP hazard rates (as generated by CPPpriorSample or CPPpostSample)

Value

See Also

CPPplotHR

Examples

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

# select a CPP prior distribution
hypars<-CPPpriorElicit(r0 = 0.1, H = 1, T00 = 50, M00 = 2)
# generate a sample of ten hazard rates
prior<-CPPpriorSample(ss = 10, hyp = hypars)

# evaluate the ten hazard rates at year multiples
CPPevalHR(time = seq(0,50), sample = prior)

# load a data set
data(earthquakes)
# generate a posterior sample
post<-CPPpostSample(hypars, times = earthquakes$ti, obs = earthquakes$ob)

# evaluate the posterior hazard rates at year multiples
CPPevalHR(time = seq(0,50), sample = post)