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(XuyX>u)=F(u+y)F(u)1F(u),0y<xFu.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

sup0y<xFuFu(y)Gξ,β(u)(y)0,uxF.\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)=1ey/β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 (Mnbn)/anHξ,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. A shared notation page is planned.

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 uxmu\ge x_m,

P(Xu>yX>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=(Xu)/uZ=(X-u)/u, conditional on X>uX>u, has survival (1+z)1.7(1+z)^{-1.7}.

Threshold uxmu\ge x_mShape ξ\xiScale β(u)\beta(u)P(Z>1X>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

Source and adaptation

Adapted from incerto-wiki, content/concepts/theorems/pickands-balkema-de-haan.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. Balkema, A. A., & de Haan, L. (1974). Residual Life Time at Great Age. The Annals of Probability, 2(5). 10.1214/aop/1176996548
  2. III, J. P. (1975). Statistical Inference Using Extreme Order Statistics. The Annals of Statistics, 3(1). 10.1214/aos/1176343003
  3. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer Berlin Heidelberg. 10.1007/978-3-642-33483-2
  4. 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