Normal-Approximation

Die Normal-Approximation ist eine Methode der Wahrscheinlichkeitsrechnung, um die Binomialverteilung für große Stichproben durch die Normalverteilung anzunähern. Hierbei handelt es sich um eine Anwendung des Satzes von Moivre-Laplace und damit auch um eine Anwendung des Zentralen Grenzwertsatzes.

Formulierung

Für eine binomialverteilte Zufallsvariable $S_n\sim \mathrm{B}\mathrm{i}\mathrm{n}(n,p)$ mit dem Erwartungswert ${\mu }_n:=𝔼[S_n]=np$ und der Standardabweichung ${\sigma }_n:=\sqrt{\mathrm{V}\mathrm{a}\mathrm{r}[S_n]}=\sqrt{np(1-p)}>0$ gilt nach dem Satz von Moivre-Laplace

$\underset{n\to \mathrm{\infty }}{\lim }(\mathrm{P}(S_n\le x)-\mathrm{\Phi }\left(\frac{x-{\mu }_n}{{\sigma }_n}\right))=0,x\in ℝ,$,

wobei $\mathrm{\Phi }$ die Wahrscheinlichkeitsverteilungsfunktion der Standardnormalverteilung bezeichnet.

Diesen asymptotischen Zusammenhang verwendet man zur Rechtfertigung der Approximation der Verteilungsfunktion von $S_n$ durch die Verteilungsfunktion der Normalverteilung $𝒩({\mu }_n,{\sigma }_n^2)$ für endliches, aber hinreichend große $n$, d. h.

$\mathrm{P}(S_n\le x)\approx \mathrm{\Phi }\left(\frac{x-{\mu }_n}{{\sigma }_n}\right),x\in ℝ.$

Mit dieser Approximation der Verteilungsfunktion erhält man für $k_1,k_2\in \{0,1,\dots ,n\}$ mit $k_1<k_2$

$\begin{array}{cc}\mathrm{P}(k_1<S_n\le k_2)& =\underset{\text{Binomialverteilung}}{\underbrace{\mathrm{P}(S_n\le k_2)-\mathrm{P}(S_n\le k_1)}}\\ & \approx \underset{\text{Normalverteilung}}{\underbrace{\mathrm{\Phi }\left(\frac{k_2-{\mu }_n}{{\sigma }_n}\right)-\mathrm{\Phi }\left(\frac{k_1-{\mu }_n}{{\sigma }_n}\right)}}.\end{array}$

Eine exakte Berechnung der Wahrscheinlichkeit würde dagegen durch

$\mathrm{P}(k_1<S_n\le k_2)=\sum _{k=k_1+1}^{k_2}(\genfrac{}{}{0ex}{}{n}{k})p^k(1-p{)}^{n-k}$

erfolgen.

Güte der Approximation

Für viele Anwendungen gilt die Näherung als hinreichend gut, falls $np(1-p)\ge 9$ gilt. Falls dies nicht gilt, so sollte zumindest $np\ge 5$ und $n(1-p)\ge 5$ gelten. Je asymmetrischer die Binomialverteilung ist, d. h. je größer die Differenz zwischen $p$ und $1-p$ ist, umso größer sollte $n$ sein.

Für $p$ nahe an 0 ist zur Näherung die Poisson-Approximation besser geeignet. Dabei wird die Binomialverteilung $\mathrm{B}\mathrm{i}\mathrm{n}(n,p)$ durch eine Poissonverteilung mit dem Parameter $np$ approximiert. Für $p$ nahe an 1 sind beide Approximationen schlecht, dann kann jedoch ${S_n}^{\prime }=n-S_n$ statt $S_n$ betrachtet werden, d. h. bei der Binomialverteilung werden Erfolge und Misserfolge vertauscht. ${S_n}^{\prime }$ ist wieder binomialverteilt mit Parametern $n$ und $1-p$ und kann daher mit der Poisson-Approximation angenähert werden.

Stetigkeitskorrektur

Mit der Stetigkeitskorrektur wird eine verbesserte Approximation von Wahrscheinlichkeiten der Binomialverteilung $\mathrm{B}\mathrm{i}\mathrm{n}(n,p)$ berechneten Wahrscheinlichkeiten durch die Wahrscheinlichkeiten aus der Normalverteilungsapproximation $𝒩(\mu ,{\sigma }^2)$ mit den Parametern $\mu =np$ und ${\sigma }^2=np(1-p)$ und der Verteilungsfunktion $x\mapsto \mathrm{\Phi }((x-\mu )/\sigma )$ angestrebt.

Die Grundidee ist, den $n+1$ Stellen $0,1,\dots ,n$ der Binomialverteilung die Wahrscheinlichkeiten der Normalverteilung in den $n+1$ Intervallen

$I_k=\{\begin{array}{cc}(-\mathrm{\infty },k+1/2)& \text{für }k=0\\ (k-1/2,k+1/2)& \text{für }k=1,\dots ,n-1\\ (k-1/2,\mathrm{\infty })& \text{für }k=n\end{array}$

zuzuordnen. Für $X\sim 𝒩(\mu ,{\sigma }^2)$ gilt dann

$P(X\in I_k)=\{\begin{array}{cc}\mathrm{\Phi }\left(\frac{k+1/2-\mu }{\sigma }\right)& \text{für }k=0\\ \mathrm{\Phi }\left(\frac{k+1/2-\mu }{\sigma }\right)-\mathrm{\Phi }\left(\frac{k-1/2-\mu }{\sigma }\right)& \text{für }k=1,\dots ,n-1\\ 1-\mathrm{\Phi }\left(\frac{k-1/2-\mu }{\sigma }\right)& \text{für }k=n\end{array}$

Für einen Wert $k\in \{0,1,\dots ,n\}$ ergeben sich mit Verwendung der Stetigkeitskorrektur die Approximationen

$P(S_n=k)\approx \mathrm{\Phi }\left(\frac{k+1/2-\mu }{\sigma }\right)-\mathrm{\Phi }\left(\frac{k-1/2-\mu }{\sigma }\right),$

$P(S_n\le k)\approx \mathrm{\Phi }\left(\frac{k+1/2-\mu }{\sigma }\right)$

und

$P(S_n\ge k)\approx 1-\mathrm{\Phi }\left(\frac{k-1/2-\mu }{\sigma }\right).$

Für zwei Werte $k_1,k_2\in \{0,1,\dots ,n\}$ mit $k_1\le k_2$ ergibt sich mit Verwendung der Stetigkeitskorrektur die Approximation

$P(k_1\le S_n\le k_2)\approx \mathrm{\Phi }\left(\frac{k_2+0,5-\mu }{\sigma }\right)-\mathrm{\Phi }\left(\frac{k_1-0,5-\mu }{\sigma }\right).$

Beispiele

Beispiel 1

Ein fairer Würfel wird 1000 Mal geworfen. Gesucht ist die Wahrscheinlichkeit, dass mindestens 100 und höchstens 150 Mal die Sechs gewürfelt wird.

• Berechnung mit Binomialverteilung

Die Zufallsvariable $S_{1000}$, die die zufällige Anzahl der gewürfelten Sechsen bei $1000$ Versuchen beschreibt, ist binomialverteilt mit den Parametern $n=1000$ und $p=1/6$, es gilt also $S_{1000}\sim \mathrm{B}\mathrm{i}\mathrm{n}(1000,1/6)$. Die gesuchte Wahrscheinlichkeit ist daher

$P(100\le S_{1000}\le 150)=\sum _{k=100}^{150}(\genfrac{}{}{0ex}{}{1000}{k}){\left(\frac{1}{6}\right)}^k{\left(\frac{5}{6}\right)}^{1000-k}\approx 0,0837.$

Mit einer Wahrscheinlichkeit von ca. 8,4 % wird also bei 1000 Versuchen zwischen 100 und 150 Mal die Sechs gewürfelt.

• Berechnung mit Normal-Approximation

Es ist $np(1-p)=1000\frac{1}{6}\frac{5}{6}=138,\overline{8}>9$. Entsprechend der Faustformel gilt die approximierte Lösung also ausreichend genau. Die beiden Parameter für die Approximation der Binomialverteilung durch eine Normalverteilung sind $\mu =np=1000/6$ und $\sigma =\sqrt{np(1-p)}=\sqrt{5000/36}.$

Die Approximation ohne Stetigkeitskorrektur, wobei die Verteilungsfunktion der Binomialverteilung durch die Verteilungsfunktion einer Normalverteilung mit denselben Parametern approximiert wird, ist

$\begin{array}{cc}P(100\le S_{1000}\le 150)& =P(S_{1000}\le 150)-P(S_{1000}\le 99)\\ & \approx \mathrm{\Phi }\left(\frac{150-\mu }{\sigma }\right)-\mathrm{\Phi }\left(\frac{99-\mu }{\sigma }\right)\\ & \approx \mathrm{\Phi }(-1,4142)-\mathrm{\Phi }(-5,74171)\approx 1-\mathrm{\Phi }(1,4142)\approx 0,0786\end{array}$

• Berechnung mit Normal-Approximation und Stetigkeitskorrektur

Die Approximation mit Stetigkeitskorrektur ist

$\begin{array}{cc}P(100\le S_{1000}\le 150)& \approx \mathrm{\Phi }\left(\frac{150+0,5-\mu }{\sigma }\right)-\mathrm{\Phi }\left(\frac{100-0,5-\mu }{\sigma }\right)\\ & \approx \mathrm{\Phi }(-1,3718)-\mathrm{\Phi }(-5,6993)\approx 1-\mathrm{\Phi }(1,3718)\approx 0,0851\end{array}$

Die Werte $\mathrm{\Phi }(z)$ können numerisch bestimmt werden oder aus Tabelle abgelesen werden, da keine explizite Stammfunktion existiert. Die approximierte Lösung ist häufig numerisch günstiger, da keine umfangreichen Berechnungen der Binomialkoeffizienten durchgeführt werden müssen.

Beispiel 2

Plot der Dichte der Normalverteilung mit μ = 12 und σ = 3 und der Binomialverteilung mit n = 48 und p = 1/4

Gegeben sei eine binomialverteilte Zufallsvariable $X$ mit den Parametern $n=48$ und $p=\frac{1}{4}$, dann hat $X$ den Erwartungswert $np=12$ und die Varianz $np(1-p)=9$. Die Binomialverteilung wird durch eine Normalverteilung $𝒩(\mu ,{\sigma }^2)$ mit dem Mittelwert $\mu =np=12$ und der Varianz ${\sigma }^2=np(1-p)=9$ approximiert.

Nun suchen wir die Antwort auf die Frage: „Wie groß ist die Wahrscheinlichkeit, dass $X$ Werte kleiner oder gleich 3 annimmt?“ bzw. „Wie groß ist die Wahrscheinlichkeit $P(0\le X\le 3)$?“ Da $3=\mu -3\sigma$ ist, handelt es sich um eine kleine Wahrscheinlichkeit im linken Verteilungsende, die exakt mit Hilfe der Binomialverteilung oder approximativ aus der Normalverteilung berechnet werden kann. Für diese Fragestellung und Parameterkonstellation ergeben sich folgende Resultate:

• Berechnung mit der Binomialverteilung

$P(0\le X\le 3)=\sum _{k=0}^3(\genfrac{}{}{0ex}{}{48}{k})\cdot (\frac{1}{4}{)}^k\cdot (\frac{3}{4}{)}^{48-k}\approx 0,0007882$

• Abschätzung mit der Normalverteilung

$\begin{array}{cc}P(0\le X\le 3)& =P(X\le 3)-P(X\le -1)\\ & \approx \mathrm{\Phi }\left(\frac{3-12}{3}\right)-\mathrm{\Phi }\left(\frac{-1-12}{3}\right)\\ & =\mathrm{\Phi }(-3)-\mathrm{\Phi }(-13/3)\\ & \approx 0,0013499-0,0000073\approx 0,001343\end{array}$

• Abschätzung mit Normalverteilung und Stetigkeitskorrektur

$\begin{array}{cc}P(0\le X\le 3)& \approx \mathrm{\Phi }\left(\frac{3+0,5-12}{3}\right)-\mathrm{\Phi }\left(\frac{0-0,5-12}{3}\right)\\ & =\mathrm{\Phi }(-8,5/3)-\mathrm{\Phi }(-12,5/3)\\ & \approx 0,002303-0.0000015\approx 0,00229\end{array}$

Die Berechnung mit der Normalverteilung ohne Stetigkeitskorrektur überschätzt die Wahrscheinlichkeit. Dies ist darauf zurückzuführen, dass die Binomialverteilung bei dieser Parameterkonstellation eine Schiefe aufweist, die dazu führt, dass die Normalverteilung im linken Verteilungsende die entsprechenden Wahrscheinlichkeiten überschätzt und im rechten Verteilungsende überschätzt. Ob dieser Approximationsfehler akzeptabel ist, hängt von der Anwendung und Fragestellung ab. Der absolute Approximationsfehler 0,001343 - 0,0007882 liegt unter 0,001 und kann damit z. B. in solchen statistischen Anwendungsfällen akzeptabel sein, bei denen Konfidenz- und Signifikanzniveaus in einer der groben Stufen $\alpha \in \{10\mathrm{\%},5\mathrm{\%},1\mathrm{\%}\}$ festgelegt werden. Anderseits ist der relative Approximationsfehler (0,001343 - 0,0007882)/0,0007882 = 70,39 %. Dies kann bei anderen Anwendungen eine unakzeptable Ungenauigkeit der zu bestimmenden Wahrscheinlichkeit sein.

Die Verwendung der Stetigkeitskorrektur verschlechtert bei dieser Konstellation die Approximation durch die Normalverteilung, da durch die Berücksichtigung zusätzlicher Wahrscheinlichkeitsmasse die Überschätzung durch die Normalverteilung weiter verstärkt wird. Das Beispiel zeigt insofern die Problematik einer standardmäßigen Anwendung der Stetigkeitskorrektur, wenn damit die Vorstellung einer gleichmäßigen Verbesserung verbunden ist.

Literatur

• Hans-Otto Georgii: Stochastik: Einführung in die Wahrscheinlichkeitstheorie und Statistik, 4. Auflage, de Gruyter, 2009, ISBN 978-3-11-021526-7, Vorlage:Doi.
• Ulrich Krengel: Einführung in die Wahrscheinlichkeitstheorie und Statistik. Vieweg, Braunschweig 1988, ISBN 978-3-528-07259-9, Vorlage:Doi.