Mathematics · Boek 4 · Bachelor Year 2

Universitaire wiskunde — Bachelor jaar 2

Universitaire wiskunde — Bachelor jaar 2 · Bachelor Year 2

21Kansrekening op aftelbare ruimten

De laatste drie hoofdstukken ontwikkelen de moderne kansrekening: kansmaten op aftelbare uitkomstenruimten, discrete toevalsveranderlijken en genererende functies. De eindige theorie van het bovenbouwvolume krijgt haar volledige infrastructuur: σ\sigma-additiviteit vervangt de eindige additiviteit, en de machinerie van de sommeerbare families uit Hoofdstuk 7 is precies wat oneindige uitkomstenruimten hanteerbaar maakt. De kernresultaten hier zijn de continuïteit van de kans langs monotone rijen gebeurtenissen en het lemma van Borel–Cantelli.

21.1 Kansruimten

Definitie 21.1 (Aftelbare kansruimte)

Zij Ω\Omega een niet-lege eindige of aftelbare verzameling (de uitkomstenruimte). Een kansmaat op Ω\Omega is een afbeelding P\P van de verzameling P(Ω)\mathcal{P}(\Omega) van alle deelverzamelingen van Ω\Omega (gebeurtenissen) naar [0,1][0, 1] zodanig dat:

  1. P(Ω)=1\P(\Omega) = 1;
  2. (σ\sigma-additiviteit) voor elke rij (An)nN(A_n)_{n\in\N} paarsgewijs disjuncte gebeurtenissen geldt

    P(nNAn)=n=0P(An).\P\Bigl(\,\bigcup_{n \in \N} A_n\Bigr) = \sum_{n=0}^{\infty} \P(A_n) .

Het paar (Ω,P)(\Omega, \P) is een (aftelbare) kansruimte.

Opmerking 21.2

Op een aftelbare Ω\Omega mogen wij alle deelverzamelingen als gebeurtenissen nemen; op overaftelbare ruimten (zoals nodig voor continue modellen in bachelorjaar 3) kan dat niet meer, en beperkt men P\P tot een geschikte collectie gebeurtenissen, een σ\sigma-algebra. Alle formules van dit hoofdstuk overleven die veralgemening woordelijk.

Propositie 21.3 (Elementaire regels)

Voor gebeurtenissen A,BA, B en een kansmaat P\P geldt: P()=0\P(\emptyset) = 0; P\P is eindig additief; P(Ac)=1P(A)\P(A^c) = 1 - \P(A); is ABA \subseteq B, dan P(A)P(B)\P(A) \leq \P(B); en

P(AB)=P(A)+P(B)P(AB).\P(A \cup B) = \P(A) + \P(B) - \P(A \cap B) .

Bewijs. De σ\sigma-additiviteit toepassen op A0=ΩA_0 = \Omega, An=A_n = \emptyset (n1n \geq 1) geeft 1=1+n1P()1 = 1 + \sum_{n\geq1}\P(\emptyset), dus P()=0\P(\emptyset) = 0; een eindige disjuncte vereniging met lege verzamelingen opvullen geeft dan de eindige additiviteit. De rest volgt als in het eindige geval (bovenbouwvolume): 1=P(A)+P(Ac)1 = \P(A) + \P(A^c) uit Ω=AAc\Omega = A \sqcup A^c; P(B)=P(A)+P(BA)P(A)\P(B) = \P(A) + \P(B \setminus A) \geq \P(A) wanneer ABA \subseteq B; en ontbinden in drie disjuncte stukken geeft

P(AB)=P(AB)+P(BA)+P(AB)=(P(A)P(AB))+(P(B)P(AB))+P(AB),\begin{align*} \P(A \cup B) &= \P(A \setminus B) + \P(B \setminus A) + \P(A \cap B)\\ &= \bigl(\P(A) - \P(A\cap B)\bigr) + \bigl(\P(B) - \P(A\cap B)\bigr) + \P(A \cap B), \end{align*}

wat de in- en uitsluiting is; de algemene versie met nn verzamelingen is Oefening 21.4.

Propositie 21.4 (Verdelingen op een aftelbare ruimte)

Een kansmaat geven op een aftelbare Ω={ω0,ω1,}\Omega = \{\omega_0, \omega_1, \dots\} komt precies neer op het geven van gewichten pi=P({ωi})0p_i = \P(\{\omega_i\}) \geq 0 met ipi=1\sum_i p_i = 1; dan is voor elke AΩA \subseteq \Omega

P(A)=ωAP({ω}),\P(A) = \sum_{\omega \in A} \P(\{\omega\}) ,

een (absoluut convergente) deelsom van de familie (pi)(p_i).

Bewijs. Gegeven P\P vormen de singletons {ω}\{\omega\}, ωA\omega \in A, een aftelbare disjuncte overdekking van AA, dus dwingt de σ\sigma-additiviteit af dat

P(A)=ωAP({ω}),\P(A) = \sum_{\omega\in A}\P(\{\omega\}),

een onvoorwaardelijke deelsom van de niet-negatieve sommeerbare familie (pi)(p_i) — herschikken is onschadelijk precies omdat de termen niet-negatief zijn (Hoofdstuk 7); in het bijzonder is ipi=P(Ω)=1\sum_ip_i = \P(\Omega) = 1. Omgekeerd, gegeven niet-negatieve gewichten met totale som 11, definieer P(A)=ωApω\P(A) = \sum_{\omega \in A}p_\omega: de familie is sommeerbaar, en de σ\sigma-additiviteit is precies de stelling over sommeren in pakketten uit Hoofdstuk 7, toegepast op de partitie van An\bigcup A_n in de AnA_n.

Voorbeeld 21.5 (Meetkundig model: wachten op de eerste kop)

Werp herhaaldelijk met een munt met kanskop p(0,1)p \in \intoo{0}{1}, en laat Ω=N{}\Omega = \N^* \cup \{\infty\} de rang van de eerste kop registreren. De natuurlijke gewichten zijn

P({k})=(1p)k1p(kN),P({})=0,\P(\{k\}) = (1 - p)^{k-1}p \quad (k \in \N^*), \qquad \P(\{\infty\}) = 0 ,

een kansmaat omdat k1(1p)k1p=p1(1p)=1\sum_{k\geq1}(1-p)^{k-1}p = \frac{p}{1 - (1-p)} = 1: met kans 11 eindigt het spel — maar de uitkomstenruimte moet toch de mogelijkheid bevatten dat het niet eindigt. De aftelbare additiviteit is wat ons toelaat te beweren dat P(het spel eindigt)=kP({k})\P(\text{het spel eindigt}) = \sum_k \P(\{k\}).

Stelling 21.6 (Monotone continuïteit)

Zij (An)(A_n) een rij gebeurtenissen.

  1. Is AnAn+1A_n \subseteq A_{n+1} voor alle nn (stijgend), dan is P(nAn)=limnP(An)\P\bigl(\bigcup_n A_n\bigr) = \lim_{n\to\infty} \P(A_n).
  2. Is AnAn+1A_n \supseteq A_{n+1} voor alle nn (dalend), dan is P(nAn)=limnP(An)\P\bigl(\bigcap_n A_n\bigr) = \lim_{n\to\infty} \P(A_n).

Bewijs. 1. Maak disjunct: zij B0=A0B_0 = A_0 en Bn=AnAn1B_n = A_n \setminus A_{n-1}. De BnB_n zijn paarsgewijs disjunct met knBk=An\bigcup_{k \leq n} B_k = A_n en nBn=nAn\bigcup_n B_n = \bigcup_n A_n. Volgens de σ\sigma-additiviteit en de eindige additiviteit is

P(nAn)=n=0P(Bn)=limNn=0NP(Bn)=limNP(AN).\P\Bigl(\bigcup_n A_n\Bigr) = \sum_{n=0}^\infty \P(B_n) = \lim_{N\to\infty}\sum_{n=0}^N \P(B_n) = \lim_{N\to\infty}\P(A_N) .

2. Ga over op de complementen: (Anc)(A_n^c) is stijgend met vereniging (An)c\bigl(\bigcap A_n\bigr)^c; pas deel 1 toe: 1P(An)=lim(1P(An))1 - \P(\bigcap A_n) = \lim (1 - \P(A_n)).

Gevolg 21.7 (Aftelbare subadditiviteit)

Voor elke rij gebeurtenissen is P(nAn)n=0P(An)\P\bigl(\bigcup_n A_n\bigr) \leq \sum_{n=0}^\infty \P(A_n).

Bewijs. De eindige subadditiviteit P(A0AN)0NP(An)\P(A_0 \cup \dots \cup A_N) \leq \sum_0^N \P(A_n) volgt met inductie uit de in- en uitsluiting (of uit de additiviteit over de disjunct gemaakte BnAnB_n \subseteq A_n). Laat NN \to \infty: het linkerlid convergeert naar P(nAn)\P(\bigcup_n A_n) wegens de monotone continuïteit toegepast op de stijgende rij CN=A0ANC_N = A_0 \cup \dots \cup A_N.

Voorbeeld 21.8 (De somgrens: ruw maar onverwoestbaar)

De subadditiviteit met eindig veel gebeurtenissen — de somgrens — ruilt nauwkeurigheid in voor universaliteit. Voor het verjaardagsprobleem met 2323 mensen geeft de kans op een botsing begrenzen door de som over de paren

P(botsing)(232)1365=2533650.693,\P(\text{botsing}) \leq \binom{23}2\cdot\frac1{365} = \frac{253}{365} \approx 0.693 ,

tegenover de werkelijke 0.5070.507: er ruim naast, omdat botsingen elkaar overlappen. Toch heeft de grens geen onafhankelijkheid nodig, geen gezamenlijke verdeling, niets dan de kansen op de paren — en daarom is in de weekendopgave en overal in Hoofdstuk 22 de somgrens het eerst getrokken wapen: blijkt zij klein te zijn, dan is de zaak beslecht zonder verdere modellering.

Voorbeeld 21.9 (Een zes komt, uiteindelijk)

Gooi eeuwig met een eerlijke dobbelsteen en zij Bn=B_n = {}“minstens één zes onder de eerste nn worpen”, een stijgende rij gebeurtenissen met P(Bn)=1(5/6)n\P(B_n) = 1 - (5/6)^n. De monotone continuïteit geeft

P(er verschijnt uiteindelijk een zes)=P(nBn)=limn(1(5/6)n)=1.\P(\text{er verschijnt uiteindelijk een zes}) = \P\Bigl(\bigcup_nB_n\Bigr) = \lim_n\bigl(1 - (5/6)^n\bigr) = 1 .

Het punt is niet de (voor de hand liggende) limiet maar de logische stap: “uiteindelijk” is een gebeurtenis over oneindig veel worpen, buiten het bereik van de eindige additiviteit, en de monotone continuïteit — dat wil zeggen de σ\sigma-additiviteit — is precies het axioma dat er een kans aan toekent. Elke uitspraak “bijna zeker” in de rest van dit boek gaat door deze ene smalle deur.

21.2 Voorwaardelijkheid en onafhankelijkheid

Definitie 21.10 (Voorwaardelijke kans)

Voor gebeurtenissen A,BA, B met P(B)>0\P(B) > 0 is de voorwaardelijke kans van AA gegeven BB

P(AB)=P(AB)P(B).\P(A \mid B) = \frac{\P(A \cap B)}{\P(B)} .

De afbeelding AP(AB)A \mapsto \P(A \mid B) is zelf een kansmaat op Ω\Omega.

Opmerking 21.11

Dat APB ⁣(A)A \mapsto \pcond BA opnieuw een kansmaat is, verdient een ogenblik: PB ⁣(Ω)=1\pcond B\Omega = 1 en de σ\sigma-additiviteit gaan door het quotiënt heen omdat doorsnijden met BB disjuncte verenigingen eerbiedigt. Het praktische gevolg: elke identiteit van dit hoofdstuk — in- en uitsluiting, monotone continuïteit, Borel–Cantelli — mag na het voorwaardelijk stellen worden toegepast, zonder nieuwe bewijzen. Kansrekenaars “werken onder PB ⁣()\pcond B{\cdot}” precies om die reden voortdurend.

Voorbeeld 21.12 (Voorwaardelijk stellen kan uniformiteit scheppen)

Gooi met twee eerlijke dobbelstenen en stel voorwaardelijk dat de som 77 is: voor elke k[ ⁣[1,6] ⁣]k \in \intint16 is

P{S=7} ⁣(X=k)=P(X=k, Y=7k)P(S=7)=1/366/36=16:\pcond{\{S = 7\}}{X = k} = \frac{\P(X = k,\ Y = 7 - k)}{\P(S = 7)} = \frac{1/36}{6/36} = \frac16 :

gegeven een som van 77 is de eerste dobbelsteen precies uniform — 77 is de enige som die met elk aantal ogen verenigbaar is, dus wist het voorwaardelijk stellen alle informatie over XX uit. Elke andere som scheeftrekt de verdeling (gegeven S=4S = 4 is de eerste dobbelsteen alleen uniform op {1,2,3}\{1, 2, 3\}). Een voorwaardelijke verdeling berekenen betekent de gezamenlijke gewichten langs de voorwaardelijke gebeurtenis hernormaliseren, meer niet.

Voorbeeld 21.13 (De tweede trekking is even goed als de eerste)

Een urne bevat 33 witte en 22 zwarte ballen; trek er twee zonder teruglegging. Iedereen is het erover eens dat P(W1)=35\P(W_1) = \frac35; wat is P(W2)\P(W_2)? Totale kans langs de eerste trekking:

P(W2)=PW1 ⁣(W2)P(W1)+PB1 ⁣(W2)P(B1)=2435+3425=1220=35:\P(W_2) = \pcond{W_1}{W_2}\,\P(W_1) + \pcond{B_1}{W_2}\,\P(B_1) = \frac24\cdot\frac35 + \frac34\cdot\frac25 = \frac{12}{20} = \frac35 :

precies P(W1)\P(W_1). Er was geen berekening nodig: wegens de symmetrie heeft elke bal evenveel kans om als tweede te worden getrokken, dus heeft de tweede trekking — onvoorwaardelijk — dezelfde verdeling als de eerste. Voorwaardelijk stellen op de eerste uitkomst verandert de kansen; haar niet kennen niet. Dit argument met verwisselbaarheid keert in het volgende hoofdstuk terug voor het steekproeven zonder teruglegging, waar het het hypergeometrische gemiddelde npnp geeft, zonder enige binomiale identiteit.

Stelling 21.14 (Samengestelde kansen, totale kans, Bayes)

  1. (Kettingregel) Is P(A1An1)>0\P(A_1 \cap \dots \cap A_{n-1}) > 0, dan

    P(A1An)=P(A1)P(A2A1)P(AnA1An1).\P(A_1 \cap \dots \cap A_n) = \P(A_1)\,\P(A_2 \mid A_1)\cdots \P(A_n \mid A_1 \cap \dots \cap A_{n-1}) .
  2. (Totale kans) Is (Bi)iI(B_i)_{i \in I} een eindige of aftelbare partitie van Ω\Omega met P(Bi)>0\P(B_i) > 0, dan is voor elke gebeurtenis AA

    P(A)=iIP(ABi)P(Bi).\P(A) = \sum_{i \in I} \P(A \mid B_i)\,\P(B_i) .
  3. (Bayes) Onder dezelfde hypothesen geldt, als bovendien P(A)>0\P(A) > 0:

    P(BjA)=P(ABj)P(Bj)iIP(ABi)P(Bi).\P(B_j \mid A) = \frac{\P(A \mid B_j)\,\P(B_j)} {\sum_{i \in I} \P(A \mid B_i)\,\P(B_i)} .

Bewijs. 1. Schrijf elke voorwaardelijke kans als een quotiënt: het rechterlid is

P(A1)P(A1A2)P(A1)P(A1A2A3)P(A1A2)P(A1An)P(A1An1),\P(A_1)\cdot\frac{\P(A_1 \cap A_2)}{\P(A_1)}\cdot \frac{\P(A_1 \cap A_2 \cap A_3)}{\P(A_1 \cap A_2)}\cdots \frac{\P(A_1 \cap \dots \cap A_n)}{\P(A_1 \cap \dots \cap A_{n-1})},

een telescoperend product: elke noemer heft de voorgaande teller op, en P(A1An)\P(A_1 \cap \dots \cap A_n) blijft over. Alle noemers zijn P(A1An1)>0\geq \P(A_1 \cap \dots \cap A_{n-1}) > 0 wegens de monotonie, dus verdwijnt er niets. (De hypothese bewaakt precies dit: voorwaardelijk stellen op een gebeurtenis met kans nul is niet gedefinieerd.) 2. De verzamelingen ABiA \cap B_i zijn paarsgewijs disjunct met vereniging AA; pas de (σ\sigma-)additiviteit en de definitie van het voorwaardelijk stellen toe. 3. Beide leden van P(BjA)P(A)=P(ABj)P(Bj)\P(B_j \mid A)\P(A) = \P(A \mid B_j)\P(B_j) zijn gelijk aan P(ABj)\P(A \cap B_j); deel door P(A)\P(A) en werk P(A)\P(A) uit met de totale kans.

Voorbeeld 21.15 (De verjaardagsbotsing, met de kettingregel)

Met nn mensen wier verjaardagen onafhankelijk en uniform over 365365 dagen verdeeld zijn, zij Dn=D_n = {}“alle nn verjaardagen verschillen”. Persoon voor persoon voorwaardelijk stellen (kettingregel) geeft

P(Dn)=k=1n1(1k365),\P(D_n) = \prod_{k=1}^{n-1}\Bigl(1 - \frac{k}{365}\Bigr),

waarbij elke nieuwe persoon de kk reeds bezette dagen moet vermijden. Voor n=23n = 23: P(D23)0.493\P(D_{23}) \approx 0.493 — een gedeelde verjaardag is al waarschijnlijker dan niet. De heuristiek die de kleinheid van 2323 verklaart: logaritmen nemen geeft lnP(Dn)k<nk365=(n2)365-\ln \P(D_n) \approx \sum_{k<n}\frac k{365} = \frac{\binom n2}{365}, en (232)=253\binom{23}2 = 253 geeft 253/3650.693ln2253/365 \approx 0.693 \approx \ln 2. Wat telt is het aantal paren, dat kwadratisch groeit: botsingsproblemen leven op de schaal n365n \sim \sqrt{365}, niet n365n \sim 365 — de verjaardagsparadox is een vierkantswortel in vermomming.

Voorbeeld 21.16 (Monty Hall, met Bayes)

Achter een van drie deuren zit uniform verdeeld een prijs. Je kiest deur 11; de presentator, die weet waar de prijs is, opent een van de andere deuren, altijd een lege (uniform kiezend wanneer hij de keuze heeft), zeg deur 33. Zij Bi=B_i = {}“prijs achter deur ii” en A=A = {}“de presentator opent deur 33”. Dan is PB1 ⁣(A)=12\pcond{B_1}{A} = \frac12, PB2 ⁣(A)=1\pcond{B_2}{A} = 1 en PB3 ⁣(A)=0\pcond{B_3}{A} = 0, dus geeft Bayes (Stelling 21.14)

P(B2A)=1131213+113+013=23:\P(B_2 \mid A) = \frac{1\cdot\frac13} {\frac12\cdot\frac13 + 1\cdot\frac13 + 0\cdot\frac13} = \frac23 :

van deur wisselen wint twee keer op drie. De berekening lokaliseert de populaire verwarring precies: de zet van de presentator is informatief (hij kon deur 22 niet openen als de prijs daar zat), en de formule van Bayes is het boekhoudkundige hulpmiddel dat die asymmetrie in de 23\frac23 omzet. Voorwaardelijk stellen op “wat werd gezien” in plaats van op “wat waar is” is de hele kunst van de formule.

Voorbeeld 21.17 (De twee weddenschappen van de Chevalier de Méré)

Twee weddenschappen uit de zeventiende eeuw, beslecht door onafhankelijkheid. Weddenschap één: minstens één zes in 44 worpen met een dobbelsteen,

P=1(56) ⁣40.518>12.\P = 1 - \Bigl(\frac56\Bigr)^{\!4} \approx 0.518 > \frac12 .

Weddenschap twee: minstens één dubbele zes in 2424 worpen met twee dobbelstenen,

P=1(3536) ⁣240.491<12.\P = 1 - \Bigl(\frac{35}{36}\Bigr)^{\!24} \approx 0.491 < \frac12 .

De Méré redeneerde dat 2424 worpen met kans 136\frac1{36} zouden moeten overeenkomen met 44 worpen met kans 16\frac16 (dezelfde verhouding 2436=46\frac{24}{36} = \frac46); het falen van deze evenredigheid — kansen op verenigingen schalen niet lineair — zou zijn brief aan Pascal hebben uitgelokt, en daarmee de geboorte van de kansrekening. De juiste vergelijking gaat via logaritmen: nn pogingen met kans pp slagen minstens eenmaal met kans 1(1p)n1enp1 - (1-p)^n \approx 1 - \eu^{-np}, dus is de eerlijke invariant npnp: hier 416=234\cdot\frac16 = \frac23 tegenover 24136=2324\cdot\frac1{36} = \frac23 — gelijk! De twee weddenschappen verschillen pas op de tweede orde in pp, en net genoeg om er één over de vijftig-procent lijn te duwen: kleine kansen vormen een gebied waar de intuïtie de exponentiële functie nodig heeft, niet de liniaal.

Opmerking 21.18 (Klassieke drogredenen bij het voorwaardelijk stellen)

Drie terugkerende verwarringen, alle zichtbaar in de voorbeelden hierboven. (i) Omkering: PB ⁣(A)\pcond BA en PA ⁣(B)\pcond AB verschillen met de factor P(A)/P(B)\P(A)/\P(B) — een test die 99%99\% nauwkeurig is op de zieken kan een positieve patiënt toch bijna zeker gezond laten wanneer de ziekte zeldzaam is (Oefening 21.3); Pziek ⁣(positief)\pcond{\text{ziek}}{ \text{positief}} noemen waar Ppositief ⁣(ziek)\pcond{\text{positief}}{ \text{ziek}} wordt bedoeld, is de drogreden van het basispercentage. (ii) Voorwaardelijk stellen op de verkeerde gebeurtenis: bij Monty Hall is de juiste voorwaardelijke gebeurtenis “de presentator opende deur 33”, niet “de prijs zit niet achter deur 33”; die twee dragen verschillende informatie, en de hele 23\frac23 hangt van het verschil af. (iii) Disjunct tegenover onafhankelijk: disjuncte gebeurtenissen met positieve kans zijn nooit onafhankelijk (P(AB)=0P(A)P(B)\P(A\cap B) = 0 \neq \P(A)\P(B)) — onafhankelijkheid is verenigbaarheid van informatie, niet afwezigheid van overlap.

Definitie 21.19 (Onafhankelijkheid)

Gebeurtenissen AA en BB heten onafhankelijk wanneer P(AB)=P(A)P(B)\P(A \cap B) = \P(A)\P(B). Een familie (Ai)iI(A_i)_{i \in I} gebeurtenissen heet (onderling) onafhankelijk wanneer voor elke eindige deelverzameling JIJ \subseteq I geldt

P(iJAi)=iJP(Ai).\P\Bigl(\bigcap_{i \in J} A_i\Bigr) = \prod_{i \in J} \P(A_i) .

Opmerking 21.20

Onderlinge onafhankelijkheid is strikt sterker dan paarsgewijze onafhankelijkheid: bij twee eerlijke muntworpen zijn de gebeurtenissen “de eerste is kop”, “de tweede is kop” en “beide stemmen overeen” paarsgewijs onafhankelijk (elk paar heeft doorsnedekans 14=1212\frac14 = \frac12\cdot\frac12), en toch heeft de drievoudige doorsnede kans 1418\frac14 \neq \frac18. Merk ook op dat als A,BA, B onafhankelijk zijn, ook A,BcA, B^c dat zijn (reken na: P(ABc)=P(A)P(AB)=P(A)(1P(B))\P(A \cap B^c) = \P(A) - \P(A\cap B) = \P(A)(1 - \P(B))), en dus ook Ac,BcA^c, B^c.

Voorbeeld 21.21 (Onafhankelijkheid afgelezen van een productstructuur)

Gooi met twee eerlijke dobbelstenen: Ω=[ ⁣[1,6] ⁣]2\Omega = \intint16^2 met uniforme gewichten. Zij A=A = {}“eerste dobbelsteen even” en B=B = {}“tweede dobbelsteen minstens 55”. Tellen: A=36=18\abs A = 3\cdot6 = 18, B=62=12\abs B = 6\cdot2 = 12, AB=32=6\abs{A\cap B} = 3\cdot2 = 6, dus

P(AB)=636=18361236=P(A)P(B):\P(A\cap B) = \frac6{36} = \frac{18}{36}\cdot\frac{12}{36} = \P(A)\,\P(B) :

onafhankelijk, en het mechanisme is zichtbaar — AA legt alleen de eerste coördinaat vast, BB alleen de tweede, en de uniforme maat op een productverzameling laat de coördinaattellingen vermenigvuldigen. Elke bewering van het type “gebeurtenissen die van disjuncte groepen worpen afhangen, zijn onafhankelijk” (massaal gebruikt in de weekendopgave) is deze berekening, met meer indices.

Voorbeeld 21.22 (Analyse van de eerste stap)

Wat is voor het meetkundige model van Voorbeeld 21.5 de kans uu dat de eerste kop op een even rang valt? Stel voorwaardelijk op de eerste worp: met kans pp is de rang 11 (oneven); met kans q=1pq = 1 - p begint het spel opnieuw met alle pariteiten omgekeerd, dus

u=p0+q(1u)u=q1+q.u = p\cdot0 + q\,(1 - u) \qquad\Longrightarrow\qquad u = \frac{q}{1 + q} .

Eén regel, geen reeks — en het stemt overeen met de rechtstreekse sommatie van Oefening 21.9, die 1u=11+q1 - u = \frac1{1+q} geeft. Deze techniek van de “eerste stap” (voorwaardelijk stellen op het eerste experiment, en een verschoven kopie van het probleem herkennen) is de kansrekenkundige vorm van een recursie, en zij is de motor achter de vergelijkingen voor de speelduur van Oefening 21.6 en de berekeningen van de eerste doorgang in de weekendopgave.

21.3 Het lemma van Borel–Cantelli

Definitie 21.23 (Limes superior van gebeurtenissen)

Voor een rij (An)(A_n) gebeurtenissen is de gebeurtenis

lim supnAn=N=0 nNAn={ωΩ:ωAn voor oneindig veel n}\limsup_n A_n = \bigcap_{N=0}^{\infty}\ \bigcup_{n \geq N} A_n = \{\omega \in \Omega : \omega \in A_n \text{ voor oneindig veel } n\}

de gebeurtenisAnA_n treedt oneindig vaak op”.

Voorbeeld 21.24 (“Oneindig vaak” en “uiteindelijk” vertaald)

Het complement van lim supnAn\limsup_nA_n is volgens de Morgan

(NnNAn) ⁣c=NnNAnc={ω:ωAn voor alle grote n},\Bigl(\bigcap_N\bigcup_{n\geq N}A_n\Bigr)^{\!c} = \bigcup_N\bigcap_{n\geq N}A_n^c = \{\omega : \omega \notin A_n \text{ voor alle grote }n\},

de gebeurtenisuiteindelijk faalt AnA_n” (geschreven lim infnAnc\liminf_nA_n^c). Dus zijn “AnA_n oneindig vaak” en “AncA_n^c uiteindelijk” complementair — dit woordenboek recht houden voorkomt de meeste ongelukken met kwantoren. Voorbeeldvertalingen voor het muntwerpen: “oneindig veel kop” is lim sup{Xn=H}\limsup\{X_n = H\}; “slechts eindig veel reeksen van 100100 maal kop” is het complement van een limsup; “de lopende frequentie convergeert naar 12\frac12” is jNnN{p^n12<1j}\bigcap_j\bigcup_N\bigcap_{n\geq N}\{\abs{\widehat p_n - \tfrac12} < \tfrac1j\} — overal aftelbare bewerkingen, zodat dit alle eerlijke gebeurtenissen zijn.

Stelling 21.25 (Borel–Cantelli)

  1. Is nP(An)<\sum_{n} \P(A_n) < \infty, dan is P(lim supnAn)=0\P\bigl(\limsup_n A_n\bigr) = 0.
  2. Zijn de gebeurtenissen AnA_n onafhankelijk en nP(An)=\sum_n \P(A_n) = \infty, dan is P(lim supnAn)=1\P\bigl(\limsup_n A_n\bigr) = 1.

Bewijs. 1. Zij CN=nNAnC_N = \bigcup_{n \geq N}A_n; de rij (CN)(C_N) is dalend met doorsnede lim supAn\limsup A_n, en volgens de aftelbare subadditiviteit (Gevolg 21.7) is

P(CN)nNP(An)N0\P(C_N) \leq \sum_{n \geq N}\P(A_n) \xrightarrow[N\to\infty]{} 0

(staart van een convergente reeks). De monotone continuïteit (Stelling 21.6) besluit: P(lim supAn)=limNP(CN)=0\P(\limsup A_n) = \lim_N \P(C_N) = 0.

2. Het volstaat aan te tonen dat P(nNAn)=1\P\bigl(\bigcup_{n\geq N}A_n\bigr) = 1 voor elke NN: hebben immers gebeurtenissen BNB_N alle kans 11, dan is

P((NBN) ⁣c)=P(NBNc)NP(BNc)=0\P\Bigl(\Bigl(\bigcap_NB_N\Bigr)^{\!c}\Bigr) = \P\Bigl(\bigcup_NB_N^c\Bigr) \leq \sum_N\P(B_N^c) = 0

wegens de aftelbare subadditiviteit (Gevolg 21.7), zodat de aftelbare doorsnede lim supAn=NnNAn\limsup A_n = \bigcap_N\bigcup_{n\geq N}A_n nog altijd kans 11 heeft. Houd NN vast en beschouw voor M>NM > N het complement:

P(n=NMAnc)=n=NM(1P(An))n=NMeP(An)=exp(n=NMP(An)),\P\Bigl(\bigcap_{n=N}^{M} A_n^c\Bigr) = \prod_{n=N}^{M}\bigl(1 - \P(A_n)\bigr) \leq \prod_{n=N}^{M} e^{-\P(A_n)} = \exp\Bigl(-\sum_{n=N}^M \P(A_n)\Bigr) ,

met de onafhankelijkheid van de complementen en de convexiteitsgrens 1xex1 - x \leq e^{-x}. Als MM \to \infty streeft de exponent naar -\infty wegens de divergentie van de reeks, dus geeft de monotone continuïteit (dalende rij) dat P(nNAnc)=0\P\bigl(\bigcap_{n \geq N}A_n^c\bigr) = 0, dat wil zeggen P(nNAn)=1\P\bigl(\bigcup_{n \geq N}A_n\bigr) = 1.

Voorbeeld 21.26 (Oneindige reeksen kop)

Werp eeuwig met een eerlijke munt, en zij AnA_n de gebeurtenis “de worpen n,n+1,,n+k1n, n+1, \dots, n + k - 1 zijn alle kop” (een reeks van kk maal kop die op tijdstip nn begint), voor vaste kk. De gebeurtenissen AjkA_{jk} (j=1,2,j = 1, 2, \dots), die van disjuncte blokken worpen afhangen, zijn onafhankelijk, elk met kans 2k2^{-k}, en j2k=\sum_j 2^{-k} = \infty: volgens Borel–Cantelli 2 zijn met kans 11 oneindig veel blokken volledig kop — elk vast patroon keert bijna zeker oneindig vaak terug. Laten wij omgekeerd de lengte van de reeks groeien, dan heeft Bn=B_n = {}“een reeks van 2log2n2\log_2 n maal kop begint in nn” de kans P(Bn)=n2\P(B_n) = n^{-2}, sommeerbaar, dus beginnen er bijna zeker slechts eindig veel van zulke lange reeksen: Borel–Cantelli ijkt precies hoe lang de langste reeksen zijn.

Voorbeeld 21.27 (De oneindige aap, gekwantificeerd)

Een aap typt onafhankelijke uniforme letters uit een alfabet van 2626 letters. Snijd het typoscript in disjuncte blokken van vier letters; de gebeurtenissen Aj=A_j = {}“blok jj spelt MATH” zijn onafhankelijk met P(Aj)=264\P(A_j) = 26^{-4}, en jP(Aj)=\sum_j\P(A_j) = \infty: volgens Borel–Cantelli 2 typt de aap bijna zeker oneindig vaak MATH — en hetzelfde geldt voor elke vaste tekst van elke lengte, met aangepaste blokken. De kwantitatieve voetnoot ontnuchtert het wonder: 264=45697626^4 = 456\,976, dus vergt de eerste MATH gemiddeld ongeveer een half miljoen aanslagen, en een toneelstuk van Shakespeare met 10510^5 tekens wacht op de orde van 2610526^{10^5} blokken — bijna zeker is een uitspraak over de horizon \infty, niet over enige horizon die een aap zal ontmoeten. Borel–Cantelli bevestigt de limiet; de grootte van de termen vertelt het verhaal op menselijke schaal.

Opmerking 21.28

In Voorbeeld 21.26 is de onderliggende uitkomstenruimte (oneindige rijen worpen) overaftelbaar, dus leeft het voorbeeld strikt genomen in het maattheoretische kader van bachelorjaar 3; de berekeningen gebruiken echter alleen de regels die in dit hoofdstuk zijn bewezen, toegepast op gebeurtenissen die door eindig veel worpen worden bepaald en op hun aftelbare combinaties. Dat is de gebruikelijke afspraak op dit niveau: de theorie wordt op aftelbare ruimten geformuleerd, en voorbeelden met oneindige spelen worden met hetzelfde gereedschap behandeld.

Opmerking 21.29 (Vooruitblik binnen dit volume)

De machinerie van dit hoofdstuk wordt door de volgende twee in haar geheel verbruikt. Indicatoren maken van gebeurtenissen toevalsveranderlijken, en de σ\sigma-additiviteit wordt de sommeerbaarheid die de verwachtingswaarde definieert (Hoofdstuk 22); Borel–Cantelli plus een sommeerbare staartgrens is precies hoe de sterke wet van de grote aantallen voor munten daar wordt bewezen. In Hoofdstuk 23 duikt de monotone continuïteit op het beslissende ogenblik weer op: de uitstervingskans van een vertakkingsproces wordt gedefinieerd als de monotone limiet limP(Zn=0)\lim\P(Z_n = 0), en de vastepuntsvergelijking waaraan zij voldoet, wordt verkregen door in die stijgende rij naar de limiet over te gaan — de laatste stelling van het boek staat op de eerste stelling van dit hoofdstuk.

Opmerking 21.30 (Methode: drie wegen naar kans één)

Uitspraken “bijna zeker” worden met drie hefbomen bewezen, in stijgende volgorde van kracht. Monotone continuïteit: geef de gebeurtenis als een stijgende vereniging (of dalende doorsnede) van gebeurtenissen met eindige horizon en berekenbare kansen (Voorbeeld 21.9). Nulverenigingen: een aftelbare vereniging van gebeurtenissen met kans nul is nul (aftelbare subadditiviteit), dus volstaat het elke slechte gebeurtenis afzonderlijk te doden — zo voegt “voor elke jj is uiteindelijk p^np<1/j\abs{\widehat p_n - p} < 1/j” zich samen tot convergentie. Borel–Cantelli: is de gebeurtenis een limsup, sommeer dan de kansen; convergentie doodt haar (zonder onafhankelijkheid nodig te hebben), en divergentie plus onafhankelijkheid bevestigt haar. De juiste hefboom kiezen is meestal het hele bewijs; de weekendopgave laat alle drie in één argument lopen.

Opmerking 21.31 (Waar dit wordt gebruikt)

De monotone continuïteit en Borel–Cantelli zijn de twee hefbomen van elke uitspraak “bijna zeker”: zij drijven de recurrentie van de toevalswandeling in de weekendopgave van dit hoofdstuk aan, de bijna-zekere kant van de wet van de grote aantallen (Hoofdstuk 22) en de analyse van het uitsterven van vertakkingsprocessen (Hoofdstuk 23). Het volume van bachelorjaar 3 bouwt de theorie opnieuw op met σ\sigma-algebra’s en de integraal van Lebesgue, waar de overaftelbare uitkomstenruimten die hier informeel worden gebruikt, volledig streng worden.

21.4 Oefeningen

Oefening 21.1

Een urne bevat nn genummerde ballen. De ballen worden één voor één zonder teruglegging getrokken. Bereken de kans dat bal nummer 11 vóór bal nummer 22 wordt getrokken. Veralgemeen: de kans dat bal 11 als eerste wordt getrokken onder de ballen 1,,k1, \dots, k.

Oplossing

Oplossing van Oefening 21.1.

Wegens de symmetrie: de trekvolgorde brengt een uniform willekeurige onderlinge volgorde op de ballen 11 en 22 voort, dus P(1 voˊoˊ2)=12\P(1 \text{ vóór } 2) = \frac12. Formeel: de posities van de ballen 11 en 22 in een trekrij verwisselen is een bijectie van de (even waarschijnlijke) uitkomsten die de gebeurtenis met haar complement verwisselt. Onder de ballen 1,,k1, \dots, k: de onderlinge volgorde van deze kk ballen is uniform over de k!k! rangschikkingen, en bal 11 staat in (k1)!(k-1)! daarvan vooraan: kans (k1)!k!=1k\frac{(k-1)!}{k!} = \frac1k.

Oefening 21.2

Toon aan dat de gewichten pk=1k(k+1)p_k = \frac{1}{k(k+1)} op Ω=N\Omega = \N^* een kansmaat definiëren, en bereken P(2N)\P(2\N^*) (even uitkomsten) als reeks; toon aan dat zij gelijk is aan 1ln21 - \ln 2. (Telescopeer 12j(2j+1)=12j12j+1\frac{1}{2j(2j+1)} = \frac{1}{2j} - \frac{1}{2j+1} en gebruik de alternerende harmonische reeks, Hoofdstuk 7.)

Oplossing

Oplossing van Oefening 21.2.

1k(k+1)=1k1k+1\frac{1}{k(k+1)} = \frac1k - \frac1{k+1}, dus telescopeert k1pk\sum_{k\geq1} p_k tot 11: een kansmaat. Even uitkomsten:

P(2N)=j=112j(2j+1)=j=1(12j12j+1)=1213+1415+\P(2\N^*) = \sum_{j=1}^{\infty}\frac{1}{2j(2j+1)} = \sum_{j=1}^{\infty}\Bigl(\frac{1}{2j} - \frac{1}{2j+1}\Bigr) = \frac12 - \frac13 + \frac14 - \frac15 + \cdots

Dit is de alternerende harmonische reeks zonder haar eerste term en met omgekeerde tekens: omdat ln2=112+1314+\ln 2 = 1 - \frac12 + \frac13 - \frac14 + \cdots (Hoofdstuk 7), is

P(2N)=(ln21)=1ln20.307.\P(2\N^*) = -\bigl(\ln 2 - 1\bigr) = 1 - \ln 2 \approx 0.307 .

Oefening 21.3

(Vals-positieven) Een ziekte treft één persoon op 1000010\,000. Een test spoort haar met kans 0.990.99 op bij de zieken, en geeft met kans 0.010.01 een vals-positief bij de gezonden. Bereken de kans ziek te zijn gegeven een positieve test, en geef commentaar.

Oplossing

Oplossing van Oefening 21.3.

Zij SS = ziek en ++ = positieve test. Bayes (Stelling 21.14) met de partitie {S,Sc}\{S, S^c\}:

P(S+)=0.99×1040.99×104+0.01×0.9999=0.0000990.000099+0.0099990.0098,\P(S \mid +) = \frac{0.99 \times 10^{-4}} {0.99 \times 10^{-4} + 0.01 \times 0.9999} = \frac{0.000099}{0.000099 + 0.009999} \approx 0.0098 ,

onder 1%1\%. Hoewel de test “99%99\% nauwkeurig” is, laat een positieve uitslag je met ongeveer 99%99\% kans gezond: de vals-positieven onder de enorme gezonde meerderheid overspoelen de echte positieven uit de piepkleine zieke minderheid. Screeningtests voor zeldzame aandoeningen moeten altijd door deze berekening met het basispercentage worden gelezen.

Oefening 21.4 ★★

Zijn A1,,AnA_1, \dots, A_n gebeurtenissen. Bewijs de formule van de in- en uitsluiting

P(i=1nAi)=J{1,,n}(1)J+1P(iJAi)\P\Bigl(\bigcup_{i=1}^n A_i\Bigr) = \sum_{\emptyset \neq J \subseteq \{1,\dots,n\}} (-1)^{\abs J + 1}\,\P\Bigl(\bigcap_{i \in J}A_i\Bigr)

door de identiteit 1i=1n(11Ai)=1Ai1 - \prod_{i=1}^n(1 - \mathbf{1}_{A_i}) = \mathbf{1}_{\bigcup A_i} over Ω\Omega te integreren (dat wil zeggen: te sommeren met gewichten P({ω})\P(\{\omega\})).

Oplossing

Oplossing van Oefening 21.4.

Puntsgewijs op Ω\Omega: ωAi\omega \in \bigcup A_i dan en slechts dan als een zekere factor 11Ai(ω)1 - \mathbf{1}_{A_i}(\omega) verdwijnt, dus is

1Ai=1i=1n(11Ai)=J{1,,n}(1)J+1iJ1Ai,\mathbf{1}_{\bigcup A_i} = 1 - \prod_{i=1}^n\bigl(1 - \mathbf{1}_{A_i}\bigr) = \sum_{\emptyset \neq J \subseteq \{1,\dots,n\}} (-1)^{\abs J + 1}\prod_{i \in J}\mathbf{1}_{A_i} ,

door het product uit te werken en de 11 over te brengen. Nu is iJ1Ai=1iJAi\prod_{i\in J}\mathbf{1}_{A_i} = \mathbf{1}_{\bigcap_{i \in J} A_i}, en sommeren tegen de gewichten P({ω})\P(\{\omega\}) — geoorloofd: eindig veel begrensde termen, elke familie sommeerbaar — maakt van elke indicator de kans op haar gebeurtenis, wat de formule geeft.

Oefening 21.5 ★★

(Het probleem van de overeenkomsten, via in- en uitsluiting) nn brieven worden uniform willekeurig in nn enveloppen gestopt, één per envelop. Toon met Oefening 21.4 aan dat de kans op geen enkele juiste overeenkomst gelijk is aan k=0n(1)kk!e1\sum_{k=0}^n \frac{(-1)^k}{k!} \to e^{-1}, en leid de kans op precies één overeenkomst af.

Oplossing

Oplossing van Oefening 21.5.

Zij AiA_i = “brief ii zit in de juiste envelop”. Voor JJ van grootte kk is P(iJAi)=(nk)!n!\P\bigl(\bigcap_{i\in J}A_i\bigr) = \frac{(n-k)!}{n!} (leg kk brieven vast, permuteer de rest). Volgens de in- en uitsluiting is

P(Ai)=k=1n(1)k+1(nk)(nk)!n!=k=1n(1)k+1k!,\P\Bigl(\bigcup A_i\Bigr) = \sum_{k=1}^n (-1)^{k+1}\binom nk \frac{(n-k)!}{n!} = \sum_{k=1}^n \frac{(-1)^{k+1}}{k!} ,

dus

P(geen overeenkomst)=1P(Ai)=k=0n(1)kk!ne10.368.\P(\text{geen overeenkomst}) = 1 - \P\Bigl(\bigcup A_i\Bigr) = \sum_{k=0}^{n}\frac{(-1)^k}{k!} \xrightarrow[n\to\infty]{} e^{-1} \approx 0.368 .

Precies één overeenkomst: een permutatie met precies één vast punt wordt bepaald door de keuze van de vaste brief (nn manieren) en een derangement (rangschikking zonder overeenkomst) van de andere n1n - 1; schrijven wij Dn1=(n1)!k=0n1(1)kk!D_{n-1} = (n-1)!\sum_{k=0}^{n-1}\frac{(-1)^k}{k!} voor het aantal derangementen (het eerste deel, met (n1)!(n-1)! geschaald), dan is

P(precies eˊeˊn overeenkomst)=nDn1n!=Dn1(n1)!=k=0n1(1)kk!ne1:\P(\text{precies één overeenkomst}) = \frac{n\,D_{n-1}}{n!} = \frac{D_{n-1}}{(n-1)!} = \sum_{k=0}^{n-1}\frac{(-1)^k}{k!} \xrightarrow[n\to\infty]{} e^{-1} :

in de limiet zijn “geen overeenkomst” en “precies één overeenkomst” even waarschijnlijk, elk met kans e1e^{-1}.

Oefening 21.6 ★★

Met een scheve munt (kanskop p(0,1)p \in \intoo{0}{1}) wordt geworpen tot er tweemaal achter elkaar kop verschijnt. Zij qnq_n de kans dat het spel langer dan nn worpen duurt. Toon, door voorwaardelijk te stellen op de eerste worp(en), aan dat qn=(1p)qn1+p(1p)qn2q_n = (1-p)\,q_{n-1} + p(1-p)\,q_{n-2} voor n2n \geq 2, en leid af dat het spel met kans 11 eindigt. (Toon qn0q_n \to 0 aan door met een meetkundige rij te vergelijken: beide nulpunten van de karakteristieke vergelijking hebben absolute waarde in (0,1)\intoo{0}{1}.)

Oplossing

Oplossing van Oefening 21.6.

Stel voorwaardelijk op het begin (kettingregel / Stelling 21.14):

  • eerste worp M (kans 1p1 - p): het spel begint opnieuw; langer duren dan nn betekent van daar af langer duren dan n1n - 1: bijdrage (1p)qn1(1-p)\,q_{n-1};
  • eerste worpen KM (kans p(1p)p(1-p)): opnieuw beginnen na twee worpen: bijdrage p(1p)qn2p(1-p)\,q_{n-2};
  • eerste worpen KK: het spel is geëindigd (binnen nn worpen, n2n \geq 2): bijdrage 00.

Bijgevolg is qn=(1p)qn1+p(1p)qn2q_n = (1-p)q_{n-1} + p(1-p)q_{n-2}. De karakteristieke vergelijking r2=(1p)r+p(1p)r^2 = (1-p)r + p(1-p) heeft nulpunten

r±=(1p)±(1p)2+4p(1p)2,r_\pm = \frac{(1-p) \pm \sqrt{(1-p)^2 + 4p(1-p)}}{2},

met r±<1\abs{r_\pm} < 1: inderdaad voldoet de veelterm χ(r)=r2(1p)rp(1p)\chi(r) = r^2 - (1-p)r - p(1-p) aan χ(1)=1(1p)p(1p)=p2>0\chi(1) = 1 - (1-p) - p(1-p) = p^2 > 0 en χ(1)=1+(1p)p(1p)>0\chi(-1) = 1 + (1-p) - p(1-p) > 0, terwijl χ(0)=p(1p)<0\chi(0) = -p(1-p) < 0: één nulpunt in (1,0)\intoo{-1}{0}, één in (0,1)\intoo{0}{1}. Dus qn=αr+n+βrn0q_n = \alpha r_+^n + \beta r_-^n \to 0. De gebeurtenissen “het spel duurt langer dan nn” dalen naar “het spel eindigt nooit”; de monotone continuïteit (Stelling 21.6) geeft P(eindigt nooit)=limqn=0\P(\text{eindigt nooit}) = \lim q_n = 0: het spel eindigt bijna zeker.

Oefening 21.7 ★★★

(Records) Trek een oneindige rij onafhankelijke uniforme rangschikkingen, in de volgende combinatorische zin: voor elke nn is de onderlinge volgorde van de eerste nn trekkingen uniform over de n!n! mogelijkheden, en Rn=R_n = {}“de nn-de trekking is een record (groter dan alle vorige)”. Aangenomen dat de gebeurtenissen RnR_n onafhankelijk zijn met P(Rn)=1/n\P(R_n) = 1/n (bewijs minstens deze laatste gelijkheid met de symmetrie), toon met Borel–Cantelli aan dat er bijna zeker oneindig veel records optreden, maar dat records op opeenvolgende tijdstippen n,n+1n, n+1 oneindig vaak optreden met kans — bereken nP(RnRn+1)\sum_n \P(R_n \cap R_{n+1}) en besluit wat Borel–Cantelli 1 geeft.

Oplossing

Oplossing van Oefening 21.7.

P(Rn)=1/n\P(R_n) = 1/n: onder de eerste nn trekkingen is elk van de nn onderlinge posities van de laatste trekking even waarschijnlijk (uniformiteit van de onderlinge volgorde), en RnR_n is de gebeurtenis dat zij de grootste is: kans 1/n1/n.

Oneindig veel records: nP(Rn)=1/n=\sum_n \P(R_n) = \sum 1/n = \infty en de RnR_n zijn onafhankelijk (aangenomen), dus geeft Borel–Cantelli 2 (Stelling 21.25) dat P(lim supRn)=1\P(\limsup R_n) = 1: records houden bijna zeker nooit op — maar zij dunnen logaritmisch uit.

Opeenvolgende records: wegens de onafhankelijkheid is

nP(RnRn+1)=n1n(n+1)<,\sum_n \P(R_n \cap R_{n+1}) = \sum_n \frac{1}{n(n+1)} < \infty ,

dus is Borel–Cantelli 1 van toepassing: bijna zeker wordt slechts eindig vaak een record onmiddellijk door een ander record gevolgd. De twee helften van het lemma werken samen: oneindig veel records, maar (bijna zeker) uiteindelijk nooit twee op rij.

Oefening 21.8 ★★★

(In de stijl van Kochen–Stone, eenvoudiger versie) Zij (An)(A_n) een rij onafhankelijke gebeurtenissen met P(An)=1n+1\P(A_n) = \frac{1}{n+1}. Toon aan dat P(lim supAn)=1\P(\limsup A_n) = 1, hoewel P(An)0\P(A_n) \to 0: “afzonderlijk zeldzaam, samen zeker”. Geef omgekeerd een rij (afhankelijke) gebeurtenissen met P(An)=\sum\P(A_n) = \infty en P(lim supAn)=0\P(\limsup A_n) = 0, wat aantoont dat de onafhankelijkheid in Borel–Cantelli 2 niet mag vervallen.

Oplossing

Oplossing van Oefening 21.8.

Eerste deel: P(An)=1n+1=\sum \P(A_n) = \sum\frac{1}{n+1} = \infty met onafhankelijkheid: Borel–Cantelli 2 geeft P(lim supAn)=1\P(\limsup A_n) = 1. Elke afzonderlijke AnA_n is steeds onwaarschijnlijker, en toch behoort bijna elke ω\omega tot oneindig veel ervan.

Tegenvoorbeeld zonder onafhankelijkheid: neem Ω=N\Omega = \N^* met de gewichten pk=1k(k+1)p_k = \frac{1}{k(k+1)} van Oefening 21.2, en An={kN:kn}A_n = \{k \in \N^* : k \geq n\}. Dan is

P(An)=kn(1k1k+1)=1n,nP(An)=,\P(A_n) = \sum_{k \geq n}\Bigl(\frac1k - \frac1{k+1}\Bigr) = \frac1n , \qquad \sum_n \P(A_n) = \infty ,

maar de AnA_n zijn dalend, dus lim supnAn=nAn=\limsup_n A_n = \bigcap_n A_n = \emptyset: P(lim supAn)=0\P(\limsup A_n) = 0. De divergentie van P(An)\sum\P(A_n) alleen waarborgt niets wanneer de gebeurtenissen zich op een krimpend deel van de ruimte opstapelen — de onafhankelijkheid is wat die samenzwering verbiedt.

Oefening 21.9

Met een munt met kanskop p(0,1)p \in \intoo01 wordt geworpen tot de eerste kop. Bereken de kans dat dit op een oneven rang gebeurt, en evalueer haar voor een eerlijke munt.

Oplossing

Oplossing van Oefening 21.9.

Met q=1pq = 1 - p valt de eerste kop op rang 2j+12j + 1 met kans q2jpq^{2j}p, dus

P(oneven rang)=j0q2jp=p1q2=11+q.\P(\text{oneven rang}) = \sum_{j\geq0}q^{2j}p = \frac{p}{1 - q^2} = \frac{1}{1 + q} .

Voor een eerlijke munt: 11+1/2=23\frac1{1 + 1/2} = \frac23. (Verstandscontrole: oneven rangen moeten waarschijnlijker zijn, want rang 11 komt als eerste — en inderdaad is 11+q>12\frac1{1+q} > \frac12 altijd.)

Oefening 21.10 ★★

Zij (An)n1(A_n)_{n\geq1} een rij onafhankelijke gebeurtenissen met P(An)=pn<1\P(A_n) = p_n < 1. Toon aan dat

P(n1Anc)=n1(1pn):=limNn=1N(1pn),\P\Bigl(\bigcap_{n\geq1}A_n^c\Bigr) = \prod_{n\geq1}(1 - p_n) := \lim_{N\to\infty}\prod_{n=1}^N(1 - p_n),

en dat deze limiet >0> 0 is dan en slechts dan als pn<\sum p_n < \infty. Breng dit in overeenstemming met Borel–Cantelli: is pn=\sum p_n = \infty, dan treedt niet alleen een zekere AnA_n bijna zeker op — er treden er oneindig veel op.

Oplossing

Oplossing van Oefening 21.10.

De gebeurtenissen BN=n=1NAncB_N = \bigcap_{n=1}^N A_n^c dalen naar nAnc\bigcap_nA_n^c, en wegens de onafhankelijkheid van de complementen is P(BN)=n=1N(1pn)\P(B_N) = \prod_{n=1}^N(1 - p_n); de monotone continuïteit (Stelling 21.6) geeft de getoonde limiet. Logaritmen nemen: (1pn)>0\prod(1 - p_n) > 0 dan en slechts dan als ln(1pn)<\sum-\ln(1 - p_n) < \infty. Is pn<\sum p_n < \infty, dan is pn0p_n \to 0 en ln(1pn)pn-\ln(1 - p_n) \sim p_n: de logaritmereeks convergeert. Is pn=\sum p_n = \infty, dan dwingt ln(1pn)pn-\ln(1 - p_n) \geq p_n divergentie af, dus is het product 00. Dit stemt overeen met Borel–Cantelli 2: voor pn=\sum p_n = \infty is niet alleen P(geen enkele An treedt op)=0\P(\text{geen enkele }A_n\text{ treedt op}) = 0, maar treden er bijna zeker oneindig veel AnA_n op.

Oefening 21.11 ★★

(De luciferdoosjes van Banach) Een roker houdt in elke zak een doosje met nn lucifers en grijpt telkens in een uniform willekeurige zak. Wanneer hij voor het eerst een doosje leeg aantreft, wat is dan de kans dat het andere doosje precies kk lucifers bevat? Toon aan dat het antwoord (2nkn)2(2nk)\binom{2n-k}{n}2^{-(2n-k)} is en ga na dat deze kansen voor n=1n = 1 tot 11 sommeren.

Oplossing

Oplossing van Oefening 21.11.

Zeg dat doosje AA het doosje is dat het eerst leeg wordt aangetroffen, terwijl het andere er kk bevat. Dat betekent: onder de eerste 2nk2n - k grepen gingen er precies nn naar AA en nkn - k naar BB (in een zekere volgorde), en greep nummer 2nk+12n - k + 1 ging opnieuw naar AA, dat leeg bleek. De grepen zijn onafhankelijke eerlijke keuzes, dus heeft deze gebeurtenis kans (2nkn)2(2nk)12\binom{2n-k}{n}2^{-(2n-k)}\cdot\frac12; verdubbelen (het lege doosje kan er een van beide zijn) geeft

P(het andere doosje bevat k)=(2nkn)2(2nk).\P(\text{het andere doosje bevat }k) = \binom{2n-k}{n}\,2^{-(2n-k)} .

Voor n=1n = 1: k=1k = 1 geeft (11)21=12\binom11 2^{-1} = \frac12 en k=0k = 0 geeft (21)22=12\binom21 2^{-2} = \frac12: samen 11, zoals het moet.

Oefening 21.12 ★★★

(σ\sigma-additiviteit is een echt axioma) (a) Toon aan dat er geen kansmaat op (N,P(N))(\N, \mathcal P(\N)) bestaat die aan alle singletons hetzelfde gewicht geeft. (b) Zij voor ANA \subseteq \N^* de grootheid d(A)=limnA[ ⁣[1,n] ⁣]nd(A) = \lim_n\frac{\abs{A\cap\intint1n}}{n} wanneer de limiet bestaat (de natuurlijke dichtheid). Toon aan dat dd eindig additief is op paren waarvoor alle drie de dichtheden bestaan, dat zij aan elk singleton dichtheid 00 en aan N\N^* dichtheid 11 geeft — en besluit dat dd niet σ\sigma-additief is. (c) Geef een verzameling zonder dichtheid. (Neem de blokken [ ⁣[22k,22k+11] ⁣]\intint{2^{2k}}{2^{2k+1}-1} afwisselend wel en niet mee.)

Oplossing

Oplossing van Oefening 21.12.

(a) Is P({n})=c\P(\{n\}) = c voor alle nn, dan dwingt de σ\sigma-additiviteit 1=nc1 = \sum_nc af: onmogelijk, of c=0c = 0 is (som 00) of c>0c > 0 (som oneindig). Er bestaat geen uniforme kans op N\N.

(b) Is AB=A \cap B = \emptyset en bestaan d(A)d(A) en d(B)d(B), dan is (AB)[ ⁣[1,n] ⁣]=A[ ⁣[1,n] ⁣]+B[ ⁣[1,n] ⁣]\abs{(A \sqcup B)\cap\intint1n} = \abs{A\cap\intint1n} + \abs{B\cap\intint1n}, dus d(AB)=d(A)+d(B)d(A \sqcup B) = d(A) + d(B): eindige additiviteit op zulke paren. Elk singleton heeft een telfunctie die uiteindelijk constant is, en dus dichtheid 00, terwijl d(N)=1d(\N^*) = 1. Was dd σ\sigma-additief, dan zou N=k{k}\N^* = \bigsqcup_k\{k\} geven dat 1=k0=01 = \sum_k 0 = 0: de dichtheid is eindig additief maar niet σ\sigma-additief — het axioma heeft inhoud.

(c) Zij A=k0[ ⁣[4k,24k1] ⁣]A = \bigcup_{k\geq0}\intint{4^k}{2\cdot4^k - 1} (blokken van 4k4^k tot 24k12\cdot4^k - 1). In n=24K1n = 2\cdot4^K - 1 is de telling kK4k434K\sum_{k\leq K}4^k \sim \frac43 4^K, wat een verhouding 23\to \frac23 geeft; in n=4K+11n = 4^{K+1} - 1 is de telling onveranderd, wat een verhouding 13\to \frac13 geeft. De verhouding oscilleert tussen de limieten 13\frac13 en 23\frac23: geen dichtheid.

21.5 Probleem: de eenvoudige toevalswandeling op Z\Z is recurrent

Vierentwintig stappen van een eenvoudige toevalswandeling; de rode stippen markeren de terugkeren naar de oorsprong. De opgave toont dat deze stippen met kans 1 nooit ophouden te verschijnen — en toch heeft de wachttijd ertussen een divergent gemiddelde.
Vierentwintig stappen van een eenvoudige toevalswandeling; de rode stippen markeren de terugkeren naar de oorsprong. De opgave toont dat deze stippen met kans 11 nooit ophouden te verschijnen — en toch heeft de wachttijd ertussen een divergent gemiddelde.

Probleem 21.1

Weekendopgave — de recurrentiestelling van Pólya op Z\Z, met onderweg het stemmenprobleem en een vleugje arcsinus

Werp eeuwig met een eerlijke munt; zij Xi=±1X_i = \pm1 de ii-de stap en Sn=X1++XnS_n = X_1 + \dots + X_n de eenvoudige toevalswandeling op Z\Z, met S0=0S_0 = 0. Zoals in Voorbeeld 21.26 worden alle gebeurtenissen hieronder door eindig veel worpen bepaald of zijn zij aftelbare combinaties van zulke gebeurtenissen, en de onafhankelijkheid van gebeurtenissen die van disjuncte blokken worpen afhangen, hoort bij het model. Wij schrijven un=P(S2n=0)u_n = \P(S_{2n} = 0) en Nn(k)N_n(k) voor het aantal ±1\pm1-paden van lengte nn van 00 naar kk.

Deel I — Paden tellen.

  1. Toon aan dat Nn(k)=(n(n+k)/2)N_n(k) = \binom{n}{(n+k)/2} wanneer n+kn + k even is en kn\abs k \leq n, en 00 anders; leid af dat P(Sn=k)=Nn(k)2n\P(S_n = k) = N_n(k)\,2^{-n}. Waarom is elk afzonderlijk pad van lengte nn even waarschijnlijk?
  2. Toon aan dat S2n+10S_{2n+1} \neq 0 en un=(2nn)4nu_n = \binom{2n}{n}4^{-n}, en bereken u1,u2,u3u_1, u_2, u_3.
  3. Bewijs un=2n12nun1u_n = \frac{2n-1}{2n}\,u_{n-1}; leid af dat (un)(u_n) daalt naar 00, en uit Voorbeeld 6.14 dat

    un1πn,zodatnun=.u_n \sim \frac{1}{\sqrt{\pi n}}, \qquad\text{zodat}\qquad \sum_n u_n = \infty .
  4. (Spiegelingsprincipe) Toon voor k1k \geq 1 aan dat de paden van lengte nn van 11 naar kk die 00 raken, in bijectie staan met de paden van 1-1 naar kk; leid af dat het aantal paden van 00 naar kk dat na tijdstip 00 boven 00 blijft, gelijk is aan Nn1(k1)Nn1(k+1)N_{n-1}(k-1) - N_{n-1}(k+1).
  5. (Stemmenstelling) Leid af dat

    P(S1>0,,Sn1>0Sn=k)=kn(k1):\P\bigl(S_1 > 0, \dots, S_{n-1} > 0 \bigm| S_n = k\bigr) = \frac kn \qquad (k \geq 1) :

    bij een telling waarin de winnaar met kk van de nn stemmen voorstaat, is de kans dat hij gedurende de hele telling voorstond gelijk aan k/nk/n. Ga dit met de hand na voor n=3n = 3, k=1k = 1.

Deel II — Terugkeer naar de oorsprong.

  1. Bewijs de sleutelidentiteit

    P(S10, S20, , S2n0)=un\P(S_1 \neq 0,\ S_2 \neq 0,\ \dots,\ S_{2n} \neq 0) = u_n

    (stel voorwaardelijk op de eerste stap, sommeer de tellingen van vraag 4 over het eindpunt, en telescopeer; eindig met 2(2n1n)=(2nn)2\binom{2n-1}{n} = \binom{2n}{n}).

  2. Leid uit de monotone continuïteit (Stelling 21.6) af dat de wandeling met kans 11 minstens eenmaal naar 00 terugkeert, en dat fn:=P(eerste terugkeer op tijdstip 2n)f_n := \P(\text{eerste terugkeer op tijdstip }2n) voldoet aan

    fn=un1un=un2n1,n1fn=1.f_n = u_{n-1} - u_n = \frac{u_n}{2n-1}, \qquad \sum_{n\geq1}f_n = 1 .
  3. Toon aan dat n2nfn=\sum_n 2n\,f_n = \infty: de terugkeer is zeker, maar de reeks die de gemiddelde wachttijd zou berekenen, divergeert (in de woordenschat van Hoofdstuk 22 heeft de terugkeertijd een oneindige verwachtingswaarde).
  4. Bewijs dat voor elke k1k \geq 1 geldt P(minstens k terugkeren naar 0)=1\P(\text{minstens } k\text{ terugkeren naar }0) = 1 (ontbind over de tijdstippen van de eerste kk terugkeren: de bijbehorende blokken worpen zijn disjunct, dus vermenigvuldigen de kansen en sommeren zij tot (nfn)k(\sum_nf_n)^k); besluit met de monotone continuïteit:

    P(Sn=0 voor oneindig veel n)=1:\P(S_n = 0 \text{ voor oneindig veel } n) = 1 :

    de eenvoudige toevalswandeling op Z\Z is recurrent.

  5. Toon aan dat de wandeling bijna zeker elke plaats kZk \in \Z bezoekt, en dus (wegens de recurrentie, herstart bij het eerste bezoek) oneindig vaak. (De tekens van de opeenvolgende uitstapjes vanuit 00 zijn onafhankelijke eerlijke munten; een positief uitstapje bezoekt 11.)

Deel III — Borel–Cantelli en de scheve wandeling.

  1. De gebeurtenissen An={S2n=0}A_n = \{S_{2n} = 0\} voldoen aan P(An)=\sum\P(A_n) = \infty; leg uit waarom Borel–Cantelli 2 er niet op van toepassing is, en wat Borel–Cantelli 1 zou geven als de reeks convergeerde. (Dat is de strategie van het hele deel.)
  2. Laat de munt nu een scheefheid p12p \neq \frac12 hebben, q=1pq = 1 - p. Toon aan dat P(S2n=0)=(2nn)(pq)n=un(4pq)n\P(S_{2n} = 0) = \binom{2n}n(pq)^n = u_n\,(4pq)^n met 4pq<14pq < 1, leid af dat nP(S2n=0)<\sum_n\P(S_{2n} = 0) < \infty, en besluit met Borel–Cantelli 1 dat de scheve wandeling bijna zeker slechts eindig vaak naar 00 terugkeert.
  3. Nog steeds voor p12p \neq \frac12: toon aan dat P(Sn=k)(nn/2)(pq)n/2(p/q)k/2\P(S_n = k) \leq \binom{n}{\floor{n/2}}\,(pq)^{n/2}\,(p/q)^{k/2} voor elke vaste kk, leid af dat elke plaats bijna zeker eindig vaak wordt bezocht, en besluit dat Sn\abs{S_n} \to \infty bijna zeker: de scheve wandeling is transiënt.
  4. Terug naar de eerlijke munt: bereken met vraag 6 de kans dat 200200 worpen geen enkele gelijkstand opleveren (Sn0S_n \neq 0 voor 1n2001 \leq n \leq 200), numeriek u1000.056u_{100} \approx 0.056. Geef commentaar op het trage verval 1/πn1/\sqrt{\pi n}: gelijkstanden zijn op de lange duur zeker, maar zeldzamer dan de intuïtie doet vermoeden.
  5. (Eerste doorgang) Zij T1T_1 het eerste tijdstip waarop de wandeling 11 raakt. Toon met het spiegelingsprincipe voor het maximum Mn=maxinSiM_n = \max_{i\leq n}S_i (bewezen in vraag 16, die niet van deze afhangt), of rechtstreeks uit vraag 7 door voorwaardelijk op de eerste stap te stellen, aan dat P(T1=2n1)=fn\P(T_1 = 2n - 1) = f_n; leid af dat P(T1<)=1\P(T_1 < \infty) = 1, terwijl de reeks (2n1)fn\sum(2n-1)f_n voor de gemiddelde tijd divergeert.

Deel IV — Maxima, laatste nulpunt, lange voorsprongen.

  1. (Spiegeling voor het maximum) Bewijs voor k1k \geq 1 dat

    P(Mnk)=2P(Sn>k)+P(Sn=k)\P(M_n \geq k) = 2\,\P(S_n > k) + \P(S_n = k)

    door het pad na zijn eerste bezoek aan niveau kk te spiegelen.

  2. Leid af dat P(M2n1)=1un\P(M_{2n} \geq 1) = 1 - u_n, dat wil zeggen P(Si0 voor alle i2n)=un\P(S_i \leq 0 \text{ voor alle } i \leq 2n) = u_n: de kans om nooit voor te staan is gelijk aan de kans om nooit op nul te zijn (vraag 6) — twee verschillende gebeurtenissen, één kans.
  3. (Laatste nulpunt) Zij L2n=max{k2n:Sk=0}L_{2n} = \max\{k \leq 2n : S_k = 0\} (even). Toon, door vraag 6 met de onafhankelijkheid van disjuncte blokken worpen te combineren, aan dat

    P(L2n=2k)=ukunk(0kn),\P(L_{2n} = 2k) = u_k\,u_{n-k} \qquad (0 \leq k \leq n),

    en leid zonder verdere berekening de binomiale identiteit k=0nukunk=1\sum_{k=0}^n u_ku_{n-k} = 1 af.

  4. Toon aan dat de verdeling van L2nL_{2n} symmetrisch is (P(L=2k)=P(L=2n2k)\P(L = 2k) = \P(L = 2n - 2k)) en, met uj1/πju_j \sim 1/\sqrt{\pi j}, dat haar uitersten haar meest waarschijnlijke waarden zijn. Zet dit voor n=5n = 5 in een tabel: P(L10=0)=u50.246\P(L_{10} = 0) = u_5 \approx 0.246 tegenover P(L10=4)=u2u30.117\P(L_{10} = 4) = u_2u_3 \approx 0.117. Interpreteer: in een lang eerlijk spel valt de laatste gelijkstand meestal heel vroeg of heel laat — lange voorsprongen zijn de regel, niet de uitzondering.
  5. Zet de vragen 16–19 in elkaar tot een alinea over het beeld van de fluctuaties van de eerlijke wandeling: de diffusieve schaal die vraag 3 suggereert, de zekerheid van de terugkeer tegenover de divergente gemiddelde wachttijd, en de volharding van voorsprongen met haar arcsinusaroma.

Deel V — De vernieuwingsidentiteit en de stelling van Pólya.

  1. Bewijs, door {S2n=0}\{S_{2n} = 0\} te partitioneren over het tijdstip van de eerste terugkeer, de vernieuwingsidentiteit

    un=k=1nfkunk(n1),en dusU(x)(1F(x))=1(0x<1),u_n = \sum_{k=1}^n f_k\,u_{n-k} \quad (n \geq 1), \qquad\text{en dus}\qquad U(x)\bigl(1 - F(x)\bigr) = 1 \quad (0 \leq x < 1),

    waarbij U(x)=n0unxnU(x) = \sum_{n\geq0}u_nx^n en F(x)=n1fnxnF(x) = \sum_{n\geq1}f_nx^n (verantwoord de stralen en het product van de reeksen met Hoofdstuk 11).

  2. Leid de recurrentietweedeling af: met x1x \to 1^- (monotone limieten van reeksen met niet-negatieve coëfficiënten) is

    nun=    nfn=1,\sum_n u_n = \infty \iff \sum_n f_n = 1 ,

    en toets haar aan de vragen 3, 7 (eerlijke wandeling) en 12 (scheve wandeling).

  3. (Dimensie 22) De eenvoudige wandeling op Z2\Z^2 neemt uniform stappen (±1,0)(\pm1, 0), (0,±1)(0, \pm1). Toon aan dat de gedraaide coördinaten Un=Xn+YnU_n = X_n + Y_n en Vn=XnYnV_n = X_n - Y_n onafhankelijke eerlijke wandelingen op Z\Z uitvoeren, leid af dat

    P(S2n(2)=(0,0))=un21πn,nun2=,\P\bigl(S^{(2)}_{2n} = (0,0)\bigr) = u_n^2 \sim \frac1{\pi n}, \qquad \sum_n u_n^2 = \infty ,

    en besluit met de vragen 21–22 (waarvan de bewijzen woordelijk overgaan) dat de wandeling op Z2\Z^2 recurrent is.

  4. (Dimensie 33) Neem voor de eenvoudige wandeling op Z3\Z^3 de lokale schatting P(S2n(3)=0)Cn3/2\P(S^{(3)}_{2n} = 0) \leq C\,n^{-3/2} aan (bewezen met de lokale limietstelling in het volume van bachelorjaar 3). Leid met Borel–Cantelli 1 af dat de wandeling op Z3\Z^3 transiënt is, en formuleer het volledige resultaat: de stelling van Pólya — de eenvoudige toevalswandeling is recurrent in de dimensies 11 en 22, en transiënt in dimensie 33 en hoger.
  5. Synthese. Som de exacte rol op van: het tellen van paden en de spiegeling; de monotone continuïteit; de onafhankelijkheid van disjuncte blokken worpen; Borel–Cantelli 1; de vernieuwingsidentiteit. Welk enkel analytisch feit (un1/πnu_n \sim 1/\sqrt{\pi n}, dus un=\sum u_n = \infty maar ook un2=\sum u_n^2 = \infty en n3/2<\sum n^{-3/2} < \infty) beslist in elke dimensie tussen recurrentie en transiëntie?
Oplossing

Oplossing van Probleem 21.1.

1. Een pad van lengte nn wordt bepaald door de verzameling van zijn opwaartse stappen; in kk eindigen betekent uu opwaartse en nun - u neerwaartse stappen met u(nu)=ku - (n - u) = k, dat wil zeggen u=n+k2u = \frac{n+k}2: mogelijk dan en slechts dan als n+kn + k even is en kn\abs k \leq n, en wel op (n(n+k)/2)\binom{n}{(n+k)/2} manieren. Elk specifiek pad is één punt van de eerlijke productmaat op nn worpen: kans 2n2^{-n}. Bijgevolg is P(Sn=k)=Nn(k)2n\P(S_n = k) = N_n(k)2^{-n}.

2. SnS_n heeft de pariteit van nn, dus S2n+10S_{2n+1} \neq 0; en un=N2n(0)4n=(2nn)4nu_n = N_{2n}(0)4^{-n} = \binom{2n}n4^{-n}. Waarden: u1=12u_1 = \frac12, u2=616=38u_2 = \frac6{16} = \frac38, u3=2064=516u_3 = \frac{20}{64} = \frac5{16}.

3. unun1=(2nn)4(2n2n1)=(2n)(2n1)4n2=2n12n<1\dfrac{u_n}{u_{n-1}} = \dfrac{\binom{2n}n}{4\binom{2n-2}{n-1}} = \dfrac{(2n)(2n-1)}{4n^2} = \dfrac{2n-1}{2n} < 1: dalend. Volgens Voorbeeld 6.14 is (2nn)4nπn\binom{2n}n \sim \frac{4^n}{\sqrt{\pi n}}, dus un1πn0u_n \sim \frac1{\sqrt{\pi n}} \to 0, en un\sum u_n divergeert door vergelijking met n1/2\sum n^{-1/2}.

4. Gegeven een pad van 11 naar kk dat 00 raakt, spiegel je zijn beginstuk (tot aan het eerste bezoek aan 00) in de horizontale as: het resultaat is een pad van 1-1 naar kk, en de bewerking is een involutie — elk pad van 1-1 naar k1k \geq 1 moet 00 kruisen, en zijn beginstuk terugspiegelen geeft het oorspronkelijke terug. Bijgevolg zijn er Nn1(k+1)N_{n-1}(k + 1) paden die raken (van 1-1 naar kk is de verplaatsing k+1k + 1). Een pad van 00 naar kk dat na tijdstip 00 boven 00 blijft, begint met een opwaartse stap en gaat dan in n1n - 1 stappen van 11 naar kk zonder 00 te raken: daarvan zijn er Nn1(k1)Nn1(k+1)N_{n-1}(k-1) - N_{n-1}(k+1).

5. Met m=n+k2m = \frac{n+k}2 en (n1m1)=mn(nm)\binom{n-1}{m-1} = \frac mn\binom nm, (n1m)=nmn(nm)\binom{n-1}{m} = \frac{n-m}n\binom nm:

Nn1(k1)Nn1(k+1)Nn(k)=(n1m1)(n1m)(nm)=m(nm)n=kn.\frac{N_{n-1}(k-1) - N_{n-1}(k+1)}{N_n(k)} = \frac{\binom{n-1}{m-1} - \binom{n-1}{m}}{\binom nm} = \frac{m - (n - m)}{n} = \frac kn .

Voor n=3n = 3, k=1k = 1: N3(1)=3N_3(1) = 3 paden (++++-, +++-+, ++-++), waarvan alleen ++++- positief blijft (+++-+ keert op tijdstip 22 naar 00 terug): één op drie, en kn=13\frac kn = \frac13.

6. Wegens de symmetrie is de kans 2P(Si>0 i2n)2\P(S_i > 0\ \forall i \leq 2n). Sommeren over het eindpunt 2k2k en vraag 4 gebruiken (met nn vervangen door 2n2n) geeft

P(Si>0 i)=22nk1(N2n1(2k1)N2n1(2k+1))=22nN2n1(1),\P(S_i > 0\ \forall i) = 2^{-2n}\sum_{k\geq1} \bigl(N_{2n-1}(2k-1) - N_{2n-1}(2k+1)\bigr) = 2^{-2n}\,N_{2n-1}(1),

een telescoperende som. Nu is N2n1(1)=(2n1n)N_{2n-1}(1) = \binom{2n-1}{n} en 2(2n1n)=(2nn)2\binom{2n-1}n = \binom{2n}n (Pascal), dus is de getoonde kans 222n(2n1n)=(2nn)4n=un2\cdot2^{-2n}\binom{2n-1}n = \binom{2n}n4^{-n} = u_n.

7. De gebeurtenissen Dn={Si0, i2n}D_n = \{S_i \neq 0,\ i \leq 2n\} dalen, met als doorsnede “nooit een terugkeer”; volgens de monotone continuïteit en vraag 6 is P(geen terugkeer)=limun=0\P(\text{geen terugkeer}) = \lim u_n = 0: de wandeling keert bijna zeker terug. Bovendien is fn=P(Dn1)P(Dn)=un1unf_n = \P(D_{n-1}) - \P(D_n) = u_{n-1} - u_n, en volgens vraag 3

un1un=un(2n2n11)=un2n1;n1fn=u0limun=1.u_{n-1} - u_n = u_n\Bigl(\frac{2n}{2n-1} - 1\Bigr) = \frac{u_n}{2n-1}; \qquad \sum_{n\geq1}f_n = u_0 - \lim u_n = 1 .

8. 2nfn=2n2n1unun2n\,f_n = \frac{2n}{2n-1}u_n \geq u_n, en un=\sum u_n = \infty (vraag 3): de reeks 2nfn\sum 2nf_n divergeert. De eerste terugkeer is zeker maar heeft geen eindige gemiddelde wachttijd — de wandeling is nulrecurrent, in de woordenschat die Hoofdstuk 22 zal leveren.

9. De gebeurtenis “minstens kk terugkeren” is de disjuncte aftelbare vereniging, over 0<n1<<nk0 < n_1 < \dots < n_k, van de gebeurtenissen “de eerste kk terugkeren gebeuren precies op de tijdstippen 2n1,,2nk2n_1, \dots, 2n_k”. Zulk een gebeurtenis is de doorsnede van kk gebeurtenissen die van de disjuncte blokken worpen [ ⁣[1,2n1] ⁣]\intint1{2n_1}, [ ⁣[2n1+1,2n2] ⁣]\intint{2n_1+1}{2n_2}, …afhangen, waarbij elk blok van een verse wandeling vraagt dat zij haar eerste terugkeer precies na het toegemeten aantal stappen maakt; wegens de onafhankelijkheid van de blokken is haar kans fn1fn2n1fnknk1f_{n_1}f_{n_2-n_1}\cdots f_{n_k-n_{k-1}}. Sommeren in pakketten (Hoofdstuk 7, alle termen niet-negatief) geeft

P(minstens k terugkeren)=(n1fn) ⁣k=1k=1.\P(\text{minstens }k\text{ terugkeren}) = \Bigl(\sum_{n\geq1}f_n\Bigr)^{\!k} = 1^k = 1 .

De gebeurtenissen dalen in kk, dus geeft de monotone continuïteit dat P(oneindig veel terugkeren)=1\P(\text{oneindig veel terugkeren}) = 1: recurrentie.

10. Volgens vraag 9 maakt de wandeling oneindig veel uitstapjes weg van 00. De eerste stap van elk uitstapje is een verse munt, onafhankelijk van alles ervoor: de kans dat de eerste mm uitstapjes alle neerwaarts beginnen, is 2m2^{-m}. Om 11 te bereiken heeft de wandeling maar één opwaarts beginnend uitstapje nodig (vanuit <0<0 moet zij door 00 voordat zij 11 bereikt, want de stappen zijn ±1\pm1), dus is P(nooit 1 raken)2m\P(\text{nooit }1\text{ raken}) \leq 2^{-m} voor elke mm: de wandeling raakt 11 bijna zeker. Ontbinden we over de (bijna zeker eindige) raaktijd, dan is de daar herstarte wandeling een verse wandeling die in 11 begint: met inductie raakt zij bijna zeker elke k1k \geq 1, en wegens de symmetrie elke k1k \leq -1. Herstarten wij ten slotte bij het eerste bezoek aan kk, dan is vraag 9 van toepassing op de verse wandeling: elke plaats wordt bijna zeker oneindig vaak bezocht.

11. De gebeurtenissen An={S2n=0}A_n = \{S_{2n} = 0\} zijn verre van onafhankelijk (in 00 zijn op tijdstip 2n2n maakt in 00 zijn op tijdstip 2n+22n + 2 veel waarschijnlijker dan un+1u_{n+1}), dus is Borel–Cantelli 2 niet beschikbaar, en het hele werk van Deel II was dan ook haar te vervangen. De andere richting heeft geen onafhankelijkheid nodig: als P(An)\sum\P(A_n) convergeert, dan levert Borel–Cantelli 1 bijna zeker eindig veel terugkeren. Die implicatie is de motor van elk bewijs van transiëntie hieronder.

12. Een terugkeer op tijdstip 2n2n vergt nn opwaartse en nn neerwaartse stappen: P(S2n=0)=(2nn)pnqn=un(4pq)n\P(S_{2n} = 0) = \binom{2n}np^nq^n = u_n(4pq)^n, en 4pq=1(pq)2<14pq = 1 - (p - q)^2 < 1 voor p12p \neq \frac12. Omdat un1u_n \leq 1, wordt de reeks P(S2n=0)\sum\P(S_{2n} = 0) gedomineerd door de meetkundige (4pq)n\sum(4pq)^n: convergent. Volgens Borel–Cantelli 1 is P(S2n=0 oneindig vaak)=0\P(S_{2n} = 0 \text{ oneindig vaak}) = 0: bijna zeker eindig veel terugkeren.

13. Voor n+kn + k even is P(Sn=k)=(nn+k2)pn+k2qnk2\P(S_n = k) = \binom{n}{\frac{n+k}2}p^{\frac{n+k}2}q^{\frac{n-k}2}; de binomiaalcoëfficiënt is hoogstens de centrale, en pn+k2qnk2=(pq)n/2(p/q)k/2p^{\frac{n+k}2}q^{\frac{n-k}2} = (pq)^{n/2}(p/q)^{k/2}, wat de gestelde grens 2n(pq)n/2(p/q)k/2=(4pq)n/2(p/q)k/2\leq 2^n(pq)^{n/2}(p/q)^{k/2} = (4pq)^{n/2}(p/q)^{k/2} geeft, sommeerbaar in nn omdat 4pq<1\sqrt{4pq} < 1. Borel–Cantelli 1: plaats kk wordt bijna zeker eindig vaak bezocht; de vereniging over kZk \in \Z van de uitzonderlijke nulgebeurtenissen is nog altijd nul (aftelbare subadditiviteit). Bijna zeker wordt elke plaats eindig vaak bezocht, dus verlaat de rij gehele getallen (Sn)(S_n) elk begrensd venster voorgoed: Sn\abs{S_n} \to \infty.

14. P(Sn0, 1n200)=u100=(200100)41001100π0.056\P(S_n \neq 0,\ 1 \leq n \leq 200) = u_{100} = \binom{200}{100}4^{-100} \approx \frac1{\sqrt{100\pi}} \approx 0.056: meer dan één kans op twintig dat 200200 eerlijke worpen nooit een gelijkstand geven. Het verval 1/πn1/\sqrt{\pi n} is tergend traag: de zekerheid van een gelijkstand (vraag 7) is verenigbaar met heel lange stukken zonder gelijkstand — een eerste voorproefje van de arcsinusverschijnselen van Deel IV.

15. Stel voorwaardelijk op de eerste stap. Is X1=+1X_1 = +1, dan is T1=1T_1 = 1, en f1=12f_1 = \frac12 klopt. Is X1=1X_1 = -1, dan moet de wandeling van 1-1 naar 11 klimmen; volgens de ontbinding in blokken splitst een eerste terugkeer naar 00 op tijdstip 2n2n als: één stap omlaag, en dan een verse wandeling die vanuit 1-1 voor het eerst 00 bereikt — gelijkwaardig: een verse wandeling die voor het eerst +1+1 bereikt — in 2n12n - 1 stappen, of de symmetrische gebeurtenis omhoog. Beide tekens dragen evenveel bij:

fn=212P(T1=2n1)=P(T1=2n1).f_n = 2\cdot\tfrac12\,\P(T_1 = 2n - 1) = \P(T_1 = 2n-1) .

Bijgevolg is P(T1<)=fn=1\P(T_1 < \infty) = \sum f_n = 1, terwijl n(2n1)fn=nun=\sum_n(2n - 1)f_n = \sum_n u_n = \infty volgens vraag 7: de wandeling bereikt 11 bijna zeker, in oneindige gemiddelde tijd.

16. Partitioneer {Mnk}\{M_n \geq k\} naar de eindwaarde Sn=mS_n = m. Voor mkm \geq k is de voorwaarde MnkM_n \geq k automatisch. Voor m<km < k spiegel je het pad na zijn eerste bezoek aan niveau kk: dit is een bijectie tussen {Mnk,Sn=m}\{M_n \geq k, S_n = m\} en {Sn=2km}\{S_n = 2k - m\} (elk pad dat in 2km>k2k - m > k eindigt, bezoekt kk; terugspiegelen is de inverse). Bijgevolg is

P(Mnk)=m>kP(Sn=m)+P(Sn=k)+m<kP(Sn=2km)=2P(Sn>k)+P(Sn=k).\P(M_n \geq k) = \sum_{m > k}\P(S_n = m) + \P(S_n = k) + \sum_{m < k}\P(S_n = 2k - m) = 2\P(S_n > k) + \P(S_n = k).

17. Op het even tijdstip 2n2n met k=1k = 1: P(S2n=1)=0\P(S_{2n} = 1) = 0 en P(S2n>1)=P(S2n2)\P(S_{2n} > 1) = \P(S_{2n} \geq 2), dus

P(M2n1)=2P(S2n2)=P(S2n2)+P(S2n2)=1un.\P(M_{2n} \geq 1) = 2\P(S_{2n} \geq 2) = \P(S_{2n} \geq 2) + \P(S_{2n} \leq -2) = 1 - u_n .

Dus P(Si0 i2n)=un\P(S_i \leq 0\ \forall i \leq 2n) = u_n: de wandeling staat in de eerste 2n2n stappen precies even vaak nooit voor als zij nooit gelijk staat (vraag 6) — twee heel verschillende gebeurtenissen, gedragen door dezelfde unu_n.

18. {L2n=2k}={S2k=0}{de wandeling van de worpen 2k+1,,2n heeft geen nulpunt}\{L_{2n} = 2k\} = \{S_{2k} = 0\} \cap \{\text{de wandeling van de worpen } 2k+1, \dots, 2n \text{ heeft geen nulpunt}\}. De twee gebeurtenissen hangen van disjuncte blokken worpen af en zijn dus onafhankelijk; de eerste heeft kans uku_k, de tweede unku_{n-k} volgens vraag 6 toegepast op de verse wandeling van 2n2k2n-2k stappen. Bijgevolg is P(L2n=2k)=ukunk\P(L_{2n} = 2k) = u_ku_{n-k}. Omdat L2nL_{2n} precies de waarden 0,2,,2n0, 2, \dots, 2n aanneemt, sommeren deze kansen tot 11: k=0nukunk=1\sum_{k=0}^nu_ku_{n-k} = 1, een binomiale identiteit geleverd door een kansrekenkundige partitie.

19. De symmetrie is onmiddellijk: ukunk=unkuku_ku_{n-k} = u_{n-k}u_k. Omdat uju_j daalt in jj, is het product ukunku_ku_{n-k} het kleinst voor centrale kk en het grootst in de uitersten k{0,n}k \in \{0, n\}, waar het gelijk is aan unu_n; kwantitatief is ukunk1πk(nk)u_ku_{n-k} \approx \frac1{\pi\sqrt{k(n-k)}} in het midden, tegenover un1πnu_n \approx \frac1{\sqrt{\pi n}} aan de randen. Voor n=5n = 5: P(L10=0)=P(L10=10)=u5=632560.246\P(L_{10} = 0) = \P(L_{10} = 10) = u_5 = \frac{63}{256} \approx 0.246, terwijl P(L10=4)=u2u3=38516=151280.117\P(L_{10} = 4) = u_2u_3 = \frac38\cdot\frac5{16} = \frac{15}{128} \approx 0.117. In een lang eerlijk spel valt de laatste gelijkstand het waarschijnlijkst vlakbij het begin of vlakbij het einde: één speler staat meestal gedurende enorme stukken voor, zonder enige scheefheid in de munt.

20. Het beeld: op tijdstip nn leeft de wandeling op de schaal n\sqrt n (de binomiale spreiding van vraag 3 — un1/πnu_n \sim 1/\sqrt{\pi n} is de hoogte van de centrale piek); zij keert met kans 11 oneindig vaak naar 00 terug (Deel II), en toch heeft de wachttijd tussen de terugkeren een divergent gemiddelde (vraag 8), en daarom kunnen afzonderlijke uitstapjes een positief aandeel van elke horizon innemen; navenant is de laatste gelijkstand van een spel met 2n2n stappen uitgesmeerd met de uiterste waarden als de waarschijnlijkste (vragen 18–19), en heeft nooit voorstaan dezelfde traag vervallende kans unu_n als nooit gelijk staan (vraag 17). Zekerheid in de limiet, volharding op elke eindige horizon: dat is de eerlijke wandeling.

21. Partitioneer {S2n=0}\{S_{2n} = 0\} (n1n \geq 1) naar de eerste terugkeertijd 2k2k, 1kn1 \leq k \leq n: het eerste blok van 2k2k worpen verwezenlijkt een eerste terugkeer, de overige 2n2k2n - 2k worpen verwezenlijken een terugkeer van een verse wandeling, en de blokken zijn onafhankelijk: un=k=1nfkunku_n = \sum_{k=1}^nf_ku_{n-k}. Beide reeksen U(x)=unxnU(x) = \sum u_nx^n en F(x)=fnxnF(x) = \sum f_nx^n hebben straal 1\geq 1 (coëfficiënten in [0,1]\intcc01), en het Cauchy-product (Hoofdstuk 11) geeft voor 0x<10 \leq x < 1

U(x)1=n1(k=1nfkunk)xn=F(x)U(x),dat wil zeggenU(x)(1F(x))=1.U(x) - 1 = \sum_{n\geq1}\Bigl(\sum_{k=1}^n f_ku_{n-k}\Bigr)x^n = F(x)\,U(x), \qquad\text{dat wil zeggen}\qquad U(x)\bigl(1 - F(x)\bigr) = 1 .

22. Als x1x \uparrow 1 stijgen U(x)U(x) en F(x)F(x) (niet-negatieve coëfficiënten); elke partiaalsom nNun\sum_{n\leq N}u_n is een limiet van nNunxnU(x)\sum_{n\leq N}u_nx^n \leq U(x), dus U(x)un(0,+]U(x) \uparrow \sum u_n \in \intoc0{+\infty}, en evenzo F(x)f=fnF(x) \uparrow f = \sum f_n. Is un=\sum u_n = \infty, dan 1F(x)=1/U(x)01 - F(x) = 1/U(x) \to 0, dus f=1f = 1. Is un=S<\sum u_n = S < \infty, dan 1f=1/S>01 - f = 1/S > 0, dus f<1f < 1. Controles: eerlijke wandeling, un=\sum u_n = \infty en f=1f = 1 (vragen 3, 7); scheve wandeling, un(4pq)n<\sum u_n(4pq)^n < \infty en navenant f=11/n0un(4pq)n<1f = 1 - 1/\sum_{n\geq0}u_n(4pq)^n < 1, in overeenstemming met de bijna-zekere eindigheid van het aantal terugkeren (vraag 12).

23. Voor de vier stappen (±1,0),(0,±1)(\pm1, 0), (0, \pm1) van de wandeling op Z2\Z^2 zijn de aangroeiingen van U=X+YU = X + Y en V=XYV = X - Y: (+,+)(+,+) voor (1,0)(1,0), (+,)(+,-) voor (0,1)(0,1), (,+)(-,+) voor (0,1)(0,-1) en (,)(-,-) voor (1,0)(-1,0) — elk tekenpaar met kans 14=1212\frac14 = \frac12\cdot\frac12: de twee coördinaatwandelingen (Un)(U_n) en (Vn)(V_n) zijn onafhankelijke eerlijke wandelingen op Z\Z. Omdat S2n(2)=(0,0)S^{(2)}_{2n} = (0,0) dan en slechts dan als U2n=0U_{2n} = 0 en V2n=0V_{2n} = 0, is

P(S2n(2)=(0,0))=un21πn,nun2=.\P\bigl(S^{(2)}_{2n} = (0,0)\bigr) = u_n^2 \sim \frac1{\pi n}, \qquad \sum_nu_n^2 = \infty .

De vernieuwingsidentiteit van vraag 21 en de tweedeling van vraag 22 gebruikten niets eendimensionaals (alleen de ontbinding over de eerste terugkeer en de onafhankelijkheid van disjuncte blokken), dus geeft un(2)=\sum u_n^{(2)} = \infty dat f(2)=1f^{(2)} = 1, en het argument van vraag 9 verscherpt het: de wandeling op Z2\Z^2 keert bijna zeker oneindig vaak naar de oorsprong terug.

24. Met de aangenomen grens P(S2n(3)=0)Cn3/2\P(S^{(3)}_{2n} = 0) \leq Cn^{-3/2} convergeert de reeks, en Borel–Cantelli 1 geeft bijna zeker eindig veel terugkeren: de wandeling op Z3\Z^3 is transiënt (en dezelfde grens met exponent d/2-d/2 handelt elke d3d \geq 3 af). Alles samen: de stelling van Pólya — de eenvoudige toevalswandeling is recurrent op Z\Z en Z2\Z^2, en transiënt op Zd\Z^d voor d3d \geq 3. Een dronken man vindt de weg naar huis; een dronken vogel misschien niet.

25. Het tellen van paden en de spiegeling brachten de exacte verdelingen voort (unu_n, de stemmenstelling, fnf_n, het maximum, het laatste nulpunt); de monotone continuïteit zette elke uitspraak over een limiet (“keert minstens eenmaal terug”, “oneindig vaak”) om in een limiet van kansen met eindige horizon; de onafhankelijkheid van disjuncte blokken dreef de vernieuwingsontbindingen aan (vragen 9, 18, 21) — zij is het aftelbare skelet van de markov-eigenschap; Borel–Cantelli 1 was het wapen voor de transiëntie (vragen 12–13, 24), zonder onafhankelijkheid nodig te hebben; en de vernieuwingsidentiteit ordende alles tot de tweedeling un=    \sum u_n = \infty \iff recurrentie. Het enige analytische ingrediënt is de lokale schatting un1/πnu_n \sim 1/\sqrt{\pi n}: haar kwadraat 1/(πn)1/(\pi n) divergeert nog altijd (dimensie 22, recurrent), terwijl n3/2n^{-3/2} convergeert (dimensie 33, transiënt) — de stelling van Pólya is uiteindelijk een uitspraak over de divergentie van nd/2\sum n^{-d/2}.