Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Generalized Extreme-Value Distribution

Statement

The generalized extreme-value distribution is the standard limit family for normalized block maxima. With location μ\mu, scale σ>0\sigma>0, shape ξ\xi, and standardized variable z=(xμ)/σz=(x-\mu)/\sigma, its CDF is

Hξ,μ,σ(x)=exp[(1+ξz)1/ξ],ξ0,H_{\xi,\mu,\sigma}(x) = \exp\left[-\left(1+\xi z\right)^{-1/\xi}\right], \qquad \xi\ne0,

where 1+ξz>01+\xi z>0. For ξ>0\xi>0, extend the CDF by zero at and below μσ/ξ\mu-\sigma/\xi; for ξ<0\xi<0, extend it by one at and above μσ/ξ\mu-\sigma/\xi. The boundary case is the Gumbel law,

H0,μ,σ(x)=exp[ez].H_{0,\mu,\sigma}(x) = \exp\left[-e^{-z}\right].

The shape parameter ξ\xi is the same extreme-value index used in threshold excesses: ξ>0\xi>0 is Frechet-type heavy-tail behavior, ξ=0\xi=0 is Gumbel-type behavior, and ξ<0\xi<0 is Weibull-type finite endpoint. Here XiX_i are observations and Mn=max1inXiM_n=\max_{1\le i\le n}X_i. The Pareto example below uses lower cutoff xm>0x_m>0 and exponent α>0\alpha>0. Gumbel attraction is broader than exponential parent tails: normal and lognormal parents are examples, despite their different tail decay.

Maxima intuition

The distribution is to block maxima what the Generalized Pareto Distribution is to threshold exceedances. For iid observations, if there are constants an>0a_n>0 and bnb_n such that (Mnbn)/an(M_n-b_n)/a_n converges to a nondegenerate law, that law is a GEV distribution up to location and scale. Existence of such a limit is an assumption; it does not hold for every parent law.

For fat-tail work, ξ>0\xi>0 is the main case. A Pareto tail with exponent α\alpha has ξ=1/α\xi=1/\alpha, so the fitted block-maximum shape should point to the same tail coordinate as a Hill or GPD threshold analysis when the modeling assumptions are reasonable.

Shape regimes

With μ=0\mu=0 and σ=1\sigma=1, the support makes the three shapes explicit.

Shape ξ\xiInterior CDFSupport
-0.3exp[(10.3x)10/3]\exp[-(1-0.3x)^{10/3}]x<10/3x<10/3; CDF is one above
0exp[ex]\exp[-e^{-x}]All real xx
0.3exp[(1+0.3x)10/3]\exp[-(1+0.3x)^{-10/3}]x>10/3x>-10/3; CDF is zero below

Every curve passes through H(0)=e1H(0)=e^{-1}. Positive shape gives a power-law right tail; negative shape gives a finite right endpoint.

Pareto block maxima example

If X1,,XnX_1,\dots,X_n are iid Pareto Type I with lower cutoff xmx_m and exponent α\alpha, then for xxmx\ge x_m,

P(Mnx)=(1(xmx)α)n.\mathbb P(M_n\le x) = \left(1-\left(\frac{x_m}{x}\right)^\alpha\right)^n.

With an=xmn1/αa_n=x_m n^{1/\alpha}, for fixed y>0y>0 and sufficiently large nn,

P(Mn/any)=(11nyα)nexp(yα),y>0.\mathbb P(M_n/a_n\le y) = \left(1-\frac{1}{ny^\alpha}\right)^n \to \exp(-y^{-\alpha}),\qquad y>0.

This is the Frechet member of the GEV family with ξ=1/α\xi=1/\alpha after a standard location-scale reparameterization.

Block-size calibration

For an exact Pareto parent with α=1.8=9/5\alpha=1.8=9/5, the limiting shape is ξ=5/9\xi=5/9. At block size n=250n=250, the normalizing scale is an=xm2505/9a_n=x_m250^{5/9}. At the normalized threshold y=1y=1,

P(M250/a2501)=(249/250)2500.367142,Φ1.8(1)=e10.367879.\mathbb P(M_{250}/a_{250}\le1) =(249/250)^{250}\approx0.367142, \qquad \Phi_{1.8}(1)=e^{-1}\approx0.367879.

This compares an exact finite-block CDF with its limit; it is not a GEV fit or an assessment of sampling uncertainty. The standard Frechet limit has GEV parameters (ξ,μ,σ)=(5/9,1,5/9)(\xi,\mu,\sigma)=(5/9,1,5/9) on the normalized scale.

Caveats

References

Source and adaptation

Adapted from incerto-wiki, content/concepts/distributions/generalized-extreme-value.md, revision 9717c9c (2026-09-13 Batch 2 import). Copyright (c) 2023 xshi19. Licensed under MIT. Links, notation, and qualifications were adapted for this site; executable figures and simulations were replaced with static calculations. No upstream execution or formal-proof verification is claimed for this adaptation.

MIT permission notice

MIT License

Copyright (c) 2023 xshi19

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the “Software”), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED “AS IS”, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

References
  1. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer Berlin Heidelberg. 10.1007/978-3-642-33483-2
  2. Resnick, S. I. (2007). Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. In Springer Series in Operations Research and Financial Engineering. Springer New York. 10.1007/978-0-387-45024-7
  3. Coles, S. (2001). An Introduction to Statistical Modeling of Extreme Values. In Springer Series in Statistics. Springer London. 10.1007/978-1-4471-3675-0