Tutorials › Statistical Theory I - Probability and Distribution Theory › Existence Threshold for Pareto Moments

Normal and Tail Families · Tutorial 187 of 994

Existence Threshold for Pareto Moments

Integrate the Pareto density to locate the moment threshold, and see exactly what that threshold implies for the mean and variance.

Intermediate 9 min read

What You'll Learn

  • State the Type I Pareto density and identify its scale and tail-index parameters.
  • Derive the raw moment formula for every order below the tail index.
  • Prove that the moment diverges at and above the threshold.
  • Determine when the mean and variance are finite and calculate their values.
  • Distinguish an infinite second moment from a variance that is not defined because the mean is infinite.

Why a Moment Threshold Matters

A distribution can have a probability density that integrates to one and still lack a finite mean or variance. For the Type I Pareto distribution, the decisive feature is its slowly decreasing density in the far right tail. Integrating that density against a power of the observation reveals a precise cutoff: the moment of order \(r\) is finite exactly when \(r\) is less than the tail-index parameter.

The cutoff is useful in practice because it identifies which familiar summaries are available. Depending on the parameter, a Pareto model can have a finite mean but no finite variance, or it can have neither. The calculation also makes an important distinction: divergence exactly at the cutoff is logarithmic, while above it the divergence is polynomial.

Definition: A random variable \(X\) has a Type I Pareto distribution with scale \(x_m>0\) and tail index \(\alpha>0\), written here as \(X\sim\operatorname{Pareto}(x_m,\alpha)\), if its density is
$$ f(x)=\alpha x_m^\alpha x^{-(\alpha+1)},\qquad x\ge x_m, $$
and \(f(x)=0\) for \(x<x_m\). The scale \(x_m\) is the lower endpoint of the support, and \(\alpha\) controls how quickly the density decreases as \(x\) increases.

This is a probability density: since \(\alpha>0\),

$$ \int_{x_m}^{\infty}\alpha x_m^\alpha x^{-(\alpha+1)}\,dx =\alpha x_m^\alpha\left[\frac{-x^{-\alpha}}{\alpha}\right]_{x_m}^{\infty} =1. $$

Because \(x_m>0\), powers \(X^r\) are well-defined for every real \(r\). In many applications, “the \(r\)th moment” refers to a nonnegative integer \(r\); the integration result below applies more generally to every real order.

Deriving the Moment Formula

For a nonnegative random variable \(Y\), its expectation is defined as the integral of \(Y\) with respect to its distribution, with the value \(+\infty\) allowed. Since \(X^r\ge 0\), this definition applies for every real \(r\). Substituting the Pareto density reduces the question of moment existence to the convergence of an improper integral.

Theorem: If \(X\sim\operatorname{Pareto}(x_m,\alpha)\), where \(x_m>0\) and \(\alpha>0\), then for every real \(r\),
$$ E[X^r]<\infty \quad\Longleftrightarrow\quad r<\alpha. $$
When \(r<\alpha\), the moment is
$$ E[X^r]=\frac{\alpha x_m^r}{\alpha-r}. $$

Proof. Using the density and the definition of expectation for a nonnegative function,

$$ E[X^r] =\int_{x_m}^{\infty}x^r\alpha x_m^\alpha x^{-(\alpha+1)}\,dx =\alpha x_m^\alpha\int_{x_m}^{\infty}x^{r-\alpha-1}\,dx. $$

First suppose \(r<\alpha\). Then \(r-\alpha<0\), and for a finite upper endpoint \(R>x_m\),

$$ \int_{x_m}^{R}x^{r-\alpha-1}\,dx =\left[\frac{x^{r-\alpha}}{r-\alpha}\right]_{x_m}^{R} =\frac{x_m^{r-\alpha}-R^{r-\alpha}}{\alpha-r}. $$

As \(R\) tends to infinity, \(R^{r-\alpha}\) tends to zero because \(r-\alpha<0\). Therefore,

$$ E[X^r] =\alpha x_m^\alpha\frac{x_m^{r-\alpha}}{\alpha-r} =\frac{\alpha x_m^r}{\alpha-r}<\infty. $$

Now suppose \(r=\alpha\). The integrand in the last integral is \(x^{-1}\), so

$$ \alpha x_m^\alpha\int_{x_m}^{R}\frac{1}{x}\,dx =\alpha x_m^\alpha\log\left(\frac{R}{x_m}\right) \longrightarrow\infty \quad\text{as }R\to\infty. $$

Finally, if \(r>\alpha\), then \(r-\alpha>0\), and

$$ \int_{x_m}^{R}x^{r-\alpha-1}\,dx =\frac{R^{r-\alpha}-x_m^{r-\alpha}}{r-\alpha} \longrightarrow\infty \quad\text{as }R\to\infty. $$

The expectation therefore diverges both at \(r=\alpha\) and for \(r>\alpha\), and is finite for all \(r<\alpha\). This proves the stated equivalence and formula. \(\square\)

The lower endpoint being strictly positive matters when considering negative orders: there is no singularity near zero because \(X\ge x_m>0\). Thus every negative-order moment is finite. For example, if \(r=-1\), the formula gives \(E[X^{-1}]=\alpha/((\alpha+1)x_m)\).

Worked Example: A Pareto Model with Finite Mean and Variance

Let \(X\sim\operatorname{Pareto}(2,3)\), so the scale is \(x_m=2\) and the tail index is \(\alpha=3\). Determine \(E[X]\), \(E[X^2]\), and \(\operatorname{Var}(X)\).

The first moment is available because \(1<3\). Substitution into the formula gives

$$ E[X]=\frac{3\cdot 2^1}{3-1}=3. $$

The second moment is also available because \(2<3\):

$$ E[X^2]=\frac{3\cdot 2^2}{3-2}=12. $$

Since the mean and second moment are finite, the variance is \(E[X^2]-(E[X])^2\):

$$ \operatorname{Var}(X)=12-3^2=3. $$

The result is consistent with the general variance formula derived below: with \(\alpha=3\) and \(x_m=2\), it gives \(3\cdot 2^2/((3-2)(3-1)^2)=3\).

Consequences for the Mean and Variance

The mean is the moment of order one, so its existence threshold is \(\alpha>1\). The variance needs more care. The standard variance \(\operatorname{Var}(X)=E[(X-E[X])^2]\) requires a finite mean; in the usual finite-variance setting, it is enough that the second moment is finite. For Pareto distributions, a finite second moment also guarantees a finite mean, because \(\alpha>2\) implies \(\alpha>1\).

Result: For \(X\sim\operatorname{Pareto}(x_m,\alpha)\):
  • If \(\alpha>1\), then \(E[X]=\dfrac{\alpha x_m}{\alpha-1}\). If \(0<\alpha\le 1\), the mean is infinite.
  • If \(\alpha>2\), then \(\operatorname{Var}(X)=\dfrac{\alpha x_m^2}{(\alpha-1)^2(\alpha-2)}\), which is finite.
  • If \(1<\alpha\le 2\), the mean is finite but the variance is infinite.
  • If \(0<\alpha\le 1\), there is no finite mean and hence no finite variance in the standard sense.

Derivation. Taking \(r=1\) in the moment theorem gives, when \(\alpha>1\),

$$ E[X]=\frac{\alpha x_m}{\alpha-1}. $$

If \(0<\alpha\le 1\), the theorem says \(E[X]=\infty\). For the variance, first take \(r=2\). When \(\alpha>2\),

$$ E[X^2]=\frac{\alpha x_m^2}{\alpha-2}. $$

Both moments are finite in this case, so subtraction is valid:

$$ \begin{aligned} \operatorname{Var}(X) &=E[X^2]-(E[X])^2\\ &=\frac{\alpha x_m^2}{\alpha-2} -\left(\frac{\alpha x_m}{\alpha-1}\right)^2\\ &=\frac{\alpha x_m^2\left((\alpha-1)^2-\alpha(\alpha-2)\right)} {(\alpha-2)(\alpha-1)^2}\\ &=\frac{\alpha x_m^2}{(\alpha-2)(\alpha-1)^2}. \end{aligned} $$

For \(1<\alpha\le2\), the mean is finite, but \(E[X^2]=\infty\). The nonnegative quantity \((X-E[X])^2\) has infinite expectation as well: for \(X\ge 2E[X]\), it is at least \(X^2/4\), and the Pareto second-moment integral over that tail diverges. Thus the variance is infinite. When \(0<\alpha\le1\), the mean itself is not finite, so the standard variance about \(E[X]\) is not defined as a finite variance.

Worked Example: Finite Mean but Infinite Variance

Let \(X\sim\operatorname{Pareto}(4,1.5)\). Determine whether its mean and variance are finite, and calculate the mean and the moment of order \(1.25\).

Since \(1<1.5\), the mean is finite:

$$ E[X]=\frac{1.5\cdot4}{1.5-1} =\frac{6}{0.5}=12. $$

The second moment is not finite because \(2\not<1.5\); hence the variance is infinite. In contrast, \(1.25<1.5\), so the moment of order \(1.25\) is finite:

$$ E[X^{1.25}] =\frac{1.5\cdot4^{1.25}}{1.5-1.25} =6\cdot4^{1.25} =6\cdot4\sqrt{2} =24\sqrt{2} \approx33.94. $$

This example separates two questions that are sometimes conflated: the mean exists because the first-order threshold is passed, but the variance fails because the second-order threshold is not.

Worked Example: Divergence at the Exact Threshold

Let \(X\sim\operatorname{Pareto}(3,2)\). Find \(E[X]\), and determine whether \(E[X^2]\) is finite.

The mean is finite because \(1<2\):

$$ E[X]=\frac{2\cdot3}{2-1}=6. $$

The second-moment order is exactly the tail index, \(r=\alpha=2\). Directly, for an upper limit \(R\),

$$ \int_3^R x^2\left(2\cdot3^2x^{-3}\right)\,dx =18\int_3^R\frac{1}{x}\,dx =18\log\left(\frac{R}{3}\right) \longrightarrow\infty. $$

Thus \(E[X^2]=\infty\), and the variance is infinite even though the mean is finite. The logarithm in the calculation shows why equality at the threshold still gives divergence.

Common Mistake: Treating the Boundary as Finite

A frequent error is to infer that \(E[X^r]\) is finite whenever \(r\le\alpha\), perhaps because the power in the density seems to balance the power \(x^r\) at equality. At \(r=\alpha\), however, the powers combine to give an integrand proportional to \(1/x\), not a constant or a decaying integrable function. The integral of \(1/x\) over \([x_m,\infty)\) diverges.

A related mistake is to call the variance finite whenever the mean is finite. The mean tests the first moment; the variance requires control of the second moment. For \(1<\alpha\le2\), the first moment exists but the second does not. Conversely, when \(\alpha\le1\), it is more precise to say that no finite mean or standard finite variance exists than to use the variance formula with an infinite mean.

Key takeaway: For a Type I Pareto random variable, compare the moment order directly with \(\alpha\). Below \(\alpha\), \(E[X^r]=\alpha x_m^r/(\alpha-r)\); at or above \(\alpha\), the moment diverges. In particular, a finite mean requires \(\alpha>1\), and a finite variance requires \(\alpha>2\).

Check Your Understanding

Use the moment threshold and formula to answer each question.

  1. For \(X\sim\operatorname{Pareto}(2,2.5)\), is \(E[X^2]\) finite? If so, calculate it.
  2. For which values of \(\alpha\) is the Type I Pareto mean finite? State the mean when it is finite.
  3. Explain why \(E[X^\alpha]\) diverges even though the density itself integrates to one.
  4. If \(1<\alpha<2\), which is finite: the mean, the second moment, or both? What does this imply about the variance?
  5. For \(X\sim\operatorname{Pareto}(5,0.8)\), is \(E[X^{-1}]\) finite? Use the formula to calculate it.