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.

Subexponentiality

Statement

Let X1,X2,X_1,X_2,\dots be iid nonnegative random variables with distribution FF and unbounded right support, so Fˉ(x)=P(X>x)>0\bar F(x)=\mathbb P(X>x)>0 for every finite xx. The distribution is subexponential when, for two independent copies,

P(X1+X2>x)2P(X>x)1,x.\frac{\mathbb P(X_1+X_2>x)}{2\mathbb P(X>x)} \to 1, \qquad x\to\infty.

Equivalently, for each fixed integer n1n\ge1, with Sn=X1++XnS_n=X_1+\cdots+X_n and Mn=max1inXiM_n=\max_{1\le i\le n}X_i,

P(Sn>x)P(Mn>x)nFˉ(x).\mathbb P(S_n>x) \sim \mathbb P(M_n>x) \sim n\bar F(x).

Equivalently,

P(Sn>x,  Mnx)=o(Fˉ(x)).\mathbb P(S_n>x,\; M_n\le x)=o(\bar F(x)).

This is the one-big-jump principle: for a fixed number of summands, a very large sum is usually caused by one summand being very large, not by all summands moving moderately together.

The maximum-tail asymptotic P(Mn>x)nFˉ(x)\mathbb P(M_n>x)\sim n\bar F(x) is not the distinctive part of the theorem: it follows for every iid distribution with an unbounded right tail and fixed nn. The subexponential content is that the sum tail has the same first-order asymptotic as the maximum tail.

Regularly varying survival tails with positive exponent, FˉRVα\bar F\in RV_{-\alpha} for α>0\alpha>0, are a standard subexponential class. We cite that theorem and use the exact Pareto family as a checkable example.

Here FF is the common CDF, nn is fixed while xx grows, and α>0\alpha>0 and xm>0x_m>0 are the exponent and cutoff in the Pareto example. The notation a(x)b(x)a(x)\sim b(x) means a(x)/b(x)1a(x)/b(x)\to1; o(b(x))o(b(x)) denotes a term negligible after division by b(x)b(x).

One-big-jump intuition

Subexponential does not mean “lighter than exponential,” and it is unrelated to the concentration-theory phrase “sub-exponential random variable” for a light-tailed ψ1\psi_1 condition. In this context it means the opposite: the tail is heavy enough that convolution barely changes the leading-order tail probability.

For two iid summands, a sufficiently large sum is about twice as likely as one summand exceeding the same threshold. For nn fixed summands, the leading factor is nn.

This is the mathematical version of a recurring Incerto warning. In thin-tail settings, a large aggregate deviation is often the result of many small coordinated deviations. In subexponential settings, the aggregate tail is usually dominated by a single large term. That changes how sums, ruin events, insurance losses, and sample moments should be read.

What we can prove directly

The statement collects several equivalent forms and one important sufficient class. We prove the elementary maximum-tail asymptotic and an exact Pareto two-summand check. The full equivalence theory and the theorem that regularly varying tails are subexponential are cited, because their proofs need more convolution-tail machinery than we develop here.

The theorem that regularly varying tails are subexponential goes back to Chistyakov’s convolution-tail criterion for sums of independent positive random variables 1964. Standard modern references with proofs are Embrechts, Klueppelberg, and Mikosch, Appendix A3 1997, and Foss, Korshunov, and Zachary, Chapter 3 2013. We do not reprove the full theorem, but the exact Pareto case gives the right calibration.

Universal maximum tail

For any distribution with unbounded right tail, Fˉ(x)0\bar F(x)\to0 as xx\to\infty. For fixed nn, iid independence gives

P(Mnx)=P(X1x,,Xnx)=F(x)n=(1Fˉ(x))n.\mathbb P(M_n\le x) =\mathbb P(X_1\le x,\dots,X_n\le x) =F(x)^n =\left(1-\bar F(x)\right)^n.

Therefore

P(Mn>x)=1(1Fˉ(x))n\mathbb P(M_n>x) =1-\left(1-\bar F(x)\right)^n

Set u=Fˉ(x)u=\bar F(x). Since u0u\to0 and nn is fixed, the binomial expansion gives

1(1u)n=nu(n2)u2+O(u3),1-(1-u)^n =nu-\binom{n}{2}u^2+O(u^3),

and therefore

1(1u)nnu=1n12u+O(u2)1.\frac{1-(1-u)^n}{nu} =1-\frac{n-1}{2}u+O(u^2) \to1.

Substituting back,

P(Mn>x)nFˉ(x).\mathbb P(M_n>x) \sim n\bar F(x).

This maximum formula only uses iid independence, fixed nn, and Fˉ(x)0\bar F(x)\to0. For nonnegative variables, Mn>xM_n>x implies Sn=X1++Xn>xS_n=X_1+\cdots+X_n>x, so it is a lower bound for the sum tail. What is special about subexponential tails is the matching upper bound:

P(Sn>x,  Mnx)=o(Fˉ(x)).\mathbb P(S_n>x,\;M_n\le x)=o(\bar F(x)).

Those are the configurations where several observations are moderately large but none exceeds xx. Subexponentiality says this residual event is negligible at the scale of one tail probability.

Pareto upper-bound check

For the Pareto tail Fˉ(x)=(xm/x)α\bar F(x)=(x_m/x)^\alpha with xxmx\ge x_m, the missing upper-bound step can be checked directly for n=2n=2. Let

Ax={X1+X2>x,  M2x}.A_x=\{X_1+X_2>x,\;M_2\le x\}.

Choose h=xm(x/xm)γh=x_m(x/x_m)^\gamma with 1/2<γ<11/2<\gamma<1. This keeps the units consistent and gives h/x0h/x\to0 and hh\to\infty.

The residual event AxA_x is the region in [0,x]2[0,x]^2 above X1+X2=xX_1+X_2=x. It is covered by two thin strips and a region where both coordinates exceed hh. The following table replaces the geometric plot.

Cover eventProbability or bound
xh<X1xx-h<X_1\le xFˉ(xh)Fˉ(x)\bar F(x-h)-\bar F(x)
xh<X2xx-h<X_2\le xFˉ(xh)Fˉ(x)\bar F(x-h)-\bar F(x)
X1>hX_1>h and X2>hX_2>hFˉ(h)2\bar F(h)^2 by independence

Outside the cover region, a point has neither coordinate in (xh,x](x-h,x] and not both coordinates above hh. Thus at least one coordinate is h\le h, and the other is xh\le x-h, so X1+X2xX_1+X_2\le x. Therefore any point in AxA_x must lie in the cover: either one summand lies in (xh,x](x-h,x], or both summands exceed hh. Hence

P(Ax)2P(xh<Xx)+P(X1>h,  X2>h).\mathbb P(A_x) \le 2\mathbb P(x-h<X\le x) +\mathbb P(X_1>h,\;X_2>h).

For the first term,

P(xh<Xx)Fˉ(x)=Fˉ(xh)Fˉ(x)Fˉ(x)=(xxh)α10,\frac{\mathbb P(x-h<X\le x)}{\bar F(x)} = \frac{\bar F(x-h)-\bar F(x)}{\bar F(x)} = \left(\frac{x}{x-h}\right)^\alpha-1 \to0,

because h/x0h/x\to0. For the second term, independence gives

P(X1>h,  X2>h)Fˉ(x)=Fˉ(h)2Fˉ(x)=(x/xm)α(12γ)0,\frac{\mathbb P(X_1>h,\;X_2>h)}{\bar F(x)} = \frac{\bar F(h)^2}{\bar F(x)} = (x/x_m)^{\alpha(1-2\gamma)} \to0,

because γ>1/2\gamma>1/2. Thus P(Ax)=o(Fˉ(x))\mathbb P(A_x)=o(\bar F(x)). Since {S2>x}\{S_2>x\} is the disjoint union of {M2>x}\{M_2>x\} and AxA_x,

P(X1+X2>x)2Fˉ(x).\mathbb P(X_1+X_2>x)\sim 2\bar F(x).

The regularly varying theorem cited above generalizes this Pareto calculation: the maximum formula is universal, while the negligible residual event is the heavy-tail property.

A useful comparison table is:

DistributionSubexponential?Regularly varying?
ParetoYesYes
LognormalYesNo
Weibull with shape 0<β<10<\beta<1YesNo
ExponentialNoNo

The Tail Class Catalog gives a broader comparison. The lognormal and stretched Weibull entries are cited examples, not consequences of the Pareto proof.

The exponential is an exact nonexample. If XiExp(λ)X_i\sim\operatorname{Exp}(\lambda), then the universal maximum formula still holds:

P(M2>x)=2eλxe2λx2eλx.\mathbb P(M_2>x)=2e^{-\lambda x}-e^{-2\lambda x} \sim2e^{-\lambda x}.

But the sum tail is much larger:

P(X1+X2>x)=eλx(1+λx),\mathbb P(X_1+X_2>x)=e^{-\lambda x}(1+\lambda x),

so

P(X1+X2>x)2P(X>x)=1+λx2,\frac{\mathbb P(X_1+X_2>x)}{2\mathbb P(X>x)} =\frac{1+\lambda x}{2}\to\infty,

not 1.

The fixed-nn assumption is load-bearing. If nn grows with xx, additional large-deviation regimes can appear.

A static bound on the residual tail

The proof supplies a finite-threshold bound. Put t=x/xmt=x/x_m and take α=1.5\alpha=1.5, γ=0.75\gamma=0.75. For sufficiently large tt with xhxmx-h\ge x_m,

0P(Ax)Fˉ(x)2[(1t1/4)3/21]+t3/4.0\le\frac{\mathbb P(A_x)}{\bar F(x)} \le 2\left[(1-t^{-1/4})^{-3/2}-1\right]+t^{-3/4}.
t=x/xmt=x/x_mUpper bound (rounded upward)
1040.343428
1080.030381
10120.003004

Together with P(M2>x)=2Fˉ(x)Fˉ(x)2\mathbb P(M_2>x)=2\bar F(x)-\bar F(x)^2, this bounds the sum tail around its leading term. It illustrates how an asymptotic statement can require large thresholds even when the model is exactly Pareto. These are conservative analytic bounds, not measured convolution probabilities.

The fixed-nn, large-threshold limit differs from the growing-sample max-to-sum ratio. A finite-mean Pareto law is subexponential even though Mn/Sn0M_n/S_n\to0 almost surely as nn\to\infty.

Caveats

References

Source and adaptation

Adapted from incerto-wiki, content/concepts/theorems/subexponentiality.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. Chistyakov, V. P. (1964). A Theorem on Sums of Independent Positive Random Variables and Its Applications to Branching Random Processes. Theory of Probability & Its Applications, 9(4), 640–648. 10.1137/1109088
  2. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer Berlin Heidelberg. 10.1007/978-3-642-33483-2
  3. Foss, S., Korshunov, D., & Zachary, S. (2013). An Introduction to Heavy-Tailed and Subexponential Distributions. In Springer Series in Operations Research and Financial Engineering. Springer New York. 10.1007/978-1-4614-7101-1
  4. Bingham, N. H., Goldie, C. M., & Teugels, J. L. (1987). Regular Variation. Cambridge University Press. 10.1017/cbo9780511721434
  5. 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