Unit 5, lesson 32

The Central Limit Theorem

What the central limit theorem says about the standardized sample mean, why it needs independent observations with finite positive variance, and how heavy tails slow or break the normal approximation.

  • Video 32
  • Sampling Distributions

Waiting times are never negative, and a few long waits give their distribution a long right tail. Suppose a model says that a single wait is exponential with mean 2 minutes, and you want the probability that the mean of 100 independent waits lies between 1.6 and 2.4 minutes. The waits themselves are far from normal, yet the mean of 100 of them has a sampling distribution that is very close to a normal curve. The central limit theorem explains when and in what sense this happens, and this lesson states it with its conditions, uses it to approximate that probability and shows where it stops working.

Key terms

Central limit theorem

Definition Let X1,X2,…X_1, X_2, \dots be independent and identically distributed random variables with mean μ\mu and finite positive variance σ2\sigma^2, and let Xˉn\bar X_n be the mean of the first nn of them. Form the standardized sample mean

Zn=Xˉn−μσ/n.Z_n = \frac{\bar X_n - \mu}{\sigma / \sqrt{n}} .

The central limit theorem says that for every real number zz

lim⁡n→∞P(Zn≤z)=Φ(z),\lim_{n \to \infty} P(Z_n \le z) = \Phi(z),

where Φ\Phi is the cumulative distribution function of the standard normal distribution N(0,1)N(0, 1). In other words, the distribution of ZnZ_n converges to N(0,1)N(0, 1). The statement is about the distribution of the mean over repeated samples of size nn. It says nothing about the shape of the observations inside one sample, and as a limit it names no sample size at which the approximation becomes good. Independence is a real condition. When observations are positively correlated, as measurements taken in related groups often are, the true standard error of the mean is larger than σ/n\sigma/\sqrt{n}, and this version of the theorem does not apply. Other versions allow observations that are not identically distributed or are weakly dependent, but they replace these conditions with others.

In plain words Averages of many independent draws from one population have a sampling distribution close to a normal curve, whatever the shape of the population, as long as its variance is finite and positive. How large the sample must be depends on that shape. Averaging changes the distribution of the mean, not the distribution of the individual observations.

Example In the exponential model with mean 2 minutes, a single wait has SD 2 minutes and a long right tail. A common measure of skewness, the standardized third moment, equals 2 for one wait and 2/n2/\sqrt{n} for the mean of nn independent waits. At n=100n = 100 it is only 0.20.2, so the sampling distribution of the mean is nearly symmetric and close to normal. The 100 waits inside one sample still form a right-skewed histogram, and taking more waits makes that histogram look more like the exponential density, not more like a normal curve.

Standardized sample mean

Definition For the mean Xˉn\bar X_n of nn independent observations with mean μ\mu and finite positive SD σ\sigma, the standardized sample mean is ZnZ_n above. It measures the deviation of Xˉn\bar X_n from μ\mu in units of the standard error

SE⁡(Xˉn)=σn.\operatorname{SE}(\bar X_n) = \frac{\sigma}{\sqrt{n}} .

Since E(Xˉn)=μE(\bar X_n) = \mu and Var⁡(Xˉn)=σ2/n\operatorname{Var}(\bar X_n) = \sigma^2/n, the standardized mean has mean 00 and variance 11 for every nn. With the sum Sn=X1+⋯+XnS_n = X_1 + \dots + X_n the same quantity is

Zn=Sn−nμσn.Z_n = \frac{S_n - n\mu}{\sigma \sqrt{n}} .

Standardizing is needed for a limit with a shape. By the law of large numbers the distribution of Xˉn\bar X_n itself piles up at μ\mu as nn grows, so its limit is a single point. The standardized mean keeps its spread at 11, which lets its shape be compared with one fixed curve. The law of large numbers says where the mean settles, and the central limit theorem describes the shape and size of its fluctuations around that point. Neither guarantees that one particular sample mean is close to μ\mu.

In plain words The standardized mean says how many standard errors the sample mean lies above or below the population mean. It is a z-score for the sample mean, computed with the SE of the mean, not with the SD of single observations.

Example For the waits with μ=2\mu = 2 and σ=2\sigma = 2 minutes, a sample of n=100n = 100 gives SE⁡(Xˉ100)=2/100=0.2\operatorname{SE}(\bar X_{100}) = 2/\sqrt{100} = 0.2 minute, and a sample mean of 2.42.4 minutes has standardized value

2.4−20.2=2.\frac{2.4 - 2}{0.2} = 2 .

With n=400n = 400 the standard error is 2/400=0.12/\sqrt{400} = 0.1 minute, and a sample mean of only 2.22.2 minutes has the same standardized value 22. The raw means crowd closer to 22 as nn grows, while the standardized values stay on one common scale.

Finite positive variance

Definition The condition 0<σ2<∞0 < \sigma^2 < \infty on the population variance. A finite variance means E(X2)<∞E(X^2) < \infty, which also makes the mean finite. The condition rules out two kinds of populations. If σ2=0\sigma^2 = 0, every observation equals μ\mu with probability 11, so Xˉn=μ\bar X_n = \mu exactly and the standardized mean is undefined. If the variance is infinite, there is no σ\sigma to standardize with, and the fluctuations of the mean need not be approximately normal at any sample size.

In plain words The population must actually vary, and its spread must be a finite number. Very heavy tails can make the variance infinite, and then this theorem gives no answer.

Example The standard Cauchy distribution has density

f(x)=1π(1+x2).f(x) = \frac{1}{\pi (1 + x^2)} .

It has no finite mean, and E(X2)E(X^2) is infinite. The mean of nn independent standard Cauchy observations has the same standard Cauchy distribution for every nn, so the probability that it is farther than 1010 from 00 is about 0.06350.0635 whether nn is 11 or 10001000. Averaging does not concentrate the mean at all. A milder case is the Pareto distribution with P(X>x)=x−1.5P(X > x) = x^{-1.5} for x≥1x \ge 1. Its mean is 33 and its variance is infinite, so the law of large numbers still applies to it, but this version of the central limit theorem does not.

Normal approximation for a sample mean

Definition Under the conditions of the theorem and for a fixed, large nn, the normal approximation treats Xˉn\bar X_n as approximately N(μ,σ2/n)N(\mu, \sigma^2/n), so that

P(a≤Xˉn≤b)≈Φ(zb)−Φ(za),\begin{aligned} &P(a \le \bar X_n \le b) \\[0.4em] &\quad \approx \Phi(z_b) - \Phi(z_a), \end{aligned}

where za=(a−μ)/(σ/n)z_a = (a - \mu)/(\sigma/\sqrt{n}) and zb=(b−μ)/(σ/n)z_b = (b - \mu)/(\sigma/\sqrt{n}) are the standardized endpoints. Its accuracy at a given nn depends on the population, especially its skewness and tails, and on the event. A central interval can be approximated very well while a far-tail probability at the same nn is off by a factor of two or more. No single sample size, such as 30, guarantees a good approximation for every population. When the population itself is normal and the observations are independent, Xˉn\bar X_n has the distribution N(μ,σ2/n)N(\mu, \sigma^2/n) exactly for every n≥1n \ge 1, so no limit is needed. When σ\sigma is unknown and replaced by the sample SD ss, the calculation changes, and later lessons treat that case.

In plain words Convert each endpoint into standard errors from μ\mu, read the area under the standard normal curve, and label the answer as an approximation.

Example For the mean of 100 waits, the endpoints 1.61.6 and 2.42.4 standardize to (1.6−2)/0.2=−2(1.6 - 2)/0.2 = -2 and (2.4−2)/0.2=2(2.4 - 2)/0.2 = 2, so

P(1.6≤Xˉ100≤2.4)≈P(−2≤Z≤2)≈0.9545.\begin{aligned} &P(1.6 \le \bar X_{100} \le 2.4) \\[0.4em] &\quad \approx P(-2 \le Z \le 2) \\[0.4em] &\quad \approx 0.9545 . \end{aligned}

In this exponential model the sum of 100 waits has a gamma distribution, and a numerical calculation from it gives the exact probability as about 0.95500.9550. The approximation is off by about 0.00050.0005 for this model, this nn and this interval, which does not establish its accuracy anywhere else. The tails show the difference. The normal curve gives 0.001350.00135 for both P(Xˉ100≥2.6)P(\bar X_{100} \ge 2.6) and P(Xˉ100≤1.4)P(\bar X_{100} \le 1.4), while the exact values are about 0.002750.00275 and 0.000430.00043. The right skew of the waits survives in the tails of the mean.

Heavy-tailed distribution

Definition Informally, a distribution is heavy-tailed when values far from its center keep much more probability than a normal distribution would give them. A common precise version for the right tail requires that, for every t>0t > 0, the product etxP(X>x)e^{tx} P(X > x) is unbounded as xx grows. Equivalently, there are no constants CC and t>0t > 0 for which

P(X>x)≤Ce−txP(X > x) \le C e^{-tx}

holds for every xx, so no exponential function bounds the tail probability. In this sense the exponential distribution is not heavy-tailed, although its right tail is heavier than a normal tail. Power-law tails, such as the Pareto tail P(X>x)=x−αP(X > x) = x^{-\alpha} for x≥1x \ge 1 with α>0\alpha > 0, are the standard example. This Pareto distribution has a finite mean exactly when α>1\alpha > 1 and a finite variance exactly when α>2\alpha > 2. Independent draws from a heavy-tailed population with finite positive variance still satisfy the central limit theorem, but the normal approximation may need a much larger sample. With infinite variance the elementary theorem does not apply.

In plain words Extreme values are rare, but not rare enough to ignore. One very large observation can move an average a long way, so the distribution of the mean can keep a long tail even in a fairly large sample.

Example The Pareto distribution with α=3\alpha = 3 has mean 1.51.5, variance 0.750.75 and P(X>10)=0.001P(X > 10) = 0.001. A normal distribution with the same mean and variance puts only about 5×10−235 \times 10^{-23} beyond 1010, which lies almost 10 SDs above its mean. The variance is finite, so the theorem applies, but slowly. In a large simulation with n=100n = 100, the standardized mean was at least 33 in about 0.0090.009 of the samples, more than six times the normal value 0.001350.00135, while the exponential waits at the same nn give about twice the normal value.

Check yourself

  1. Take 400 independent waits from the same exponential population. Find the standard error of the mean, the interval that extends two standard errors on each side of 22 minutes, and the approximate probability that the sample mean falls inside it. Are the individual waits now closer to normal?
Show the answer

The standard error is 2/400=0.12/\sqrt{400} = 0.1 minute, so two standard errors give a margin of 0.20.2 minute and the interval runs from 1.81.8 to 2.22.2 minutes. Its endpoints standardize to −2-2 and 22, so the normal approximation again gives about 0.95450.9545. For this model the exact gamma calculation gives about 0.95460.9546.

The individual waits are still exponential with mean 2 minutes and SD 2 minutes. The theorem describes the sampling distribution of the mean over repeated samples. A larger sample does not change how each wait is produced.

  1. Someone argues that with n=400n = 400 the waits are approximately normal with mean 2 and SD 2 minutes, so the chance that a single wait exceeds 6 minutes is P(Z>2)≈0.023P(Z > 2) \approx 0.023. What is wrong, and what is the correct probability?
Show the answer

The central limit theorem is about the mean of the waits, not about a single wait, which stays exponential with mean 2 minutes. For that distribution

P(X>6)=e−6/2=e−3≈0.0498,\begin{aligned} P(X > 6) &= e^{-6/2} = e^{-3} \\[0.4em] &\approx 0.0498, \end{aligned}

more than twice the claimed value. The normal model would also give a single wait a probability of about 0.1590.159 of being negative, which is impossible for a waiting time.

  1. A study records 100 waits, but they come from 10 tour groups of 10 people, and everyone in a group waits exactly as long as the group. The 10 group waits are independent and exponential with mean 2 minutes. Find the actual standard error of the mean of the 100 recorded waits, and decide whether P(1.6≤Xˉ≤2.4)P(1.6 \le \bar X \le 2.4) is still about 0.950.95.
Show the answer

Each group value appears 10 times, so the mean of the 100 recorded waits equals the mean of the 10 independent group waits. Its standard error is

210≈0.632 minute,\frac{2}{\sqrt{10}} \approx 0.632 \text{ minute},

about 3.163.16 times the value 0.20.2 that treats the 100 waits as independent. The sum of the 10 group waits has a gamma distribution, and a numerical calculation gives

P(1.6≤Xˉ≤2.4)≈0.474,\begin{aligned} &P(1.6 \le \bar X \le 2.4) \\[0.4em] &\quad \approx 0.474, \end{aligned}

about half of the 0.95450.9545 that the independence assumption would claim. One hundred rows carry only 10 independent pieces of information here, and the independence condition of the theorem fails.

How this lesson connects