Books Probability and Statistical Inference

Approximations for Discrete Distributions

The central limit theorem, Theorem 5.10, applies not only to continuous random variables but also to discrete ones.

As an example, consider π‘ŒβˆΌBin⁑(𝑛,𝑝)π‘ŒβˆΌ\operatorname{Bin}(𝑛,𝑝). Recall that this distribution arises as the sum of 𝑛𝑛 independent Bernoulli trials: if π‘‹π‘–βˆΌBer⁑(𝑝){𝑋}_{𝑖}∼\operatorname{Ber}(𝑝) are i.i.d.i.i.d. and π‘Œ=βˆ‘π‘–=1π‘›π‘‹π‘–π‘Œ={\sum }_{𝑖=1}^{𝑛}{𝑋}_{𝑖}, then π‘ŒβˆΌBin⁑(𝑛,𝑝)π‘ŒβˆΌ\operatorname{Bin}(𝑛,𝑝). For large 𝑛𝑛, the central limit theorem suggests approximating π‘Œπ‘Œ by π‘ŠβˆΌπ’©(𝑛𝑝,𝑛𝑝(1βˆ’π‘))π‘ŠβˆΌπ’©(𝑛𝑝,𝑛𝑝(1-𝑝)). Because π‘Œπ‘Œ is discrete and π‘Šπ‘Š is continuous, we need one more convention before replacing binomial probabilities by normal ones.

For each π‘˜π‘˜ in the support of π‘Œπ‘Œ, namely {0,1,…,𝑛}\{0,1,\ldots ,𝑛\}, draw a rectangle of height π‘“π‘Œ(π‘˜){𝑓}_{π‘Œ}(π‘˜) and width 11 centered at π‘˜π‘˜. Equivalently, define a histogram density hh by

h(𝑦)=π‘“π‘Œ(π‘˜)βˆ€π‘˜βˆ’1/2≀𝑦<π‘˜+1/2.h(𝑦)={𝑓}_{π‘Œ}(π‘˜)\quad \forall π‘˜-1/2\leq 𝑦<π‘˜+1/2.

This probability histogram represents the binomial masses as a continuous density on the real line. The point probability P(π‘Œ=π‘˜)β„™(π‘Œ=π‘˜) is the area of the rectangle centered at π‘˜π‘˜, which we approximate by

P(π‘Œ=π‘˜)β‰ˆP(π‘˜βˆ’1/2β‰€π‘Š<π‘˜+1/2).β„™(π‘Œ=π‘˜)\approx β„™(π‘˜-1/2\leq π‘Š<π‘˜+1/2).

More generally,

P(π‘˜β‰€π‘Œβ‰€π‘š)β‰ˆP(π‘˜βˆ’1/2β‰€π‘Š<π‘š+1/2).β„™(π‘˜\leq π‘Œ\leq π‘š)\approx β„™(π‘˜-1/2\leq π‘Š<π‘š+1/2).

The half-unit extension on each side is called the continuity correction.

Figure 5.11 Probability histogram for π‘ŒβˆΌBin⁑(10,1/2)π‘ŒβˆΌ\operatorname{Bin}(10,1/2) with normal approximation 𝒩(5,5/2)𝒩(5,5/2): each unit-width rectangle has area P(π‘Œ=π‘˜)β„™(π‘Œ=π‘˜); the shaded band shows the continuity correction P(π‘˜βˆ’1/2β‰€π‘Š<π‘˜+1/2)β„™(π‘˜-1/2\leq π‘Š<π‘˜+1/2) at π‘˜=5π‘˜=5.

Figure 5.11 overlays the histogram with the approximating normal density. The dashed curve follows the binomial masses closely, and replacing the area of one rectangle by the area under the normal curve over the same interval is the idea behind the continuity correction.