Concentration Inequalities, the Law of Large Numbers, and the Central Limit Theorem

Markov and Chebyshev inequalities as distribution-free tail bounds, and Chebyshev confidence intervals for a population proportion with sample-size calculations (Hajek §2.9). The Gaussian-approximation heuristic, the law of large numbers and the central limit theorem are pending.

· 6 min read · 1,030 words

Two inequalities that bound how much probability mass can sit far from the mean, using nothing but the mean (Markov) or the mean and variance (Chebyshev), plus the two limit theorems they enable: the law of large numbers (the sample average converges to the true mean — provable directly via Chebyshev) and the central limit theorem (properly normalized, the sample average converges in distribution to a standard normal). Chebyshev alone gives a cheap, distribution-free confidence interval for a sample average; the CLT gives a tighter one once its approximation is justified. Also covers the informal, Gaussian-approximation-heuristic version of the CLT (sums of many independent, comparably-sized random effects look Gaussian, which is why the normal approximation shows up everywhere in practice) ahead of — and as motivation for — the formal statement.

Markov inequality

TheoremMarkov inequality (Hajek §2.9)

If Y≥0Y\ge0 and c>0c > 0, P{Y≥c}≤E[Y]c.P\{Y\ge c\}\le\frac{E[Y]}c. Equality holds iff pY(0)+pY(c)=1p_Y(0) + p_Y(c) = 1.

ProofBound each value from below

Each possible value uiu_i is at least 0, and at least cc when ui≥cu_i\ge c. So E[Y]=∑iuipY(ui)≥∑i: ui≥cc pY(ui)=cP{Y≥c}.E[Y] = \sum_iu_ip_Y(u_i)\ge\sum_{i:\,u_i\ge c}c\,p_Y(u_i) = cP\{Y\ge c\}. (This is the discrete case. For continuous YY, replace the sum by an integral.)

∎

ExampleBalls in buckets (Hajek)

200 balls sit in 100 buckets in some unknown way. At most how many buckets can hold 5 or more? Pick a bucket at random and let YY be its ball count, so E[Y]=2E[Y] = 2. Markov gives P{Y≥5}≤25P\{Y\ge5\}\le\frac25, so at most 40 buckets. This is tight: put 5 balls in each of 40 buckets.

IntuitionOnly the mean is known

If too much mass sat far out, the mean would have to be larger. Markov is the best bound available from the mean alone, as the tight example shows.

Chebyshev inequality

TheoremChebyshev inequality

If XX has mean μ\mu and variance σ2\sigma^2, then for d>0d > 0 and a>0a > 0,

P{∣X−μ∣≥d}≤σ2d2,P{∣X−μ∣≥aσ}≤1a2.(2.12–2.13)P\{|X - \mu|\ge d\}\le\frac{\sigma^2}{d^2},\qquad P\{|X - \mu|\ge a\sigma\} \le\frac1{a^2}. \tag{2.12–2.13}
ProofMarkov applied to the squared deviation

Apply Markov to Y=(X−μ)2Y = (X - \mu)^2 with c=d2c = d^2: P{∣X−μ∣≥d}=P{Y≥d2}≤E[Y]/d2=σ2/d2P\{|X - \mu|\ge d\} = P\{Y\ge d^2\}\le E[Y]/d^2 = \sigma^2/d^2. Then take d=aσd = a\sigma.

∎

ImportantDistribution-free, therefore conservative

Chebyshev holds for every distribution with the given mean and variance, so for any particular one it is usually loose. At 2σ2\sigma it guarantees 75%, while a normal puts about 95% within 2σ2\sigma (continuous-random-variables-and-distributions).

Confidence intervals from concentration bounds

DefinitionPoint and interval estimators

To estimate a population fraction pp, sample nn at random and let XX be the number with the property, so p^=X/n\hat p = X/n is a point estimator. If the population is much larger than nn, XX is approximately binomial(n,p)(n,p), with σ=np(1−p)\sigma = \sqrt{np(1-p)}.

TheoremA Chebyshev confidence interval

For any a>0a > 0,

P{p∈(p^−a2n, p^+a2n)}≥1−1a2.(2.15)P\Big\{p\in\Big(\hat p - \frac a{2\sqrt n},\ \hat p + \frac a{2\sqrt n}\Big) \Big\}\ge1 - \frac1{a^2}. \tag{2.15}
ProofChebyshev, then remove the unknown pp

Step 1: Chebyshev on XX. P{∣X−np∣<aσ}≥1−1/a2P\{|X - np| < a\sigma\}\ge1 - 1/a^2. Divide by nn: P{∣p^−p∣<ap(1−p)/n}≥1−1/a2P\{|\hat p - p| < a\sqrt{p(1-p)/n}\}\ge1 - 1/a^2.

Step 2: The interval still depends on pp. Since p(1−p)≤14p(1-p)\le\frac14 (discrete-random-variables-and-distributions), replacing p(1−p)p(1-p) by 14\frac14 widens the interval to half-width a2n\frac a{2\sqrt n}, which can only raise the probability.

∎

FormulaConfidence levels

aa2510
1−1/a21 - 1/a^275%96%99%

WarningThe interval is random, not pp

“96% confidence” is a statement before sampling: the random interval will cover the fixed, unknown pp with probability at least 96%. After you see the data and compute a particular interval, it is not true that pp lies in it with probability 96%.

ExampleHow many phone numbers to test? (Hajek)

To estimate the fraction of busy numbers within 0.1 with 96% confidence, take a=5a = 5. The half-width 2.5n≤0.1\frac{2.5}{\sqrt n}\le0.1 needs n≥(2.5/0.1)2=625n\ge(2.5/0.1)^2 = 625.

RemarkChebyshev versus the normal approximation

Chebyshev intervals are conservative. Intervals from the normal approximation are narrower for the same confidence (continuous-random-variables-and-distributions). The formal justification is the central limit theorem, still to be written below.


Still to write

NoteTODO — The central limit theorem and the Gaussian approximation (informal heuristic)

Absorbed from gaussian-distribution-and-clt. Needs Hajek’s informal treatment; source text not yet in sources/. The DeMoivre–Laplace special case is already in continuous-random-variables-and-distributions.

NoteTODO — Law of large numbers

Needs Ross Ch. 8 and the corresponding Hajek material. Source text not yet in sources/.

NoteTODO — Central limit theorem (formal statement)

Needs Ross Ch. 8 and the corresponding Hajek material. Source text not yet in sources/.


Quick reference

ObjectStatementNote
MarkovP{Y≥c}≤E[Y]/cP\{Y\ge c\}\le E[Y]/c for Y≥0Y\ge0tight: mass only at 0 and cc
ChebyshevP{∣X−μ∣≥aσ}≤1/a2P\{\lvert X-\mu\rvert\ge a\sigma\}\le1/a^2Markov on (X−μ)2(X-\mu)^2
Proportion CIp^±a2n\hat p\pm\frac a{2\sqrt n} with confidence ≥1−1/a2\ge1 - 1/a^2uses p(1−p)≤14p(1-p)\le\frac14
Sample sizen≥(a/(2ϵ))2n\ge(a/(2\epsilon))^2 for half-width ϵ\epsilon625 for ϵ=.1\epsilon = .1, 96%

Facts to keep to hand

  • Tail bounds from moments alone. Markov needs only the mean, Chebyshev the mean and variance. Both hold for every distribution, which makes them robust and conservative.
  • Confidence belongs to the procedure, not the interval. The random interval covers the fixed pp with the stated probability before the data are seen.
  • Precision costs 1/ϵ21/\epsilon^2 samples. Halving the half-width quadruples the sample size, because the spread of p^\hat p shrinks like 1/n1/\sqrt n.

References

← All notes