evt0 (version 1.1.5)

otherPORT.q: 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"))

Value

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

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").

Author

B G Manjunath

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).

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")