Calculates the MLEs of the parameters (\(\mu\), \(\sigma\), \(\xi\)), and joint asymptotic covariance matrix of these MLEs over a range of thresholds as supplied by the user.
.Joint_MLE_NHPP(x, u = NULL, k, q1, q2 = 1, par, M)
a list with components
mle matrix of MLEs above the supplied thresholds; columns are (\(\mu\), \(\sigma\), \(\xi\))
Cov.all joint asymptotic covariance matrix of all MLEs
Cov.mu joint asymptotic covariance matrix of MLEs for \(\mu\)
Cov.sig joint asymptotic covariance matrix of MLEs for \(\sigma\)
Cov.xi joint asymptotic covariance matrix of MLEs for \(\xi\)
vector of data
optional vector of thresholds. If not supplied, then k thresholds between quantiles (q1, q2) will be used
number of thresholds to consider if u
not supplied
lower quantile to consider for threshold
upper quantile to consider for threshold. Default to 1
starting values for the optimization
number of superpositions or 'blocks' / 'years' the process corresponds to. It affects the estimation of \(mu\) and \(sigma\), but these can be changed post-hoc to correspond to any number)
Jennifer L. Wadsworth