angcdf {extremis}

Empirical-Likelihood Based Inference for the Angular Measure

Description

This function computes empirical-likelihood based estimators for the angular distribution function of a bivariate extreme value distribution.

Usage

angcdf(Y, tau = 0.95, method = "euclidean", raw = TRUE)

Arguments

Y

data frame with two columns from which the estimate is to be computed.

tau

value used to threshold the data; by default it is set as the 0.95 quantile of the pseudo-radius tau = 0.95.

method

a character string setting the method to be used. By default method = "euclidean", the other option being method = "empirical". See details.

raw

logical; if TRUE, Y will be converted to unit Fréchet scale. If FALSE, Y will be understood as already in unit Fréchet scale.

Details

method = "euclidean" implements the maximum Euclidean likelihood spectral distribution function as introduced by de Carvalho et al (2013). method = "empirical" implements the maximum Empirical likelihood spectral distribution function as introduced by Einmahl and Segers (2009).

Value

H

angular distribution function.

w

pseudo-angles.

Y

data.

The plot method depicts the empirical likelihood-based angular distribution function.

Author(s)

Miguel de Carvalho

References

de Carvalho, M., Oumow, B., Segers, J. and Warchol, M. (2013) A Euclidean likelihood estimator for bivariate tail dependence. Communications in Statistics—Theory and Methods, 42, 1176–1192.

Einmahl, J. H. J., and Segers, J. (2009) Maximum empirical likelihood estimation of the spectral measure of an extreme-value distribution. The Annals of Statistics, 37, 2953–2989.

Examples

## de Carvalho et al (2013, Fig. 7)
data(beatenberg)
attach(beatenberg)
fit <- angcdf(beatenberg, tau = 0.98, raw = FALSE)
plot(fit)
rug(fit$w)