Generalized extreme value distribution

From Wikipedia, the free encyclopedia

Generalized extreme value
Probability density function
Cumulative distribution function
Parameters	$\mu \in [-\infty,\infty] \,$ location (real) $\sigma \in (0,\infty] \,$ scale (real) $\xi\in [-\infty,\infty] \,$ shape (real)
Support	$x>\mu-\sigma/\xi\,\;(\xi > 0)$ $x<\mu-\sigma/\xi\,\;(\xi < 0)$ $x \in [-\infty,\infty]\,\;(\xi = 0)$
Probability density function (pdf)	$\frac{1}{\sigma}(1\!+\!\xi z)^{-1/\xi-1}e^{-(1\!+\!\xi z)^{-1/\xi}}$ where $z=\frac{x-\mu}{\sigma}$
Cumulative distribution function (cdf)	$e^{-(1+\xi z)^{-1/\xi}}$
Mean	$\mu-\frac{\sigma}{\xi}+\frac{\sigma}{\xi}g_1$ where $g k = Γ(1 - k ξ)$
Median	$\frac{\ln^{-\xi}(2)-1}{\xi}$
Mode	$\frac{(1+\xi)^{-\xi}-1}{\xi}$
Variance	$\frac{\sigma^2}{\xi^2}(g_2-g_1^2)$
Skewness	$\frac{-g_3+3g_1g_2-2g_1^3}{(g_2-g_1^2)^{3/2}}$
Excess kurtosis	$\frac{g_4-4g_1g_3+6g_2g_1^2-3g_1^4}{(g_2-g_1^2)^{2}}$
Entropy
Moment-generating function (mgf)
Characteristic function

In probability theory and statistics, the generalized extreme value distribution (GEV) is a family of continuous probability distributions developed within extreme value theory to combine the Gumbel, Fréchet and Weibull families also known as type I, II and III extreme value distributions. Its importance arises from the fact that it is the limit distribution of the maxima of a sequence of independent and identically distributed random variables. Because of this, the GEV is used as an approximation to model the maxima of long (finite) sequences of random variables.

1 Specification
2 Mean, standard deviation, mode, skewness and kurtosis excess
3 Link to Fréchet, Weibull and Gumbel families
4 Extremal types theorem
5 References

[edit] Specification

The generalized extreme value distribution has cumulative distribution function

$F(x;\mu,\sigma,\xi) = \exp\left\{-\left[1+\xi\left(\frac{x-\mu}{\sigma}\right)\right]^{-1/\xi}\right\}$

for $1 + ξ(x - μ) / σ > 0$ , where $\mu\in\mathbb R$ is the location parameter, $σ > 0$ the scale parameter and $\xi\in\mathbb R$ the shape parameter.

The density function is, consequently

$f(x;\mu,\sigma,\xi) = \frac{1}{\sigma}\left[1+\xi\left(\frac{x-\mu}{\sigma}\right)\right]^{-1/\xi-1}$

$\exp\left\{-\left[1+\xi\left(\frac{x-\mu}{\sigma}\right)\right]^{-1/\xi}\right\}$

again, for $1 + ξ(x - μ) / σ > 0$ .

[edit] Mean, standard deviation, mode, skewness and kurtosis excess

$\operatorname{E}(X) = \mu-\frac{\sigma}{\xi}+\frac{\sigma}{\xi}g_1$

$\operatorname{Var}(X) = \frac{\sigma^2}{\xi^2}(g_2-g_1^2)$

$\operatorname{Mode}(X) = \mu+\frac{\sigma}{\xi}[(1+\xi)^{-\xi}-1]$

The skewness is

$\operatorname{skewness}(X) = \frac{-g_3+3g_1g_2-2g_1^3}{(g_2-g_1^2)^{3/2}}$

The kurtosis is:

$\operatorname{kurtosis}(X) = \frac{g_4-4g_1g_3+6g_2g_1^2-3g_1^4}{(g_2-g_1^2)^{2}}$

where $g k = Γ(1 - k ξ)$ , k=1,2,3,4, and $Γ(t)$ is Gamma function.

[edit] Link to Fréchet, Weibull and Gumbel families

The factual accuracy of this article is disputed.
Please see the relevant discussion on the talk page.

The shape parameter $ξ$ governs the tail behaviour of the distribution. The sub-families defined by $\xi\to 0$ , $ξ > 0$ and $ξ < 0$ correspond, respectively, to the Gumbel, Fréchet and Weibull families, whose cumulative distribution functions are displayed below.

Gumbel or type I extreme value distribution

$F(x;\mu,\sigma)=e^{-e^{-(x-\mu)/\sigma}}\;\;\; for\;\; x\in\mathbb R$

Fréchet or type II extreme value distribution

$F(x;\mu,\sigma,\alpha)=\begin{cases} 0 & x\leq \mu \\ e^{-((x-\mu)/\sigma)^{-\alpha}} & x>\mu \end{cases}$

Weibull or type III extreme value distribution

$F(x;\mu,\sigma,\alpha)=\begin{cases} e^{-(-(x-\mu)/\sigma)^{\alpha}} & x<\mu \\ 1 & x\geq \mu \end{cases}$

where $σ > 0$ and $α > 0$ .

Remark I: For reliability issues the Weibull distribution by itself is used with the variable $t = μ - x$ , which gives a strictly positive support - in contrast to the use in extreme value theory.

Remark II: Note the differences in the ranges of interest for the three extreme value distributions: Gumbel is unlimited, Fréchet has a lower limit, while Weibull has an upper limit.

One can link the type I to types II and III the following way: if the cumulative distribution function of some random variable $X$ is of type II: $F (x;0,σ,α)$ , then the cumulative distribution function of $ln X$ is of type I, namely $F (x;lnσ,1 / α)$ . Similarly, if the cumulative distribution function of $X$ is of type III: $F (x;0,σ,α)$ , the cumulative distribution function of $ln X$ is of type I: $F (x;lnσ, - 1 / α)$ .

[edit] Extremal types theorem

Credit for the extremal types theorem (or convergence to types theorem) is given to Gnedenko (1948), previous versions were stated by Fisher and Tippett in 1928 and Fréchet in 1927.

Let $X_1,X_2\ldots$ be a sequence of independent and identically distributed random variables, let $M_n=\max\{X_1,\ldots,X_n\}$ . If two sequences of real numbers $a n, b n$ exist such that $a n > 0$ and

$\lim_{n \to \infty}P\left(\frac{M_n-b_n}{a_n}\leq x\right) = F(x)$

then if $F$ is a non degenerate distribution function, it belongs to either the Gumbel, the Fréchet or the Weibull family.

Clearly, the theorem can be reformulated saying that $F$ is a member of the GEV family.

It is worth noting that the result, which is stated for maxima, can be applied to minima by taking the sequence $- X n$ instead of the sequence $X n$ .

For the practical application this theorem means: For samples taken from a well behaving, arbitrary distribution $X$ the resulting extreme value distribution $M n$ can be approximated and parametrised with the extreme value distribution with the appropriate support.

Thus the role of extremal types theorem for maxima is similar to that of central limit theorem for averages. The latter states that the limit distribution of arithmetic mean of a sequence $X n$ of random variable is the normal distribution no matter what the distribution of the $X n$ , The extremal types theorem is similar in scope where maxima is substituted for average and GEV distribution is substituted for normal distribution.

[edit] References

Leadbetter, M.R., Lindgreen, G. and Rootzén, H. (1983). Extremes and related properties of random sequences and processes. Springer-Verlag. ISBN 0-387-90731-9.
Resnick, S.I. (1987). Extreme values, regular variation and point processes. Springer-Verlag. ISBN 0-387-96481-9.
Coles, Stuart (2001). An Introduction to Statistical Modeling of Extreme Values,. Springer-Verlag. ISBN 1-85233-459-2.

	Probability distributions [ view • talk • edit ]
	Univariate	Multivariate
Discrete:	Benford • Bernoulli • binomial • Boltzmann • categorical • compound Poisson • degenerate • Gauss-Kuzmin • geometric • hypergeometric • logarithmic • negative binomial • parabolic fractal • Poisson • Rademacher • Skellam • uniform • Yule-Simon • zeta • Zipf • Zipf-Mandelbrot	Ewens • multinomial • multivariate Polya
Continuous:	Beta • Beta prime • Cauchy • chi-square • Dirac delta function • Erlang • exponential • exponential power • F • fading • Fisher's z • Fisher-Tippett • Gamma • generalized extreme value • generalized hyperbolic • generalized inverse Gaussian • Half-Logistic • Hotelling's T-square • hyperbolic secant • hyper-exponential • hypoexponential • inverse chi-square (scaled inverse chi-square)• inverse Gaussian • inverse gamma (scaled inverse gamma) • Kumaraswamy • Landau • Laplace • Lévy • Lévy skew alpha-stable • logistic • log-normal • Maxwell-Boltzmann • Maxwell speed • normal (Gaussian) • normal inverse Gaussian • Pareto • Pearson • polar • raised cosine • Rayleigh • relativistic Breit-Wigner • Rice • shifted Gompertz • Student's t • triangular • type-1 Gumbel • type-2 Gumbel • uniform • Variance-Gamma • Voigt • von Mises • Weibull • Wigner semicircle • Wilks' lambda	Dirichlet • inverse-Wishart • Kent • matrix normal • multivariate normal • multivariate Student • von Mises-Fisher • Wigner quasi • Wishart
Miscellaneous:	Cantor • conditional • exponential family • infinitely divisible • location-scale family • marginal • maximum entropy • phase-type • posterior • prior • quasi • sampling • singular