Skip to content

Limit Theorems | Mathematics - Wyatt's Notes

Theorem 4.1 (Weak Law of Large Numbers). Let X1,X2,X_1, X_2, \ldots be i.i.d. With E[Xi]=μE[X_i] = \mu and Var(Xi)=σ2<\mathrm{Var}(X_i) = \sigma^2 < \infty. Then for every ε>0\varepsilon > 0:

limnP(1ni=1nXiμε)=0\lim_{n \to \infty} P\left(\left|\frac{1}{n}\sum_{i=1}^{n} X_i - \mu\right| \geq \varepsilon\right) = 0

Proof. Let Sn=1ni=1nXiS_n = \frac{1}{n}\sum_{i=1}^{n} X_i. Then E[Sn]=μE[S_n] = \mu and Var(Sn)=σ2/n\mathrm{Var}(S_n) = \sigma^2/n. By Chebyshev”s inequality:

P(Snμε)Var(Sn)ε2=σ2nε20as nP(|S_n - \mu| \geq \varepsilon) \leq \frac{\mathrm{Var}(S_n)}{\varepsilon^2} = \frac{\sigma^2}{n\varepsilon^2} \to 0 \quad \mathrm{as\ } n \to \infty

\blacksquare

Theorem 4.2 (Strong Law of Large Numbers). Under the same conditions:

P(limn1ni=1nXi=μ)=1P\left(\lim_{n \to \infty} \frac{1}{n}\sum_{i=1}^{n} X_i = \mu\right) = 1

The sample mean converges to the population mean almost surely.

Theorem 4.3 (Central Limit Theorem). Let X1,X2,X_1, X_2, \ldots be i.i.d. With E[Xi]=μE[X_i] = \mu and Var(Xi)=σ2(0,)\mathrm{Var}(X_i) = \sigma^2 \in (0, \infty). Then

SnnμσndN(0,1)\frac{S_n - n\mu}{\sigma\sqrt{n}} \xrightarrow{d} N(0, 1)

Where Sn=i=1nXiS_n = \sum_{i=1}^{n} X_i and d\xrightarrow{d} denotes convergence in distribution.

Equivalently, for large nn:

P(Snnμσnz)Φ(z)P\left(\frac{S_n - n\mu}{\sigma\sqrt{n}} \leq z\right) \approx \Phi(z)

Where Φ\Phi is the CDF of the standard normal.

Proof (using characteristic functions). Let φX(t)=E[eitX]\varphi_X(t) = E[e^{itX}] be the characteristic function of X1X_1. The characteristic function of (Snnμ)/(σn)(S_n - n\mu)/(\sigma\sqrt{n}) is:

φn(t)=[φX(tσn)]neitnμ/σ\varphi_n(t) = \left[\varphi_X\left(\frac{t}{\sigma\sqrt{n}}\right)\right]^n \cdot e^{-it\sqrt{n}\mu/\sigma}

Expanding φX\varphi_X around 0: φX(s)=1+iμs(σ2+μ2)s22+o(s2)\varphi_X(s) = 1 + i\mu s - \frac{(\sigma^2 + \mu^2)s^2}{2} + o(s^2). Substituting s=t/(σn)s = t/(\sigma\sqrt{n}):

φn(t)=[1+iμtσn(σ2+μ2)t22σ2n+o(1n)]neitnμ/σ\varphi_n(t) = \left[1 + \frac{i\mu t}{\sigma\sqrt{n}} - \frac{(\sigma^2 + \mu^2)t^2}{2\sigma^2 n} + o\left(\frac{1}{n}\right)\right]^n \cdot e^{-it\sqrt{n}\mu/\sigma}

Using limn(1+an/n)n=eliman\lim_{n \to \infty}(1 + a_n/n)^n = e^{\lim a_n}:

limnφn(t)=exp(iμtσ(σ2+μ2)t22σ2)exp(iμtσ)=et2/2\lim_{n \to \infty} \varphi_n(t) = \exp\left(\frac{i\mu t}{\sigma} - \frac{(\sigma^2 + \mu^2)t^2}{2\sigma^2}\right) \cdot \exp\left(-\frac{i\mu t}{\sigma}\right) = e^{-t^2/2}

This is the characteristic function of N(0,1)N(0, 1). By Levy’s continuity theorem, the convergence in distribution follows. \blacksquare

Problem. A fair die is rolled 100 times. Approximate the probability that the sum exceeds 370.

Solution

Let XiX_i be the value of the ii-th roll. Then E[Xi]=7/2=3.5E[X_i] = 7/2 = 3.5 and Var(Xi)=35/122.917\mathrm{Var}(X_i) = 35/12 \approx 2.917.

S100=i=1100XiS_{100} = \sum_{i=1}^{100} X_i. By the CLT:

S10035010035/12N(0,1)\frac{S_{100} - 350}{\sqrt{100 \cdot 35/12}} \approx N(0, 1)

P(S100>370)=P(Z>370350291.7)P(Z>1.17)0.121P(S_{100} > 370) = P\left(Z > \frac{370 - 350}{\sqrt{291.7}}\right) \approx P(Z > 1.17) \approx 0.121

\blacksquare

Worked Example: Sample Mean Distribution

Solution. A population has mean 50 and standard deviation 10. Find the probability that the mean of a sample of 64 observations exceeds 52.

By the CLT, XˉN(50,100/64)=N(50,1.5625)\bar{X} \approx N(50, 100/64) = N(50, 1.5625).

P(Xˉ>52)=P(Z>52501.5625)=P(Z>1.6)0.0548P(\bar{X} > 52) = P\left(Z > \frac{52 - 50}{\sqrt{1.5625}}\right) = P(Z > 1.6) \approx 0.0548

\blacksquare

4.4 Intuition: Why Do the LLN and CLT Work?

Section titled “4.4 Intuition: Why Do the LLN and CLT Work?”

The Law of Large Numbers says that the sample mean converges to the true mean as the sample size grows. The intuition is simple: random fluctuations in individual observations cancel out when you average many of them. Positive deviations and negative deviations are equally likely, and their effects average to zero. The weak law guarantees convergence in probability (most samples give a mean close to the true mean), while the strong law guarantees almost sure convergence (the sample mean is eventually and permanently close).

The Central Limit Theorem is deeper: it says that the shape of the distribution of a sum converges to a normal distribution, regardless of the original distribution. The reason is that each individual observation contributes a small amount to the sum, and the combined effect of many small independent contributions is approximately normal. This is analogous to how a random walk in many dimensions ends up approximately Gaussian. The CLT explains why the normal distribution appears so frequently in nature: any quantity that is the sum of many small independent effects will be approximately normally distributed. The Berry-Esseen theorem quantifies how fast this convergence happens, bounding the error by a constant times ρ/(σ3n)\rho/(\sigma^3 \sqrt{n}).

  • The CLT does not apply to small samples. The CLT is an asymptotic result. For small nn ( n<30n < 30), the normal approximation can be poor unless the underlying distribution is already close to normal. Use the Berry—Esseen theorem for finite-sample bounds.
  • Independence is critical for the LLN and CLT. If the XiX_i are dependent, the sample mean may not converge to the population mean, or the convergence rate may differ. For stationary sequences with weak dependence, versions of these theorems still hold, but the …/1-number-and-algebra/3_proof-and-logics are more involved.
  • Convergence in distribution is weaker than convergence in probability. The CLT gives convergence in distribution of the standardised sum, not convergence of the sum itself. The LLN gives the latter (convergence in probability).
  • Berry-Esseen bound: For i.i.d. variables with E[X3]=ρ<E[|X|^3] = \rho < \infty, the rate of CLT convergence is bounded by Fn(z)Φ(z)Cρ/(σ3n)|F_n(z) - \Phi(z)| \leq C\rho/(\sigma^3\sqrt{n}) where C<0.4748C < 0.4748.
  • Lindeberg-Feller CLT: Generalises the CLT to independent (but not identically distributed) random variables. Requires the Lindeberg condition: no single variable dominates the sum.
  • CLT for sample proportions: If YiBernoulli(p)Y_i \sim \text{Bernoulli}(p), then p^=YˉN(p,p(1p)/n)\hat{p} = \bar{Y} \approx N(p, p(1-p)/n) for large nn. This is the basis of confidence intervals for proportions.
  • Delta method: If n(Xˉμ)dN(0,σ2)\sqrt{n}(\bar{X} - \mu) \xrightarrow{d} N(0, \sigma^2), then n(g(Xˉ)g(μ))dN(0,[g(μ)]2σ2)\sqrt{n}(g(\bar{X}) - g(\mu)) \xrightarrow{d} N(0, [g'(\mu)]^2\sigma^2) for differentiable gg. Extends the CLT to nonlinear functions of the mean.

4.6 Worked Example: Confidence Intervals via CLT

Section titled “4.6 Worked Example: Confidence Intervals via CLT”

A poll surveys 1000 voters and finds 540 support a candidate. Construct a 95% confidence interval for the true proportion pp.

p^=540/1000=0.54\hat{p} = 540/1000 = 0.54. By the CLT, p^N(p,p(1p)/n)\hat{p} \approx N(p, p(1-p)/n). For a 95% CI:

p^±z0.025p^(1p^)n=0.54±1.960.54×0.461000=0.54±0.031\hat{p} \pm z_{0.025}\sqrt{\frac{\hat{p}(1-\hat{p})}{n}} = 0.54 \pm 1.96\sqrt{\frac{0.54 \times 0.46}{1000}} = 0.54 \pm 0.031

The 95% CI is (0.509,0.571)(0.509, 0.571). Since the interval includes 0.5, we cannot reject the hypothesis that the race is tied at the 5% significance level.

4.7 Worked Example: Poisson Approximation via CLT

Section titled “4.7 Worked Example: Poisson Approximation via CLT”

Events occur at rate λ=50\lambda = 50 per hour. Approximate P(S60)P(S \leq 60) where SS is the total count in one hour.

The Poisson distribution with λ=50\lambda = 50 has mean 50 and variance 50. By the CLT:

P(S60)P(Z605050)=P(Z1.414)0.9214P(S \leq 60) \approx P\left(Z \leq \frac{60 - 50}{\sqrt{50}}\right) = P(Z \leq 1.414) \approx 0.9214

The exact Poisson probability is 0.9278. The normal approximation is accurate to within 1%.

A fair die is rolled 60 times. Approximate the probability that the total is between 200 and 240.

Each roll has E[Xi]=3.5E[X_i] = 3.5 and Var(Xi)=35/12\mathrm{Var}(X_i) = 35/12. By the CLT, the total S60S_{60} satisfies:

P(200<S60<240)=P(20021060×35/12<Z<24021060×35/12)P(200 < S_{60} < 240) = P\left(\frac{200 - 210}{\sqrt{60 \times 35/12}} < Z < \frac{240 - 210}{\sqrt{60 \times 35/12}}\right)

=P(1013.23<Z<3013.23)=P(0.756<Z<2.268)= P\left(\frac{-10}{13.23} < Z < \frac{30}{13.23}\right) = P(-0.756 < Z < 2.268)

Φ(2.268)Φ(0.756)=0.98830.2248=0.7635\approx \Phi(2.268) - \Phi(-0.756) = 0.9883 - 0.2248 = 0.7635

flowchart TD
A[4_Limit Theorems] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
D --> G[Real-world usage]

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Ensure you have mastered the prerequisite material before attempting this advanced content.