- Generalized extreme value distribution
Probability distribution
name =Generalized extreme value
type =density
pdf_
cdf_
parameters = location (real)
scale (real)
shape (real)
support =
wherez=frac{x-mu}{sigma}
cdf =e^{-(1+xi z)^{-1/xi
mean =mu-frac{sigma}{xi}+frac{sigma}{xi}g_1
whereg_k=Gamma(1-kxi)
median =mu + sigma frac{ln^{-xi}(2)-1}{xi}
mode =mu + sigma 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
g_k =g_k=Gamma(1-kxi)
kurtosis =frac{g_4-4g_1g_3+6g_2g_1^2-3g_1^4}{(g_2-g_1^2)^{2
entropy =
mgf =
char =In
probability theory andstatistics , the generalized extreme value distribution (GEV) is a family of continuousprobability distribution s developed withinextreme 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.pecification
The generalized extreme value distribution has cumulative distribution function
:
F(x;mu,sigma,xi) = expleft{-left [1+xileft(frac{x-mu}{sigma} ight) ight] ^{-1/xi} ight} for
1+xi(x-mu)/sigma>0 , wheremuinmathbb R is the location parameter,sigma>0 the scale parameter andxiinmathbb R the shape parameter.The density function is, consequently
:
f(x;mu,sigma,xi) = frac{1}{sigma}left [1+xileft(frac{x-mu}{sigma} ight) ight] ^{-1/xi-1} :expleft{-left [1+xileft(frac{x-mu}{sigma} ight) ight] ^{-1/xi} ight} again, for
1+xi(x-mu)/sigma>0 .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 excess is::
operatorname{kurtosis excess}(X) = frac{g_4-4g_1g_3+6g_2g_1^2-3g_1^4}{(g_2-g_1^2)^{2 where
g_k=Gamma(1-kxi) , k=1,2,3,4, andGamma(t) is Gamma function.Link to Fréchet, Weibull and Gumbel families
The shape parameter
xi governs the tail behaviour of the distribution. The sub-families defined byxi o 0 ,xi>0 andxi<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;; xinmathbb R
* Fréchet or type II extreme value distribution:F(x;mu,sigma,alpha)=egin{cases} 0 & xleq mu \ e^{-((x-mu)/sigma)^{-alpha & x>mu end{cases}
* Reversed Weibull or type III extreme value distribution:F(x;mu,sigma,alpha)=egin{cases} e^{-(-(x-mu)/sigma)^{alpha & x where sigma>0 andalpha>0 .Remark I: The theory here relates to maxima and the distribution being discussed is an extreme value distribution for maxima. A Generalised Extreme Value distribution for minima can be obtained, for example by substituting (-x) for x in the distribution function and this yields a separate family of distributions.
Remark II: The ordinary Weibull distribution arises in reliability applications and is obtained from the distribution here by using the variable
t = mu - x , which gives a strictly positive support - in contrast to the use in the extreme value theory here. This arises because the Weibull distribution is used in cases that deal with the minimum rather than the maximum. The distribution here has an addition parameter compared to the usual form of the Weibull distribution and, in addition, is reversed so that the distribution has an upper bound rather than a lower bound. Importantly, in applications of the GEV, the upper bound is unknown and so must be estimated while when applying the Weibull distribution the lower bound is known to be zero.Remark III: Note the differences in the ranges of interest for the three extreme value distributions: Gumbel is unlimited, Fréchet has a lower limit, while the reversed 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,sigma,alpha) , then the cumulative distribution function ofln X is of type I, namelyF(x;ln sigma, 1/alpha) . Similarly, if the cumulative distribution function ofX is of type III:F(x;0,sigma,alpha) , the cumulative distribution function ofln X is of type I:F(x;-ln sigma, 1/alpha) .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_2ldots be a sequence of independent and identically distributed random variables, letM_n=max{X_1,ldots,X_n} . If two sequences of real numbersa_n, b_n exist such thata_n>0 and:lim_{n o infty}Pleft(frac{M_n-b_n}{a_n}leq x ight) = F(x) then ifF 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 sequenceX_n .For the practical application this theorem means: For samples taken from a well behaving, arbitrary distribution
X the resulting extreme value distributionM_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 sequenceX_n of random variable is thenormal distribution no matter what the distribution of theX_n , The extremal types theorem is similar in scope where maxima is substituted for average and GEV distribution is substituted for normal distribution.References
* Embrechts, P., C. Klüppelberg, and T. Mikosch (1997) Modelling extremal events for insurance and finance. Berlin: Spring Verlag
*
*
*
Wikimedia Foundation. 2010.