evir (version 1.7-4)

gpd.q: Add Quantile Estimates to plot.gpd

Description

Calculates quantile estimates and confidence intervals for high quantiles above the threshold in a GPD analysis, and adds a graphical representation to an existing plot.

Usage

gpd.q(x, pp, ci.type = c("likelihood", "wald"), ci.p = 0.95,
    like.num = 50)

Arguments

x

a list object returned by plot.gpd or tailplot

pp

the desired probability for quantile estimate (e.g. 0.99 for the 99th percentile)

ci.type

method for calculating a confidence interval: "likelihood" or "wald"

ci.p

probability for confidence interval (must be less than 0.999)

like.num

number of times to evaluate profile likelihood

Details

The GPD approximation in the tail is used to estimate quantile. The "wald" method uses the observed Fisher information matrix to calculate confidence interval. The "likelihood" method reparametrizes the likelihood in terms of the unknown quantile and uses profile likelihood arguments to construct a confidence interval.

See Also

gpd, plot.gpd, gpd.sfall, tailplot

Examples

# NOT RUN {
data(danish)
# }
# NOT RUN {
out <- gpd(danish, 10)
# }
# NOT RUN {
tp <- tailplot(out)
# }
# NOT RUN {
gpd.q(tp, 0.999)
# }
# NOT RUN {
# Estimates 99.9th percentile of Danish fire losses 
# }