Unit 8, lesson 53

Markov and Chebyshev Inequalities

How Markov's inequality bounds an upper tail of a nonnegative variable by its mean, how Chebyshev's inequality bounds two-sided deviations by the variance, when these bounds are sharp, and why a bound is not an estimate.

  • Video 53
  • Averages and Probability Bounds

A delay is never negative and averages 22 minutes, and nothing else about its distribution is known. Could the probability of a delay of 1010 minutes or more be 0.50.5? It could not, and the reason takes one line. This lesson derives Markov’s and Chebyshev’s inequalities, which turn a mean or a variance into a guaranteed ceiling on a tail probability, states the conditions each one needs, and shows that some distributions reach these ceilings exactly.

Key terms

Probability bound

Definition An upper bound for the probability of an event AA is a number bb with P(A)≤bP(A) \le b that is proved from partial information about the model, such as the possible values, the mean or the variance, and therefore holds for every distribution with that information. A lower bound P(A)≥bP(A) \ge b is defined the same way. A tail probability is the probability of an event such as {X≥a}\{X \ge a\} or {∣X−μ∣≥d}\{\lvert X - \mu \rvert \ge d\}, in which a variable is large or far from its mean. Every probability is at most 11, so an upper bound b≥1b \ge 1 is true but carries no information, and the useful statement is P(A)≤min⁡(1,b)P(A) \le \min(1, b).

In plain words A bound is a ceiling, not an estimate. It promises that the probability is no larger, while the actual probability can be far smaller, even 00. Distributions that share the same partial information can have very different tail probabilities under the same ceiling.

Example Here are three nonnegative delay models, each with mean 22 minutes. The first always gives exactly 22 minutes. The second gives 00 with probability 0.90.9 and 2020 with probability 0.10.1. The third gives 00 with probability 0.80.8 and 1010 with probability 0.20.2. Their probabilities of a delay of at least 1010 minutes are 00, 0.10.1 and 0.20.2. The next card shows that this probability is at most 0.20.2 for every nonnegative variable with mean 22. For the threshold 11 the same inequality gives the ratio 2/1=22/1 = 2, which only says that the probability is at most 11.

Markov’s inequality

Definition Let XX be a random variable with P(X≥0)=1P(X \ge 0) = 1 and finite expectation, and let a>0a > 0. Then

P(X≥a)≤E[X]a.P(X \ge a) \le \frac{E[X]}{a}.

To prove it, let II be the indicator of the event {X≥a}\{X \ge a\}. With probability 11 the inequality X≥aIX \ge aI holds, because X≥aX \ge a where I=1I = 1 and X≥0X \ge 0 where I=0I = 0. Taking expectations gives E[X]≥a E[I]=a P(X≥a)E[X] \ge a\,E[I] = a\,P(X \ge a), and dividing by a>0a > 0 finishes the proof. For a discrete XX the same argument reads

E[X]=∑xx pX(x)≥∑x≥ax pX(x)≥a∑x≥apX(x)=a P(X≥a).\begin{aligned} E[X] &= \sum_x x \, p_X(x) \\[0.4em] &\ge \sum_{x \ge a} x \, p_X(x) \\[0.4em] &\ge a \sum_{x \ge a} p_X(x) \\[0.4em] &= a\,P(X \ge a). \end{aligned}

The first inequality drops the terms with x<ax < a, which are nonnegative because X≥0X \ge 0, and the second replaces each remaining value by the smaller or equal number aa. Since {X>a}\{X > a\} is contained in {X≥a}\{X \ge a\}, the same bound holds for P(X>a)P(X > a).

In plain words If more than E[X]/aE[X]/a of the probability sat at aa or above, those values alone would contribute more than E[X]E[X] to the mean. Nonnegativity matters because negative values could pull the mean back down and hide a heavy upper tail. The threshold has to be positive so that dividing by it is possible and keeps the direction of the inequality.

Example A nonnegative delay with mean 22 minutes has P(X≥10)≤2/10=0.2P(X \ge 10) \le 2/10 = 0.2, and the minutes cancel. Without nonnegativity the argument fails. If X=−10X = -10 or X=10X = 10 with probability 12\tfrac12 each, the mean is 00, yet P(X≥10)=12P(X \ge 10) = \tfrac12, larger than the value 0/10=00/10 = 0 the formula would give. The negative values canceled the positive ones in the mean.

Two-sided deviation event

Definition Let YY have finite mean μ\mu and let d>0d > 0. The two-sided deviation event at distance dd is {∣Y−μ∣≥d}\{\lvert Y - \mu \rvert \ge d\}. It is the union of the lower tail {Y≤μ−d}\{Y \le \mu - d\} and the upper tail {Y≥μ+d}\{Y \ge \mu + d\}, and both boundary points μ−d\mu - d and μ+d\mu + d belong to it. Its complement is the event {∣Y−μ∣<d}\{\lvert Y - \mu \rvert < d\}, which is {μ−d<Y<μ+d}\{\mu - d < Y < \mu + d\} with strict inequalities. The strict event {∣Y−μ∣>d}\{\lvert Y - \mu \rvert > d\} is contained in the two-sided event and has smaller probability when the boundary points carry probability. Because both sides are nonnegative, ∣Y−μ∣≥d\lvert Y - \mu \rvert \ge d holds exactly when (Y−μ)2≥d2(Y - \mu)^2 \ge d^2. A one-sided event such as {Y≥μ+d}\{Y \ge \mu + d\} is contained in the two-sided event, so its probability is no larger.

In plain words A two-sided event counts misses in either direction. For a discrete variable, whose values can sit exactly on a boundary, it pays to write out both tails, check which endpoints are included and take complements with care.

Example Let YY be the roll of a fair die, so μ=3.5\mu = 3.5, and take d=2.5d = 2.5. The event ∣Y−3.5∣≥2.5\lvert Y - 3.5 \rvert \ge 2.5 is {Y≤1}∪{Y≥6}={1,6}\{Y \le 1\} \cup \{Y \ge 6\} = \{1, 6\}, with probability 13\tfrac13. Its complement {1<Y<6}={2,3,4,5}\{1 < Y < 6\} = \{2, 3, 4, 5\} has probability 23\tfrac23. The strict event ∣Y−3.5∣>2.5\lvert Y - 3.5 \rvert > 2.5 contains no possible roll and has probability 00.

Chebyshev’s inequality

Definition Let YY have finite mean μ\mu and finite variance σ2\sigma^2, and let d>0d > 0. Then

P(∣Y−μ∣≥d)≤σ2d2,P(\lvert Y - \mu \rvert \ge d) \le \frac{\sigma^2}{d^2},

and taking complements,

P(∣Y−μ∣<d)≥1−σ2d2.P(\lvert Y - \mu \rvert < d) \ge 1 - \frac{\sigma^2}{d^2}.

To prove it, note that (Y−μ)2(Y - \mu)^2 is nonnegative with expectation σ2\sigma^2, that d2>0d^2 > 0, and that {∣Y−μ∣≥d}={(Y−μ)2≥d2}\{\lvert Y - \mu \rvert \ge d\} = \{(Y - \mu)^2 \ge d^2\}. Markov’s inequality applied to (Y−μ)2(Y - \mu)^2 with threshold d2d^2 gives the bound.

In plain words A small variance leaves little room for probability far from the mean, because probability at distance dd or more adds at least d2d^2 times its size to the variance. The variance goes in the numerator, not the standard deviation, and the squared units cancel. The inequality needs no symmetry, no bell shape and no continuity.

Example Let YY have mean 5050 units and variance 44 square units. Then P(∣Y−50∣≥6)≤4/36=1/9P(\lvert Y - 50 \rvert \ge 6) \le 4/36 = 1/9, so P(44<Y<56)≥8/9P(44 < Y < 56) \ge 8/9. The one-sided event {Y≥56}\{Y \ge 56\} is part of the two-sided event, so its probability is also at most 1/91/9. Putting the standard deviation 22 in the numerator would give 2/36=1/182/36 = 1/18, which is not a valid bound, since the sharp example below has mean 5050, variance 44 and deviation probability exactly 1/91/9. For the fair die of the previous card, σ2=35/12\sigma^2 = 35/12 and the bound is (35/12)/2.52=7/15≈0.467(35/12)/2.5^2 = 7/15 \approx 0.467, above the actual probability 13\tfrac13.

Bound in standard deviation units

Definition Let YY have finite mean μ\mu and finite standard deviation σ>0\sigma > 0, and let k>0k > 0. Taking d=kσd = k\sigma in Chebyshev’s inequality gives

P(∣Y−μ∣≥kσ)≤1k2.P(\lvert Y - \mu \rvert \ge k\sigma) \le \frac{1}{k^2}.

Equivalently, at least 1−1/k21 - 1/k^2 of the probability lies strictly within kk standard deviations of the mean. The bound is below 11 exactly when k>1k > 1, and for k≤1k \le 1 it says nothing. If σ=0\sigma = 0, then P(Y=μ)=1P(Y = \mu) = 1 and the deviation event has probability 00 for every d>0d > 0.

In plain words Measure the distance from the mean in standard deviations and the ceiling depends on that count alone. For every distribution with finite positive variance, at most a quarter of the probability lies 22 or more standard deviations from the mean, and at most a ninth lies 33 or more away. These are worst cases, not the percentages that hold for a normal distribution.

Example At k=2k = 2 the bound is 14\tfrac14 and at k=3k = 3 it is 19\tfrac19. For a normal distribution the actual probabilities are about 0.0460.046 and 0.00270.0027, far below these ceilings. For YY with mean 5050 and standard deviation 22, the distance 66 is k=3k = 3 standard deviations, which matches the bound 19\tfrac19 found before. The distance 11 is k=12k = \tfrac12, and the bound 44 gives no information.

Sharp bound

Definition An inequality that holds for every distribution in a class is sharp for that class if some distribution in the class attains it with equality. No smaller number can then serve as a bound based on the same information. Markov’s inequality is sharp. For a mean m≥0m \ge 0 and a threshold a>0a > 0 with a≥ma \ge m, the variable equal to aa with probability m/am/a and to 00 otherwise is nonnegative, has mean mm and has P(X≥a)=m/aP(X \ge a) = m/a. Chebyshev’s inequality is sharp too. For any μ\mu, any σ≥0\sigma \ge 0 and any d>0d > 0 with d≥σd \ge \sigma, put probability σ2/(2d2)\sigma^2/(2d^2) at each of μ−d\mu - d and μ+d\mu + d and the remaining 1−σ2/d21 - \sigma^2/d^2 at μ\mu. This distribution has mean μ\mu, variance σ2\sigma^2 and

P(∣Y−μ∣≥d)=σ2d2.P(\lvert Y - \mu \rvert \ge d) = \frac{\sigma^2}{d^2}.

When a<ma < m or d<σd < \sigma the ratio exceeds 11, and the capped value 11 is attained by the constant mm or by the values μ±σ\mu \pm \sigma with probability 12\tfrac12 each.

In plain words Sharp means the worst case really happens, so nothing better can be promised from the same information. It does not mean the bound is close for the distribution at hand. A better answer for a particular model needs more information about it.

Example The delay equal to 00 with probability 0.80.8 and 1010 with probability 0.20.2 is nonnegative with mean 22 and has P(X≥10)=0.2P(X \ge 10) = 0.2, exactly Markov’s bound. The variable equal to 4444, 5050 or 5656 with probabilities 118\tfrac{1}{18}, 89\tfrac89 and 118\tfrac{1}{18} has mean 5050 and variance 36⋅218=436 \cdot \tfrac{2}{18} = 4, and P(∣Y−50∣≥6)=19P(\lvert Y - 50 \rvert \ge 6) = \tfrac19, exactly Chebyshev’s bound. In both examples all the tail probability sits exactly on the boundary, so the strict events {X>10}\{X > 10\} and {∣Y−50∣>6}\{\lvert Y - 50 \rvert > 6\} have probability 00.

Check yourself

  1. A cost CC is nonnegative with mean 1212. Bound P(C≥40)P(C \ge 40). What does Markov’s inequality say about P(C≥6)P(C \ge 6)?
Show the answer

Markov’s inequality gives P(C≥40)≤12/40=0.3P(C \ge 40) \le 12/40 = 0.3. This is a ceiling, and the cost equal to 4040 with probability 0.30.3 and 00 otherwise reaches it. For the threshold 66 the ratio is 12/6=212/6 = 2, which is true but says no more than P(C≥6)≤1P(C \ge 6) \le 1. Indeed a cost that always equals 1212 has P(C≥6)=1P(C \ge 6) = 1.

  1. A measurement MM has mean 100100 and variance 99. Bound P(∣M−100∣≥6)P(\lvert M - 100 \rvert \ge 6) and give a lower bound for P(94<M<106)P(94 < M < 106). A classmate answers 3/363/36 for the first part. What went wrong?
Show the answer

Chebyshev’s inequality gives P(∣M−100∣≥6)≤9/36=14P(\lvert M - 100 \rvert \ge 6) \le 9/36 = \tfrac14. The standard deviation is 33, so the distance 66 is k=2k = 2 standard deviations and 1/k2=141/k^2 = \tfrac14 agrees. The complement gives P(94<M<106)≥34P(94 < M < 106) \ge \tfrac34, with strict inequalities because the bounded event includes 9494 and 106106. The classmate put the standard deviation 33 in the numerator instead of the variance 99. The value 3/36=1123/36 = \tfrac{1}{12} is not a valid bound, since the measurement equal to 9494 or 106106 with probability 18\tfrac18 each and to 100100 with probability 34\tfrac34 has mean 100100, variance 36⋅14=936 \cdot \tfrac14 = 9 and deviation probability exactly 14\tfrac14.

  1. A daily profit GG, in hundreds of dollars, can be negative and has mean 22. Does Markov’s inequality give P(G≥10)≤0.2P(G \ge 10) \le 0.2? If in addition G≥−6G \ge -6 on every outcome, what bound follows, and can it be improved with only this information?
Show the answer

No, because GG is not nonnegative. The profit equal to 1010 or −6-6 with probability 12\tfrac12 each has mean 22 but P(G≥10)=12P(G \ge 10) = \tfrac12, larger than 0.20.2. With G≥−6G \ge -6, the shifted variable H=G+6H = G + 6 is nonnegative with mean 88, and

P(G≥10)=P(H≥16)≤816=12.\begin{aligned} P(G \ge 10) &= P(H \ge 16) \\[0.4em] &\le \frac{8}{16} = \frac12 . \end{aligned}

The profit just described satisfies G≥−6G \ge -6, has mean 22 and attains 12\tfrac12, so this bound is sharp and cannot be lowered with only these facts.

How this lesson connects