evir (version 1.7-4)

gpd.sfall: Add Expected Shortfall Estimates to a GPD Plot

Description

Calculates expected shortfall (tail conditional expectation) 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.sfall(x, pp, ci.p = 0.95, like.num = 50)

Arguments

x

a list object returned by plot.gpd or tailplot

pp

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

ci.p

probability for confidence interval (must be less than 0.999)

like.num

number of times to evaluate profile likelihood

Details

Expected shortfall is the expected size of the loss, given that a particular quantile of the loss distribution is exceeded. The GPD approximation in the tail is used to estimate expected shortfall. The likelihood is reparametrised in terms of the unknown expected shortfall and profile likelihood arguments are used to construct a confidence interval.

See Also

gpd, plot.gpd, tailplot, gpd.q

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 
# }
# NOT RUN {
gpd.sfall(tp, 0.999)
# }
# NOT RUN {
# Estimates associated expected shortfall for Danish fire losses 
# }