Extreme value theory

From Wikipedia, the free encyclopedia
Jump to: navigation, search
Extreme value theory is used to model the risk of extreme, rare events, such as the 1755 Lisbon earthquake.

Extreme value analysis (EVA) is a branch of statistics dealing with the extreme deviations from the median of probability distributions. It seeks to assess, from a given ordered sample of a given random variable, the probability of events that are more extreme than any observed prior. Extreme value analysis is widely used in many disciplines, ranging from structural engineering, finance, earth sciences, traffic prediction, geological engineering, etc. For example, EVA might be used in the field of hydrology to estimate the value an unusually large flooding event, such as the 100-year flood. Similarly, for the design of a breakwater, a coastal engineer would seek to estimate the 50-year wave and design the structure accordingly.

Contents

[edit] Data sampling

Two approaches exist to fit the tail of a sample empirical cumulative distribution function (ECDF)to one of the three possible distribution functions. The first method relies on approximating a distribution from a so-called block maxima (minima) series. In operational statistics situations, it is customary and convenient to apply a sampling method that consists in extracting the annual maxima. In doing so, a so-called "Annual Maxima Series" (AMS) is generated. The second method relies on sampling points from the data set that exceeds a certain threshold (falls below a certain floor). This method is generally referred to as the "Point Over Threshold" method (POT).:

  1. Basic theory approach as described in the book by Burry (1975). In general this conforms to the first theorem in extreme value theory (Fisher and Tippett, 1928; Gnedenko, 1943).
  2. Most common at this moment is the tail-fitting approach based on the second theorem in extreme value theory (Pickands, 1975; Balkema and de Haan, 1974).

The difference between the two theorems is due to the nature of the data generation. For Theorem I the data are generated in full range, while in Theorem II data is only generated when it surpasses a certain threshold, called Peak Over Threshold models (POT). The POT approach has been developed largely in the insurance business, where only losses (pay outs) above a certain threshold are accessible to the insurance company. Strangely, this approach is often used for cases where Theorem I applies, which creates problems with the basic model assumptions.

Extreme value distributions are the limiting distributions for the minimum or the maximum of a very large collection of independent random variables from the same arbitrary distribution. Emil Julius Gumbel (1958) showed that for any well-behaved initial distribution (i.e., F(x) is continuous and has an inverse), only a few models are needed, depending on whether you are interested in the maximum or the minimum, and also if the observations are bounded above or below.

[edit] Applications

Applications of extreme value theory include predicting the probability distribution of:

[edit] History

The field of extreme value theory was pioneered by Leonard Tippett (1902–1985). Tippett was employed by the British Cotton Industry Research Association, where he worked to make cotton thread stronger. In his studies, he realized that the strength of a thread was controlled by the strength of its weakest fibers. With the help of R. A. Fisher, Tippet obtained three asymptotic limits describing the distributions of extremes. The German mathematician and anti-Nazi activist Emil Julius Gumbel codified this theory in his 1958 book Statistics of Extremes, including the Gumbel distributions that bear his name.

A summary of historically important publications relating to extreme values theory can be found on the article List of publications in statistics.

[edit] Univariate theory

[edit] Classical extreme value theory and models

Let X_1, \dots, X_n be a sequence of independent and identically distributed variables with distribution function F and let M_n =\max(X_1,\dots,X_n) denote the maximum.

In theory, the exact distribution of the maximum can be derived:


\begin{align}
\Pr(M_n \leq z) & = \Pr(X_1 \leq z, \dots, X_n \leq z) \\
& = \Pr(X_1 \leq z) \cdots \Pr(X_n \leq z) = (F(z))^n
\end{align}

In practice, we might not have the distribution function F but the Fisher–Tippett–Gnedenko theorem provides the following asymptotic result

If there exist sequences of constants {an > 0} and {bn} such that

 \Pr\{(M_n-b_n)/a_n \leq z\} \rightarrow G(z) as n \rightarrow \infty and G is a non-degenerate distribution then G belongs to one of the following families:
 G(z)=\exp\left\{-\exp\left(-\left(\frac{z-b}{a}\right)\right)\right\}\text{ for }z\in\mathbb R.
 G(z)=\begin{cases} 0 & z\leq b \\ \exp\left\{\left(-\left(\frac{z-b}{a}\right)\right)^{-\alpha}\right\} & z>b. \end{cases}
 G(z)=\begin{cases} \exp\left\{-\left( -\left( \frac{z-b}{a} \right) \right)^\alpha\right\} & z<b \\ 1 & z\geq b. \end{cases}

where α > 0.

[edit] Models for exceedances

[edit] See also

[edit] Citations

  1. ^ Alvardo (1998, p.68.)
  2. ^ "Ultimate 100m World Records Through Extreme-Value Theory", CentER Discussion Paper, Tilburg University, 57, 2009, http://arno.uvt.nl/show.cgi?fid=95436, retrieved 2009-08-12 

[edit] References

[edit] External links

Personal tools
Namespaces
Variants
Actions
Navigation
Interaction
Toolbox
Print/export
Languages