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.

Tail Class Catalog

Main facts

We catalog common right-tail examples by two asymptotic properties: subexponentiality and regular variation. The classifications are cited results; the appendix proves the elementary convolution decomposition. Standard references are Bingham, Goldie, and Teugels for regular variation 1987, and Embrechts, Klueppelberg, and Mikosch 1997 plus Foss, Korshunov, and Zachary 2013, 2nd ed. for subexponential tails.

Let Fˉ(x)=P(X>x)\bar F(x)=\mathbb P(X>x) be the survival function. In the table, “regularly varying” means FˉRVα\bar F\in RV_{-\alpha} for some α>0\alpha>0, so Fˉ(tx)/Fˉ(x)tα\bar F(tx)/\bar F(x)\to t^{-\alpha} for fixed t>0t>0. Here FF is the CDF, xx is the threshold, and α\alpha is a positive tail exponent. All examples below have nonnegative support, with positive scale parameters, positive gamma shape and rate, and nonzero lognormal log-variance. Shared notation is planned.

DistributionSubexponential?Regularly varying?Diagnostic tail behavior
Pareto, Fˉ(x)=(xm/x)α\bar F(x)=(x_m/x)^\alphaYesYesExact power tail.
LognormalYesNoHeavier than any exponential, lighter than any power.
Weibull, Fˉ(x)=exp(xβ)\bar F(x)=\exp(-x^\beta) with 0<β<10<\beta<1YesNoStretched exponential tail.
ExponentialNoNoMemoryless light tail.
GammaNoNoExponential tail with a polynomial factor.

For the lognormal row, “heavier than any exponential, lighter than any power” means the survival function satisfies

ecxFˉ(x)for every c>0,xpFˉ(x)0for every p>0.e^{cx}\bar F(x)\to\infty \quad\text{for every }c>0, \qquad x^p\bar F(x)\to0 \quad\text{for every }p>0.

The inclusion direction to remember is:

FˉRVα, α>0F is subexponential,\bar F\in RV_{-\alpha},\ \alpha>0 \quad\Longrightarrow\quad F\text{ is subexponential},

for distributions supported on [0,)[0,\infty), or under the standard corresponding right-tail assumptions. The converse fails: lognormal and stretched-Weibull tails are subexponential without being regularly varying.

Reusable diagnostics

For a nonnegative continuous law with unbounded support, compare

At(x)=Fˉ(tx)Fˉ(x),R(x)=P(X1+X2>x)2Fˉ(x),A_t(x)=\frac{\bar F(tx)}{\bar F(x)},\qquad R(x)=\frac{\mathbb P(X_1+X_2>x)}{2\bar F(x)},

where X1,X2X_1,X_2 are iid and t>1t>1 is fixed. The first ratio checks multiplicative tail scaling; the second is the defining two-summand subexponential ratio. These are population quantities. Finite evaluations of them do not certify an asymptotic class from data.

Distribution and parametersMultiplier ratio At(x)A_t(x)Limit of R(x)R(x)
Pareto Type I, xxmx\ge x_mtαt^{-\alpha}1
Lognormal, logXN(0,1)\log X\sim N(0,1)Φˉ(logx+logt)/Φˉ(logx)0\bar\Phi(\log x+\log t)/\bar\Phi(\log x)\to01
Weibull, 0<β<10<\beta<1exp[(tβ1)xβ]0\exp[-(t^\beta-1)x^\beta]\to01
Exponential, rate 1e(t1)x0e^{-(t-1)x}\to0\infty
Gamma, shape 2, rate 1e(t1)x(1+tx)/(1+x)0e^{-(t-1)x}(1+tx)/(1+x)\to0\infty

Here Φˉ\bar\Phi is the standard normal survival function. For the lognormal row, the normal-tail asymptotic Φˉ(z)ϕ(z)/z\bar\Phi(z)\sim\phi(z)/z, with ϕ\phi the standard normal density, yields

At(x)logxlogx+logtexp[(logt)logx(logt)22]0.A_t(x)\sim \frac{\log x}{\log x+\log t} \exp\left[-(\log t)\log x-\frac{(\log t)^2}{2}\right] \longrightarrow0.

The asymptotic for Φˉ\bar\Phi follows from the complementary-error-function expansion in NIST DLMF, Section 7.12. Thus the multiplier diagnostic separates power tails from all four other examples, but cannot distinguish subexponential lognormal and stretched-Weibull tails from the two light-tailed examples. The classifications of lognormal and stretched Weibull are cited results; a ratio table is not a proof of their subexponentiality.

Exact light-tail checks

Exponential and integer-shape gamma sums give elementary checks of RR. With shape–rate notation, Gamma(k,1)\operatorname{Gamma}(k,1) has survival exj=0k1xj/j!e^{-x}\sum_{j=0}^{k-1}x^j/j! for x0x\ge0. Independence adds the shapes. Consequently,

RExp(1)(x)=1+x2,RGamma(2,1)(x)=1+x+x2/2+x3/62(1+x).R_{\operatorname{Exp}(1)}(x)=\frac{1+x}{2}, \qquad R_{\operatorname{Gamma}(2,1)}(x) =\frac{1+x+x^2/2+x^3/6}{2(1+x)}.
Threshold xxExponential sum ratioGamma shape 2 sum ratio
21.519/181.05619/18\approx1.056
105.5683/6610.348683/66\approx10.348
10050.5515303/606850.335515303/606\approx850.335

Both ratios diverge, although the gamma ratio happens to be near one at x=2x=2. A single moderate-threshold observation is insufficient to infer a limiting tail class. The Cramér condition provides the exponential-moment contrast for these two light-tailed laws.

Caveats

Appendix: convolution calculation

This appendix proves the elementary two-summand convolution decomposition used by the subexponential diagnostic. We work with independent continuous variables, then specialize to an iid lower-bounded distribution. Discrete, dependent, and two-sided variants require separate assumptions.

For independent continuous random variables X1X_1 and X2X_2 with densities f1f_1, f2f_2 and survival functions Fˉ1\bar F_1, Fˉ2\bar F_2,

P(X1+X2>x)=f1(y)Fˉ2(xy)dy.\mathbb P(X_1+X_2>x) = \int_{-\infty}^{\infty} f_1(y)\bar F_2(x-y)\,dy.

To see the iid reduction, take a common density ff, survival function Fˉ\bar F, and lower support endpoint aa, with x>2ax>2a. Split the event {X1+X2>x}\{X_1+X_2>x\} into three cases, ignoring probability-zero boundary points: first X2x/2X_2\le x/2, second X1x/2X_1\le x/2, and third both variables exceed x/2x/2. In the first case, conditioning on X2=yX_2=y with ayx/2a\le y\le x/2 leaves the requirement X1>xyX_1>x-y, so this part contributes ax/2f(y)Fˉ(xy)dy\int_a^{x/2} f(y)\bar F(x-y)\,dy. The second case contributes the same quantity by iid symmetry. The remaining upper-right square has probability Fˉ(x/2)2\bar F(x/2)^2 by independence.

For a=0a=0, the partition can be recorded as disjoint regions rather than a contour plot. Boundaries have zero probability under the continuous-law assumption.

Region within {X1+X2>x}\{X_1+X_2>x\}Probability contribution
0X2x/20\le X_2\le x/2, X1>xX2X_1>x-X_20x/2f(y)Fˉ(xy)dy\int_0^{x/2}f(y)\bar F(x-y)\,dy
0X1x/20\le X_1\le x/2, X2>xX1X_2>x-X_1The same integral by iid symmetry
X1>x/2X_1>x/2, X2>x/2X_2>x/2Fˉ(x/2)2\bar F(x/2)^2

For iid variables with lower support endpoint aa, symmetry gives the more useful half-line formula

P(X1+X2>x)=2ax/2f(y)Fˉ(xy)dy+Fˉ(x/2)2.\mathbb P(X_1+X_2>x) = 2\int_a^{x/2} f(y)\bar F(x-y)\,dy+\bar F(x/2)^2.

Dividing by 2Fˉ(x)2\bar F(x) yields a normalized integral:

R(x)=ax/2f(y)Fˉ(xy)Fˉ(x)dy+Fˉ(x/2)22Fˉ(x)=ax/2exp{logf(y)+logFˉ(xy)logFˉ(x)}dy+12exp{2logFˉ(x/2)logFˉ(x)}.\begin{aligned} R(x) &=\int_a^{x/2}\frac{f(y)\bar F(x-y)}{\bar F(x)}\,dy +\frac{\bar F(x/2)^2}{2\bar F(x)}\\ &=\int_a^{x/2} \exp\{\log f(y)+\log\bar F(x-y)-\log\bar F(x)\}\,dy\\ &\quad+\frac12\exp\{2\log\bar F(x/2)-\log\bar F(x)\}. \end{aligned}

The logarithmic form applies where the factors are positive, with zero terms understood by limits. It avoids forming tiny unnormalized probabilities before taking their ratio. The decomposition is an ordinary conditioning argument; no numerical integration or executable package helper is needed for it.

References

Source and adaptation

Adapted from incerto-wiki, content/concepts/distributions/tail-class-catalog.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. Bingham, N. H., Goldie, C. M., & Teugels, J. L. (1987). Regular Variation. Cambridge University Press. 10.1017/cbo9780511721434
  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