otherPORT.q: Other peaks over random threshold high quantile estimate

In evt0: Mean of Order P, Peaks over Random Threshold Hill and High Quantile Estimates

Other peaks over random threshold high quantile estimate

Description

This function computes peaks over random threshold (PORT) high quantile or value-at-risk (VaR) based on moment (MO), generalized Hill (GH) and mixed moment (MM) extreme value index (EVI) estimates.

Usage

otherPORT.q(x, k, q1, q2, method = c("MO", "GH", "MM"))

Arguments

x	Data vector.
k	a vector of number of upper order statistics.
q1	quantile for PORT.
q2	quantile level.
method	Method used, moment estimate("MO", default), generalized Hill ("GH") and mixed moment ("MM").

Details

The computation of high quantile estimate is based on the method by Weissman (1978) and the EVI estimators are given in Dekkers et al. (1989), Beirlant et al. (1996) and Fraga Alves et al. (2009).

Value

a k dimensional vector of PORT EVI and high quantil estimates.

Author(s)

B G Manjunath [e-mail]

References

Araujo Santos, P., Fraga Alves, M.I. and Gomes, M.I. (2006). Peaks over random threshold methodology for tail index and quantile estimation. Revstat, 4(3), 227–247.

Dekkers, A., Einmahl, J. and L. de Haan. (1989). A moment estimator for the index of an extreme-value distribution. Ann. Statist., 17, 1833– 1855.

Beirlant, J., Vynckier, P. and Teugels, J. (1996). Excess functions and estimation of the extreme-value index. Bernoulli, 2, 293–318.

Fraga Alves, M.I., Gomes, M.I., de Haan, L. and Neves, C. (2009). The mixed moment estimator and location invariant alternatives. Extremes, 12, 149–185.

Weissman, I. (1978). Estimation of parameters and large quantiles based on the k largest observations. J. Amer. Statist. Assoc., 73, 812– 815.

See Also

other.EVI

Examples

# generate random samples               
x = rfrechet(50000, loc = 0, scale = 1,shape = 1/0.5)

# estimate PORT EVI and high quantile at level q2
otherPORT.q(x,c(500,5000),0.1,0.5,"MO")

evt0 documentation built on April 22, 2023, 1:15 a.m.