Open main menu
Home
Random
Recent changes
Special pages
Community portal
Preferences
About Wikipedia
Disclaimers
Incubator escapee wiki
Search
User menu
Talk
Dark mode
Contributions
Create account
Log in
Editing
Generalized extreme value distribution
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
{{Short description|Family of probability distributions}} {{More citations needed|date=May 2020}} {{Probability distribution | type = density | notation = <math>\textrm{GEV}(\mu, \sigma, \xi)</math> | parameters = <math>\mu \in \mathbb{R}</math> ([[location parameter|location]]) <br/> <math>\sigma > 0</math> ([[scale parameter|scale]]) <br/> <math>\xi \in \mathbb{R}</math> ([[shape parameter|shape]]) | support = <math>x \in \big[ \mu - \tfrac{\sigma}{\xi}, +\infty\big)</math> when <math>\xi > 0</math> <br/> <math>x \in (-\infty, \infty)</math> when <math>\xi = 0</math> <br/> <math>x \in \big(-\infty, \mu - \tfrac{\sigma}{\xi}\big]</math> when <math>\xi < 0</math> | pdf = <math>\tfrac{1}{\sigma} \, t(x)^{\xi + 1}\,\mathrm{e}^{-t(x)}</math> <br/> where <math>t(x) = \begin{cases} \left[ 1 + \xi \left( \tfrac{x - \mu}{\sigma} \right) \right]^{-1/\xi} & \text{if } \xi \neq 0 \\ \exp\left(-\tfrac{x - \mu}{\sigma} \right) & \text{if } \xi = 0 \end{cases}</math> | cdf = <math>\mathrm{e}^{-t(x)}</math> for <math>x</math> in the support (''see above'') | mean = <math>\begin{cases} \mu + \dfrac{\sigma(g_1 - 1)}{\xi} & \text{if } \xi \neq 0, \xi < 1 \\ \mu + \sigma\gamma & \text{if } \xi = 0 \\ \infty & \text{if } \xi \geq 1 \end{cases}</math> <br/> where <math>g_k = \Gamma(1 - k\xi)</math> (''see'' [[Gamma function]]) <br/> and <math>\gamma</math> is [[Eulerâs constant]] | median = <math>\begin{cases}\mu + \sigma \, \dfrac{(\ln2)^{-\xi} - 1}{\xi} & \text{if } \xi \neq 0 \\ \mu - \sigma \ln\ln2 & \text{if } \xi = 0 \end{cases}</math> | mode = <math>\begin{cases}\mu + \sigma\, \dfrac{( 1 + \xi)^{-\xi} - 1}{\xi} & \text{if } \xi \neq 0 \\ \mu & \text{if } \xi = 0 \end{cases}</math> | variance = <math>\begin{cases} \sigma^2 \, \dfrac{g_2 - g_1^2}{\xi^2} & \text{if } \xi \neq 0 \text{ and } \xi < \tfrac{1}{2} \\ \sigma^2 \, \frac{\pi^2}{6} & \text{if } \xi = 0 \\ \infty & \text{if } \xi \geq \tfrac{1}{2} \end{cases}</math> | skewness = <math>\begin{cases} \sgn(\xi) \, \dfrac{g_3 - 3 g_2 g_1 + 2 g_1^3}{(g_2 - g_1^2)^{3/2}} & \text{if } \xi \neq 0 \text{ and } \xi < \tfrac{1}{3} \\ \dfrac{12 \sqrt{6} \, \zeta(3)}{\pi^3} & \text{if } \xi = 0 \end{cases}</math> <br> where <math>\sgn(x)</math> is the [[sign function]] <br/> and <math>\zeta(x)</math> is the [[Riemann zeta function]] | kurtosis = <math>\begin{cases} \dfrac{g_4 - 4 g_3 g_1 + 6 g_1^2 g_2 - 3 g_1^4}{(g_2 - g_1^2)^2} & \text{if } \xi \neq 0 \text{ and } \xi < \tfrac{1}{4} \\ \tfrac{12}{5} & \text{if } \xi = 0\end{cases}</math> | entropy = <math>\ln(\sigma) + \gamma\xi + \gamma + 1</math> | mgf = ''see'' {{harvp|Muraleedharan|Guedes Soares|Lucas|2011}}<ref name=Muraldhrn-Soares-Lucas-2011/> | char = ''see'' {{harvp|Muraleedharan|Guedes Soares|Lucas|2011}}<ref name=Muraldhrn-Soares-Lucas-2011/> }} In [[probability theory]] and [[statistics]], the '''generalized extreme value''' ('''GEV''') '''distribution'''<ref> {{cite web |last=Weisstein |first=Eric W. |title=Extreme value distribution |website=mathworld.wolfram.com |url=https://mathworld.wolfram.com/ExtremeValueDistribution.html |access-date=2021-08-06 |lang=en }} </ref> is a family of continuous [[probability distribution]]s developed within [[extreme value theory]] to combine the [[Gumbel distribution|Gumbel]], [[FrĂ©chet distribution|FrĂ©chet]] and [[Weibull distribution|Weibull]] families also known as type I, II and III extreme value distributions. By the [[FisherâTippettâGnedenko theorem|extreme value theorem]] the GEV distribution is the only possible limit distribution of properly normalized maxima of a sequence of independent and identically distributed random variables.<ref name=haan/> Note that a limit distribution needs to exist, which requires regularity conditions on the tail of the distribution. Despite this, the GEV distribution is often used as an approximation to model the maxima of long (finite) sequences of random variables. In some fields of application the generalized extreme value distribution is known as the '''FisherâTippett distribution''', named after [[Ronald Fisher|R.A. Fisher]] and [[L. H. C. Tippett|L.H.C. Tippett]] who recognised three different forms outlined below. However usage of this name is sometimes restricted to mean the special case of the [[Gumbel distribution]]. The origin of the common functional form for all three distributions dates back to at least {{harvp|Jenkinson|1955}},<ref> {{cite journal |last=Jenkinson |first=Arthur F. |title=The frequency distribution of the annual maximum (or minimum) values of meteorological elements |journal=[[Quarterly Journal of the Royal Meteorological Society]] |pages=158â171 |year=1955 |volume=81 |issue=348 |doi=10.1002/qj.49708134804 |bibcode=1955QJRMS..81..158J }} </ref> though allegedly<ref name=haan> {{cite book |last1=Haan |first1=Laurens |last2=Ferreira |first2=Ana |year=2007 |title=Extreme Value Theory: An introduction |publisher=Springer }} </ref> it could also have been given by {{harvp|von Mises|1936}}.<ref name=mises> {{cite journal |last=von Mises |first=R. |year=1936 |title=La distribution de la plus grande de {{mvar|n}} valeurs |journal=Rev. Math. Union Interbalcanique |volume=1 |pages=141â160 }} </ref> ==Specification== Using the standardized variable <math>s = \tfrac{x - \mu}{\sigma}</math>, where <math>\mu</math>, the location parameter, can be any real number, and <math>\sigma > 0</math> is the scale parameter; the cumulative distribution function of the GEV distribution is then :<math display="block"> F(s; \xi) = \begin{cases} \exp (-\mathrm{e}^{-s}) & \text{for } \xi = 0 , \\ \exp \bigl( \! - ( 1 + \xi s)^{-1/\xi} \bigr) & \text{for } \xi \neq 0 \text{ and } \xi s > -1 , \\ 0 & \text{for } \xi > 0 \text{ and } s \le -\tfrac{1}{\xi} , \\ 1 & \text{for } \xi < 0 \text{ and } s \ge \tfrac{1}{|\xi|} , \end{cases}</math> where <math>\xi</math>, the shape parameter, can be any real number. Thus, for <math>\xi > 0</math>, the expression is valid for <math>s > -\tfrac{1}{\xi}</math>, while for <math>\xi < 0</math> it is valid for <math>s < - \tfrac{1}{\xi}</math>. In the first case, <math>-\tfrac{1}{\xi}</math> is the negative, lower end-point, where <math>F</math> is {{math|0}}; in the second case, <math>-\tfrac{1}{\xi}</math> is the positive, upper end-point, where <math>F</math> is 1. For <math>\xi = 0</math>, the second expression is formally undefined and is replaced with the first expression, which is the result of taking the limit of the second, as <math>\xi \to 0</math> in which case <math>s</math> can be any real number. In the special case of <math>x = \mu</math>, we have <math>s = 0</math>, so <math>F(0; \xi) = \mathrm{e}^{-1} \approx 0.368</math> regardless of the values of <math>\xi</math> and <math>\sigma</math>. The probability density function of the standardized distribution is :<math display="block">f(s; \xi) = \begin{cases} \mathrm{e}^{-s} \exp (-\mathrm{e}^{-s} ) & \text{for } \xi = 0 , \\ (1 + \xi s)^{-( 1 + 1/\xi )} \exp\bigl(\! -( 1 + \xi s )^{-1/\xi} \bigr) & \text{for } \xi \neq 0 \text{ and } \xi s > -1 , \\ 0 & \text{otherwise;} \end{cases}</math> again valid for <math>s > -\tfrac{1}{\xi}</math> in the case <math>\xi > 0</math>, and for <math>s < -\tfrac{1}{\xi}</math> in the case <math>\xi < 0</math>. The density is zero outside of the relevant range. In the case <math>\xi = 0</math>, the density is positive on the whole real line. Since the cumulative distribution function is invertible, the quantile function for the GEV distribution has an explicit expression, namely :<math display="block">Q(p; \mu, \sigma, \xi) = \begin{cases} \mu - \sigma \ln ( -\ln p ) & \text{for } \xi = 0 \text{ and } p \in (0, 1) , \\ \mu + \dfrac{\sigma}{\xi} \big( ( -\ln p)^{-\xi} - 1 \big) & \text{for } \xi > 0 \text{ and } p \in [0, 1) , \text{ or } \xi < 0 \text{ and } p \in (0, 1] ; \end{cases}</math> and therefore the quantile density function <math>q = \tfrac{\mathrm{d}Q}{\mathrm{d}p}</math> is :<math>q(p; \sigma, \xi) = \frac{\sigma}{(-\ln p)^{\xi + 1} \,p} \qquad \text{for } p \in (0, 1) ,</math> valid for <math>\sigma > 0</math> and for any real <math>\xi</math>. [[File:GevDensity 2.svg|Example of probability density functions for distributions of the GEV family.]] <ref name=norton/> ==Summary statistics== Using <math>\ g_k \equiv \Gamma(1 - k\ \xi) ~</math> for <math>~ k \in\{\ 1, 2, 3, 4\ \}\ ,</math> where <math>\ \Gamma(\cdot)\ </math> is the [[gamma function]], some simple statistics of the distribution are given by:{{citation needed|date=May 2011}} :<math>\ \operatorname{\mathbb E}\left( X \right) = \mu + \bigl(\ g_1 - 1\ \bigr)\frac{\ \sigma\ }{ \xi } \quad </math> for <math>\quad \xi < 1\ ,</math> :<math>\ \operatorname{Var}\left( X \right) = \bigl(\ g_2 - g_1^2\ \bigr) \frac{~ \sigma^2\ }{~ \xi^2\ }\ ,</math> :<math>\ \operatorname{Mode}\left( X \right) = \mu + \Bigl(\ \left( 1 + \xi \right)^{-\xi} - 1\ \Bigr) \frac{\ \sigma\ }{ \xi } ~.</math> The [[skewness]] is :<math> \ \operatorname{skewness}\left( X \right) = \begin{cases} \frac{\ \bigl(\ g_3\ -\ 3\ g_2\ g_1\ +\ 2\ g_1^3\ \bigr)\ }{ \bigl(\ g_2 - g_1^2\ \bigr)^{3/2} } \cdot \sgn(\xi) \quad ~~ \mathsf{ for }\quad \xi \ne 0\ , \\ \\ \quad \frac{\ 12\sqrt{6\ }\ \zeta(3)\ }{~ \pi^3 } \quad \approx \quad 1.14 \quad\qquad \mathsf{ for }\quad \xi = 0 ~. \end{cases} </math> The excess [[kurtosis]] is: :<math>\ \operatorname{kurtosis\ excess}\left( X \right) = \frac{\ \left(\ g_4 - 4\ g_3\ g_1 + 6\ g_2\ g_1^2 - 3\ g_1^4\ \right)\ }{ \left(\ g_2 - g_1^2\ \right)^2 } - 3 ~.</math> <!-- Unsourced image removed: [[Image:gevDensity.jpeg||Example of density functions for distributions of the GEV family.]] --> == Link to FrĂ©chet, Weibull, and Gumbel families == The shape parameter <math>\ \xi\ </math> governs the tail behavior of the distribution. The sub-[[locationâscale family|families]] defined by three cases: <math>\ \xi = 0\ ,</math> <math>\ \xi > 0\ ,</math> and <math>\ \xi < 0\ ;</math> these correspond, respectively, to the [[Gumbel distribution|''Gumbel'']], [[FrĂ©chet distribution|''FrĂ©chet'']], and [[Weibull distribution|''Weibull'']] [[locationâscale family|families]], whose cumulative distribution functions are displayed below. * Type I or ''[[Gumbel distribution|'''Gumbel''']]'' extreme value distribution, case <math>~ \xi = 0\ , \quad</math> for all <math> \quad x \in \Bigl(\ -\infty\ ,\ +\infty\ \Bigr)\ :</math> :<math> F(\ x;\ \mu,\ \sigma,\ 0\ ) = \exp \left( - \exp \left( -\frac{\ x - \mu\ }{\sigma} \right) \right) ~. </math> * Type II or ''[[FrĂ©chet distribution|'''FrĂ©chet''']]'' extreme value distribution, case <math>~ \xi > 0\ , \quad </math> for all <math> \quad x \in \left(\ \mu - \tfrac{\sigma}{\ \xi\ }\ ,\ +\infty\ \right)\ :</math> :Let <math>\quad \alpha \equiv \tfrac{\ 1\ }{ \xi } > 0 \quad </math> and <math> \quad y \equiv 1 + \tfrac{\xi}{\sigma} (x-\mu)\ ;</math> :<math> F(\ x;\ \mu,\ \sigma,\ \xi\ ) = \begin{cases} 0 & y \leq 0 \quad \mathsf{~ or\ equiv. ~} \quad x \leq \mu - \tfrac{\sigma}{\ \xi\ } \\ \exp\left( -\frac{1}{~ y^\alpha\ } \right) & y > 0 \quad \mathsf{~ or\ equiv. ~} \quad x > \mu - \tfrac{\sigma}{\ \xi\ } ~. \end{cases}</math> * Type III or ''reversed [[Weibull distribution|'''Weibull''']]'' extreme value distribution, case <math>~ \xi < 0\ , \quad </math> for all <math> \quad x \in \left( -\infty\ ,\ \mu + \tfrac{ \sigma }{\ |\ \xi\ |\ }\ \right)\ :</math> :Let <math> \quad \alpha \equiv - \tfrac{1}{\ \xi\ } > 0 \quad </math> and <math> \quad y \equiv 1 - \tfrac{\ |\ \xi\ |\ }{\sigma} (x - \mu)\ ;</math> :<math> F(\ x;\ \mu,\ \sigma,\ \xi\ ) = \begin{cases} \exp\left( -y^{\alpha} \right) & y > 0 \quad \mathsf{~ or\ equiv. ~} \quad x < \mu + \tfrac{ \sigma }{\ |\ \xi\ |\ } \\ 1 & y \leq 0 \quad \mathsf{~ or\ equiv. ~} \quad x \geq \mu + \tfrac{ \sigma }{\ |\ \xi\ |\ } ~. \end{cases}</math> The subsections below remark on properties of these distributions. === Modification for minima rather than maxima === The theory here relates to data maxima and the distribution being discussed is an extreme value distribution for maxima. A generalised extreme value distribution for data minima can be obtained, for example by substituting <math>\ -x\;</math> for <math>\;x\;</math> in the distribution function, and subtracting the cumulative distribution from one: That is, replace <math>\ F(x)\ </math> with {{nobr|<math>\ 1 - F(-x)\ </math> .}} Doing so yields yet another family of distributions. === Alternative convention for the Weibull distribution === The ordinary Weibull distribution arises in reliability applications and is obtained from the distribution here by using the variable <math>\ t = \mu - x\ ,</math> which gives a strictly positive support, in contrast to the use in the formulation of extreme value theory here. This arises because the ordinary Weibull distribution is used for cases that deal with data ''minima'' rather than data maxima. 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, whereas when applying the ordinary Weibull distribution in reliability applications the lower bound is usually known to be zero. === Ranges of the distributions === Note the differences in the ranges of interest for the three extreme value distributions: [[Gumbel distribution|Gumbel]] is unlimited, [[FrĂ©chet distribution|FrĂ©chet]] has a lower limit, while the reversed [[Weibull distribution|Weibull]] has an upper limit. More precisely, [[Extreme value theory#Univariate theory|univariate extreme value theory]] describes which of the three is the limiting law according to the initial law {{mvar| X }} and in particular depending on the original distribution's tail. === Distribution of log variables === One can link the type I to types II and III in the following way: If the cumulative distribution function of some random variable <math>\ X\ </math> is of type II, and with the positive numbers as support, i.e. <math>\ F(\ x;\ 0,\ \sigma,\ \alpha\ )\ ,</math> then the cumulative distribution function of <math>\ln X</math> is of type I, namely <math>\ F(\ x;\ \ln \sigma,\ \tfrac{1}{\ \alpha\ },\ 0\ ) ~.</math> Similarly, if the cumulative distribution function of <math>\ X\ </math> is of type III, and with the negative numbers as support, i.e. <math>\ F(\ x;\ 0,\ \sigma,\ -\alpha\ )\ ,</math> then the cumulative distribution function of <math>\ \ln (-X)\ </math> is of type I, namely <math>\ F(\ x;\ -\ln \sigma,\ \tfrac{\ 1\ }{\alpha},\ 0\ ) ~.</math> ==Link to logit models (logistic regression)== [[Multinomial logit]] models, and certain other types of [[logistic regression]], can be phrased as [[latent variable]] models with [[error variable]]s distributed as [[Gumbel distribution]]s (type I generalized extreme value distributions). This phrasing is common in the theory of [[discrete choice]] models, which include [[logit model]]s, [[probit model]]s, and various extensions of them, and derives from the fact that the difference of two type-I GEV-distributed variables follows a [[logistic distribution]], of which the [[logit function]] is the [[quantile function]]. The type-I GEV distribution thus plays the same role in these logit models as the [[normal distribution]] does in the corresponding probit models. ==Properties== The [[cumulative distribution function]] of the generalized extreme value distribution solves the [[stability postulate]] equation.{{Citation needed|date=May 2011}} The generalized extreme value distribution is a special case of a max-stable distribution, and is a transformation of a min-stable distribution. ==Applications== * The GEV distribution is widely used in the treatment of "tail risks" in fields ranging from insurance to finance. In the latter case, it has been considered as a means of assessing various financial risks via metrics such as [[value at risk]].<ref> {{cite report |last=Moscadelli |first=Marco |date=30 July 2004 |title=The modelling of operational risk: Experience with the analysis of the data collected by the Basel Committee |ssrn=557214 |doi=10.2139/ssrn.557214 |type=non-peer reviewed article |url=http://www.unalmed.edu.co/~ndgirald/Archivos%20Lectura/Archivos%20curso%20Riesgo%20Operativo/moscadelli%202004.pdf |via=Archivos curso Riesgo Operativo de N.D. Girald (unalmed.edu.co) }} </ref><ref> {{cite journal |last1 = GuĂ©gan |first1 = D. |last2 = Hassani |first2 = B.K. |year = 2014 |title = A mathematical resurgence of risk management: An extreme modeling of expert opinions |journal = Frontiers in Finance and Economics |volume = 11 |issue = 1 |pages = 25â45 |ssrn = 2558747 }} </ref> [[File:GEV Surinam.png|thumb|300px|Fitted GEV probability distribution to monthly maximum one-day rainfalls in October, Surinam<ref name=CumFreq-ref> {{cite web |title=CumFreq for probability distribution fitting |website=waterlog.info |url=https://www.waterlog.info/cumfreq.htm }} See also [[CumFreq]]. </ref>]] * However, the resulting shape parameters have been found to lie in the range leading to undefined means and variances, which underlines the fact that reliable data analysis is often impossible.<ref> {{cite web |first=Kjersti |last=Aas |date=23 January 2008 |title={{grey|[no title cited]}} |type=lecture |publisher=[[Norwegian University of Science and Technology|Norges teknisk-naturvitenskapelige universitet]] |place=Trondheim, NO |website=citeseerx.ist.psu.edu |citeseerx=10.1.1.523.6456 |url=http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.523.6456&rep=rep1&type=pdf |format=PDF |url-status=dead }} </ref>{{full citation|date=October 2024}} * In [[hydrology]] the GEV distribution is applied to extreme events such as annual maximum one-day rainfalls and river discharges.<ref> {{cite journal |last1=Liu |first1=Xin |last2=Wang |first2=Yu |date=September 2022 |title=Quantifying annual occurrence probability of rainfall-induced landslide at a specific slope |journal=Computers and Geotechnics |volume=149 |page=104877 |doi=10.1016/j.compgeo.2022.104877 |doi-access=free |bibcode=2022CGeot.14904877L |s2cid=250232752 |lang=en }} </ref> The blue picture, made with [[CumFreq]],<ref name=CumFreq-ref/> illustrates an example of fitting the GEV distribution to ranked annually maximum one-day rainfalls showing also the 90% [[confidence belt]] based on the [[binomial distribution]]. The rainfall data are represented by [[plotting position]]s as part of the [[cumulative frequency analysis]]. ===Example for Normally distributed variables=== Let <math>\ \left\{\ X_i\ \big|\ 1 \le i \le n\ \right\}\ </math> be [[i.i.d.]] [[Normal Distribution|normally distributed]] random variables with mean {{math|0}} and variance {{math|1}}. The [[FisherâTippettâGnedenko theorem]]<ref> {{cite book |last1=David |first1=Herbert A. |title=Order Statistics |last2=Nagaraja |first2=Haikady N. |publisher=John Wiley & Sons |year=2004 |page=299 }} </ref> tells us that <math>\ \max \{\ X_i\ \big|\ 1 \le i \le n\ \} \sim GEV(\mu_n, \sigma_n, 0)\ ,</math> where <math> \begin{align} \mu_n &= \Phi^{-1}\left( 1 - \frac{\ 1\ }{ n } \right) \\ \sigma_n &= \Phi^{-1}\left( 1 - \frac{ 1 }{\ n\ \mathrm{e}\ } \right)- \Phi^{-1}\left(1-\frac{\ 1\ }{ n } \right) ~. \end{align} </math> This allow us to estimate e.g. the mean of <math>\ \max \{\ X_i\ \big|\ 1 \le i \le n\ \}\ </math> from the mean of the GEV distribution: <math> \begin{align} \operatorname{\mathbb E}\left\{\ \max\left\{\ X_i\ \big|\ 1 \le i \le n\ \right\}\ \right\} & \approx \mu_n + \gamma_{\mathsf E}\ \sigma_n \\ &= (1 - \gamma_{\mathsf E})\ \Phi^{-1}\left( 1 - \frac{\ 1\ }{ n } \right) + \gamma_{\mathsf E}\ \Phi^{-1}\left( 1 - \frac{1}{\ e\ n\ } \right) \\ &= \sqrt{\log \left(\frac{ n^2 }{\ 2 \pi\ \log \left(\frac{n^2}{2\pi} \right)\ }\right) ~}\ \cdot\ \left(1 + \frac{ \gamma }{\ \log n\ } + \mathcal{o} \left(\frac{ 1 }{\ \log n\ } \right) \right)\ , \end{align} </math> where <math>\ \gamma_{\mathsf E}\ </math> is the [[EulerâMascheroni constant]]. ==Related distributions== # If <math>\ X \sim \textrm{GEV}(\mu,\,\sigma,\,\xi)\ </math> then <math>\ m X + b \sim \textrm{GEV}(m \mu+b,\ |m|\sigma,\ \xi)\ </math> # If <math>\ X \sim \textrm{Gumbel}(\mu,\ \sigma)\ </math> ([[Gumbel distribution]]) then <math>\ X \sim \textrm{GEV}(\mu,\,\sigma,\,0)\ </math> # If <math>\ X \sim \textrm{Weibull}(\sigma,\,\mu)\ </math> ([[Weibull distribution]]) then <math>\ \mu\left(1-\sigma\log \tfrac{X}{\sigma} \right) \sim \textrm{GEV}(\mu,\,\sigma,\,0)\ </math> # If <math>\ X \sim \textrm{GEV}(\mu,\,\sigma,\,0)\ </math> then <math>\ \sigma \exp (-\tfrac{X-\mu}{\mu \sigma} ) \sim \textrm{Weibull}(\sigma,\,\mu)\ </math> ([[Weibull distribution]]) # If <math>\ X \sim \textrm{Exponential}(1)\ </math> ([[Exponential distribution]]) then <math>\ \mu - \sigma \log X \sim \textrm{GEV}(\mu,\,\sigma,\,0)\ </math> # If <math>\ X \sim \mathrm{Gumbel}(\alpha_X, \beta)\ </math> and <math>\ Y \sim \mathrm{Gumbel}(\alpha_Y, \beta)\ </math> then <math>\ X-Y \sim \mathrm{Logistic}(\alpha_X-\alpha_Y,\beta)\ </math> (see [[Logistic distribution]]). # If <math>\ X\ </math> and <math>\ Y \sim \mathrm{Gumbel}(\alpha, \beta)\ </math> then <math>\ X+Y \nsim \mathrm{Logistic}(2 \alpha,\beta)\ </math> (The sum is ''not'' a logistic distribution). :: Note that <math>\ \operatorname{\mathbb E}\{\ X + Y\ \} = 2\alpha+2\beta\gamma \neq 2\alpha = \operatorname{\mathbb E}\left\{\ \operatorname{Logistic}(2 \alpha,\beta)\ \right\} ~.</math> ===Proofs=== 4. Let <math>\ X \sim \textrm{ Weibull }(\sigma,\,\mu)\ ,</math> then the cumulative distribution of <math>\ g(x) = \mu\left(1-\sigma\log\frac{X}{\sigma} \right)\ </math> is: :<math> \begin{align} \operatorname{\mathbb P}\left\{\ \mu \left(1-\sigma\log\frac{\ X\ }{ \sigma } \right) < x\ \right\} &= \operatorname{\mathbb P}\left\{\ \log\frac{X}{\sigma} > \frac{1 - x/\mu}{\sigma}\ \right\} \\ {} \\ & \mathsf{\ Since\ the\ logarithm\ is\ always\ increasing:\ } \\ {} \\ &= \operatorname{\mathbb P}\left\{\ X > \sigma \exp\left[ \frac{1 - x/\mu}{\sigma} \right]\ \right\} \\ &= \exp\left( - \left(\cancel{\sigma} \exp\left[ \frac{1 - x/\mu}{\sigma} \right] \cdot \cancel{\frac{1}{\sigma}} \right)^\mu \right) \\ &= \exp\left( - \left( \exp\left[ \frac{\cancelto{\mu}{1} - x/\cancel{\mu}}{\sigma} \right] \right)^\cancel{\mu} \right) \\ &= \exp\left( - \exp\left[ \frac{\mu - x}{\sigma} \right] \right) \\ &= \exp\left( - \exp\left[ - s \right] \right), \quad s = \frac{x - \mu}{\sigma}\ , \end{align} </math> :which is the cdf for <math>\sim \textrm{GEV}(\mu,\,\sigma,\,0) ~.</math> 5. Let <math>\ X \sim \textrm{Exponential}(1)\ ,</math> then the cumulative distribution of <math>\ g(X) = \mu - \sigma \log X\ </math> is: :<math> \begin{align} \operatorname{\mathbb P}\left\{\ \mu - \sigma \log X < x\ \right\} &= \operatorname{\mathbb P}\left\{\ \log X > \frac{\mu - x}{\sigma}\ \right\} \\ {} \\ & \mathsf{\ Since\ the\ logarithm\ is\ always\ increasing:\ } \\ {} \\ &= \operatorname{\mathbb P}\left\{\ X > \exp\left( \frac{\ \mu - x\ }{ \sigma } \right)\ \right\} \\ &= \exp\left[- \exp\left( \frac{\ \mu - x\ }{ \sigma } \right) \right] \\ &= \exp\left[- \exp(-s) \right]\ , \quad ~\mathsf{ where }~ \quad s \equiv \frac{x - \mu}{\sigma}\ ; \end{align} </math> :which is the cumulative distribution of <math>\ \operatorname{GEV}(\mu, \sigma, 0) ~.</math> ==See also== * [[Extreme value theory#Univariate theory|Extreme value theory (univariate theory)]] * [[FisherâTippettâGnedenko theorem]] * [[Generalized Pareto distribution]] * [[German tank problem]], opposite question of population maximum given sample maximum * [[PickandsâBalkemaâDe Haan theorem]] ==References== {{reflist|refs= <ref name=Muraldhrn-Soares-Lucas-2011> {{cite book |last1=Muraleedharan |first1=G. |last2=Guedes Soares |first2=C. |last3=Lucas |first3=ClĂĄudia |year=2011 |section=Characteristic and moment generating functions of generalised extreme value distribution (GEV) |editor-first=Linda L. |editor-last=Wright |title=Sea Level Rise, Coastal Engineering, Shorelines, and Tides |at=Chapter 14, pp. 269â276 |publisher=Nova Science Publishers |isbn=978-1-61728-655-1 }} </ref> <ref name=norton>{{cite journal |last1 = Norton |first1 = Matthew |last2 = Khokhlov |first2 = Valentyn |last3 = Uryasev |first3 = Stan |year = 2019 |title = Calculating CVaR and bPOE for common probability distributions with application to portfolio optimization and density estimation |journal = Annals of Operations Research |volume = 299 |issue = 1â2 |pages = 1281â1315 |publisher = Springer |doi = 10.1007/s10479-019-03373-1 |arxiv = 1811.11301 |s2cid = 254231768 |url = http://uryasev.ams.stonybrook.edu/wp-content/uploads/2019/10/Norton2019_CVaR_bPOE.pdf |access-date = 2023-02-27 |archive-date = 2023-03-31 |archive-url = https://web.archive.org/web/20230331230821/http://uryasev.ams.stonybrook.edu/wp-content/uploads/2019/10/Norton2019_CVaR_bPOE.pdf |url-status = dead }}</ref> }} ==Further reading== {{refbegin|25em|small=yes}} * {{cite book |last1=Embrechts |first1=Paul |last2=KlĂŒppelberg |first2=Claudia |author2-link= Claudia KlĂŒppelberg |last3=Mikosch |first3=Thomas |title=Modelling Extremal Events for Insurance and Finance |year=1997 |location=Berlin, DE |publisher=Springer Verlag |isbn=9783540609315 |url=https://books.google.com/books?id=BXOI2pICfJUC |via=Google books }} * {{cite book | last1=Leadbetter | first1=M.R. | last2=Lindgren | first2=G. | last3=RootzĂ©n | first3=H. | year=1983 | title=Extremes and Related Properties of Random Sequences and Processes | publisher=Springer-Verlag | isbn = 0-387-90731-9 }} * {{cite book | last=Resnick | first=S.I. | year=1987 | title=Extreme Values, Regular Variation, and Point Processes | publisher=Springer-Verlag | isbn = 0-387-96481-9 }} * {{cite book | last=Coles | first=Stuart | year=2001 | title=An Introduction to Statistical Modeling of Extreme Values | url = https://books.google.com/books?id=2nugUEaKqFEC&pg=PP1 | publisher=Springer-Verlag | isbn = 1-85233-459-2 }} {{refend}} {{ProbDistributions|continuous-variable}} [[Category:Continuous distributions]] [[Category:Extreme value data]] [[Category:Location-scale family probability distributions]] [[Category:Stability (probability)]]
Edit summary
(Briefly describe your changes)
By publishing changes, you agree to the
Terms of Use
, and you irrevocably agree to release your contribution under the
CC BY-SA 4.0 License
and the
GFDL
. You agree that a hyperlink or URL is sufficient attribution under the Creative Commons license.
Cancel
Editing help
(opens in new window)
Pages transcluded onto the current version of this page
(
help
)
:
Template:Citation needed
(
edit
)
Template:Cite book
(
edit
)
Template:Cite journal
(
edit
)
Template:Cite report
(
edit
)
Template:Cite web
(
edit
)
Template:Full citation
(
edit
)
Template:Harvp
(
edit
)
Template:Math
(
edit
)
Template:More citations needed
(
edit
)
Template:Mvar
(
edit
)
Template:Nobr
(
edit
)
Template:ProbDistributions
(
edit
)
Template:Probability distribution
(
edit
)
Template:Refbegin
(
edit
)
Template:Refend
(
edit
)
Template:Reflist
(
edit
)
Template:Short description
(
edit
)