Tutorials › Statistical Theory I - Probability and Distribution Theory › Worked Tail Comparisons for Discrete and Continuous Laws

Stochastic Orders and Comparisons · Tutorial 949 of 994

Worked Tail Comparisons for Discrete and Continuous Laws

Work through explicit tail and expectation calculations to determine when these discrete and continuous distributions are stochastically ordered.

Advanced 10 min read

What You'll Learn

  • Calculate CDFs and survival probabilities for Bernoulli, binomial, and exponential laws.
  • Compare Bernoulli distributions by their success probabilities.
  • Compare binomial tails and increasing-function expectations when the success probability changes.
  • Identify why equal means alone cannot establish stochastic order.
  • Compare exponential rates, including the equality of survival probabilities at zero.
  • Compute expectations of selected increasing functions directly from the distributions.

What a Tail Comparison Tells Us

A comparison of two distributions is often most useful when stated in terms of thresholds: how likely is each random variable to exceed a given value? For random variables \(X\) and \(Y\), the usual stochastic order \(X\le_{\mathrm{st}}Y\) means that \(P(X>t)\le P(Y>t)\) for every real \(t\). Equivalently, their CDFs satisfy \(F_X(t)\ge F_Y(t)\) for every \(t\). This tutorial develops that comparison by calculation for Bernoulli, binomial, and exponential laws.

A tail comparison also predicts the direction of expectations of increasing functions, provided the expectations exist. The examples below make that connection concrete. They also show why checking a mean, a single threshold, or a single endpoint is not enough to establish an order.

Definition: The survival probability of \(X\) at threshold \(t\) is \(P(X>t)\). For the usual stochastic order, \(X\le_{\mathrm{st}}Y\) exactly when \(P(X>t)\le P(Y>t)\) for all real \(t\), or, equivalently, \(F_X(t)\ge F_Y(t)\) for all real \(t\). A strict inequality at some thresholds does not replace the requirement that the inequality hold at every threshold.

Bernoulli Distributions: One Tail Determines the Comparison

If \(X\sim\operatorname{Bernoulli}(p)\), then \(P(X=1)=p\) and \(P(X=0)=1-p\), where \(0\le p\le1\). Its CDF and survival probability are

$$ F_X(t)= \begin{cases} 0,&t<0,\\ 1-p,&0\le t<1,\\ 1,&t\ge1, \end{cases} \qquad P(X>t)= \begin{cases} 1,&t<0,\\ p,&0\le t<1,\\ 0,&t\ge1. \end{cases} $$

Thus, for two Bernoulli variables with success probabilities \(p_1\) and \(p_2\), the only interval on which their survival probabilities can differ is \(0\le t<1\). It follows directly that \(\operatorname{Bernoulli}(p_1)\le_{\mathrm{st}}\operatorname{Bernoulli}(p_2)\) exactly when \(p_1\le p_2\).

Worked Example: Comparing Two Bernoulli Laws

Let \(X\sim\operatorname{Bernoulli}(0.2)\) and \(Y\sim\operatorname{Bernoulli}(0.7)\). For \(0\le t<1\), \(P(X>t)=0.2\) and \(P(Y>t)=0.7\). For \(t<0\), both survival probabilities equal 1; for \(t\ge1\), both equal 0. Therefore \(X\le_{\mathrm{st}}Y\).

For an increasing function \(g\) on \(\{0,1\}\), monotonicity means \(g(1)\ge g(0)\). Its expectation under a Bernoulli law is

$$ E[g(X)]=(1-p)g(0)+pg(1) =g(0)+p\{g(1)-g(0)\}. $$

For example, take \(g(x)=1+2x\), which is increasing on the support. Then \(E[g(X)]=0.8(1)+0.2(3)=1.4\), whereas \(E[g(Y)]=0.3(1)+0.7(3)=2.4\). The threshold comparison and this increasing-function expectation both point in the same direction.

Binomial Distributions: Comparing Every Integer Threshold

For \(X\sim\operatorname{Binomial}(n,p)\), the probability mass at \(j\in\{0,\ldots,n\}\) is \(\binom{n}{j}p^j(1-p)^{n-j}\). Since \(X\) is integer-valued, its survival probability changes only at integer thresholds: for integer \(k\), \(P(X\ge k)\) is the tail at the threshold just below \(k\). Comparing these tails for \(k=1,\ldots,n\), together with the trivial tails outside the support, suffices to compare two binomial laws with the same \(n\).

The following derivative calculation verifies that, for fixed \(n\), every such tail is nondecreasing in \(p\). For \(1\le k\le n\), write

$$ S_k(p)=P(X\ge k)=\sum_{j=k}^{n}\binom{n}{j}p^j(1-p)^{n-j}. $$

For \(0<p<1\), differentiate term by term; the sum is finite. Using \(j\binom{n}{j}=n\binom{n-1}{j-1}\) and \((n-j)\binom{n}{j}=n\binom{n-1}{j}\), the positive terms cancel the negative terms from adjacent indices, leaving

$$ \begin{aligned} S_k'(p) &=n\sum_{j=k}^{n}\binom{n-1}{j-1}p^{j-1}(1-p)^{n-j}\\ &\quad-n\sum_{j=k}^{n-1}\binom{n-1}{j}p^j(1-p)^{n-j-1}\\ &=n\binom{n-1}{k-1}p^{k-1}(1-p)^{n-k}\ge0. \end{aligned} $$

Continuity of \(S_k(p)\), which is a polynomial, extends the nondecreasing property to \(p=0\) and \(p=1\). Therefore, with \(n\) fixed and \(p_1\le p_2\), \(\operatorname{Binomial}(n,p_1)\le_{\mathrm{st}}\operatorname{Binomial}(n,p_2)\). For \(0<p_1<p_2<1\), the derivative is positive for every \(1\le k\le n\), so each of these integer-threshold tails is strictly larger at \(p_2\).

Worked Example: Binomial Tail Tables and Increasing Expectations

Compare \(X\sim\operatorname{Binomial}(4,0.25)\) and \(Y\sim\operatorname{Binomial}(4,0.5)\). The probability masses, listed in order from zero to four, are

$$ P(X=j)=\frac{(81,108,54,12,1)_j}{256}, \qquad P(Y=j)=\frac{(1,4,6,4,1)_j}{16}. $$

For instance, \(P(X=2)=\binom42(0.25)^2(0.75)^2=54/256\), and \(P(Y=2)=\binom42(0.5)^4=6/16\). Summing the masses from each threshold upward gives

$$ \begin{array}{c|cccc} k&1&2&3&4\\ \hline P(X\ge k)&175/256&67/256&13/256&1/256\\ P(Y\ge k)&15/16&11/16&5/16&1/16 \end{array} $$

Each tail for \(Y\) is at least the corresponding tail for \(X\): for example, at \(k=2\), \(11/16=176/256>67/256\), and at \(k=4\), \(1/16=16/256>1/256\). Together with the equal trivial tails outside the support, this verifies \(X\le_{\mathrm{st}}Y\).

Here are two increasing-function expectations. First, \(g(x)=x\) is increasing, and \(E[X]=4(0.25)=1\), while \(E[Y]=4(0.5)=2\). Second, \(g(x)=x^2\) is increasing on \(\{0,1,2,3,4\}\). Since a binomial variable has variance \(np(1-p)\), \(E[X^2]=\operatorname{Var}(X)+(E[X])^2=4(0.25)(0.75)+1^2=1.75\), while \(E[Y^2]=4(0.5)(0.5)+2^2=5\).

The qualification “with \(n\) fixed” matters. Equal means alone do not establish a stochastic order. For example, let \(U\sim\operatorname{Binomial}(1,0.5)\) and \(V\sim\operatorname{Binomial}(2,0.25)\). Both have mean \(0.5\), but

$$ P(U\ge1)=\frac12>\frac7{16}=P(V\ge1), \qquad P(U\ge2)=0<\frac1{16}=P(V\ge2). $$

The tail inequalities go in opposite directions, so neither law is stochastically larger than the other. Each calculation is necessary: checking only one threshold, or only the mean, would miss this failure of ordering.

Increasing-Function Expectations from Integer Tails

For a nonnegative integer-valued random variable, tail probabilities also give a direct formula for expectations. Let \(g\) be nondecreasing on the nonnegative integers, with finite real values, and define \(\Delta g(k)=g(k)-g(k-1)\ge0\) for \(k\ge1\). For each nonnegative integer \(m\), telescoping gives \(g(m)=g(0)+\sum_{k=1}^{m}\Delta g(k)\). Applying this identity to \(X\) yields

$$ g(X)=g(0)+\sum_{k=1}^{\infty}\Delta g(k)\mathbf{1}_{\{X\ge k\}}. $$

Every summand after \(g(0)\) is nonnegative. Tonelli's theorem therefore permits interchanging the expectation and the infinite sum, even if the resulting expectation is infinite:

$$ E[g(X)]=g(0)+\sum_{k=1}^{\infty}\Delta g(k)P(X\ge k). $$

In particular, if \(X\le_{\mathrm{st}}Y\) and both are nonnegative integer-valued, each tail term for \(X\) is no larger than the corresponding term for \(Y\). The formula then gives \(E[g(X)]\le E[g(Y)]\) for every such nondecreasing \(g\), with the expectations understood as extended expectations when necessary. For a finite-support variable, the sum has only finitely many nonzero terms.

Exponential Distributions: Rate Changes and Continuous Tails

Let \(X\sim\operatorname{Exponential}(\lambda)\), where \(\lambda>0\) is the rate. Its density is \(f_X(x)=\lambda e^{-\lambda x}\) for \(x\ge0\), and zero otherwise. Integrating the density gives, for \(t\ge0\),

$$ F_X(t)=1-e^{-\lambda t}, \qquad P(X>t)=e^{-\lambda t}. $$

For \(t<0\), \(F_X(t)=0\) and \(P(X>t)=1\), because \(X\) is nonnegative. If \(0<\lambda_1<\lambda_2\), then \(e^{-\lambda_2t}\le e^{-\lambda_1t}\) for every \(t\ge0\). Equality holds at \(t=0\), since both survival probabilities are 1 there; the inequality is strict for every \(t>0\). For \(t<0\), both probabilities equal 1. Thus the larger rate has no larger survival probability at every nonnegative threshold, and a strictly smaller survival probability at every positive threshold. In stochastic-order notation, \(\operatorname{Exponential}(\lambda_2)\le_{\mathrm{st}}\operatorname{Exponential}(\lambda_1)\).

Worked Example: Exponential Tails and Bounded Increasing Functions

Let \(X\sim\operatorname{Exponential}(2)\) and \(Y\sim\operatorname{Exponential}(1)\). At \(t=0\), both survival probabilities are \(e^0=1\). At \(t=1\), \(P(X>1)=e^{-2}\approx0.1353\), while \(P(Y>1)=e^{-1}\approx0.3679\). More generally, for every \(t>0\), \(e^{-2t}<e^{-t}\), and for \(t<0\) both survival probabilities are 1. Hence \(X\le_{\mathrm{st}}Y\).

To compare expectations of an increasing function, take \(g(x)=\min(x,1)\), which is bounded and nondecreasing for \(x\ge0\). Direct integration, split at 1, gives for a general rate \(\lambda>0\)

$$ \begin{aligned} E[\min(X,1)] &=\int_0^1 x\lambda e^{-\lambda x}\,dx+\int_1^\infty \lambda e^{-\lambda x}\,dx\\ &=\left[-xe^{-\lambda x}-\frac{e^{-\lambda x}}{\lambda}\right]_0^1+e^{-\lambda}\\ &=\frac{1-e^{-\lambda}}{\lambda}. \end{aligned} $$

Consequently, \(E[\min(X,1)]=(1-e^{-2})/2\approx0.4323\), while \(E[\min(Y,1)]=1-e^{-1}\approx0.6321\). The expectation for the larger-rate distribution is smaller, as its tail comparison predicts.

Common Mistakes and a Practical Check

  • Calling the exponential inequality strict at zero. At \(t=0\), every positive-rate exponential law has survival probability 1. The comparison for rates \(\lambda_2>\lambda_1\) is non-strict at \(t=0\) and strict only for \(t>0\).
  • Reversing the direction when comparing CDFs and tails. If \(X\le_{\mathrm{st}}Y\), then \(P(X>t)\le P(Y>t)\), but \(F_X(t)\ge F_Y(t)\). The CDF inequality points in the opposite direction.
  • Inferring an order from means or one threshold. Equal means, or a favorable comparison at one cutoff, says nothing about all other cutoffs. The binomial pair with equal means above has tails crossing at \(k=1\) and \(k=2\).
  • Using a function that is not increasing on the relevant support. The expectation comparison concerns increasing functions. For example, \(x^2\) is increasing on nonnegative values, but not on the whole real line.

For a discrete law, list the possible integer thresholds and calculate each tail; for an exponential law, use its survival formula and check the boundary \(t=0\) separately. After verifying the comparison at every threshold, expectations of selected increasing functions provide useful numerical checks, not substitutes for the tail calculation.

Check Your Understanding

Use the support and formulas for each distribution to justify each answer.

  1. For \(A\sim\operatorname{Bernoulli}(0.35)\) and \(B\sim\operatorname{Bernoulli}(0.6)\), compare their survival probabilities for \(t<0\), \(0\le t<1\), and \(t\ge1\).
  2. For fixed \(n\), why is the binomial tail \(P(X\ge k)\) nondecreasing in \(p\)? Identify the derivative formula and its range of validity.
  3. Compare \(P(X\ge2)\) for \(X\sim\operatorname{Binomial}(3,0.2)\) and \(Y\sim\operatorname{Binomial}(3,0.6)\), showing the calculation for both.
  4. For exponential rates \(\lambda_1<\lambda_2\), state where their survival probabilities are equal and where the higher-rate law has strictly smaller survival.
  5. Explain why equal means do not imply usual stochastic order, using a pair of distributions whose tails cross at two thresholds.