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.

Frechet Distribution and Frechet-Type Limits

Statement

The Frechet distribution is the positive-shape extreme-value law for normalized maxima with heavy right tails. In the standard one-parameter form with shape α>0\alpha>0, its CDF is

Φα(y)={0,y≤0,exp⁡(−y−α),y>0.\Phi_\alpha(y) = \begin{cases} 0, & y\le0,\\ \exp(-y^{-\alpha}), & y>0. \end{cases}

For y>0y>0, the density is

ϕα(y)=αy−α−1exp⁡(−y−α).\phi_\alpha(y) = \alpha y^{-\alpha-1}\exp(-y^{-\alpha}).

It is the ξ>0\xi>0 member of the generalized extreme-value distribution. With ξ=1/α\xi=1/\alpha, location μ=1\mu=1, and scale σ=ξ\sigma=\xi,

Hξ,1,ξ(y)=exp⁡[−y−1/ξ]=Φα(y),y>0.H_{\xi,1,\xi}(y) = \exp\left[-y^{-1/\xi}\right] = \Phi_\alpha(y), \qquad y>0.

The phrase Frechet-type refers to distributions whose normalized maxima converge to this law. In the usual right-tail setting, this is the regularly varying maximum-domain-of-attraction case: a Pareto-type survival tail with exponent α\alpha has extreme-value index ξ=1/α>0\xi=1/\alpha>0. Here FF is the parent CDF, Fˉ=1−F\bar F=1-F its survival function, XiX_i its iid draws, and xm>0x_m>0 the lower cutoff in the exact Pareto example. Φα\Phi_\alpha and ϕα\phi_\alpha denote the Frechet CDF and density.

Fréchet density, CDF and log-log survival for three exponents, with the support boundary at zero and curved survival approaching a power tail.

Figure 1:Unit-scale Fréchet laws with α=0.8,1.5,3\alpha=0.8,1.5,3, using the same exponent colors as the Pareto card. The boundary is zero. Survival approaches a power tail at large positive thresholds; it is not an exact power over the whole support. The last panel’s lower visible probability is 10-8.

Maxima interpretation

Frechet-type behavior is a statement about extremes, not necessarily about the full parent distribution. If X1,…,XnX_1,\dots,X_n are iid and the right tail is regularly varying with exponent α>0\alpha>0, then an appropriate positive normalization ana_n puts the maximum

Mn=max⁡1≤i≤nXiM_n=\max_{1\le i\le n} X_i

on the Frechet scale. The full regular-variation domain-of-attraction theorem is cited in the references; the exact Pareto calculation below shows the core mechanism.

For Pareto Type I with lower cutoff xmx_m and exponent α\alpha, take an=xmn1/αa_n=x_m n^{1/\alpha}. Then for fixed y>0y>0 and all large enough nn,

P(Mn/an≤y)=(1−1nyα)n→exp⁡(−y−α)=Φα(y).\mathbb P(M_n/a_n\le y) = \left(1-\frac{1}{ny^\alpha}\right)^n \to \exp(-y^{-\alpha}) = \Phi_\alpha(y).

The same shape parameter is therefore seen in two coordinates:

αfor the Pareto-type survival exponent,ξ=1αfor the extreme-value index.\alpha \quad \text{for the Pareto-type survival exponent}, \qquad \xi=\frac1\alpha \quad \text{for the extreme-value index}.

Tail behavior

Although the Frechet law is a limit law for maxima, it is itself heavy-tailed. Its survival function satisfies

Φˉα(y)=1−exp⁡(−y−α)∼y−α,y→∞.\bar\Phi_\alpha(y) = 1-\exp(-y^{-\alpha}) \sim y^{-\alpha}, \qquad y\to\infty.

Thus the Frechet survival tail is regularly varying with index −α-\alpha. This is why the ξ>0\xi>0 GEV regime is also the Pareto-type regime: the block-maximum limit and the parent survival tail share the same reciprocal coordinate ξ=1/α\xi=1/\alpha.

Finite-block calculation

For any α>0\alpha>0, the exact Pareto probability at y=1y=1 is (1−1/n)n(1-1/n)^n. The limit is e−1e^{-1}, independently of α\alpha because the normalization already includes n1/αn^{1/\alpha}.

Block size nnExact P(Mn/an≤1)\mathbb P(M_n/a_n\le1)Frechet limit
100.3486780.367879
1000.3660320.367879
6000.3675730.367879

For fixed y>0y>0, set t=y−αt=y^{-\alpha}. Expanding the logarithm gives

nlog⁡(1−t/n)=−t−t22n+O(n−2).n\log(1-t/n)=-t-\frac{t^2}{2n}+O(n^{-2}).

This explains the finite-block correction for an exact Pareto parent at a fixed yy. It is not a uniform error bound over all thresholds or a rate for arbitrary regularly varying parents. No simulated or fitted CDF is used here.

Exact CDFs of normalized Pareto maxima for block sizes 1, 4, 16 and 64 approach the Fréchet CDF. A second panel displays the finite-block CDF difference.

Figure 2:Exact population CDFs for iid Pareto parents with α=1.5\alpha=1.5 and normalization an=xmn1/αa_n=x_m n^{1/\alpha}; no maxima are simulated or fitted. Dotted lines mark the finite-block lower boundaries n−1/αn^{-1/\alpha}. The second panel plots the limit CDF minus each finite-block CDF, making small corrections visible without displacing overlapping curves.

Use frechet and the existing pareto_normalized_maximum_cdf.

from xmath.incerto.distributions import frechet
from xmath.incerto.glossary import pareto_normalized_maximum_cdf

for n in (1, 4, 16, 64):
    finite = pareto_normalized_maximum_cdf([1], alpha=1.5, n=n)[0]
    print(f"n={n:d}: cdf(1)={finite:.4g}, limit={frechet(1.5).cdf(1):.4g}")

Output

n=1: cdf(1)=0, limit=0.3679
n=4: cdf(1)=0.3164, limit=0.3679
n=16: cdf(1)=0.3561, limit=0.3679
n=64: cdf(1)=0.365, limit=0.3679

Regenerate both figures with uv run python demos/incerto/distribution_cards.py.

Caveats

References

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