extremalindex {extRemes} | R Documentation |
Estimates the extremal index based on the intervals estimator due to Ferro and Segers (2003).
extremalindex(xdat, u)
xdat |
n X 1 numeric vector of the data. |
u |
User defined threshold. May be a single number or a numeric n X 1 vector. |
The extremal index is a useful indicator of how much clustering of exceedances of a threshold occurs in the limit of the distribution. If theta = 1, then the data is independent and if theta < 1, then there is some dependency (clustering) in the limit.
There are many possible estimators of the extremal index. The one used here is the intervals estimator described in Ferro and Segers (2003). It is unbiased in the mean and can be used to estimate the number of clusters, which is also done by this function.
A list with components:
theta |
The intervals estimate of the extremal index (theta). |
msg |
A message stating whether any interexceedance times were above 2 and possibly other pertinent information. |
C |
The estimated number of clusters (handles ties). |
run.length |
The estimated run length for runs declustering. |
Eric Gilleland
Ferro, C.A.T. and Segers, J. (2003): Inference for clusters of extreme values, J.R. Statist. Soc. B, 65 (2): 545-556.
Gilleland, Eric and Katz, Richard W. Tutorial for the 'Extremes Toolkit: Weather and Climate Applications of Extreme Value Statistics.' http://www.assessment.ucar.edu/toolkit, 2005.
dclust
data( Tphap) temp <- -Tphap[,"MinT"] look <- extremalindex( temp, -70) look # See Gilleland et al. (2005) for more.