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.

Mean Excess Function

Statement

For a random variable XX and a threshold uu with P(X>u)>0\mathbb P(X>u)>0 and finite conditional expectation above uu, the mean excess function is

e(u)=E[XuX>u].e(u)=\mathbb E[X-u\mid X>u].

It measures the expected overshoot beyond uu after the threshold has already been crossed. For a Pareto Type I tail with exponent α>1\alpha>1 and uxmu\ge x_m,

e(u)=uα1.e(u)=\frac{u}{\alpha-1}.

For a generalized Pareto variable with location zero, shape ξ\xi, scale β>0\beta>0, and threshold uu in its support with β+ξu>0\beta+\xi u>0, the mean excess function is linear when ξ<1\xi<1:

e(u)=β+ξu1ξ.e(u)=\frac{\beta+\xi u}{1-\xi}.

The exponential distribution is the boundary case ξ=0\xi=0, where e(u)e(u) is constant. Here uu is the threshold, ξ\xi the GPD shape, and β\beta its scale at threshold zero.

Diagnostic intuition

The mean excess function asks what remains after an event is already large. Thin-tail intuition often expects the remaining excess to be tame once a high threshold has been crossed. A Pareto tail says the opposite: the expected additional excess grows in proportion to the threshold itself.

This makes the mean excess plot a practical threshold diagnostic. A roughly flat plot suggests exponential-like exceedances. A roughly increasing linear plot suggests heavy-tail generalized Pareto behavior. A downward line suggests a finite endpoint. Strong curvature usually says the chosen threshold range is mixing body and tail behavior, or that the model class is too simple.

Pareto derivation and GPD sketch

We derive the Pareto mean-excess formula from the conditional survival function. The generalized Pareto linear formula follows from standard GPD threshold stability and the GPD mean formula. The GPD page derives the stability identity; its moment formulas are cited rather than reproved here.

For Pareto Type I,

P(X>x)=(xmx)α,xxm.\mathbb P(X>x)=\left(\frac{x_m}{x}\right)^\alpha, \qquad x\ge x_m.

For uxmu\ge x_m and α>1\alpha>1,

e(u)=E[XuX>u]=0P(Xu>yX>u)dy.e(u) = \mathbb E[X-u\mid X>u] = \int_0^\infty \mathbb P(X-u>y\mid X>u)\,dy.

Taking the ratio Fˉ(u+y)/Fˉ(u)\bar F(u+y)/\bar F(u) for y0y\ge0 gives

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

Therefore

e(u)=0(1+yu)αdy=u1zαdz=uα1.e(u) = \int_0^\infty \left(1+\frac{y}{u}\right)^{-\alpha}\,dy = u\int_1^\infty z^{-\alpha}\,dz = \frac{u}{\alpha-1}.

For a generalized Pareto distribution, the threshold-stability property gives another generalized Pareto distribution above threshold uu with updated scale β+ξu\beta+\xi u, provided uu is in the support and β+ξu>0\beta+\xi u>0. Its mean exists only for ξ<1\xi<1, and equals (β+ξu)/(1ξ)(\beta+\xi u)/(1-\xi).

Mean-excess diagnostic

Exact theoretical shapes provide a baseline for the empirical diagnostic. Always report how many observations exceed each threshold.

The theoretical shapes can be compared directly:

DistributionThreshold rangeMean excess e(u)e(u)
Pareto, xm=1x_m=1, α=1.6\alpha=1.6u1u\ge1u/0.6u/0.6
Exponential with mean 1u0u\ge01
Uniform on [0,5][0,5]0u<50\le u<5(5u)/2(5-u)/2

For observations x1,,xnx_1,\ldots,x_n, let k(u)=i1{xi>u}k(u)=\sum_i\mathbf1_{\{x_i>u\}}. The empirical diagnostic is

e^n(u)=i:xi>u(xiu)k(u),k(u)>0.\widehat e_n(u)=\frac{\sum_{i:x_i>u}(x_i-u)}{k(u)},\qquad k(u)>0.

For the fixed observations (1,2,2,4)(1,2,2,4) at u=2u=2, only 4 exceeds the threshold, so k(2)=1k(2)=1 and e^n(2)=2\widehat e_n(2)=2. At u=4u=4 the mean excess estimate is undefined because the exceedance count is zero. The counting example uses the same strict inequality convention.

What to notice. Pareto mean excess rises with the threshold, exponential mean excess is flat, and bounded-tail mean excess slopes downward toward the endpoint. Empirical values are diagnostics; high thresholds aim to reduce tail-model bias but leave fewer exceedances and more noise.

Caveats

References

Source and adaptation

Adapted from incerto-wiki, content/concepts/theorems/mean-excess-function.md, revision 9717c9c (2026-09-13 import). Copyright (c) 2023 xshi19. Licensed under MIT. Links and notation 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. Davison, A. C., & Smith, R. L. (1990). Models for Exceedances Over High Thresholds. Journal of the Royal Statistical Society Series B: Statistical Methodology, 52(3), 393–425. 10.1111/j.2517-6161.1990.tb01796.x
  2. 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
  3. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer Berlin Heidelberg. 10.1007/978-3-642-33483-2