Tutorials › Statistical Theory I - Probability and Distribution Theory › Pareto Type I: Support, Quantiles, and Tail Probabilities

Normal and Tail Families · Tutorial 186 of 994

Pareto Type I: Support, Quantiles, and Tail Probabilities

Learn to translate Pareto Type I parameters into distribution and tail calculations, including quantile-based risk thresholds.

Intermediate 10 min read

What You'll Learn

  • Define the Pareto Type I distribution by its minimum value and positive tail index.
  • Derive its cumulative distribution function and probability density.
  • Invert the CDF to find quantiles and upper-tail risk thresholds.
  • Calculate exceedance probabilities at and above the support boundary.
  • Use the survival function to evaluate conditional exceedance probabilities.
  • Check parameterization, endpoint, and percentage-versus-probability conventions.

Introduction

The Pareto Type I distribution describes positive quantities whose possible values begin at a fixed minimum and whose upper tail decreases according to a power law. It is useful when the question is not only how likely a particular value is, but how likely an observation is to exceed a high threshold. Examples include modeling large losses, unusually large claims, or other positive measurements with a long upper tail.

Two parameters determine the law: a minimum value \(x_m\), which sets the support and measurement scale, and a positive tail index \(\alpha\), which controls how quickly exceedance probabilities decline. This tutorial derives the distribution function, density, quantiles, and tail probabilities, then applies them to threshold-risk calculations. It does not require calculating moments.

Definition and Support

Definition: A random variable \(X\) has a Pareto Type I distribution with minimum value \(x_m>0\) and tail index \(\alpha>0\) if its cumulative distribution function is $$ F(x)= \begin{cases} 0, & x<x_m,\\ 1-\left(\dfrac{x_m}{x}\right)^\alpha, & x\ge x_m. \end{cases} $$ We write \(X\sim\operatorname{ParetoI}(x_m,\alpha)\). Its support is \([x_m,\infty)\).

The support says that values below \(x_m\) are impossible, while values at or above it are possible. The distribution has no point mass at the boundary: substituting \(x=x_m\) gives \(F(x_m)=0\). The parameter \(x_m\) has the same units as \(X\); the index \(\alpha\) is dimensionless. Larger \(\alpha\) means faster tail decay when \(x_m\) is held fixed.

Deriving the Density and Survival Function

For \(x>x_m\), differentiate the CDF. On the region below \(x_m\), the density is zero. The result is the Pareto Type I density.

Density: For \(X\sim\operatorname{ParetoI}(x_m,\alpha)\), $$ f(x)= \begin{cases} \alpha x_m^\alpha x^{-(\alpha+1)}, & x\ge x_m,\\ 0, & x<x_m. \end{cases} $$

The density formula is consistent with the CDF because, for \(x>x_m\), $$ \frac{d}{dx}\left[1-x_m^\alpha x^{-\alpha}\right] =\alpha x_m^\alpha x^{-(\alpha+1)}. $$ At the single endpoint \(x_m\), the value assigned to a density does not affect probabilities; the displayed convention gives the right-hand formula there. To check that the density integrates to one, use \(\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} =x_m^\alpha x_m^{-\alpha} =1. $$

The survival function, or upper-tail probability, is \(P(X>x)=1-F(x)\). Therefore,

$$ P(X>x)= \begin{cases} 1, & x<x_m,\\ \left(\dfrac{x_m}{x}\right)^\alpha, & x\ge x_m. \end{cases} $$

At \(x=x_m\), this probability is one: \(X\) is at least \(x_m\), and it has no atom at \(x_m\). Since the distribution is continuous, \(P(X>x)=P(X\ge x)\) for every fixed \(x\). This equality is useful in threshold calculations, but it would not hold automatically for a distribution with point masses.

Worked Example: CDF, Density, and an Interval Probability

Suppose \(X\sim\operatorname{ParetoI}(2,3)\), with \(X\) measured in the same units as \(x_m=2\). Find \(F(5)\), the density at \(5\), and \(P(3<X\le 5)\).

The CDF at 5 is $$ F(5)=1-\left(\frac{2}{5}\right)^3 =1-\frac{8}{125} =\frac{117}{125}=0.936. $$

The density at 5 is $$ f(5)=3(2)^3(5)^{-4} =\frac{24}{625} =0.0384. $$ A density is not itself a probability at a point; here it gives probability per unit of \(X\) near 5.

For the interval probability, subtract the two CDF values: $$ P(3<X\le 5)=F(5)-F(3) =\left(1-\frac{8}{125}\right)-\left(1-\frac{8}{27}\right) =\frac{8}{27}-\frac{8}{125} =\frac{784}{3375}\approx 0.2323. $$ As a check, the upper-tail probabilities at 3 and 5 are approximately \(0.2963\) and \(0.0640\); their difference is approximately \(0.2323\), as required.

Quantiles and Upper-Tail Risk Thresholds

A quantile reverses the CDF: it gives the value below which a specified fraction of the distribution lies. For \(0<p<1\), define the \(p\)-quantile by \(q_p=\inf\{x:F(x)\ge p\}\). Because the Pareto CDF is continuous and strictly increasing on \([x_m,\infty)\), this quantile is the unique value satisfying \(F(q_p)=p\).

Quantile formula: If \(X\sim\operatorname{ParetoI}(x_m,\alpha)\) and \(0<p<1\), then $$ q_p=x_m(1-p)^{-1/\alpha}. $$

To obtain the formula, solve the CDF equation: $$ p=1-\left(\frac{x_m}{q_p}\right)^\alpha \quad\Longrightarrow\quad \left(\frac{x_m}{q_p}\right)^\alpha=1-p \quad\Longrightarrow\quad q_p=x_m(1-p)^{-1/\alpha}. $$ The solution is at least \(x_m\), since \(0<1-p<1\) implies \((1-p)^{-1/\alpha}>1\).

At the boundary, the quantiles approach \(x_m\) as \(p\) decreases to zero, and increase without bound as \(p\) approaches one. Thus it is convenient to use \(q_0=x_m\) as the lower endpoint by continuous extension; there is no finite \(q_1\). For a high confidence level \(p\), the formula also gives a risk threshold: the probability of exceeding \(q_p\) is \(1-p\). For instance, \(q_{0.99}\) is the threshold exceeded with probability \(0.01\).

Worked Example: A 95th-Percentile Threshold

Let \(X\sim\operatorname{ParetoI}(2,3)\). Find the 95th percentile and verify its exceedance probability.

Using \(p=0.95\), $$ q_{0.95}=2(1-0.95)^{-1/3} =2(0.05)^{-1/3} =2(20)^{1/3} \approx 2(2.7144) \approx 5.4288. $$

The formula for the survival probability verifies the result: $$ P(X>q_{0.95}) =\left(\frac{2}{2(20)^{1/3}}\right)^3 =\left(\frac{1}{(20)^{1/3}}\right)^3 =\frac{1}{20} =0.05. $$ Thus about 5% of the distribution lies above this threshold. The 95th percentile is a value, not a probability: its units are the same as those of \(X\).

Exceedance Probabilities and Conditional Tails

For a threshold \(t\), the probability of exceeding it is particularly simple when \(t\) is within the support. If \(t\ge x_m\), it is \((x_m/t)^\alpha\); if \(t<x_m\), exceedance is certain. This makes it possible to answer either of two common questions directly: “What is the chance of exceeding this value?” and “How high must the threshold be to make exceedance this unlikely?”

The same survival formula gives a conditional exceedance probability. For \(t\ge u\ge x_m\), divide the probability of exceeding \(t\) by the probability of exceeding \(u\): $$ P(X>t\mid X>u) =\frac{P(X>t)}{P(X>u)} =\frac{(x_m/t)^\alpha}{(x_m/u)^\alpha} =\left(\frac{u}{t}\right)^\alpha. $$ The denominator is positive because \(P(X>u)=(x_m/u)^\alpha>0\). This calculation describes the chance of going past a higher threshold among observations already known to exceed \(u\).

Worked Example: Rare-Event Threshold and Conditional Exceedance

A loss \(X\), measured in thousands of dollars, follows \(\operatorname{ParetoI}(500,2.5)\). Find (a) the probability that a loss exceeds 1,000, (b) the threshold exceeded with probability 0.01, and (c) the chance that a loss exceeds 2,000 given that it already exceeds 1,000.

(a) Since 1,000 is above the minimum, $$ P(X>1000)=\left(\frac{500}{1000}\right)^{2.5} =(0.5)^{2.5} =(0.5)^2\sqrt{0.5} =0.25(0.7071068) \approx 0.1768. $$ So the exceedance probability is about 17.68%.

(b) A threshold exceeded with probability 0.01 is the 99th percentile. Using the quantile formula, $$ q_{0.99}=500(1-0.99)^{-1/2.5} =500(0.01)^{-0.4} =500(100)^{0.4} \approx 500(6.3096) \approx 3154.8. $$ The threshold is approximately 3,154.8 thousand dollars. Substitution checks the tail probability: $$ P(X>3154.8)\approx\left(\frac{500}{3154.8}\right)^{2.5}\approx 0.01. $$

(c) With \(t=2000\) and \(u=1000\), $$ P(X>2000\mid X>1000) =\left(\frac{1000}{2000}\right)^{2.5} =(0.5)^{2.5} \approx 0.1768. $$ This conditional probability is not \(P(X>2000)\): it is calculated among losses already known to exceed 1,000.

How the Parameters Shape Tail Calculations

The formulas clarify the separate roles of the two parameters. Multiplying \(x_m\) by a positive constant multiplies every quantile and threshold by that same constant; it changes the scale without changing the tail index. Holding the threshold as a multiple of the minimum, say \(t=cx_m\) for \(c\ge1\), gives $$ P(X>cx_m)=c^{-\alpha}. $$ For fixed \(c>1\), increasing \(\alpha\) decreases this probability because \(c^{-\alpha}\) becomes smaller. At a fixed percentile \(p\), the multiplier from the minimum is \((1-p)^{-1/\alpha}\), which likewise decreases as \(\alpha\) increases.

These statements concern tail probabilities and quantiles only. They do not imply that all summaries of the distribution exist or are finite; moment existence requires a separate calculation. Keeping those questions distinct prevents a valid threshold calculation from being mistaken for a statement about an average or another moment.

Common Mistakes

  • Using the wrong tail expression. The CDF is \(1-(x_m/x)^\alpha\), while the exceedance probability is \((x_m/x)^\alpha\) for \(x\ge x_m\). Check whether the question asks for probability below or above the threshold.
  • Reversing the quantile exponent. Solving the CDF gives \(q_p=x_m(1-p)^{-1/\alpha}\), not \(x_m(1-p)^{1/\alpha}\). A high percentile must exceed \(x_m\), which the correct expression guarantees.
  • Forgetting the support case. The power formula for the survival probability applies when \(t\ge x_m\). For \(t<x_m\), \(P(X>t)=1\), since all possible values exceed \(t\).
  • Confusing percent and probability. A 1% exceedance probability is \(0.01\), so its threshold is \(q_{0.99}\), not \(q_{0.01}\).
  • Confusing conditional and unconditional risk. \(P(X>t\mid X>u)\) is the fraction exceeding \(t\) among observations already exceeding \(u\); it is not the unconditional probability \(P(X>t)\).
Key takeaway: For a Pareto Type I variable, the essential calculations are governed by \(P(X>x)=(x_m/x)^\alpha\) for \(x\ge x_m\) and \(q_p=x_m(1-p)^{-1/\alpha}\). One formula moves from a threshold to its risk; the other moves from a risk level to its threshold.

Check Your Understanding

Use the support, CDF, quantile, and survival formulas to answer these questions.

  1. For \(X\sim\operatorname{ParetoI}(4,2)\), find \(F(8)\) and \(P(X>8)\). Verify that the two probabilities sum to one.
  2. For \(X\sim\operatorname{ParetoI}(3,1.5)\), write the formula for the 90th percentile. What exceedance probability does that threshold have?
  3. If \(X\sim\operatorname{ParetoI}(10,3)\), what is \(P(X>5)\), and why does the answer differ from applying the power formula at 5?
  4. For \(t\ge u\ge x_m\), derive \(P(X>t\mid X>u)\) from the survival function.
  5. For a fixed \(c>1\), explain how \(P(X>cx_m)\) changes when \(\alpha\) increases.