Tutorials › Stochastic Processes › Normal Approximations for Sums and Sample Means

Convergence of Random Variables · Tutorial 61 of 994

Normal Approximations for Sums and Sample Means

Standardize a count, sum, or sample mean and use the standard normal CDF to obtain a probability approximation.

Intermediate 10 min read

What You'll Learn

  • Standardize sums and sample means using their means and variances.
  • Approximate one-sided probabilities with the standard normal CDF.
  • Apply the central limit theorem to binomial counts and iid sums.
  • Approximate interval probabilities for sample means.
  • Recognize when the approximation is inapplicable or can be misleading.

From a Central Limit Theorem to a Probability

The central limit theorem describes the distribution of a standardized sum when the number of independent observations becomes large. A practical use of that result is to estimate probabilities that would otherwise require a complicated sum or convolution. The calculation has two parts: express the event in standard-deviation units, then evaluate the standard normal cumulative distribution function.

Let \(X_1,X_2,\ldots\) be iid random variables on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), with mean \(\mu\) and finite, strictly positive variance \(\sigma^2\). Define \(S_n=\sum_{i=1}^n X_i\) and \(\bar X_n=S_n/n\). Independence gives \(\mathbb{E}S_n=n\mu\) and \(\operatorname{Var}(S_n)=n\sigma^2\). Consequently, \(\mathbb{E}\bar X_n=\mu\) and \(\operatorname{Var}(\bar X_n)=\sigma^2/n\).

Definition: The standard normal CDF is \(\Phi(z)=\mathbb{P}(Z\le z)\), where \(Z\) has the normal distribution with mean \(0\) and variance \(1\). Equivalently, \(\Phi(z)=\int_{-\infty}^{z}(2\pi)^{-1/2}e^{-u^2/2}\,du\). A normal approximation replaces the distribution of a standardized sum or mean by this standard normal distribution.

Standardizing means subtracting the expectation and dividing by the standard deviation. The exact identities are

$$ \frac{S_n-n\mu}{\sigma\sqrt n} =\frac{\bar X_n-\mu}{\sigma/\sqrt n}. $$

Thus a threshold \(c\) for the sum corresponds to the standardized threshold \((c-n\mu)/(\sigma\sqrt n)\). A threshold \(a\) for the sample mean corresponds to \((a-\mu)/(\sigma/\sqrt n)\). These identities are algebra, not approximations; the approximation enters when the standardized variable is treated as standard normal.

Theorem (Normal approximation for sums and means): Suppose \(X_1,X_2,\ldots\) are iid with finite mean \(\mu\) and variance \(0<\sigma^2<\infty\). For each fixed real \(z\), as \(n\to\infty\), \[ \mathbb{P}\!\left(\frac{S_n-n\mu}{\sigma\sqrt n}\le z\right)\longrightarrow \Phi(z). \] Consequently, for fixed real \(a\le b\), \[ \mathbb{P}\!\left(a\le\frac{S_n-n\mu}{\sigma\sqrt n}\le b\right) \longrightarrow \Phi(b)-\Phi(a). \] For a particular finite \(n\), using the corresponding standard normal probability gives an approximation, not an identity.

Proof. By the central limit theorem for standardized sums, the standardized variables \(W_n=(S_n-n\mu)/(\sigma\sqrt n)\) converge in distribution to a standard normal variable \(Z\). Convergence in distribution implies convergence of CDFs at every continuity point of the limiting CDF. Since \(\Phi\) is continuous at every real number, \(\mathbb{P}(W_n\le z)\to\Phi(z)\) for each fixed \(z\). For the interval statement, write \(\mathbb{P}(a\le W_n\le b)=\mathbb{P}(W_n\le b)-\mathbb{P}(W_n<a)\). For any \(\epsilon>0\), \(\mathbb{P}(W_n\le a-\epsilon)\le\mathbb{P}(W_n<a)\le\mathbb{P}(W_n\le a)\). Taking lower and upper limits in \(n\) bounds the lower and upper limits of \(\mathbb{P}(W_n<a)\) between \(\Phi(a-\epsilon)\) and \(\Phi(a)\). Letting \(\epsilon\) decrease to zero and using continuity gives \(\mathbb{P}(W_n<a)\to\Phi(a)\). The two CDF limits give the stated interval limit. \(\square\)

The theorem is an asymptotic statement: it does not give an error bound for a chosen sample size. It also requires a nonzero finite variance. If \(\sigma^2=0\), each \(X_i=\mu\) almost surely, so the standardized expression divides by zero and is not defined; the sum and mean are instead constant. If the variance is infinite, this version of the central limit theorem does not apply.

Three Standardization Calculations

Worked Example: A Binomial Count

Let \(X\) be the number of successes in \(100\) independent trials, each with success probability \(p=\tfrac25\). Approximate \(\mathbb{P}(X\le45)\) by standardizing the count.

Write \(X=\sum_{i=1}^{100}I_i\), where \(I_i\) is \(1\) on a success and \(0\) otherwise. The indicators are iid, with mean \(p=\tfrac25\) and variance \(p(1-p)=\tfrac25\cdot\tfrac35=\tfrac6{25}\). Therefore the mean of \(X\) is \(100\cdot\tfrac25=40\), and its variance is \(100\cdot\tfrac6{25}=24\). The standard deviation is \(\sqrt{24}=2\sqrt6\). Each trial has \(0<p<1\), so the variance is positive, and the iid finite-variance assumptions hold.

For the event \(X\le45\), the standardized endpoint is

$$ z=\frac{45-40}{\sqrt{24}}=\frac{5}{2\sqrt6}\approx1.0206. $$

Using the standard normal CDF table, \(\Phi(1.0206)\approx0.8462\), so the normal approximation is \(\mathbb{P}(X\le45)\approx0.8462\). This is a CDF approximation at the event's stated endpoint; it is not an exact binomial probability.

Check: The endpoint \(45\) is five counts above the mean \(40\), and the standard deviation is \(\sqrt{24}\), so the standardized endpoint is positive and a little greater than \(1\). The resulting CDF value, about \(0.8462\), is greater than \(\Phi(1)\approx0.8413\), as expected for a larger argument.

Worked Example: A Sum of Non-Bernoulli Observations

Suppose \(X_1,\ldots,X_{36}\) are iid, each equal to \(-1\) or \(5\) with probability \(\tfrac12\) each. Approximate \(\mathbb{P}(S_{36}\le81)\).

For one observation, \(\mu=\tfrac12(-1)+\tfrac12(5)=2\). Its variance is \(\tfrac12(-1-2)^2+\tfrac12(5-2)^2=\tfrac12(9)+\tfrac12(9)=9\), so \(\sigma=3\). Thus \(\mathbb{E}S_{36}=36\cdot2=72\), \(\operatorname{Var}(S_{36})=36\cdot9=324\), and the standard deviation of the sum is \(\sqrt{324}=18\). The observations are independent and have finite, positive variance.

The event threshold standardizes to \(z=(81-72)/18=9/18=\tfrac12\). The normal approximation is therefore \(\mathbb{P}(S_{36}\le81)\approx\Phi(\tfrac12)\approx0.6915\), using the standard normal CDF table.

Check: Standardizing the sample mean gives the same threshold. The event is \(\bar X_{36}\le81/36=9/4\), and its standardized endpoint is \((\tfrac94-2)/(3/\sqrt{36})=(\tfrac14)/(\tfrac12)=\tfrac12\). This agrees with the calculation for the sum.

Worked Example: An Interval for a Sample Mean

Let \(X_1,\ldots,X_{25}\) be iid, with mean \(50\) and variance \(100\). Approximate \(\mathbb{P}(48\le\bar X_{25}\le53)\).

The standard deviation of one observation is \(\sqrt{100}=10\), so the standard deviation of the sample mean is \(10/\sqrt{25}=2\). The lower endpoint standardizes to \((48-50)/2=-1\), and the upper endpoint standardizes to \((53-50)/2=3/2\). The central limit theorem applies because the observations are iid with finite, positive variance. The interval approximation is

$$ \mathbb{P}(48\le\bar X_{25}\le53) \approx\Phi(3/2)-\Phi(-1) \approx0.9332-0.1587 =0.7745. $$

The CDF values are standard normal table values at the two computed endpoints. The interval is inclusive; the limiting normal distribution assigns probability zero to each endpoint, so inclusion or exclusion of an endpoint does not change the limiting interval probability.

Check: The standardized interval is \([-1,3/2]\), which contains zero, the standard normal mean. Its approximate probability is positive and exceeds one half, consistent with an interval that extends one standard deviation below and one and a half above the mean.

Reading the Approximation Correctly

For a lower-tail probability, standardize the event's cutoff and read \(\Phi(z)\). For an upper-tail probability, use the complement \(1-\Phi(z)\), taking care to express the original event with the correct inequality. For an interval, subtract the lower-end CDF value from the upper-end CDF value. These operations approximate probabilities under the standardized distribution; they do not change the original event or make its finite-sample distribution continuous.

For a binomial variable, the event \(X\le k\) has the exact standardized form \((X-np)/\sqrt{np(1-p)}\le(k-np)/\sqrt{np(1-p)}\). The theorem motivates evaluating \(\Phi\) at that standardized cutoff for large \(n\). A discrete distribution can have jumps, so a normal CDF does not reproduce its finite-sample probabilities exactly. How to adjust cutoffs for discreteness, and how to assess approximation quality, are separate questions.

Common mistake: using the standard deviation of one observation for the sample mean. In the sample-mean example, an incorrect calculation might divide the endpoint differences by \(10\), obtaining \((48-50)/10=-0.2\) and \((53-50)/10=0.3\). This would give \(\Phi(0.3)-\Phi(-0.2)\approx0.6179-0.4207=0.1972\). The error is that \(10\) is the standard deviation of one observation, not of the mean of \(25\) observations. The correct standard deviation is \(10/\sqrt{25}=2\), giving endpoints \(-1\) and \(3/2\) and the approximation \(0.7745\). Dividing by the wrong scale changes the event's location in standard-deviation units.

Finally, the theorem does not guarantee that every finite-sample approximation is accurate. The assumptions explain when the asymptotic result applies; accuracy at a particular \(n\) can depend on the underlying distribution and the event being estimated. Treat the displayed normal probability as an approximation, and retain enough digits in the standardized endpoint to avoid unnecessary rounding error before consulting a CDF table.

Check Your Understanding

For each question, write the standardization before evaluating or interpreting a normal CDF.

  1. If \(S_n\) is a sum of iid variables with mean \(\mu\) and variance \(\sigma^2>0\), what is its standard deviation, and what is the standardized endpoint for \(S_n\le c\)?
  2. A binomial count has \(n=64\) and \(p=\tfrac12\). Find its mean, variance, and standardized endpoint for the event \(X\le36\).
  3. If an iid population has mean \(10\), variance \(16\), and sample size \(n=16\), what is the standard deviation of the sample mean?
  4. Why does the normal approximation theorem not apply in its stated form when the observation variance is zero?
  5. For a standardized interval \(a\le W_n\le b\), which two standard normal CDF values determine the limiting probability?