Books Probability and Statistical Inference

The Central Limit Theorem

Let 𝑋1,…,𝑋𝑛{𝑋}_{1},\ldots ,{𝑋}_{𝑛} be a random sample of size 𝑛𝑛 from a distribution with mean ΞΌ\mu and variance Οƒ2{\sigma }^{2}. By Theorem 5.4, the sample mean 𝑋𝑋 has mean ΞΌ\mu and variance Οƒ2/𝑛{\sigma }^{2}/𝑛. As π‘›β†’βˆžπ‘›\to \infty , the variance of 𝑋𝑋 approaches 00, so 𝑋𝑋 degenerates (or converges in distribution) to a constant random variable with value ΞΌ\mu . In a later section, we will introduce the concept of convergence in distribution.

Let

π‘Š=π‘‹βˆ’ΞΌΟƒ/𝑛.π‘Š=\frac{𝑋-\mu }{\sigma /\sqrt{𝑛}}.

The mean of π‘Šπ‘Š is given by

E(π‘Š)=1Οƒ/𝑛E(π‘‹βˆ’ΞΌ)=0.\mathbb{E}(π‘Š)=\frac{1}{\sigma /\sqrt{𝑛}}\mathbb{E}(𝑋-\mu )=0.

And the variance is

Var⁑(π‘Š)=E[(π‘‹βˆ’ΞΌΟƒ/𝑛)2]=E[(π‘‹βˆ’ΞΌ)2]Οƒ2/𝑛=Var⁑(𝑋)Οƒ2/𝑛=Οƒ2/𝑛σ2/𝑛=1.\operatorname{Var}(π‘Š)=\mathbb{E}\left[{\left(\frac{𝑋-\mu }{\sigma /\sqrt{𝑛}}\right)}^{2}\right]=\frac{\mathbb{E}[{(𝑋-\mu )}^{2}]}{{\sigma }^{2}/𝑛}=\frac{\operatorname{Var}(𝑋)}{{\sigma }^{2}/𝑛}=\frac{{\sigma }^{2}/𝑛}{{\sigma }^{2}/𝑛}=1.

Though the numerator π‘‹βˆ’ΞΌπ‘‹-\mu of π‘Šπ‘Š degenerates to zero, the denominator Οƒ/𝑛\sigma /\sqrt{𝑛} spreads out the probability density over a wider range to prevent the degeneration.

We are interested in the distribution of π‘Šπ‘Š when 𝑛𝑛 is large. If the random sample is from a normal distribution, then we know that π‘‹βˆΌπ’©(ΞΌ,Οƒ2/𝑛)π‘‹βˆΌπ’©(\mu ,{\sigma }^{2}/𝑛) and hence π‘ŠβˆΌπ’©(0,1)π‘ŠβˆΌπ’©(0,1). But for a general distribution, π‘Šπ‘Š needs not to be a standard normal random variable. However, as 𝑛𝑛 increases, the distribution of π‘Šπ‘Š approaches the standard normal distribution. That is, the asymptotic behavior of π‘Šπ‘Š is like a standard normal random variable. This is the content of the Central Limit Theorem.

Theorem 5.10Central Limit Theorem

Let 𝑋1,…,𝑋𝑛{𝑋}_{1},\ldots ,{𝑋}_{𝑛} be a random sample of size 𝑛𝑛 from a distribution with mean ΞΌ\mu and variance Οƒ2{\sigma }^{2}. Then the distribution of

π‘Š=π‘‹βˆ’ΞΌΟƒ/𝑛=βˆ‘π‘–=1π‘›π‘‹π‘–βˆ’π‘›ΞΌπ‘›Οƒπ‘Š=\frac{𝑋-\mu }{\sigma /\sqrt{𝑛}}=\frac{{\sum }_{𝑖=1}^{𝑛}{𝑋}_{𝑖}-𝑛\mu }{\sqrt{𝑛}\sigma }

is 𝒩(0,1)𝒩(0,1) in limit as π‘›β†’βˆžπ‘›\to \infty .

In practice, then, for large 𝑛𝑛 we can treat π‘Šπ‘Š as approximately 𝒩(0,1)𝒩(0,1), or equivalently 𝑋𝑋 as approximately 𝒩(ΞΌ,Οƒ2/𝑛)𝒩(\mu ,{\sigma }^{2}/𝑛) and βˆ‘π‘–=1𝑛𝑋𝑖{\sum }_{𝑖=1}^{𝑛}{𝑋}_{𝑖} as approximately 𝒩(𝑛μ,𝑛σ2)𝒩(𝑛\mu ,𝑛{\sigma }^{2}). Most often it is the sum that we need, so we usually work with βˆ‘π‘–=1𝑛𝑋𝑖{\sum }_{𝑖=1}^{𝑛}{𝑋}_{𝑖}.

Example 5.8

Let 𝑋1,𝑋2,𝑋3,𝑋4βˆΌπ’°(0,1){𝑋}_{1},{𝑋}_{2},{𝑋}_{3},{𝑋}_{4}βˆΌπ’°(0,1) be i.i.d. Each 𝑋𝑖{𝑋}_{𝑖} has mean ΞΌ=1/2\mu =1/2 and variance Οƒ2=1/12{\sigma }^{2}=1/12, so a sum of 𝑛𝑛 such variables has mean 𝑛μ𝑛\mu and variance 𝑛σ2𝑛{\sigma }^{2}. The central limit theorem suggests that for large 𝑛𝑛 the sum is well approximated by 𝒩(𝑛μ,𝑛σ2)𝒩(𝑛\mu ,𝑛{\sigma }^{2}). Here we compare the exact PDF of a sum with this normal approximation for 𝑛=2𝑛=2 and 𝑛=4𝑛=4.

For π‘Œ=𝑋1+𝑋2π‘Œ={𝑋}_{1}+{𝑋}_{2} we have 𝑛=2𝑛=2, so the approximating density is 𝒩(1,1/6)𝒩(1,1/6). From Exercise 5.6, the exact PDF of π‘Œπ‘Œ is

π‘“π‘Œ(𝑦)={𝑦,0≀𝑦<12βˆ’π‘¦,1≀𝑦≀20,otherwise{𝑓}_{π‘Œ}(𝑦)=\begin{cases}𝑦,\quad & 0\leq 𝑦<1 \\ 2-𝑦,\quad & 1\leq 𝑦\leq 2 \\ 0,\quad & \text{otherwise}\end{cases}

For π‘Š=𝑋1+𝑋2+𝑋3+𝑋4π‘Š={𝑋}_{1}+{𝑋}_{2}+{𝑋}_{3}+{𝑋}_{4} we have 𝑛=4𝑛=4, so the approximating density is 𝒩(2,1/3)𝒩(2,1/3). Write π‘Œ1=𝑋1+𝑋2{π‘Œ}_{1}={𝑋}_{1}+{𝑋}_{2} and π‘Œ2=𝑋3+𝑋4{π‘Œ}_{2}={𝑋}_{3}+{𝑋}_{4}. Then π‘Œ1{π‘Œ}_{1} and π‘Œ2{π‘Œ}_{2} are independent and each has the same distribution as π‘Œπ‘Œ, and π‘Š=π‘Œ1+π‘Œ2π‘Š={π‘Œ}_{1}+{π‘Œ}_{2}. By Proposition 5.1,

π‘“π‘Š(𝑀)=(π‘“π‘Œβˆ—π‘“π‘Œ)(𝑀)={𝑀36,0≀𝑀<1βˆ’3𝑀3+12𝑀2βˆ’12𝑀+46,1≀𝑀<23𝑀3βˆ’24𝑀2+60π‘€βˆ’446,2≀𝑀<3βˆ’π‘€3+12𝑀2βˆ’48𝑀+646,3≀𝑀<40,otherwise{𝑓}_{π‘Š}(𝑀)=({𝑓}_{π‘Œ}\ast {𝑓}_{π‘Œ})(𝑀)=\begin{cases}\frac{{𝑀}^{3}}{6},\quad & 0\leq 𝑀<1 \\ \frac{-3{𝑀}^{3}+12{𝑀}^{2}-12𝑀+4}{6},\quad & 1\leq 𝑀<2 \\ \frac{3{𝑀}^{3}-24{𝑀}^{2}+60𝑀-44}{6},\quad & 2\leq 𝑀<3 \\ \frac{-{𝑀}^{3}+12{𝑀}^{2}-48𝑀+64}{6},\quad & 3\leq 𝑀<4 \\ 0,\quad & \text{otherwise}\end{cases}
Figure 5.10 Sums of i.i.d.
𝒰(0,1)𝒰(0,1) variables compared with normal approximations from the central limit theorem: (a) π‘Œ=𝑋1+𝑋2π‘Œ={𝑋}_{1}+{𝑋}_{2} with 𝒩(1,1/6)𝒩(1,1/6); (b) π‘Š=𝑋1+𝑋2+𝑋3+𝑋4π‘Š={𝑋}_{1}+{𝑋}_{2}+{𝑋}_{3}+{𝑋}_{4} with 𝒩(2,1/3)𝒩(2,1/3).

Figure 5.10 overlays each exact density with its normal approximation. The fit at 𝑛=2𝑛=2 is only moderate, but at 𝑛=4𝑛=4 the two curves are already nearly the same, even though the theorem concerns only the limit as π‘›β†’βˆžπ‘›\to \infty .