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.

Pickands-Balkema-de Haan Theorem

Statement

Let FF be a distribution function with finite or infinite right endpoint xF=sup⁡{x:F(x)<1}x_F=\sup\{x:F(x)<1\}. For a high threshold u<xFu<x_F, define the conditional excess distribution

Fu(y)=P(X−u≤y∣X>u)=F(u+y)−F(u)1−F(u),0≤y<xF−u.F_u(y)=\mathbb P(X-u\le y\mid X>u) =\frac{F(u+y)-F(u)}{1-F(u)}, \qquad 0\le y<x_F-u.

If FF is in the maximum domain of attraction of an extreme-value law with extreme-value index ξ\xi, then there is a positive scale function β(u)\beta(u) such that

sup⁡0≤y<xF−u∣Fu(y)−Gξ,β(u)(y)∣→0,u↑xF.\sup_{0\le y<x_F-u} \left|F_u(y)-G_{\xi,\beta(u)}(y)\right| \to 0, \qquad u\uparrow x_F.

The comparison is over the excess support, using the full GPD CDF: for ξ<0\xi<0, set Gξ,β(y)=1G_{\xi,\beta}(y)=1 at and beyond −β/ξ-\beta/\xi. The generalized Pareto term has scale β>0\beta>0 and is

Gξ,β(y)=1−(1+ξyβ)−1/ξ,ξ≠0,G_{\xi,\beta}(y) = 1-\left(1+\frac{\xi y}{\beta}\right)^{-1/\xi}, \qquad \xi\ne0,

on the support where 1+ξy/β>01+\xi y/\beta>0, and

G0,β(y)=1−e−y/βG_{0,\beta}(y)=1-e^{-y/\beta}

for ξ=0\xi=0. The theorem is the asymptotic justification for the peaks-over-threshold model: sufficiently high threshold exceedances are modeled by a generalized Pareto distribution.

XX denotes a draw from FF, uu the threshold, xFx_F its right endpoint, ξ\xi the extreme-value index, and β(u)>0\beta(u)>0 the excess scale. The maximum-domain-of-attraction assumption means that for iid draws, some an>0a_n>0 and bnb_n give (Mn−bn)/an⇒Hξ,0,1(M_n-b_n)/a_n\Rightarrow H_{\xi,0,1}. Here MnM_n is the sample maximum and ⇒\Rightarrow is convergence in distribution. The GEV page defines HH. See Shared notation.

Threshold intuition

Block-maxima theory says that properly normalized maxima have only a small number of possible limiting shapes, collected in the Generalized Extreme-Value Distribution. The Pickands-Balkema-de Haan theorem says the matching threshold view has the same discipline: once we condition on being far enough into the tail, the remaining excess has an approximately Generalized Pareto Distribution.

For fat-tail work, the case ξ>0\xi>0 is the main bridge. It corresponds to a regularly varying right tail with exponent α=1/ξ\alpha=1/\xi. The theorem does not say every high observation is generated by a clean Pareto law. It says the excess distribution over moving high thresholds has a universal limit shape under the same domain-of-attraction assumptions used in extreme-value theory.

Exact Pareto excess special case

The full Pickands-Balkema-de Haan theorem is cited here, not proved. The calculation below proves the exact Pareto special case: once the threshold is above the lower cutoff, the excess distribution is already generalized Pareto, not merely asymptotically close to it.

The full theorem is a classical result of Balkema and de Haan and, independently, Pickands. Its proof uses the equivalence between convergence of normalized maxima and convergence of normalized threshold excesses. We record the statement and use exact Pareto algebra as a checkable special case.

If XX is Pareto Type I with lower cutoff xmx_m and tail exponent α\alpha, then for u≥xmu\ge x_m,

P(X−u>y∣X>u)=P(X>u+y)P(X>u)=(1+yu)−α.\mathbb P(X-u>y\mid X>u) = \frac{\mathbb P(X>u+y)}{\mathbb P(X>u)} = \left(1+\frac{y}{u}\right)^{-\alpha}.

This is exactly a generalized Pareto survival function with

ξ=1α,β(u)=uα.\xi=\frac1\alpha, \qquad \beta(u)=\frac{u}{\alpha}.

Thus exact Pareto tails do not merely approach the generalized Pareto form; they have it at every threshold above xmx_m.

Exact threshold calibration

For α=1.7=17/10\alpha=1.7=17/10, the population excess parameters are ξ=10/17\xi=10/17 and β(u)=10u/17\beta(u)=10u/17. The scaled excess Z=(X−u)/uZ=(X-u)/u, conditional on X>uX>u, has survival (1+z)−1.7(1+z)^{-1.7}.

Threshold u≥xmu\ge x_mShape ξ\xiScale β(u)\beta(u)P(Z>1∣X>u)\mathbb P(Z>1\mid X>u)
xmx_m10/1710/1710xm/1710x_m/172-1.7
2xm2x_m10/1710/1720xm/1720x_m/172-1.7
4xm4x_m10/1710/1740xm/1740x_m/172-1.7

These are exact population identities, not fitted estimates. In finite data, estimated shape and scale fluctuate, and real tails may approach the GPD family slowly. Threshold selection examines that modeling choice; the theorem supplies no finite-sample cutoff.

Caveats

References

References
  1. Balkema, A. A., & de Haan, L. (1974). Residual Life Time at Great Age. The Annals of Probability, 2(5), 792–804. 10.1214/aop/1176996548
  2. Pickands, J. (1975). Statistical Inference Using Extreme Order Statistics. The Annals of Statistics, 3(1), 119–131. 10.1214/aos/1176343003
  3. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer. 10.1007/978-3-642-33483-2
  4. Coles, S. (2001). An Introduction to Statistical Modeling of Extreme Values. Springer. 10.1007/978-1-4471-3675-0