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.

Pareto Moment Existence

Statement

Let XX have a Pareto Type I distribution with lower cutoff xm>0x_m>0 and tail exponent α>0\alpha>0:

Fˉ(x)=P(X>x)=(xmx)α,x≥xm.\bar F(x)=\mathbb P(X>x)=\left(\frac{x_m}{x}\right)^\alpha, \qquad x\ge x_m.

For any moment order p>0p>0,

E[Xp]<∞⟺p<α.\mathbb E[X^p] < \infty \quad\Longleftrightarrow\quad p<\alpha.

When p<αp<\alpha,

E[Xp]=αxmpα−p.\mathbb E[X^p]=\frac{\alpha x_m^p}{\alpha-p}.

At p=αp=\alpha, the truncated moment diverges logarithmically. For p>αp>\alpha, it diverges as a power of the upper cutoff. In particular, the Pareto mean exists only for α>1\alpha>1, and the variance exists only for α>2\alpha>2.

Here Fˉ\bar F denotes survival probability and E\mathbb E denotes expectation; pp is a positive moment order.

Proof by direct integration

We prove the Pareto Type I moment boundary and the moment formula by direct integration. The same calculation also explains the logarithmic boundary and power-divergent regimes in the truncated-moment section. The broader regularly varying moment test is handled by Karamata’s Theorem.

For x≥xmx\ge x_m, the Pareto density is

f(x)=αxmαx−(α+1).f(x)=\alpha x_m^\alpha x^{-(\alpha+1)}.

For p>0p>0,

E[Xp]=∫xm∞xpαxmαx−(α+1) dx=αxmα∫xm∞xp−α−1 dx.\mathbb E[X^p] = \int_{x_m}^{\infty}x^p\alpha x_m^\alpha x^{-(\alpha+1)}\,dx = \alpha x_m^\alpha \int_{x_m}^{\infty}x^{p-\alpha-1}\,dx.

The integral converges exactly when p−α−1<−1p-\alpha-1<-1, equivalently p<αp<\alpha. Evaluating the convergent case gives

αxmα⋅xmp−αα−p=αxmpα−p.\alpha x_m^\alpha \cdot\frac{x_m^{p-\alpha}}{\alpha-p} = \frac{\alpha x_m^p}{\alpha-p}.

Adding α+1\alpha+1 to both sides proves the algebraic equivalence p−α−1<−1⟺p<αp-\alpha-1<-1 \Longleftrightarrow p<\alpha.

This exact Pareto calculation is the simplest instance of the Karamata moment test for regularly varying tails.

Truncated-moment behavior

For a finite upper cutoff b≥xmb\ge x_m, define

Mp(b)=E[Xp1{X≤b}].M_p(b)=\mathbb E[X^p\mathbf 1_{\{X\le b\}}].

Direct integration gives

Mp(b)={αxmαα−p(xmp−α−bp−α),p<α,αxmαlog⁡(b/xm),p=α,αxmαp−α(bp−α−xmp−α),p>α.M_p(b)= \begin{cases} \dfrac{\alpha x_m^\alpha}{\alpha-p} \left(x_m^{p-\alpha}-b^{p-\alpha}\right), & p<\alpha,\\[1.1em] \alpha x_m^\alpha\log(b/x_m), & p=\alpha,\\[0.8em] \dfrac{\alpha x_m^\alpha}{p-\alpha} \left(b^{p-\alpha}-x_m^{p-\alpha}\right), & p>\alpha. \end{cases}

Thus the same theorem has three operational regimes: convergence to a finite moment, logarithmic boundary growth, and power growth.

For α=1.5\alpha=1.5 and xm=1x_m=1, the three regimes reduce to:

Moment order ppTruncated moment Mp(b)M_p(b), b≥1b\ge1Limit as b→∞b\to\infty
13(1−b−1/2)3(1-b^{-1/2})3
1.51.5log⁡b1.5\log bInfinite, logarithmic growth
23(b−1)3(\sqrt b-1)Infinite, power growth

For p<αp<\alpha, the truncated moment levels off. At p=αp=\alpha, it keeps growing, but only logarithmically. For p>αp>\alpha, the upper tail contributes a visible power-law rise.

Products and moment integrals

Vertical panels show independent Pareto product survivals, three truncated moment regimes on log-log axes, and a logarithmic integral divided by its Karamata approximation.

Top: exact independent Pareto products with exponents (1,2)(1,2) and (1,1)(1,1), compared with one Pareto of exponent 1; all minima are one. Middle: the three exact truncated moments in the table above, with α=1.5\alpha=1.5, on log-log axes. The zero value at b=1b=1 is omitted from this panel; the positive curves have no linear-to-log scale transition. Bottom: ∫1xtρlog⁡t dt\int_1^x t^\rho\log t\,dt divided by its Karamata leading term, for ρ=0,1\rho=0,1. These deterministic calculations illustrate a slowly varying correction and the strict moment boundary; no observations were sampled.

Use pareto_truncated_moment.

from xmath.incerto.power_laws import pareto_truncated_moment

for order in (1, 1.5, 2):
    moments = pareto_truncated_moment([10, 100, 10000], alpha=1.5, order=order)
    print(f"order={order:.4g}: [" + ", ".join(f"{v:.4g}" for v in moments) + "]")

Output

order=1: [2.051, 2.7, 2.97]
order=1.5: [3.454, 6.908, 13.82]
order=2: [6.487, 27, 297]

These are indicator-truncated population moments; no sample is censored.

Figure regeneration: uv run python demos/incerto/power_laws.py. Each panel uses 241 cutoffs from 1 to 104 (the middle omits the zero endpoint, leaving 240 points; the last begins at 1.01 to avoid a zero leading term). The Pareto product proof and Karamata calculation own the other two formulas. The cutoff is the indicator truncation in Mp(b)M_p(b), not conditioning on X≤bX\le b and not the capped moment Emin⁡(X,b)p\mathbb E\min(X,b)^p.

The chapter’s transformation discussion (§5.2.3, p. 100) and squared Student example (§5.6, pp. 103–105) use different moment orders Taleb, 2025. The moment-estimation note already explains why an existing second moment need not give a finite variance to the empirical average of squares.

Examples

The theorem is exact, but sample estimates can still look misleading. With α\alpha close to a boundary, a finite run can appear calm until a new large observation changes the empirical moment.

For α=1.2\alpha=1.2 and xm=1x_m=1, the population mean is 6 while the second raw moment is infinite. Any finite sample of finite observations still has a finite sample second moment. Its finiteness does not settle the population question.

Caveats

References

References
  1. Taleb, N. N. (2025). Statistical Consequences of Fat Tails (3rd ed.). https://arxiv.org/abs/2001.10488v4
  2. Bingham, N. H., Goldie, C. M., & Teugels, J. L. (1987). Regular Variation. Cambridge University Press. 10.1017/CBO9780511721434
  3. Resnick, S. I. (2007). Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. Springer. 10.1007/978-0-387-45024-7