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 B i n (n, p)$ mit dem Erwartungswert $μ_{n} := 𝔼 [S_{n}] = n p$ und der Standardabweichung $σ_{n} := \sqrt{V a r [S_{n}]} = \sqrt{n p (1 - p)} > 0$ gilt nach dem Satz von Moivre-Laplace

lim_{n \to \infty} (P (S_{n} \leq x) - Φ (\frac{x - μ_{n}}{σ_{n}})) = 0, x \in ℝ,

,

wobei $Φ$ die Wahrscheinlichkeitsverteilungsfunktion der Standardnormalverteilung bezeichnet.

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

P (S_{n} \leq x) \approx Φ (\frac{x - μ_{n}}{σ_{n}}), 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{matrix} P (k_{1} < S_{n} \leq k_{2}) & = \underset{Binomialverteilung}{\underset{⏟}{P (S_{n} \leq k_{2}) - P (S_{n} \leq k_{1})}} \\ \approx \underset{Normalverteilung}{\underset{⏟}{Φ (\frac{k_{2} - μ_{n}}{σ_{n}}) - Φ (\frac{k_{1} - μ_{n}}{σ_{n}})}} . \end{matrix}

Eine exakte Berechnung der Wahrscheinlichkeit würde dagegen durch

P (k_{1} < S_{n} \leq k_{2}) = \sum_{k = k_{1} + 1}^{k_{2}} (\binom{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 $n p (1 - p) \geq 9$ gilt.^[1]^[2]^[3] Falls dies nicht gilt, so sollte zumindest $n p \geq 5$ und $n (1 - p) \geq 5$ gelten.^[4]^[5] 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 $B i n (n, p)$ durch eine Poissonverteilung mit dem Parameter $n p$ approximiert. Für $p$ nahe an 1 sind beide Approximationen schlecht, dann kann jedoch ${S_{n}}^{'} = n - S_{n}$ statt $S_{n}$ betrachtet werden, d. h. bei der Binomialverteilung werden Erfolge und Misserfolge vertauscht. ${S_{n}}^{'}$ 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 $B i n (n, p)$ berechneten Wahrscheinlichkeiten durch die Wahrscheinlichkeiten aus der Normalverteilungsapproximation $𝒩 (μ, σ^{2})$ mit den Parametern $μ = n p$ und $σ^{2} = n p (1 - p)$ und der Verteilungsfunktion $x \mapsto Φ ((x - μ) / σ)$ 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{matrix} (- \infty, k + 1 / 2) & für k = 0 \\ (k - 1 / 2, k + 1 / 2) & für k = 1, \dots, n - 1 \\ (k - 1 / 2, \infty) & für k = n \end{matrix}

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

P (X \in I_{k}) = {\begin{matrix} Φ (\frac{k + 1 / 2 - μ}{σ}) & für k = 0 \\ Φ (\frac{k + 1 / 2 - μ}{σ}) - Φ (\frac{k - 1 / 2 - μ}{σ}) & für k = 1, \dots, n - 1 \\ 1 - Φ (\frac{k - 1 / 2 - μ}{σ}) & für k = n \end{matrix}

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

P (S_{n} = k) \approx Φ (\frac{k + 1 / 2 - μ}{σ}) - Φ (\frac{k - 1 / 2 - μ}{σ}),

P (S_{n} \leq k) \approx Φ (\frac{k + 1 / 2 - μ}{σ})

und

P (S_{n} \geq k) \approx 1 - Φ (\frac{k - 1 / 2 - μ}{σ}) .

Für zwei Werte $k_{1}, k_{2} \in {0, 1, \dots, n}$ mit $k_{1} \leq k_{2}$ ergibt sich mit Verwendung der Stetigkeitskorrektur die Approximation

P (k_{1} \leq S_{n} \leq k_{2}) \approx Φ (\frac{k_{2} + 1 / 2 - μ}{σ}) - Φ (\frac{k_{1} - 1 / 2 - μ}{σ}) .

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 B i n (1000, 1 / 6)

.^[6] Die gesuchte Wahrscheinlichkeit ist daher

P (100 \leq S_{1000} \leq 150) = \sum_{k = 100}^{150} (\binom{1000}{k}) {(\frac{1}{6})}^{k} {(\frac{5}{6})}^{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

n p (1 - p) = 1000 \frac{1}{6} \frac{5}{6} = 138, \bar{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

μ = n p = 1000 / 6

und

σ = \sqrt{n p (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{matrix} P (100 \leq S_{1000} \leq 150) & = P (S_{1000} \leq 150) - P (S_{1000} \leq 99) \\ \approx Φ (\frac{150 - μ}{σ}) - Φ (\frac{99 - μ}{σ}) \\ \approx Φ (- 1,4142) - Φ (- 5,74171) \approx 1 - Φ (1,4142) \approx 0,0786 \end{matrix}

Berechnung mit Normal-Approximation und Stetigkeitskorrektur

Die Approximation mit Stetigkeitskorrektur ist

\begin{matrix} P (100 \leq S_{1000} \leq 150) & \approx Φ (\frac{150 + 0,5 - μ}{σ}) - Φ (\frac{100 - 0,5 - μ}{σ}) \\ \approx Φ (- 1,3718) - Φ (- 5,6993) \approx 1 - Φ (1,3718) \approx 0,0851 \end{matrix}

Die Werte $Φ (z)$ können numerisch bestimmt werden oder aus der Standardnormalverteilungstabelle 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

Gegeben sei eine binomialverteilte Zufallsvariable $X$ mit den Parametern $n = 48$ und $p = \frac{1}{4}$ , dann hat $X$ den Erwartungswert $n p = 12$ und die Varianz $n p (1 - p) = 9$ . Die Binomialverteilung wird durch eine Normalverteilung $𝒩 (μ, σ^{2})$ mit dem Mittelwert $μ = n p = 12$ und der Varianz $σ^{2} = n p (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 \leq X \leq 3)$ ?“ Da $3 = μ - 3 σ$ 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 \leq X \leq 3) = \sum_{k = 0}^{3} (\binom{48}{k}) \cdot (\frac{1}{4})^{k} \cdot (\frac{3}{4})^{48 - k} \approx 0,0007882

Abschätzung mit der Normalverteilung

\begin{matrix} P (0 \leq X \leq 3) & = P (X \leq 3) - P (X \leq - 1) \\ \approx Φ (\frac{3 - 12}{3}) - Φ (\frac{- 1 - 12}{3}) \\ = Φ (- 3) - Φ (- 13 / 3) \\ \approx 0,0013499 - 0,0000073 \approx 0,001343 \end{matrix}

Abschätzung mit Normalverteilung und Stetigkeitskorrektur

\begin{matrix} P (0 \leq X \leq 3) & \approx Φ (\frac{3 + 0,5 - 12}{3}) - Φ (\frac{0 - 0,5 - 12}{3}) \\ = Φ (- 8,5 / 3) - Φ (- 12,5 / 3) \\ \approx 0,002303 - 0.0000015 \approx 0,00229 \end{matrix}

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 $α \in {10 %, 5 %, 1 %}$ 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, doi:10.1515/9783110215274.
Ulrich Krengel: Einführung in die Wahrscheinlichkeitstheorie und Statistik. Vieweg, Braunschweig 1988, ISBN 978-3-528-07259-9, doi:10.1007/978-3-322-96418-2.

Einzelnachweise

↑ Horst Rinne: Taschenbuch der Statistik. 4. Auflage. Harri Deutsch, Frankfurt am Main 2008, ISBN 978-3-8171-1827-4, S. 304.
↑ Michael Sachs: Wahrscheinlichkeitsrechnung und Statistik für Ingenieurstudenten an Fachhochschulen. Fachbuchverlag Leipzig, München 2003, ISBN 3-446-22202-2, S. 129–130
↑ Christian Hassold, Sven Knoth, Detlef Steuer; Formelsammlung Statistik I & II. Beschreibende Statistik - Wahrscheinlichkeitsrechnung - Schließende Statistik; Hamburg 2010, S. 25 (Memento vom 9. Februar 2016 im Internet Archive), zuletzt abgerufen am 9. Februar 2016.
↑ K.Zirkelbach, W.Schmid; Kommentierte Formelsammlung Statistik I und II. Deskriptive Statistik - Wahrscheinlichkeitsrechnung; Frankfurt(Oder) 2008, S. 29.
↑ Formelsammlung zur Vorlesung Statistik I/II für Statistiker, Mathematiker und Informatiker (WS 08/09); LMU München 2008, S. 23, zuletzt abgerufen am 9. Februar 2016.
↑ Etwas "exakter": Zur Modellierung definiert man den Wahrscheinlichkeitsraum $(Ω, Σ, P)$ mit der Ergebnismenge $Ω := {0, \dots, 1000}$ der Anzahl der gewürfelten Sechsen. Die σ-Algebra ist dann kanonisch die Potenzmenge der Ergebnismenge $Σ := 𝒫 (Ω)$ und die Wahrscheinlichkeitsverteilung die Binomialverteilung $P ({k}) := B_{n, p} ({k})$ , wobei $n = 1000$ ist und $p = 1 / 6$ .

[1] Horst Rinne: Taschenbuch der Statistik. 4. Auflage. Harri Deutsch, Frankfurt am Main 2008, ISBN 978-3-8171-1827-4, S. 304.

[SACHS-2] Michael Sachs: Wahrscheinlichkeitsrechnung und Statistik für Ingenieurstudenten an Fachhochschulen. Fachbuchverlag Leipzig, München 2003, ISBN 3-446-22202-2, S. 129–130

[3] Christian Hassold, Sven Knoth, Detlef Steuer; Formelsammlung Statistik I & II. Beschreibende Statistik - Wahrscheinlichkeitsrechnung - Schließende Statistik; Hamburg 2010, S. 25 (Memento vom 9. Februar 2016 im Internet Archive), zuletzt abgerufen am 9. Februar 2016.

[4] K.Zirkelbach, W.Schmid; Kommentierte Formelsammlung Statistik I und II. Deskriptive Statistik - Wahrscheinlichkeitsrechnung; Frankfurt(Oder) 2008, S. 29.

[5] Formelsammlung zur Vorlesung Statistik I/II für Statistiker, Mathematiker und Informatiker (WS 08/09); LMU München 2008, S. 23, zuletzt abgerufen am 9. Februar 2016.

[6] Etwas "exakter": Zur Modellierung definiert man den Wahrscheinlichkeitsraum $(Ω, Σ, P)$ mit der Ergebnismenge $Ω := {0, \dots, 1000}$ der Anzahl der gewürfelten Sechsen. Die σ-Algebra ist dann kanonisch die Potenzmenge der Ergebnismenge $Σ := 𝒫 (Ω)$ und die Wahrscheinlichkeitsverteilung die Binomialverteilung $P ({k}) := B_{n, p} ({k})$ , wobei $n = 1000$ ist und $p = 1 / 6$ .

[1]

[2]

[3]

[4]

[5]

[6]