Featured image of post Warum kleine Zufälle am Anfang zu massiver Ungleichheit führen: Das Pólya-Urnenmodell und verstärkte stochastische Prozesse

Warum kleine Zufälle am Anfang zu massiver Ungleichheit führen: Das Pólya-Urnenmodell und verstärkte stochastische Prozesse

Eine mathematische Analyse des Pólya-Urnenmodells: Austauschbarkeit, Beta-Binomialverteilung, Doobscher Martingalkonvergenzsatz zur Beta(r,b)-Verteilung, Vergleich mit dem starken Gesetz der großen Zahlen, Satz von de Finetti, Pfadabhängigkeit nach W. Brian Arthur, Preferential Attachment und Dirichlet-Prozesse.

Einleitung: Warum kristallisieren sich winzige Anfangszufälle zu unumkehrbaren Asymmetrien?

Zwei Videos oder Standards gleicher Qualität starten gleichzeitig: Durch einen winzigen Zufallsvorsprung gewinnt eines durch positives Feedback die Marktdominanz (Matthäus-Effekt), während das andere verschwindet. Das starke Gesetz der großen Zahlen (SLLN) lehrt für unabhängige (i.i.d.) Münzwürfe das Gegenteil: Anfangsschwankungen waschen sich für \(n \to \infty\) aus und der Anteil konvergiert fast sicher gegen die Konstante \(p\).

Reale Prozesse sind oft verstärkte stochastische Prozesse, bei denen das Eintreten eines Ereignisses die Wahrscheinlichkeit seiner Wiederholung erhöht. 1923 führten George Pólya und Florian Eggenberger das Pólya-Urnenmodell ein, das zwei Paradoxa zeigt:

  1. Erstes Paradoxon (Austauschbarkeit): Obwohl sich der Urneninhalt nach jedem Zug ändert, bleibt \(\mathbb{P}(I_n = 1) = \frac{r}{r+b}\) für alle \(n\) konstant, und die Wahrscheinlichkeit einer Folge hängt nur von der Anzahl roter und blauer Kugeln ab, nicht von ihrer Reihenfolge (Exchangeability).
  2. Zweites Paradoxon (Zufälliger Grenzwert): Für \(n \to \infty\) konvergiert der rote Anteil \(X_n\) fast sicher gegen eine Zufallsvariable \(X_\infty \sim \mathrm{Beta}(r, b)\). Für \(r=1, b=1\) ist \(X_\infty \sim \mathrm{Uniform}(0, 1)\) gleichverteilt auf \([0, 1]\)!

Kapitel 1: Drei Urnen-Universen — Erschöpfung, Gedächtnislosigkeit und Pólya-Verstärkung

Bei \(n = 0\) enthalte eine Urne \(R_0 = r\) rote und \(B_0 = b\) blaue Kugeln (\(N_0 = r + b\)). In jedem Schritt wird eine Kugel gezogen und mit \(c\) Kugeln derselben Farbe zurückgelegt: 1. Ohne Zurücklegen (\(c = -1\)) (negatives Feedback, hypergeometrisch), 2. Mit Zurücklegen (\(c = 0\)) (i.i.d. \(\mathrm{Binomial}(n, \frac{r}{r+b})\)), 3. Pólya-Urne (\(c = +1\)) (positives Feedback: gezogene Kugel plus \(1\) Kugel gleicher Farbe).

  flowchart TD
    Start["Start (n = 0): Rot r, Blau b, X₀ = r/(r+b)"]
    Draw["Schritt n+1: 1 Kugel ziehen, P(Iₙ₊₁ = 1 | Fₙ) = Xₙ"]
    BranchW["1. Ohne Zurücklegen (c = -1)"]
    BranchB["2. Mit Zurücklegen (c = 0)"]
    BranchP["3. Polya-Urne (c = +1)"]
    LimW["Endet bei n = r+b mit r/(r+b)"]
    LimB["SLLN: Xₙ → r/(r+b) f.s."]
    LimP["Doob: Xₙ → X_∞ ~ Beta(r, b) f.s."]
    Start --> Draw
    Draw --> BranchW --> LimW
    Draw --> BranchB --> LimB
    Draw --> BranchP --> LimP

Mit \(I_n \in \{0, 1\}\) und \(S_n = \sum_{k=1}^n I_k\) gilt \(R_n = r + S_n\) und \(N_n = r + b + n\):

\[X_n := \frac{R_n}{N_n} = \frac{r + S_n}{r + b + n}, \qquad \mathbb{P}(I_{n+1} = 1 \mid \mathcal{F}_n) = X_n, \qquad |X_{n+1} - X_n| \le \frac{1}{r + b + n + 1}\]
Eigenschaft1. Ohne Zurücklegen (\(c = -1\))2. Mit Zurücklegen (\(c = 0\))3. Pólya-Urne (\(c = +1\))
RückkopplungNegativ (Erschöpfung)Keine (gedächtnislos)Positiv (Selbstverstärkung)
AbhängigkeitNegativ korreliertUnabhängig (i.i.d.)Unendlich austauschbar
Verteilung \(S_n\)\(\mathrm{Hypergeom}(N_0, r, n)\)\(\mathrm{Binomial}(n, \frac{r}{r+b})\)\(\mathrm{BetaBinom}(n, r, b)\)
Grenzwert \(X_\infty\)Konstante \(\frac{r}{r+b}\)Konstante \(\frac{r}{r+b}\) f.s.Zufallsvariable \(\mathrm{Beta}(r, b)\) f.s.

Kapitel 2: Austauschbarkeit und Beta-Binomialverteilung

Für jede Folge \((i_1, \dots, i_n) \in \{0, 1\}^n\) mit \(k = \sum_{j=1}^n i_j\) roten Zügen ist das Nennerprodukt stets \((r+b)^{(n)}\) und das Zählerprodukt kommutativ stets \(r^{(k)} b^{(n-k)}\):

\[\mathbb{P}(I_1 = i_1, \dots, I_n = i_n) = \frac{r^{(k)} \, b^{(n-k)}}{(r+b)^{(n)}}, \qquad \mathbb{P}(S_n = k) = \binom{n}{k} \frac{\mathrm{B}(r+k, \, b+n-k)}{\mathrm{B}(r, \, b)}\]

Für \(r = 1, b = 1\) kürzen sich alle Fakultäten zur diskreten Gleichverteilung \(\mathbb{P}(S_n = k) = \frac{1}{n + 1}\) (\(k = 0, 1, \dots, n\)).


Kapitel 3: Doobscher Martingalkonvergenzsatz und \(\mathrm{Beta}(r,b)\)-Grenzwert

Wegen \(\mathbb{E}[R_{n+1} \mid \mathcal{F}_n] = R_n + X_n = (N_n + 1)X_n\) ist \(X_n\) ein beschränktes Martingal auf \([0, 1]\):

\[\mathbb{E}[X_{n+1} \mid \mathcal{F}_n] = \frac{R_n + \mathbb{E}[I_{n+1} \mid \mathcal{F}_n]}{N_n + 1} = \frac{(N_n + 1) X_n}{N_n + 1} = X_n\]

Nach dem Doobschen Martingalkonvergenzsatz gilt \(X_n \xrightarrow{\text{f.s.}} X_\infty\). Mit \(\frac{\Gamma(z+\alpha)}{\Gamma(z+\beta)} \sim z^{\alpha-\beta}\) für \(k = \lfloor np \rfloor\) folgt die Dichte der Betaverteilung \(\mathrm{Beta}(r, b)\):

\[\boxed{f_{X_\infty}(p) = \frac{\Gamma(r+b)}{\Gamma(r)\Gamma(b)} \, p^{r-1} (1-p)^{b-1} \quad (0 < p < 1), \qquad \mathbb{E}[X_\infty] = \frac{r}{r+b}, \qquad \mathrm{Var}(X_\infty) = \frac{rb}{(r+b)^2(r+b+1)} > 0}\]

Kapitel 4: Satz von de Finetti, Phasenübergänge und Anwendungen

  1. Darstellungssatz von de Finetti (1931): Jede unendlich austauschbare Folge ist eine Mischung von i.i.d. Bernoulli-Versuchen bezüglich des Prior-Maßes \(\mathrm{Beta}(r,b)\): \(\mathbb{P}(I_1 = i_1, \dots, I_n = i_n) = \int_0^1 p^k (1-p)^{n-k} \, \frac{\Gamma(r+b)}{\Gamma(r)\Gamma(b)} p^{r-1}(1-p)^{b-1} \, dp\)
  2. Nichtlineare Urnen (\(f(x) = \frac{x^\gamma}{x^\gamma + (1-x)^\gamma}\)): \(X_{n+1} - X_n = \frac{1}{N_n+1}[(f(X_n)-X_n) + \xi_{n+1}]\) liefert drei Phasen:
  flowchart TD
    Root["Nichtlineare Urne: f(x) = x^γ / (x^γ + (1-x)^γ)"]
    Sub["1. Sublinear (0 < γ < 1)"]
    Lin["2. Linear (γ = 1): g(x) ≡ 0"]
    Sup["3. Superlinear (γ > 1)"]
    ResSub["Koexistenz: Xₙ → 1/2 f.s."]
    ResLin["Pfadabhängig: Xₙ → Beta(r, b) f.s."]
    ResSup["Winner-Take-All: P(X_∞ ∈ {0, 1}) = 1"]
    Root --> Sub --> ResSub
    Root --> Lin --> ResLin
    Root --> Sup --> ResSup
  1. Anwendungen & Grenzen: W. Brian Arthurs Lock-In (1989) erklärt Monopole bei \(\gamma > 1\) (QWERTY, VHS); Preferential Attachment (Yule–Simon / Barabási–Albert) erzeugt \(P(k) \sim k^{-(2 + \frac{1}{1-\alpha})}\); neue Farben mit Gewicht \(\theta > 0\) aus \(H_0\) führen zum Dirichlet-Prozess \(\mathrm{DP}(\theta, H_0)\) (CRP); MusicLab (Salganik et al. 2006) zeigt vier Grenzen: Fitness \(\eta_i\), Sättigung, Gedächtnis und Netzwerk.
comments powered by Disqus