Warum funktioniert die CLT Arbeit für

Wir wissen also, dass eine Summe von Poisson mit dem $n$Parameter $\lambda$ selbst eine Poisson mit $n\lambda$ . So hypothetisch könnte man nehmen $x \sim poisson(\lambda = 1)$ und sagen , dass es tatsächlich ist $\sum_1^n x_i \sim poisson(\lambda = 1)$ , wobei jedes $x_i$ ist: $x_i \sim poisson(\lambda = 1/n)$ und nehmen Sie ein großes n, um CLT zum Laufen zu bringen.

Das funktioniert (offensichtlich) nicht. Ich nehme an, das hat etwas damit zu tun, wie CLT für Zufallsvariablen "schneller" arbeitet, die "näher" am Normalen liegen, und je kleiner Lambda ist, desto mehr erhalten wir eine Zufallsvariable, die meistens 0 ist und selten etwas anderes variiert.

Was ich jedoch erklärte, ist meine Intuition. Gibt es eine formellere Erklärung, warum dies der Fall ist?

Vielen Dank!

Ich stimme @whuber zu, dass die Wurzel der Verwirrung darin zu liegen scheint, die asymptotische Summierung in CLT durch eine Art Teilung in Ihrem Argument zu ersetzen. In CLT erhalten wir die feste Verteilung dann ziehen n Zahlen x i aus ihr und berechnen die Summe ˉ x n = $f(x,\lambda)$$n$$x_i$ . Wenn wirnweitererhöhen, passiert etwas Interessantes: $\bar x_n=\frac{1}{n}\sum_{i=1}^nx_i$$n$ wobeiμ,σ2 derMittelwert und die Varianz der Verteilungf(x) sind.

$$\sqrt n (\bar x_n-\mu)\rightarrow\mathcal{N}(0,\sigma^2)$$ $\mu,\sigma^2$$f(x)$

Was Sie mit Poisson vorschlagen, ist etwas rückwärts: Anstatt die Variablen aus einer festen Verteilung zu summieren , möchten Sie die feste Verteilung in sich ständig ändernde Teile aufteilen . Mit anderen Worten, Sie nehmen eine Variable x aus einer festen Verteilung f ( x , λ ) und teilen sie dann in x i, so dass n ∑ i = 1 x i ≡ x$x$$f(x,\lambda)$$x_i$

$$\sum_{i=1}^nx_i\equiv x$$

Was sagt CLT über diesen Prozess? Nichts. Beachten Sie, wie wir in CLT jemals und seine sichänderndeVerteilungfn(x), die zu einerfestenVerteilungN(0,σ2)konvergiert$\sqrt n(\bar x_n-\mu)$$f_n(x)$$\mathcal{N}(0,\sigma^2)$

In Ihrem Setup ändern sich weder die Summe noch deren Verteilung f ( x , λ ) ! Sie sind repariert. Sie ändern sich nicht, sie laufen zu nichts zusammen. CLT hat also nichts zu sagen.$x$$f(x,\lambda)$

Außerdem sagt CLT nichts über die Anzahl der Elemente in der Summe aus. Sie können eine Summe von 1000 Variablen aus Poisson (0,001) haben, und CLT sagt nichts über die Summe aus. Es heißt nur, wenn Sie N weiter erhöhen, sieht diese Summe irgendwann aus wie eine Normalverteilung . In der Tat, wenn N = 1.000.000, erhalten Sie die enge Annäherung an die Normalverteilung.$\frac{1}{N}\sum_{i=1}^N x_i, x_i\sim Poisson(0.001)$

Ihre Intuition stimmt nur in Bezug auf die Anzahl der Elemente in der Summe, dh wenn sich die Anfangsverteilung von der normalen Verteilung unterscheidet, müssen Sie mehr Elemente summieren, um zur normalen Verteilung zu gelangen. Je mehr formale (aber informell noch) Art und Weise , indem man die charakteristische Funktion von Poisson würde: Wenn Sie & lgr; > > 1 , erhalten Sie mit der Taylor - Entwicklung (WRT t ) des verschachtelten Exponenten: ≈ exp ( i λ t - λ / 2 t

$$\exp(\lambda (\exp(it)-1))$$ $\lambda>>1$$t$ Dies ist die charakteristische Funktion der Normalverteilung N ( λ , λ 2 $$\approx\exp(i\lambda t-\lambda/2t^2)$$ $\mathcal{N}(\lambda,\lambda^2)$

Ihre Intuition wird jedoch nicht richtig angewendet: Wenn Sie die Summe in CLT durch eine Art Division verschieben, wird die Sache durcheinander gebracht und CLT wird nicht anwendbar.

Das Problem bei Ihrem Beispiel ist, dass Sie zulassen, dass sich die Parameter ändern, wenn sich ändert. Die CLT sagt Ihnen, dass für eine feste Verteilung mit einem endlichen Mittelwert und sd, als n → ∞ ,$n$$n \rightarrow \infty$

,$\frac {\sum x - \mu} {\sqrt n} \rightarrow_d N(0, \sigma)$

wobei und σ aus dem Mittelwert und sd der Verteilung von x stammen .$\mu$$\sigma$$x$

Natürlich sind für verschiedene Verteilungen (z. B. höhere Versetzungen) größere erforderlich, bevor die aus diesem Theorem abgeleitete Approximation vernünftig wird. In Ihrem Beispiel für λ m = 1 / m , ein n > > m erforderlich ist , bevor die normale Annäherung sinnvoll ist.$n$$\lambda_m = 1/m$$n >> m$

BEARBEITEN

Es wird diskutiert, wie die CLT nicht für Summen gilt, sondern für standardisierte Summen (dh nicht∑xi). Theoretisch ist dies natürlich richtig: Die nicht standardisierte Summe wird in den meisten Fällen eine undefinierte Verteilung haben.$\sum x_i / \sqrt n$$\sum x_i$

In der Praxis können Sie die vom CLT begründete Annäherung jedoch durchaus auf Summen anwenden! Wenn durch eine normale CDF für großes n angenähert werden kann , dann kann F ∑ x sicherlich auch, da die Multiplikation mit einem Skalar die Normalität bewahrt. Und das können Sie gleich in diesem Problem sehen: Denken Sie daran, dass, wenn X i ∼ P o i s ( λ ) , dann Y = ∑ n i = 1 X i ∼ P o i s ( n λ $F_{\bar x}$$n$$F_{\sum x}$$X_i \sim Pois(\lambda)$$Y = \sum_{i = 1}^n X_i \sim Pois(n\lambda)$. Und wir alle haben in unserem oberen Teilungswahrscheinlichkeitskurs gelernt, dass für großes die CDF eines P o i s ( λ ) recht gut durch eine Normale mit μ = λ , σ 2 = λ angenähert werden kann . Also für jeden festen λ , können wir die CDF von annähernd Y ~ P o i s ( n λ ) recht gut mit Φ ( y - n $\lambda$$Pois(\lambda)$$\mu = \lambda$$\sigma^2 = \lambda$ $\lambda$$Y \sim Pois(n\lambda)$für ein ausreichend großesn,wenn& lgr;>0(eine Näherung kann trivial angewendet werden, wenn& lgr;=0 ist, aber nicht die Berechnung der CDF, wie ich sie geschrieben habe).$\Phi( \frac{y - n\lambda}{\sqrt{n\lambda} })$$n$$\lambda > 0$$\lambda = 0$

Während die CLT nicht ohne weiteres für Summen gilt, gilt die auf der CLT basierende Annäherung mit Sicherheit. Ich glaube, darauf bezog sich das OP, als es darum ging, die CLT auf die Summe anzuwenden.

Ich behaupte, die Frage ist interessanter, wenn man allgemeiner darüber nachdenkt, wie die Verteilung des übergeordneten Poisson von abhängt , beispielsweise mit den Parametern λ n und λ n = 1 als Sonderfall. Ich halte es für vernünftig zu fragen, warum und wie wir das verstehen können, dass ein zentraler Grenzwertsatz für die Summe S n = ∑ n i = 1 X i , n nicht gilt . Schließlich ist es üblich, eine CLT auch bei Problemen anzuwenden, bei denen die Verteilungen der Komponenten der Summe von n abhängen$n$$\lambda_n$$\lambda_n = 1$$S_n = \sum_{i=1}^n X_{i,n}$$n$. Es ist auch üblich, Poisson-Verteilungen als Verteilung einer Summe von Poisson-Variablen zu zerlegen und dann eine CLT anzuwenden.

$X_{i, n}$$n$$S_n$$n$$S_n \sim Poi(n)$$Poi(\lambda_n)$

$S_n \sim Poi(1)$$n$$S_n$

$s_n^2 = \mathrm{Var(S_n)}$$\forall \epsilon >0$

$$\frac{1}{s_n^2}\sum_{i=1}^n\mathbb E[X_{i,n} - 1/n]^2I(\vert X_{i,n} - 1/n \vert >\epsilon s_n) \to 0,n\to\infty$$

Nun, für den vorliegenden Fall, stirbt die Varianz der Terme in der Summe so schnell ab $n$$s_n = 1$$n$$n$$X_{i,n}$

$$n\mathbb E[X_{1,n} - 1/n]^2I(\vert X_{1,n} - 1/n \vert >\epsilon) \to 0.$$

$\epsilon$$n$

$$\begin{align} n\mathbb E[X_{1,n} - 1/n]^2I(\vert X_{1,n} - 1/n \vert >\epsilon) &>n\epsilon^2P(X_{1,n}>0) \\ &=\epsilon^2n[1 - e^{-1/n}] \\ &= \epsilon^2n[1-(1 - 1/n + o(1/n))] \\ &= \epsilon^2 + o(1), \end{align}$$

$S_n$$n$$n$