Mathematics · Boek 5 · Bachelor Year 3

Universitaire wiskunde — Bachelor jaar 3

Universitaire wiskunde — Bachelor jaar 3 · Bachelor Year 3

22Kansrekening: fundamenten en de wet van de grote aantallen

Het volume van bachelorjaar 2 bouwde de kansrekening op aftelbare ruimten; de maattheorie neemt nu elke beperking weg. Een kansruimte is een maatruimte met totale massa 11, toevalsveranderlijken zijn meetbare afbeeldingen, de verwachtingswaarde is de lebesgue-integraal — en meteen is het hele analytische arsenaal (Hoofdstukken 9, 10 en 11) op het toeval van toepassing. Dit hoofdstuk installeert het woordenboek, construeert oneindige rijen onafhankelijke toevalsveranderlijken (op [0,1]\intcc01, uit binaire cijfers: het toeval zit verstopt in de lebesguemaat), bewijst de lemma’s van Borel–Cantelli en de nul-een-wet van Kolmogorov, sorteert de convergentiewijzen uit, en bewijst de wet van de grote aantallen — de stelling die frequenties naar kansen laat convergeren en de statistiek mogelijk maakt. De weekendopgave geeft het bewijs van Etemadi van de sterke wet in haar definitieve L1L^1-vorm.

22.1 Het woordenboek

Definitie 22.1

Een kansruimte is een maatruimte (Ω,A,P)(\Omega, \mathcal A, \P) met P(Ω)=1\P(\Omega) = 1; de elementen van A\mathcal A heten gebeurtenissen, en een eigenschap geldt bijna zeker (b.z.) als haar gebeurtenis kans 11 heeft. Een toevalsveranderlijke is een meetbare afbeelding X ⁣:ΩRX \colon \Omega \to \R (of Rd\R^d: een toevalsvector); haar verdeling is de beeldkansmaat PX=XP\P_X = X_*\P op R\R (Oefening 11.9), vastgelegd door de verdelingsfunctie FX(t)=P(Xt)F_X(t) = \P(X \leq t) (Oefening 9.3). XX heeft dichtheid ff als PX=f ⁣dλ\P_X = f\,\dd\lambda; ze heet discreet als PX\P_X een aftelbare combinatie van diracmassa’s is. De verwachtingswaarde is

E[X]=ΩX ⁣dP(X0 of XL1(P)),\E[X] = \int_\Omega X\,\dd\P \qquad (X \geq 0 \text{ of } X \in L^1(\P)),

en de overdrachtsformule (Oefening 11.9) berekent haar in de verdeling: E[g(X)]=Rg ⁣dPX\E[g(X)] = \int_\R g\,\dd\P_X=g(xk)pk= \sum g(x_k)p_k in het discrete geval, =g(x)f(x) ⁣dx= \int g(x)f(x)\dd x in het geval met een dichtheid: de formules uit het volume van bachelorjaar 2, nu stellingen van één theorie. De variantie is V(X)=E[(XEX)2]=E[X2](EX)2\V(X) = \E[(X - \E X)^2] = \E[X^2] - (\E X)^2 voor XL2X \in L^2.

Voorbeeld 22.2

De standaardverdelingen en hun opmerkelijke transformaties: bernoulli B(p)\mathcal B(p), binomiaal B(n,p)\mathcal B(n, p), meetkundig, poisson P(λ)\mathcal P(\lambda) (discreet: de tabellen uit het volume van bachelorjaar 2 blijven geldig); uniform op [0,1]\intcc01 (de lebesguemaat zelf); exponentieel E(λ)\mathcal E(\lambda) (dichtheid λeλx1x>0\lambda\eu^{-\lambda x}\mathbf 1_{x>0}); de gaussische verdeling N(m,σ2)\mathcal N(m, \sigma^2) met dichtheid 1σ2πexp((xm)22σ2)\frac1{\sigma\sqrt{2\pi}}\exp\bigl(-\frac{(x - m)^2}{2\sigma^2}\bigr) — een kansdichtheid volgens Probleem 10.1, met gemiddelde mm en variantie σ2\sigma^2 (de gaussische momenten, Oefening 11.10).

Propositie 22.3 (Markov en Chebyshev)

Voor X0X \geq 0 en a>0a > 0 is P(Xa)EXa\P(X \geq a) \leq \frac{\E X}{a}; voor XL2X \in L^2 is P(XEXa)V(X)a2\P\bigl(\abs{X - \E X} \geq a\bigr) \leq \frac{\V(X)}{a^2}.

Bewijs. Oefening 10.5(a); Chebyshev is Markov toegepast op (XEX)2(X - \E X)^2.

22.2 Onafhankelijkheid

Definitie 22.4

Deel-σ\sigma-algebra’s A1,,AnA\mathcal A_1, \dots, \mathcal A_n \subseteq \mathcal A heten onafhankelijk als P(A1An)=P(Ai)\P(A_1\cap\dots\cap A_n) = \prod\P(A_i) voor alle AiAiA_i \in \mathcal A_i; gebeurtenissen zijn onafhankelijk als de σ\sigma-algebra’s {,Ai,Aic,Ω}\{\varnothing, A_i, A_i^c, \Omega\} dat zijn; en toevalsveranderlijken X1,,XnX_1, \dots, X_n als de σ\sigma-algebra’s σ(Xi)=Xi1(B(R))\sigma(X_i) = X_i^{-1}(\mathcal B(\R)) dat zijn. Een oneindige familie is onafhankelijk als elke eindige deelfamilie dat is.

Stelling 22.5

X1,,XnX_1, \dots, X_n zijn onafhankelijk dan en slechts dan als de verdeling van de vector (X1,,Xn)(X_1, \dots, X_n) de productmaat PX1PXn\P_{X_1}\otimes\cdots\otimes\P_{X_n} is. In dat geval is voor gi0g_i \geq 0 (of zodanig dat de producten integreerbaar zijn)

E[igi(Xi)]=iE[gi(Xi)],\E\Bigl[\prod_ig_i(X_i)\Bigr] = \prod_i\E[g_i(X_i)],

in het bijzonder E[XY]=EXEY\E[XY] = \E X\,\E Y en V(X1++Xn)=V(Xi)\V(X_1 + \dots + X_n) = \sum\V(X_i) voor onafhankelijke veranderlijken in L2L^2.

Bewijs. Zijn de XiX_i onafhankelijk, dan vallen de twee kansmaten P(X1,,Xn)\P_{(X_1,\dots,X_n)} en PXi\bigotimes\P_{X_i} samen op alle producten B1××BnB_1\times\dots\times B_n van borelverzamelingen — een π\pi-systeem dat B(Rn)\mathcal B(\R^n) voortbrengt (Propositie 11.2(b)) — en dus overal (Stelling 9.7). Omgekeerd ontbindt een productverdeling alle gebeurtenissen iXi1(Bi)\bigcap_iX_i^{-1}(B_i): onafhankelijkheid. De formule voor de verwachtingswaarde is dan Tonelli of Fubini (Stelling 11.5) via de overdrachtsformule; E[XY]=EXEY\E[XY] = \E X\E Y is het geval gi=idg_i = \mathrm{id}, en het kwadraat ontwikkelen geeft de additiviteit van de varianties (de kruistermen E[(XiEXi)(XjEXj)]=0\E[(X_i - \E X_i)(X_j - \E X_j)] = 0).

Stelling 22.6 (Bestaan van onafhankelijke rijen)

Op ([0,1],L,λ)\bigl(\intcc01, \mathcal L, \lambda\bigr) bestaat er een rij (Un)n1(U_n)_{n\geq1} onafhankelijke toevalsveranderlijken, elk uniform op [0,1]\intcc01. Bijgevolg bestaan er voor elke voorgeschreven verdelingen (μn)(\mu_n) op R\R onafhankelijke (Xn)(X_n) met PXn=μn\P_{X_n} = \mu_n.

Bewijs. Cijfers. Zij voor ω[0,1]\omega \in \intcc01 (bk(ω))(b_k(\omega)) haar binaire cijfers (ω=bk2k\omega = \sum b_k2^{-k}; kies de ontwikkeling die niet op louter enen eindigt — de dubbelzinnigheid betreft slechts een aftelbare, dus verwaarloosbare verzameling). Elke bkb_k is een toevalsveranderlijke ({bk=1}\{b_k = 1\} is een eindige vereniging van dyadische intervallen) en de vector (b1,,bm)(b_1, \dots, b_m) neemt elke waarde in {0,1}m\{0,1\}^m aan op een dyadisch interval van lengte 2m2^{-m}: de bkb_k zijn onafhankelijk en bernoulli(12)(\frac12).

Hergroeperen. Splits N\N^* in oneindig veel disjuncte oneindige verzamelingen (In)(I_n) (bijvoorbeeld met priemmachten, of met diagonalen); zij (kjn)j(k^n_j)_j een opsomming van InI_n en stel

Un=j1bkjn2j.U_n = \sum_{j\geq1} b_{k^n_j}\,2^{-j} .

Elke UnU_n is uniform: haar binaire cijfers zijn onafhankelijke eerlijke bits, dus P(Un[l2m,(l+1)2m))=2m\P(U_n \in [l2^{-m}, (l+1)2^{-m})) = 2^{-m} voor elk dyadisch interval, en de dyadische intervallen leggen de verdeling vast (Stelling 9.7). De UnU_n zijn onafhankelijk: het zijn functies van disjuncte blokken van de onafhankelijke familie (bk)(b_k) — formeel hangen de gebeurtenissen {UnDn}\{U_n \in D_n\} voor dyadische DnD_n af van eindig veel cijfers uit disjuncte verzamelingen, en ontbinden ze; het argument met het π\pi-systeem tilt dat op tot alle borelverzamelingen.

Willekeurige verdelingen. Zij Gn(u)=inf{t:Fμn(t)u}G_n(u) = \inf\{t : F_{\mu_n}(t) \geq u\} (de kwantielfunctie van de verdelingsfunctie FμnF_{\mu_n}); de sleutelgelijkwaardigheid Gn(u)t    uFμn(t)G_n(u) \leq t \iff u \leq F_{\mu_n}(t) (rechtscontinuïteit van FF, monotonie) toont dat Xn=Gn(Un)X_n = G_n(U_n) meetbaar is met P(Xnt)=P(UnFμn(t))=Fμn(t)\P(X_n \leq t) = \P(U_n \leq F_{\mu_n}(t)) = F_{\mu_n}(t): dus verdeling μn\mu_n; de onafhankelijkheid wordt geërfd (functies van onafhankelijke veranderlijken, Oefening 22.3).

Voorbeeld 22.7 (Het verjaardagsprobleem, eerlijk)

Onder nn mensen met onafhankelijke, uniform over N=365N = 365 dagen verdeelde verjaardagen is de kans dat alle verjaardagen verschillen

pn=k=1n1(1kN),p_n = \prod_{k=1}^{n-1}\Bigl(1 - \frac kN\Bigr),

door herhaald te conditioneren (of rechtstreeks: de gunstige N(N1)(Nn+1)N(N-1)\cdots(N - n + 1) op het totaal NnN^n, een telargument dat de productformule van de onafhankelijkheid streng maakt). Logaritmen nemen en ln(1x)=x+O(x2)-\ln(1 - x) = x + O(x^2) gebruiken geeft

lnpn=n(n1)2N+O(n3N2),duspnen2/2N.\ln p_n = -\frac{n(n-1)}{2N} + O\Bigl(\frac{n^3}{N^2}\Bigr), \qquad\text{dus}\qquad p_n \approx \eu^{-n^2/2N} .

Het kantelpunt pn=12p_n = \frac12 ligt bij n2Nln21.18Nn \approx \sqrt{2N\ln2} \approx 1.18\sqrt N: voor N=365N = 365 is n=23n = 23 (p23=0.4927p_{23} = 0.4927). Twee moralen. Ten eerste duiken botsingen tussen nn voorwerpen in NN dozen op op de schaal nNn \sim \sqrt N, niet nNn \sim N — de verjaardagsschaling die de botsingen van hashfuncties en de kostprijs N\sqrt N van verjaardagsaanvallen in de cryptografie regeert. Ten tweede is de berekening een sjabloon: de (n2)\binom n2 gebeurtenissen “dit paar botst” zijn niet onafhankelijk, en toch gedraagt het antwoord zich alsof ze het wel waren (e(n2)/N\eu^{-\binom n2/N} is precies de heuristiek met onafhankelijke paren) — een eerste geval van de poissonbenadering die in de weekendopgave van Hoofdstuk 23 streng wordt gemaakt (de ongelijkheid van Le Cam).

22.3 Borel–Cantelli en de nul-een-wet

Stelling 22.8 (Borel–Cantelli)

Zij (An)(A_n) gebeurtenissen en lim supAn=NnNAn\limsup A_n = \bigcap_N \bigcup_{n\geq N}A_n (“AnA_n treedt oneindig vaak op”).

  1. Is P(An)<\sum\P(A_n) < \infty, dan is P(lim supAn)=0\P(\limsup A_n) = 0.
  2. Is P(An)=\sum\P(A_n) = \infty en zijn de AnA_n onafhankelijk, dan is P(lim supAn)=1\P(\limsup A_n) = 1.

Bewijs. (1) is Oefening 9.4. (2): voor NMN \leq M geeft de onafhankelijkheid van de complementen (Oefening 22.3)

P(n=NMAnc)=n=NM(1P(An))exp(n=NMP(An))M0\P\Bigl(\bigcap_{n=N}^{M}A_n^c\Bigr) = \prod_{n=N}^M\bigl(1 - \P(A_n)\bigr) \leq \exp\Bigl(-\sum_{n=N}^M\P(A_n)\Bigr) \xrightarrow[M \to \infty]{} 0

(1xex1 - x \leq \eu^{-x}; de reeks divergeert). Dus is P(nNAn)=1\P\bigl(\bigcup_{n\geq N}A_n\bigr) = 1 voor elke NN, en heeft de dalende doorsnede over NN nog steeds kans 11 (continuïteit van boven, Propositie 9.6).

Stelling 22.9 (Nul-een-wet van Kolmogorov)

Zij (Xn)(X_n) onafhankelijk en T=Nσ(XN,XN+1,)\mathcal T = \bigcap_N\sigma(X_N, X_{N+1}, \dots) de staart-σ\sigma-algebra (gebeurtenissen die ongevoelig zijn voor elk eindig aantal van de XnX_n: de convergentie van Xn\sum X_n, van Snn\frac{S_n}n, waarden van lim sup\limsup’s, …). Dan heeft elke TTT \in \mathcal T P(T){0,1}\P(T) \in \{0, 1\}.

Bewijs. Leg NN vast. De σ\sigma-algebra’s σ(X1,,XN)\sigma(X_1, \dots, X_N) en σ(XN+1,)\sigma(X_{N+1}, \dots) zijn onafhankelijk: gebeurtenissen die van disjuncte blokken afhangen, ontbinden op de voortbrengende π\pi-systemen (cilinders iN{XiBi}\bigcap_{i\leq N}\{X_i \in B_i\} respectievelijk eindige voorwaarden op latere veranderlijken), en Dynkin (Stelling 9.4, tweemaal toegepast, telkens aan één kant) breidt de ontbinding uit. Een staartgebeurtenis TT ligt voor elke NN in σ(XN+1,)\sigma(X_{N+1}, \dots): dus is TT onafhankelijk van elke σ(X1,,XN)\sigma(X_1, \dots, X_N), en dus van de σ\sigma-algebra die ze voortbrengen, σ(X1,X2,)\sigma(X_1, X_2, \dots) (nog eens Dynkin: de vereniging van de σ(X1,,XN)\sigma(X_1,\dots,X_N) is een π\pi-systeem dat haar voortbrengt). Maar Tσ(X1,X2,)T \in \sigma(X_1, X_2, \dots) eveneens: TT is onafhankelijk van zichzelf, dus P(T)=P(TT)=P(T)2\P(T) = \P(T\cap T) = \P(T)^2: P(T){0,1}\P(T) \in \{0, 1\}.

22.4 Convergentiewijzen

Definitie 22.10

XnXX_n \to X bijna zeker als P(XnX)=1\P(X_n \to X) = 1; in kans als P(XnXε)0\P(\abs{X_n - X} \geq \varepsilon) \to 0 voor elke ε>0\varepsilon > 0; in LpL^p als EXnXp0\E\abs{X_n - X}^p \to 0.

Propositie 22.11

(a) b.z. convergentie impliceert convergentie in kans; (b) LpL^p-convergentie impliceert convergentie in kans; (c) convergentie in kans impliceert b.z. convergentie langs een deelrij; (d) geen enkele andere implicatie geldt in het algemeen.

Bewijs. (a) P(XnXε)P(supmnXmXε)P(lim sup{XmXε})=0\P(\abs{X_n - X} \geq \varepsilon) \leq \P\bigl(\sup_{m\geq n}\abs{X_m - X} \geq \varepsilon\bigr) \downarrow \P\bigl(\limsup\{\abs{X_m - X} \geq \varepsilon\}\bigr) = 0 bij b.z. convergentie (continuïteit van boven; de limsup-gebeurtenis sluit convergentie uit). (b) Markov: P(XnXε)εpEXnXp\P(\abs{X_n - X} \geq \varepsilon) \leq \varepsilon^{-p}\,\E\abs{X_n - X}^p. (c) Kies nkn_k met P(XnkX2k)2k\P(\abs{X_{n_k} - X} \geq 2^{-k}) \leq 2^{-k}; Borel–Cantelli (1) maakt dat uiteindelijk b.z. XnkX<2k\abs{X_{n_k} - X} < 2^{-k}. (d) De schrijfmachine (Oefening 12.3) op ([0,1],λ)(\intcc01, \lambda) convergeert in L1L^1 en in kans maar nergens puntsgewijs; n1(0,1/n)0n\mathbf 1_{\intoo0{1/n}} \to 0 b.z. maar niet in L1L^1; de details en de overige tegenvoorbeelden staan in Oefening 22.6.

22.5 De wet van de grote aantallen

Overal zijn (Xn)(X_n) onafhankelijk met dezelfde verdeling (i.i.d.), en Sn=X1++XnS_n = X_1 + \dots + X_n.

Stelling 22.12 (Zwakke wet van de grote aantallen)

Is X1L2X_1 \in L^2, met m=EX1m = \E X_1, dan is

P(Snnmε)V(X1)nε2n0:\P\Bigl(\Bigl|\frac{S_n}{n} - m\Bigr| \geq \varepsilon\Bigr) \leq \frac{\V(X_1)}{n\,\varepsilon^2} \xrightarrow[n\to\infty]{} 0 :

Snnm\frac{S_n}n \to m in kans (en in L2L^2).

Bewijs. ESnn=m\E\frac{S_n}n = m en V(Snn)=nV(X1)n2\V\bigl(\frac{S_n}n\bigr) = \frac{n\V(X_1)}{n^2} (Stelling 22.5); Chebyshev.

Stelling 22.13 (Sterke wet van de grote aantallen)

Is X1L1X_1 \in L^1, dan is

Snnnb.z.E[X1].\frac{S_n}{n} \xrightarrow[n\to\infty]{\text{b.z.}} \E[X_1].

We bewijzen dat hier onder de sterkere hypothese X1L4X_1 \in L^4; het algemene geval (L1L^1: het bewijs van Etemadi) is de weekendopgave.

Bewijs onder EX14<\E X_1^4 < \infty. Na centreren (XiXimX_i \mapsto X_i - m) mogen we m=0m = 0 aannemen. Ontwikkel:

E[Sn4]=i,j,k,lE[XiXjXkXl]=nE[X14]+3n(n1)(E[X12])2Cn2,\E[S_n^4] = \sum_{i,j,k,l}\E[X_iX_jX_kX_l] = n\,\E[X_1^4] + 3n(n-1)\,\bigl(\E[X_1^2]\bigr)^2 \leq C\,n^2 ,

want de onafhankelijkheid en het centreren doden elke term met een geïsoleerde factor (E[XiXjXkXl]=E[Xi]E[]=0\E[X_iX_jX_kX_l] = \E[X_i]\E[\cdots] = 0 tenzij de indices paren vormen: de enige overlevenden zijn de nn termen met i=j=k=li=j=k=l en de 3n(n1)3n(n-1) termen met twee verschillende paren). Markov:

P(Snnε)=P(Sn4n4ε4)Cn2n4ε4=Cε4n2,\P\Bigl(\Bigl|\frac{S_n}n\Bigr| \geq \varepsilon\Bigr) = \P\bigl(S_n^4 \geq n^4\varepsilon^4\bigr) \leq \frac{Cn^2}{n^4\varepsilon^4} = \frac{C}{\varepsilon^4n^2},

sommeerbaar: Borel–Cantelli (1) geeft voor elke rationale ε\varepsilon dat uiteindelijk b.z. Sn/n<ε\abs{S_n/n} < \varepsilon; doorsnijden over εQ+\varepsilon \in \Q_+^* (aftelbaar veel gebeurtenissen met kans 11) geeft Sn/n0S_n/n \to 0 b.z.

Voorbeeld 22.14 (Wat de sterke wet oplevert)

(a) Frequenties: voor i.i.d. muntworpen convergeert de waargenomen frequentie van kop b.z. naar pp — de empirische rechtvaardiging van de kansrekening zelf. (b) Monte carlo: voor gL1([0,1])g \in L^1(\intcc01) en (Un)(U_n) i.i.d. uniform (Stelling 22.6) is 1nkng(Uk)01g\frac1n\sum_{k\leq n}g(U_k) \to \int_0^1g b.z.: integralen door te bemonsteren, in elke dimensie, met de van de dimensie onafhankelijke snelheid n1/2\sim n^{-1/2} die in Hoofdstuk 23 precies wordt gemaakt. (c) Normale getallen: bijna elk reëel getal heeft in zijn binaire ontwikkeling asymptotische frequentie 12\frac12 aan enen (pas de sterke wet toe op de cijferveranderlijken van Stelling 22.6) — de stelling van Borel, een uitspraak over de meest alledaagse getallen bewezen met maat: Probleem 22.1 maakt haar in alle bases af.

Methode 22.15

De werkvolgorde voor asymptotische uitspraken over toevallige rijen: (1) Is de gebeurtenis een staartgebeurtenis? Dan is haar kans 00 of 11 (Stelling 22.9) en hoef je enkel te beslissen welke van beide. (2) Om b.z. uitspraken te bewijzen: Borel–Cantelli — sommeerbare kansen voor de “slechte” gebeurtenissen, via afschattingen van het type Markov of Chebyshev op de momenten die bestaan; de onafhankelijkheid is enkel voor de omgekeerde richting nodig. (3) Deelrij plus insluiting: bewijs de convergentie langs een hanteerbare deelrij en beheers de oscillatie ertussen met monotonie of maximaalongelijkheden — het skelet van het bewijs van Etemadi. (4) Voor limieten in verdeling: wacht op Hoofdstuk 23.

22.6 Oefeningen

Oefening 22.1

(a) Zij XX met een continue strikt stijgende verdelingsfunctie FF. Toon aan dat F(X)F(X) uniform is op [0,1]\intcc01, en dat G(U)FG(U) \sim F voor UU uniform en G=F1G = F^{-1}: simulatie door inversie. (b) Bereken de verdelingsfunctie en de dichtheid van X2X^2 voor XX uniform op [1,1]\intcc{-1}1, en van 1λlnU-\frac1\lambda\ln U voor UU uniform op (0,1)\intoo01.

Oplossing

Oplossing van Oefening 22.1.

(a) Voor u(0,1)u \in \intoo01: P(F(X)u)=P(XF1(u))=F(F1(u))=u\P(F(X) \leq u) = \P(X \leq F^{-1}(u)) = F(F^{-1}(u)) = u (de continuïteit en de strikte monotonie maken FF een bijectie op (0,1)\intoo01 met {F(X)u}={XF1(u)}\{F(X) \leq u\} = \{X \leq F^{-1}(u)\}): dus is F(X)F(X) uniform. Omgekeerd is P(G(U)t)=P(UF(t))=F(t)\P(G(U) \leq t) = \P(U \leq F(t)) = F(t): om een verdeling te simuleren, pas je de inverse verdelingsfunctie toe op een uniforme steekproef.

(b) Y=X2Y = X^2 met XX uniform op [1,1]\intcc{-1}1: voor t[0,1]t \in \intcc01 is FY(t)=P(tXt)=tF_Y(t) = \P(-\sqrt t \leq X \leq \sqrt t) = \sqrt t: dichtheid 12t1(0,1)\frac1{2\sqrt t}\mathbf 1_{\intoo01}. En P(1λlnUt)=P(Ueλt)=1eλt\P\bigl(-\frac1\lambda\ln U \leq t\bigr) = \P(U \geq \eu^{-\lambda t}) = 1 - \eu^{-\lambda t}: de exponentiële verdeling E(λ)\mathcal E(\lambda) — de inversie in actie.

Oefening 22.2

(a) Bereken het gemiddelde en de variantie van de poissonverdeling P(λ)\mathcal P(\lambda) en van de meetkundige verdeling met de overdrachtsformule. (b) Toon aan dat een positieve toevalsveranderlijke TT met P(T>t)>0\P(T > t) > 0 voor alle tt de geheugenloosheid P(T>t+sT>t)=P(T>s)\P(T > t + s \mid T > t) = \P(T > s) voor alle s,t0s, t \geq 0 bezit dan en slechts dan als TT exponentieel is. (De overlevingsfunctie voldoet aan de functionaalvergelijking van Cauchy; de monotonie vervangt de continuïteit.)

Oplossing

Oplossing van Oefening 22.2.

(a) Poisson: EX=k0keλλkk!=λ\E X = \sum_{k\geq0}k\,\eu^{-\lambda} \frac{\lambda^k}{k!} = \lambda en E[X(X1)]=λ2\E[X(X-1)] = \lambda^2, dus V=λ2+λλ2=λ\V = \lambda^2 + \lambda - \lambda^2 = \lambda. Meetkundig (P(X=k)=p(1p)k1\P(X = k) = p(1-p)^{k-1}): EX=1p\E X = \frac1p en V=1pp2\V = \frac{1-p}{p^2} (differentieer de meetkundige reeks tweemaal).

(b) G(t)=P(T>t)G(t) = \P(T > t) is niet-stijgend met G(0+)G(0^+)\dots; G ⁣:[0,)(0,1]G \colon \intco0\infty \to \intoc01, en de geheugenloosheid luidt G(t+s)=G(t)G(s)G(t + s) = G(t)G(s). Dan is G(nt)=G(t)nG(n t) = G(t)^n en G(t/n)=G(t)1/nG(t/n) = G(t)^{1/n}: dus G(q)=G(1)qG(q) = G(1)^q voor rationale q0q \geq 0; met G(1)=eλG(1) = \eu^{-\lambda} ((0,1)\in \intoo01: G(1)=1G(1) = 1 zou G1G \equiv 1 afdwingen, onmogelijk voor een eindige toevalsveranderlijke, en G(1)=0G(1) = 0 is door de hypothese uitgesloten) en een willekeurige tt tussen rationale getallen inklemmen (monotonie): G(t)=eλtG(t) = \eu^{-\lambda t} — de exponentiële verdeling. De omkering is een berekening.

Oefening 22.3 ★★

(a) Toon aan dat als X1,,XnX_1, \dots, X_n onafhankelijk zijn en de fif_i borelfuncties, de fi(Xi)f_i(X_i) onafhankelijk zijn. (b) Toon aan dat gebeurtenissen A1,,AnA_1, \dots, A_n onafhankelijk zijn dan en slechts dan als hun complementen dat zijn, en dan en slechts dan als de indicatoren 1Ai\mathbf 1_{A_i} onafhankelijke toevalsveranderlijken zijn. (c) (Paarsgewijs is zwakker) Twee eerlijke munten: A=A = de eerste is kop, B=B = de tweede is kop, C=C = de twee stemmen overeen. Toon aan dat A,B,CA, B, C paarsgewijs onafhankelijk zijn maar niet onafhankelijk.

Oplossing

Oplossing van Oefening 22.3.

(a) σ(fi(Xi))=fi(Xi)1(B)Xi1(B)=σ(Xi)\sigma(f_i(X_i)) = f_i(X_i)^{-1}(\mathcal B) \subseteq X_i^{-1}(\mathcal B) = \sigma(X_i) (fif_i is borel), en deel-σ\sigma-algebra’s van onafhankelijke σ\sigma-algebra’s zijn onafhankelijk (de definiërende identiteit geldt a fortiori).

(b) σ(Ai)={,Ai,Aic,Ω}=σ(Aic)=σ(1Ai)\sigma(A_i) = \{\varnothing, A_i, A_i^c, \Omega\} = \sigma(A_i^c) = \sigma(\mathbf 1_{A_i}): alle drie de uitspraken beweren de onafhankelijkheid van dezelfde σ\sigma-algebra’s. (Dat de ontbinding over de AiA_i zich naar de complementen voortplant, is het argument met het λ\lambda-systeem binnen de gelijkwaardigheid in Definitie 22.4 — of rechtstreeks in- en uitsluiting.)

(c) P(A)=P(B)=P(C)=12\P(A) = \P(B) = \P(C) = \frac12; en op paren is AB=AC=BCA\cap B = A\cap C = B\cap C: elke doorsnede is “twee keer kop” of analoog, met kans 14\frac14: paarsgewijs onafhankelijk. Maar P(ABC)=P(KK)=1418\P(A\cap B\cap C) = \P(\text{KK}) = \frac14 \neq \frac18: niet onafhankelijkCC ligt door AA en BB vast.

Oefening 22.4 ★★

(a) (De oneindige aap) Een i.i.d. rij uniforme toetsaanslagen op een eindig alfabet bevat b.z. elke eindige tekst oneindig vaak: bewijs dat met Borel–Cantelli (2) op disjuncte blokken. (b) (Reeksen) Zij voor i.i.d. eerlijke bits RnR_n de lengte van de reeks enen die op plaats nn begint. Toon aan dat b.z. Rn(1+ε)log2nR_n \geq (1+\varepsilon)\log_2n slechts eindig vaak, en Rnlog2nR_n \geq \log_2 n oneindig vaak (beide helften van Borel–Cantelli; ga voor de tweede over op disjuncte blokken om onafhankelijkheid te winnen): de langste reeks in de eerste nn cijfers groeit als log2n\log_2n.

Oplossing

Oplossing van Oefening 22.4.

(a) Zij de tekst TT van lengte LL en q=aLq = a^{-L} (aa de grootte van het alfabet). De gebeurtenissen Ek={E_k = \{de plaatsen kL+1,,(k+1)LkL+1, \dots, (k+1)L spellen T}T\} zijn onafhankelijk (disjuncte blokken i.i.d. letters), elk met kans q>0q > 0: dus P(Ek)=\sum\P(E_k) = \infty, en Borel–Cantelli (2) geeft b.z. oneindig veel voorkomens.

(b) Boven: P(Rn(1+ε)log2n)2(1+ε)log2n=n(1+ε)\P\bigl(R_n \geq (1+\varepsilon)\log_2n\bigr) \leq 2^{-(1+\varepsilon)\log_2n} = n^{-(1+\varepsilon)}, sommeerbaar: volgens Borel–Cantelli (1) zijn er b.z. slechts eindig veel zulke nn. Onder: pak disjuncte blokken — het jj-de van lengte j=log2sj\ell_j = \lceil\log_2s_j\rceil beginnend in sj=i<jis_j = \sum_{i<j}\ell_i; de gebeurtenissen “blok jj bestaat enkel uit enen” zijn onafhankelijk met kans 2j1sj1jlog2j2^{-\ell_j} \asymp \frac1{s_j} \asymp \frac1{j\log_2 j}, waarvan de som divergeert: Borel–Cantelli (2) geeft oneindig veel blokken met enkel enen, dat wil zeggen Rsjlog2sjR_{s_j} \geq \log_2 s_j oneindig vaak. Samen: de maximale reekslengte in de eerste nn cijfers is b.z. (1+o(1))log2n(1 + o(1))\log_2n.

Oefening 22.5 ★★

Zij (Xn)(X_n) onafhankelijk. (a) Toon aan dat de convergentiestraal van Xnzn\sum X_n z^n b.z. een constante is (mogelijk 00 of \infty). (b) Toon aan dat P(Xn convergeert){0,1}\P(\sum X_n \text{ convergeert}) \in \{0, 1\} en P(Sn/nm){0,1}\P(S_n/n \to m) \in \{0,1\}. (c) Geef een gebeurtenis over (Xn)(X_n) die geen staartgebeurtenis is, en ga na dat de nul-een-wet er kan falen.

Oplossing

Oplossing van Oefening 22.5.

(a) R=(lim supXn1/n)1R = \bigl(\limsup\abs{X_n}^{1/n}\bigr)^{-1} blijft ongewijzigd als eindig veel XnX_n worden aangepast: voor elke NN is RR dus σ(XN,XN+1,)\sigma(X_N, X_{N+1}, \dots)-meetbaar, dat wil zeggen staartmeetbaar. Dan heeft elke gebeurtenis {Rc}\{R \leq c\} kans 00 of 11 (Stelling 22.9), zodat de verdelingsfunctie van RR enkel de waarden 00 en 11 aanneemt: ze springt in één enkel punt c0[0,+]c_0 \in \intcc0{+\infty}, en R=c0R = c_0 b.z.

(b) De convergentie van Xn\sum X_n en van Snn\frac{S_n}n is ongevoelig voor het wijzigen van eindig veel termen (voor de tweede: de gewijzigde termen dragen O(1/n)0O(1/n) \to 0 bij): staartgebeurtenissen; nul-een-wet.

(c) {X1>0}\{X_1 > 0\} hangt van X1X_1 af: voor i.i.d. tekens (P(X1=±1)=12\P(X_1 = \pm1) = \frac12) is haar kans 12{0,1}\frac12 \notin \{0,1\} — geen tegenspraak, want het is geen staartgebeurtenis.

Oefening 22.6 ★★

Geef op ([0,1],λ)(\intcc01, \lambda) — met bewijzen — toevalsveranderlijken zodat: (a) Xn0X_n \to 0 in kans en in elke LpL^p, maar nergens b.z.; (b) Xn0X_n \to 0 b.z. maar in geen enkele LpL^p; (c) Xn0X_n \to 0 in L1L^1 maar niet in L2L^2; (d) en toon aan: convergeert XnXX_n \to X in kans met XnYL1\abs{X_n} \leq Y \in L^1, dan is XnXX_n \to X in L1L^1 (deelrijen + gedomineerde convergentie + de truc met de deeldeelrij).

Oplossing

Oplossing van Oefening 22.6.

Werk op ([0,1],λ)(\intcc01, \lambda). (a) De schrijfmachine 1In\mathbf 1_{I_n} (Oefening 12.3): Xnpp=λ(In)0\norm{X_n}_p^p = \lambda(I_n) \to 0 (alle p<p < \infty), en dus ook in kans; in elke ω\omega keren de waarden 00 en 11 allebei terug: nergens puntsgewijze convergentie. (b) Xn=n1(0,1/n)0X_n = n\mathbf 1_{\intoo0{1/n}} \to 0 buiten 00, maar Xnpn11/p1\norm{X_n}_p \geq n^{1 - 1/p} \geq 1. (c) Xn=n1(0,1/n)X_n = \sqrt n\,\mathbf 1_{\intoo0{1/n}}: EXn=n1/20\E\abs{X_n} = n^{-1/2} \to 0 en EXn2=1\E X_n^2 = 1. (d) Haal uit elke deelrij (convergentie in kans) een verdere deelrij die b.z. convergeert (Propositie 22.11(c)); de gedomineerde convergentie geeft er L1L^1-convergentie langs, met dezelfde limiet XX. Elke deelrij van de getallenrij EXnX\E\abs{X_n - X} heeft dus een deeldeelrij die naar 00 gaat: de hele rij gaat naar 00.

Oefening 22.7 ★★

Een opiniepeiling schat een onbekende verhouding pp met de empirische frequentie p^n\hat p_n van nn onafhankelijke trekkingen. (a) Chebyshev: toon aan dat P(p^npε)14nε2\P(\abs{\hat p_n - p} \geq \varepsilon) \leq \frac1{4n\varepsilon^2} (gebruik p(1p)14p(1-p) \leq \frac14). (b) Hoeveel trekkingen waarborgen met deze grens een fout 3%\leq 3\% met kans 95%\geq 95\%? (Het ware antwoord, via Hoofdstuk 23, is ongeveer 10701070: Chebyshev is eerlijk maar grof.)

Oplossing

Oplossing van Oefening 22.7.

(a) p^n=Snn\hat p_n = \frac{S_n}n met SnS_n binomiaal: V(p^n)=p(1p)n14n\V(\hat p_n) = \frac{p(1-p)}n \leq \frac1{4n}, en Chebyshev (Propositie 22.3) geeft de grens. (b) Los 14n(0.03)20.05\frac1{4n(0.03)^2} \leq 0.05 op: n140.00090.055556n \geq \frac{1}{4\cdot0.0009\cdot0.05} \approx 5556. De centrale limietstelling zal n1070n \approx 1070 rechtvaardigen voor dezelfde waarborg: Chebyshev betaalt haar algemeenheid met een factor 5\approx 5.

Oefening 22.8 ★★★

(Bernstein) Definieer voor fC([0,1])f \in \mathcal C(\intcc01) de bernsteinveelterm Bnf(x)=k=0n(nk)xk(1x)nkf(kn)B_nf(x) = \sum_{k=0}^n\binom nkx^k(1-x)^{n-k}f\bigl(\frac kn\bigr). (a) Herken Bnf(x)=E[f(Snn)]B_nf(x) = \E\bigl[f\bigl(\frac {S_n}n\bigr)\bigr] voor SnS_n binomiaal B(n,x)\mathcal B(n, x). (b) Bewijs dat BnffB_nf \to f uniform op [0,1]\intcc01: splits op {Snnxδ}\{\abs{\frac{S_n}n - x} \leq \delta\} en haar complement, met de uniforme continuïteit en Chebyshev met de uniforme grens V(Snn)14n\V(\frac{S_n}n) \leq \frac1{4n}. (c) Besluit: een tweede, kanstheoretisch bewijs van de benaderingsstelling van Weierstrass (Gevolg 7.16), met de expliciete snelheid Bnff32ωf(n1/2)\norm{B_nf - f}_\infty \leq \frac32\,\omega_f(n^{-1/2}) voor de continuïteitsmodulus ωf\omega_f — bewijs minstens de vorm O(ωf(n1/2))O(\omega_f(n^{-1/2})).

Oplossing

Oplossing van Oefening 22.8.

(a) Is SnB(n,x)S_n \sim \mathcal B(n, x), dan geeft de overdrachtsformule E[f(Snn)]=k(nk)xk(1x)nkf(kn)=Bnf(x)\E\bigl[f(\frac{S_n}n)\bigr] = \sum_k\binom nkx^k(1-x)^{n-k}f(\frac kn) = B_nf(x).

(b)–(c) Zij ω=ωf\omega = \omega_f de continuïteitsmodulus (f(u)f(v)ω(uv)\abs{f(u) - f(v)} \leq \omega(\abs{u - v}), en ω(cδ)(1+c)ω(δ)\omega(c \delta) \leq (1 + c)\,\omega(\delta) door stappen aaneen te schakelen). Dan is voor elke δ>0\delta > 0

f(u)f(x)(1+(ux)2δ2)ω(δ)\abs{f(u) - f(x)} \leq \Bigl(1 + \frac{(u - x)^2}{\delta^2}\Bigr)\omega(\delta)

(is uxδ\abs{u - x} \leq \delta, dan duidelijk; anders is ω(ux)(1+uxδ)ω(δ)(1+(ux)2δ2)ω(δ)\omega(\abs{u-x}) \leq (1 + \frac{\abs{u-x}}\delta) \omega(\delta) \leq (1 + \frac{(u-x)^2}{\delta^2}) \omega(\delta)). Neem verwachtingswaarden in u=Snnu = \frac{S_n}n:

Bnf(x)f(x)(1+V(Sn/n)δ2)ω(δ)(1+14nδ2)ω(δ);\abs{B_nf(x) - f(x)} \leq \Bigl(1 + \frac{\V(S_n/n)}{\delta^2}\Bigr)\omega(\delta) \leq \Bigl(1 + \frac{1}{4n\delta^2}\Bigr)\omega(\delta) ;

met δ=n1/2\delta = n^{-1/2}: Bnff54ω(n1/2)32ω(n1/2)0\norm{B_nf - f}_\infty \leq \frac54\,\omega\bigl(n^{-1/2}\bigr) \leq \frac32\,\omega\bigl(n^{-1/2}\bigr) \to 0 (uniforme continuïteit op de compacte verzameling): een kanstheoretische stelling van Weierstrass, met een expliciete en uniforme snelheid.

Oefening 22.9 ★★★

(Plaatjesverzamelaar) Kaarten van nn soorten worden uniform met teruglegging getrokken; zij TnT_n het aantal trekkingen tot alle soorten gezien zijn. (a) Schrijf Tn=k=1nτkT_n = \sum_{k=1}^{n}\tau_k met τk\tau_k meetkundig met parameter nk+1n\frac{n - k + 1}n en de τk\tau_k onafhankelijk, en leid af dat ETn=nHnnlnn\E T_n = n\,H_n \sim n\ln n (HnH_n het harmonische getal) en V(Tn)π26n2\V(T_n) \leq \frac{\pi^2}6n^2. (b) Chebyshev: Tnnlnn1\frac{T_n}{n\ln n} \to 1 in kans. (c) Verfijn met Borel–Cantelli: toon rechtstreeks aan dat P(Tn>βnlnn)n1β\P(T_n > \beta n\ln n) \leq n^{1 - \beta} voor β>1\beta > 1 (de somgrens op de gebeurtenis dat een soort na βnlnn\beta n\ln n trekkingen wordt gemist, met 1xex1 - x \leq \eu^{-x}), en leid af dat langs n=2mn = 2^m b.z. uiteindelijk TnβnlnnT_n \leq \beta n\ln n voor elke β>2\beta > 2.

Oplossing

Oplossing van Oefening 22.9.

(a) Zijn k1k - 1 soorten verzameld, dan is elke trekking nieuw met kans pk=nk+1np_k = \frac{n-k+1}n: τk\tau_k is meetkundig (pk)(p_k), en de τk\tau_k zijn onafhankelijk (dat zijn de trekkingen). Sommen: ETn=knnk+1=nHnnlnn\E T_n = \sum_k\frac n{n-k+1} = nH_n \sim n\ln n; en V(Tn)=1pkpk2n2j=1n1j2π26n2\V(T_n) = \sum\frac{1 - p_k}{p_k^2} \leq n^2\sum_{j=1}^n\frac1{j^2} \leq \frac{\pi^2}6n^2.

(b) Chebyshev: P(TnnHnεnlnn)π2n2/6ε2n2ln2n0\P\bigl(\abs{T_n - nH_n} \geq \varepsilon n\ln n\bigr) \leq \frac{\pi^2n^2/6}{\varepsilon^2n^2\ln^2n} \to 0, en nHnnlnn1\frac{nH_n}{n\ln n} \to 1: dus Tnnlnn1\frac{T_n}{n\ln n} \to 1 in kans.

(c) Somgrens: Tn>tT_n > t betekent dat een soort na t\lceil t\rceil trekkingen ongezien is, dus P(Tn>t)n(11n)tnet/n\P(T_n > t) \leq n(1 - \frac1n)^{t} \leq n\,\eu^{-t/n}; bij t=βnlnnt = \beta n\ln n geeft dat n1β\leq n^{1 - \beta}. Voor β>1\beta > 1 is m2m(1β)<\sum_m 2^{m(1-\beta)} < \infty: Borel–Cantelli geeft langs n=2mn = 2^m b.z. uiteindelijk TnβnlnnT_n \leq \beta n\ln n — in het bijzonder voor elke β>2\beta > 2 zoals gevraagd (langs de deelrij werkt elke β>1\beta > 1).

Oefening 22.10 ★★

Met de constructie met cijfers (Stelling 22.6): (a) ga met een rechtstreekse berekening na dat U=b2k2kU = \sum b_{2k}2^{-k} (de cijfers met even index van een uniforme ω\omega) uniform is en onafhankelijk van V=b2k12kV = \sum b_{2k-1}2^{-k}; (b) leid een meetbare bijectie op nulverzamelingen na af tussen [0,1]\intcc01 en [0,1]2\intcc01^2 die de maat behoudt, en becommentarieer: één uniform toevalsgetal bevat er twee (en aftelbaar veel) onafhankelijke — vergelijk met de peanokromme (Probleem 6.1), die wel surjectiviteit haalde maar geen maatbehoud of injectiviteit.

Oplossing

Oplossing van Oefening 22.10.

(a) De cijfers met even index (b2k)k(b_{2k})_k zijn i.i.d. eerlijke bits (een deelfamilie van de onafhankelijke cijferfamilie), dus geeft U=kb2k2kU = \sum_kb_{2k}2^{-k} aan elk dyadisch interval de juiste kans (als in Stelling 22.6): uniform; en evenzo VV; en (U,V)(U, V) hangen van disjuncte cijferblokken af: onafhankelijk (ontbinding op dyadische rechthoeken, dan Dynkin).

(b) Φ(ω)=(U(ω),V(ω))\Phi(\omega) = (U(\omega), V(\omega)) is meetbaar met Φλ=λλ=λ2\Phi_*\lambda = \lambda\otimes\lambda = \lambda_2 (overeenstemming op dyadische rechthoeken plus uniciteit). De cijfers afwisselen definieert een inverse buiten de (verwaarloosbare) verzameling dyadische rationale getallen in een van beide factoren: een maatbehoudende bijectie tussen deelverzamelingen van volle maat van [0,1]\intcc01 en [0,1]2\intcc01^2. Contrast met Peano (Probleem 6.1): de continuïteit dwong surjectiviteit af zonder injectiviteit; de continuïteit inruilen voor louter meetbaarheid koopt een maatisomorfisme — de dimensie is onzichtbaar voor de maattheorie en zichtbaar voor de topologie.

Oefening 22.11 ★★

(Records) Zij (Xn)n1(X_n)_{n\geq1} i.i.d. met continue verdelingsfunctie, en zeg dat er op tijdstip nn een record optreedt als Xn>max(X1,,Xn1)X_n > \max(X_1, \dots, X_{n-1}) (tijdstip 11 is een record). Zij RnR_n de recordindicator. (a) Toon aan dat P(Rn=1)=1n\P(R_n = 1) = \frac1n (wegens de symmetrie is elk van de n!n! ordeningen van X1,,XnX_1, \dots, X_n even waarschijnlijk, en gelijke waarden hebben kans 00). (b) Toon aan dat de RnR_n onafhankelijk zijn (tel de ordeningen die met voorgeschreven recordposities verenigbaar zijn, of beargumenteer dat de onderlinge volgorde van X1,,Xn1X_1, \dots, X_{n-1} onafhankelijk is van de rang van XnX_n ertussen). (c) Leid met Borel–Cantelli (Stelling 22.8, beide helften) af dat er b.z. oneindig veel records optreden, maar dat records op opeenvolgende tijdstippen n,n+1n, n+1 oneindig vaak optreden met kans — beslis welke! — en bereken nP(Rn=1,Rn+1=1)\sum_n\P(R_n = 1, R_{n+1} = 1).

Oplossing

Oplossing van Oefening 22.11.

(a) De continuïteit van de verdeling maakt gelijke waarden tot nulgebeurtenissen (als in de argumenten met ordestatistieken uit het hoofdstuk), en de n!n! onderlinge ordeningen van (X1,,Xn)(X_1, \dots, X_n) zijn verwisselbaar, dus even waarschijnlijk. Rn=1R_n = 1 betekent dat het maximum op de laatste plaats zit: kans (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n.

(b) Leg nn vast en conditioneer op de onderlinge ordening van X1,,Xn1X_1, \dots, X_{n-1}: XnX_n in een van de nn mogelijke rangplaatsen invoegen is uniform en onafhankelijk van die ordening (de verwisselbaarheid van het nn-tal). Dus is RnR_n (de gebeurtenis “XnX_n neemt de bovenste plaats in”) onafhankelijk van de hele recordgeschiedenis (R1,,Rn1)(R_1, \dots, R_{n-1}), die een functie is van de onderlinge ordening van de eerste n1n - 1 veranderlijken. Inductie geeft de volledige onafhankelijkheid met P(Rn=1)=1n\P(R_n = 1) = \frac1n.

(c) P(Rn=1)=1n=\sum\P(R_n = 1) = \sum\frac1n = \infty met onafhankelijkheid: de tweede helft van Borel–Cantelli geeft b.z. oneindig veel records (records houden nooit op — maar ze worden logaritmisch zeldzamer: E[#recordsn]=Hnlnn\E[\#\text{records} \leq n] = H_n \approx \ln n). Opeenvolgende records: P(Rn=Rn+1=1)=1n(n+1)\P(R_n = R_{n+1} = 1) = \frac1{n(n+1)} (onafhankelijkheid), en

n1n(n+1)=n(1n1n+1)=1<:\sum_n\frac1{n(n+1)} = \sum_n\Bigl(\frac1n - \frac1{n+1}\Bigr) = 1 < \infty :

de eerste helft van Borel–Cantelli is van toepassing — er treden b.z. slechts eindig veel paren opeenvolgende records op.

Oefening 22.12 ★★

(De langste reeks kop) Werp oneindig vaak met een eerlijke munt, en zij LnL_n de lengte van de langste reeks opeenvolgende koppen binnen de eerste nn worpen. (a) Toon aan dat voor elke ε>0\varepsilon > 0 b.z. uiteindelijk Ln(1+ε)log2nL_n \leq (1 + \varepsilon)\log_2n (de kans dat er onder de eerste nn worpen een reeks van lengte \ell begint, is hoogstens n2n2^{-\ell}; Borel–Cantelli langs n=2kn = 2^k). (b) Toon aan dat b.z. uiteindelijk Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n (hak de eerste nn worpen in n/\lfloor n/\ell\rfloor disjuncte blokken van lengte =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil; de blokken zijn onafhankelijk, elk volledig kop met kans 22^{-\ell}, en de kans dat geen enkel blok volledig kop is, is hoogstens exp(n2/)\exp(-n2^{-\ell}/\ell); sommeer opnieuw langs n=2kn = 2^k). (c) Besluit Lnlog2n1\frac{L_n}{\log_2n} \to 1 b.z.: in een miljoen eerlijke worpen mag je een reeks van ongeveer 2020 koppen verwachten — en een gegevensverzameling zonder zo’n reeks is waarschijnlijk verzonnen.

Oplossing

Oplossing van Oefening 22.12.

(a) Een reeks van lengte \ell die op plaats ini \leq n begint, heeft kans 22^{-\ell}; somgrens: P(Ln)n2\P(L_n \geq \ell) \leq n2^{-\ell}. Met n=(1+ε)log2n\ell_n = (1 + \varepsilon)\log_2n: P(Lnn)nε\P(L_n \geq \ell_n) \leq n^{-\varepsilon}. Langs n=2kn = 2^k is k2kε<\sum_k2^{-k\varepsilon} < \infty, dus b.z. uiteindelijk L2k<(1+ε)kL_{2^k} < (1+\varepsilon)k (Borel–Cantelli); voor algemene nn kies je 2k1<n2k2^{k-1} < n \leq 2^k en gebruik je de monotonie van LnL_n plus log22k1log2n\log_22^{k-1} \leq \log_2n: LnL2k<(1+ε)k(1+ε)kk1log2nL_n \leq L_{2^k} < (1 + \varepsilon)k \leq (1 + \varepsilon)\frac{k}{k-1} \log_2n, en de extra factor wordt opgeslorpt door ε\varepsilon licht te vergroten.

(b) Met =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil en m=n/m = \lfloor n/\ell\rfloor disjuncte blokken: de blokken zijn onafhankelijk, elk volledig kop met kans 2n(1ε)/22^{-\ell} \geq n^{-(1-\varepsilon)}/2, dus is

P(Ln<)(12)mexp(m2)exp(cnεlog2n)\P(L_n < \ell) \leq \bigl(1 - 2^{-\ell}\bigr)^{m} \leq \exp\bigl(-m2^{-\ell}\bigr) \leq \exp\Bigl(-c\,\frac{n^{\varepsilon}}{\log_2n}\Bigr)

voor een constante c>0c > 0 en grote nn. Die kansen zijn sommeerbaar langs n=2kn = 2^k (zelfs over alle nn): Borel–Cantelli geeft b.z. uiteindelijk Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n (de monotonie vult de gaten tussen de 2k2^k op als in (a), zonder schade).

(c) Beide grenzen langs een rij ε=1j\varepsilon = \frac1j, met de aftelbaar vele gebeurtenissen van volle maat doorsneden: Lnlog2n1\frac{L_n}{\log_2n} \to 1 b.z. Voor n=106n = 10^6 is log2n19.9\log_2n \approx 19.9 — een reeks van 20\approx 20 koppen is geen verdachte anomalie maar een wiskundige zekerheid, en haar afwezigheid is een aanwijzing dat een mens “toeval” heeft nagebootst (mensen durven zelden meer dan 55 of 66 koppen op rij op te schrijven).

22.7 Probleem: het bewijs van Etemadi van de sterke wet

Probleem 22.1

Weekendopgave — de sterke wet van de grote aantallen voor i.i.d. integreerbare veranderlijken

De sterke wet van Kolmogorov — SnnEX1\frac{S_n}n \to \E X_1 b.z. voor i.i.d. XnL1X_n \in L^1 — had lang enkel ingewikkelde bewijzen; in 1981 vond N. Etemadi er een van opvallende zuinigheid, dat niets voorbij dit hoofdstuk gebruikt (en de onafhankelijkheid zelfs verzwakt tot paarsgewijze onafhankelijkheid). Wij volgen het. Zij (Xn)(X_n) paarsgewijs onafhankelijk, identiek verdeeld en integreerbaar; m=EX1m = \E X_1 en Sn=X1++XnS_n = X_1 + \dots + X_n.

Deel I — Herleidingen.

  1. Toon aan dat het volstaat Xn0X_n \geq 0 te behandelen (splits Xn=Xn+XnX_n = X_n^+ - X_n^-: ga na dat de twee helften opnieuw paarsgewijs onafhankelijk, identiek verdeeld en integreerbaar zijn). Neem voortaan Xn0X_n \geq 0 aan.
  2. (Afknotting) Zij Yn=Xn1XnnY_n = X_n\,\mathbf 1_{X_n \leq n} en Sn=Y1++YnS_n^* = Y_1 + \dots + Y_n. Toon aan dat

    n1P(XnYn)=n1P(X1>n)E[X1]<\sum_{n\geq1}\P(X_n \neq Y_n) = \sum_{n\geq1}\P(X_1 > n) \leq \E[X_1] < \infty

    (Oefening 11.3), en leid met Borel–Cantelli af dat SnSnn0\frac{S_n - S_n^*}{n} \to 0 b.z.: het volstaat dus Snnm\frac{S^*_n}n \to m b.z. te bewijzen.

  3. Toon aan dat EYn=E[X11X1n]m\E Y_n = \E\bigl[X_1\mathbf 1_{X_1\leq n}\bigr] \to m (monotone convergentie), en dus 1nknEYkm\frac1n\sum_{k\leq n}\E Y_k \to m (Cesàro): het volstaat dus SnESnn0\frac{S_n^* - \E S_n^*}{n} \to 0 b.z. te bewijzen.

Deel II — De variantieafschatting.

  1. Toon aan dat

    V(Yn)E[Yn2]=E[X121X1n]\V(Y_n) \leq \E[Y_n^2] = \E\bigl[X_1^2\,\mathbf 1_{X_1 \leq n}\bigr]

    en, met de laagjesformule (Propositie 11.8), de sleutelgrens

    n1V(Yn)n2n11n2E[X121X1n]CE[X1]<\sum_{n\geq1}\frac{\V(Y_n)}{n^2} \leq \sum_{n\geq1}\frac1{n^2}\, \E\bigl[X_1^2\mathbf 1_{X_1\leq n}\bigr] \leq C\,\E[X_1] < \infty

    (verwissel de som en de verwachtingswaarde — Tonelli voor reeksen — en schat nx1n22max(x,1)\sum_{n \geq x}\frac1{n^2} \leq \frac2{\max(x,1)} af voor de inwendige schatting x2nxn22xx^2\sum_{n\geq x}n^{-2} \leq 2x).

Deel III — Convergentie langs meetkundige deelrijen. Leg α>1\alpha > 1 vast en zij kj=αjk_j = \lfloor\alpha^j\rfloor.

  1. Toon met de paarsgewijze onafhankelijkheid (de varianties tellen op, Stelling 22.5 — ga na dat de additiviteit van de varianties enkel paarsgewijze onafhankelijkheid vergt) en Chebyshev aan dat voor elke ε>0\varepsilon > 0

    j1P(SkjESkjkjε)1ε2j11kj2nkjV(Yn)=1ε2n1V(Yn)j:kjn1kj2.\sum_{j\geq1}\P\Bigl(\Bigl| \frac{S^*_{k_j} - \E S^*_{k_j}}{k_j}\Bigr| \geq \varepsilon\Bigr) \leq \frac1{\varepsilon^2}\sum_{j\geq1}\frac1{k_j^2} \sum_{n\leq k_j}\V(Y_n) = \frac1{\varepsilon^2}\sum_{n\geq1}\V(Y_n) \sum_{j\,:\,k_j\geq n}\frac1{k_j^2} .
  2. Toon aan dat j:kjnkj2Cαn2\sum_{j : k_j \geq n}k_j^{-2} \leq \frac{C_\alpha}{n^2} (meetkundige reeksen; pas op met het gehele deel: zorg voor iets als kjαj2k_j \geq \frac{\alpha^j}2 zodra αj2\alpha^j \geq 2), en besluit met vraag 4 en Borel–Cantelli:

    SkjESkjkjjb.z.0,en dusSkjkjm b.z.\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \xrightarrow[j\to\infty]{\text{b.z.}} 0, \qquad\text{en dus}\qquad \frac{S^*_{k_j}}{k_j} \to m \ \text{b.z.}

Deel IV — Insluiting en besluit.

  1. Toon voor kjnkj+1k_j \leq n \leq k_{j+1} met de monotonie van SnS^*_n (niet-negatieve termen!) aan dat

    kjkj+1Skjkj    Snn    kj+1kjSkj+1kj+1,\frac{k_j}{k_{j+1}}\,\frac{S^*_{k_j}}{k_j} \;\leq\; \frac{S^*_n}{n} \;\leq\; \frac{k_{j+1}}{k_j}\,\frac{S^*_{k_{j+1}}}{k_{j+1}},

    en leid b.z. af:

    mαlim infSnnlim supSnnαm.\frac m\alpha \leq \liminf\frac{S^*_n}n \leq \limsup\frac{S^*_n}n \leq \alpha\,m .
  2. Laat α1\alpha \downarrow 1 gaan langs een rij en besluit Snnm\frac{S_n^*}n \to m b.z., en dus (Deel I) de sterke wet van de grote aantallen:

     Snnnb.z.E[X1]. \boxed{\ \frac{S_n}{n} \xrightarrow[n\to\infty]{\text{b.z.}} \E[X_1].\ }
  3. Waar precies volstond de paarsgewijze (in plaats van de volledige) onafhankelijkheid? Som de drie plaatsen op waar hypothesen van het onafhankelijkheidstype werden ingeroepen.

Deel V — Dividenden.

  1. (De normale getallen van Borel) Toon aan dat λ\lambda-bijna elke x[0,1]x \in \intcc01 normaal is in elke basis b2b \geq 2: elk cijfer 0,,b10, \dots, b-1 verschijnt met asymptotische frequentie 1b\frac1b (leg bb en een cijfer vast, pas de sterke wet toe op de indicatorveranderlijken — verantwoord dat de cijfers in basis bb van een uniforme veranderlijke i.i.d. uniform op {0,,b1}\{0,\dots,b-1\} zijn, als in Stelling 22.6 — en snijd daarna de aftelbaar vele gebeurtenissen met kans één). Geef één expliciet niet-normaal getal, en bedenk: de stelling beweert de normaliteit van bijna alle getallen, en toch blijft de normaliteit van 2\sqrt2 of π\pi bewijzen een open probleem.
  2. (Monte carlo, gewaarborgd) Verantwoord de methode van Voorbeeld 22.14(b) volledig voor gL1([0,1]d)g \in L^1(\intcc01^d): construeer de i.i.d. uniforme steekproef op [0,1]d\intcc01^d met Stelling 22.6 en Oefening 22.10, en formuleer wat de sterke wet oplevert.

Deel VI — Wat volledige onafhankelijkheid oplevert: maximaalongelijkheden en toevallige reeksen. Etemadi geeft enkel paarsgewijze onafhankelijkheid uit; de resterende delen buiten de volledige (onderlinge) versie uit. Zij (Zn)(Z_n) onafhankelijke gecentreerde veranderlijken in L2L^2 en Sk=Z1++ZkS_k = Z_1 + \dots + Z_k (een nieuwe notatie, los van de XnX_n hierboven).

  1. (Maximaalongelijkheid van Kolmogorov) Bewijs voor ε>0\varepsilon > 0 dat

    P(max1knSkε)    1ε2k=1nV(Zk):\P\Bigl(\max_{1\leq k\leq n}\abs{S_k} \geq \varepsilon\Bigr) \;\leq\; \frac1{\varepsilon^2}\sum_{k=1}^n\V(Z_k) :

    de prijs van Chebyshev koopt het maximum (deel de gebeurtenis op naar de eerste index kk met Skε\abs{S_k} \geq \varepsilon; schrijf op dat stuk Sn2Sk2+2Sk(SnSk)S_n^2 \geq S_k^2 + 2S_k(S_n - S_k) en gebruik de onafhankelijkheid van de coalities (Z1,,Zk)(Z_1, \dots, Z_k) en (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n), Stelling 22.5). Wijs de stap aan waar paarsgewijze onafhankelijkheid niet langer zou volstaan.

  2. (Eenreeksstelling van Khinchin–Kolmogorov) Leid af: is nV(Zn)<\sum_n\V(Z_n) < \infty, dan convergeert nZn\sum_nZ_n bijna zeker (toon aan dat de partiële sommen b.z. een cauchyrij vormen: laat mm \to \infty gaan in de maximaalongelijkheid toegepast op ZN+1,,ZN+mZ_{N+1}, \dots, Z_{N+m}, en daarna NN \to \infty).
  3. (Rademacherreeksen) Zij (εn)(\varepsilon_n) i.i.d. tekens met P(εn=±1)=12\P(\varepsilon_n = \pm1) = \frac12 (Stelling 22.6), en zij (xn)(x_n) een rij reële getallen. Toon aan dat nxnεn\sum_nx_n\varepsilon_n b.z. convergeert zodra nxn2<\sum_nx_n^2 < \infty; toon ook aan dat, wat (xn)(x_n) ook is, de kans dat nxnεn\sum_nx_n\varepsilon_n convergeert 00 of 11 is (Stelling 22.9).
  4. De omkering, elementair. Stel Tn=knxkεkT_n = \sum_{k\leq n}x_k\varepsilon_k en sn2=knxk2s_n^2 = \sum_{k\leq n}x_k^2, en neem aan dat sns_n \to \infty. (a) Bewijs de ongelijkheid van Paley–Zygmund: voor Z0Z \geq 0 met EZ2<\E Z^2 < \infty en 0<θ<10 < \theta < 1 is

    P(Z>θEZ)    (1θ)2(EZ)2EZ2\P\bigl(Z > \theta\,\E Z\bigr) \;\geq\; (1 - \theta)^2\,\frac{(\E Z)^2}{\E Z^2}

    (splits EZ\E Z op het niveau θEZ\theta\E Z en pas Cauchy–Schwarz toe op het bovenste stuk). (b) Toon aan dat ETn43sn4\E T_n^4 \leq 3s_n^4. (c) Leid P(Tn>sn2)316\P\bigl(\abs{T_n} > \frac{s_n}2\bigr) \geq \frac3{16} af en besluit dat nxnεn\sum_nx_n\varepsilon_n b.z. divergeert; en dus de dichotomie

    nxnεn convergeert b.z.    nxn2<.\sum_nx_n\varepsilon_n\ \text{convergeert b.z.} \iff \sum_nx_n^2 < \infty .
  5. (Toevallige harmonische reeks) Besluit dat nεnns\sum_n\frac{\varepsilon_n}{n^s} b.z. convergeert dan en slechts dan als s>12s > \frac12. Voor 12<s1\frac12 < s \leq 1 convergeert de reeks b.z. terwijl nns=\sum_nn^{-s} = \infty: toevallige tekens brengen opheffing ter sterkte van een vierkantswortel voort — vergelijk met de alternerende reeks n(1)nns\sum_n\frac{(-1)^n}{n^s}, die voor elke s>0s > 0 convergeert.

Deel VII — Concentratie: de ongelijkheid van Hoeffding. De sterke wet zegt dat Snnm\frac{S_n}n \to m; concentratieongelijkheden zeggen hoe onwaarschijnlijk een afwijking is bij elke vaste nn.

  1. (Lemma van Hoeffding) (a) Toon aan dat coshλeλ2/2\cosh\lambda \leq \eu^{\lambda^2/2} voor alle λR\lambda \in \R, door de twee reeksen term voor term te vergelijken. (b) Zij ZZ gecentreerd met aZba \leq Z \leq b en a<ba < b. Toon aan dat

    EeλZexp(λ2(ba)28)\E\,\eu^{\lambda Z} \leq \exp\Bigl(\frac{\lambda^2(b - a)^2}8\Bigr)

    (begrens eλz\eu^{\lambda z} op [a,b]\intcc ab met haar koorde, neem verwachtingswaarden, en bestudeer φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t) met p=abap = \frac{-a}{b-a} en t=λ(ba)t = \lambda(b - a): toon aan dat φ(0)=φ(0)=0\varphi(0) = \varphi'(0) = 0 en φ14\varphi'' \leq \frac14).

  2. (Ongelijkheid van Hoeffding) Zij X1,,XnX_1, \dots, X_n onafhankelijk met aiXibia_i \leq X_i \leq b_i en Sn=X1++XnS_n = X_1 + \dots + X_n. Bewijs voor t>0t > 0

    P(SnESnt)exp(2t2i=1n(biai)2),\P\bigl(S_n - \E S_n \geq t\bigr) \leq \exp\Bigl(\frac{-2t^2}{\sum_{i=1}^n(b_i - a_i)^2}\Bigr),

    en dezelfde grens voor de onderste staart (exponentiële Chebyshev: begrens Eeλ(SnESn)\E\,\eu^{\lambda(S_n - \E S_n)} met de onafhankelijkheid en vraag 17, en optimaliseer daarna over λ>0\lambda > 0).

  3. (De sterke wet, begrensd geval, met een snelheid) Zij de XiX_i i.i.d. met waarden in [a,b]\intcc ab en m=EX1m = \E X_1. Toon aan dat

    P(Snnmε)2exp(2nε2(ba)2)\P\Bigl(\Bigl|\frac{S_n}n - m\Bigr| \geq \varepsilon\Bigr) \leq 2\exp\Bigl(\frac{-2n\varepsilon^2}{(b - a)^2}\Bigr)

    en vind Snnm\frac{S_n}n \to m b.z. terug met Borel–Cantelli: een tweede bewijs van de sterke wet voor begrensde veranderlijken — geen afknotting, een exponentiële snelheid bij elke eindige nn, maar begrensde termen en volledige onafhankelijkheid. Vergelijk de hypothesen met die van Etemadi.

  4. (Monte carlo, gewaarborgd bij vaste nn) Zij g ⁣:[0,1]d[0,1]g \colon \intcc01^d \to \intcc01 meetbaar en (Uk)(U_k) de i.i.d. uniforme steekproef uit vraag 11. Toon voor gegeven ε,δ>0\varepsilon, \delta > 0 aan dat

    nlog(2/δ)2ε2    P(1nk=1ng(Uk)g ⁣dλdε)δ,n \geq \frac{\log(2/\delta)}{2\varepsilon^2} \implies \P\Bigl(\Bigl|\frac1n\sum_{k=1}^ng(U_k) - \int g\,\dd\lambda_d\Bigr| \geq \varepsilon\Bigr) \leq \delta,

    en evalueer de drempel voor ε=δ=102\varepsilon = \delta = 10^{-2}. De grens bevat dd niet: vergelijk met vraag 11 en met deterministische roosters.

Deel VIII — Hoe groot is een toevalswandeling? Op weg naar de geïtereerde logaritme. Zij Sn=ε1++εnS_n = \varepsilon_1 + \dots + \varepsilon_n de eenvoudige toevalswandeling, opgebouwd uit i.i.d. eerlijke tekens.

  1. (Subgaussische staarten) Toon aan dat EeλSn=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} en leid voor x>0x > 0 af dat

    P(Snx)ex2/(2n),P(Snx)2ex2/(2n).\P(S_n \geq x) \leq \eu^{-x^2/(2n)}, \qquad \P(\abs{S_n} \geq x) \leq 2\,\eu^{-x^2/(2n)} .
  2. Leid met Borel–Cantelli af dat

    lim supnSn2nlogn1b.z.\limsup_{n\to\infty}\frac{\abs{S_n}} {\sqrt{2n\log n}} \leq 1 \quad\text{b.z.}

    (sommeer voor η>0\eta > 0 de staartgrenzen in x=(1+η)2nlognx = (1 + \eta)\sqrt{2n\log n}, en snijd daarna over η=1p\eta = \frac1p). In het bijzonder leeft de wandeling op de schaal n\sqrt n van de centrale limietstelling, op een logaritmische factor na — ver onder de grove grens Snn\abs{S_n} \leq n.

  3. Toon langs de verdubbelende deelrij nj=2jn_j = 2^j aan dat

    lim supjSnj2njloglognj1b.z.,\limsup_{j\to\infty}\frac{S_{n_j}} {\sqrt{2n_j\log\log n_j}} \leq 1 \quad\text{b.z.},

    en bedenk: de wet van de geïtereerde logaritme (Khinchin; Hartman–Wintner voor algemene gecentreerde termen in L2L^2) zegt dat

    lim supnSn2nloglogn=1b.z.\limsup_{n\to\infty}\frac{S_n} {\sqrt{2n\log\log n}} = 1 \quad\text{b.z.}

    Leg precies uit wat de zonet bewezen schatting langs een deelrij scheidt van de bovenste helft van die uitspraak (men moet maxnjnnj+1Sn\max_{n_j \leq n \leq n_{j+1}}S_n binnen elk blok beheersen, wat een maximaalongelijkheid op de exponentiële schaal vergt) en ga kwantitatief na dat de ongelijkheid van vraag 12 daarvoor te zwak is. De onderste helft berust op het tweede lemma van Borel–Cantelli toegepast op onafhankelijke blokken; beide helften zijn eerlijk materiaal van bachelorjaar 3 voor een aparte cursus kansrekening.

  4. (Uniforme afwijking over een eindige klasse) Zij A1,,ANA_1, \dots, A_N gebeurtenissen in een herhaalbaar experiment, en schat elke kans met haar empirische frequentie p^i\hat p_i over nn i.i.d. herhalingen. Toon door de ongelijkheid van Hoeffding met een somgrens te combineren aan dat

    P(maxiNp^iP(Ai)>ε)    2Ne2nε2,\P\Bigl(\max_{i\leq N}\,\abs{\hat p_i - \P(A_i)} > \varepsilon\Bigr) \;\leq\; 2N\,\eu^{-2n\varepsilon^2},

    en leid de regel voor de steekproefgrootte af: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2} waarborgt dat alle NN schattingen tegelijk ε\varepsilon-nauwkeurig zijn met kans 1δ\geq 1 - \delta. Bereken nn voor N=106N = 10^6, ε=0.01\varepsilon = 0.01 en δ=0.05\delta = 0.05: de logaritmische prijs van de uniformiteit.

  5. (Het toevallige harmonische venster) Toon door de twee helften van de theorie van de toevallige reeksen te combineren aan dat voor i.i.d. tekens (εn)(\varepsilon_n) de reeks nεnnα\sum_n\frac{\varepsilon_n}{n^\alpha} b.z. convergeert als α>12\alpha > \frac12 en b.z. divergeert als α12\alpha \leq \frac12; contrasteer met de absolute convergentie (die α>1\alpha > 1 vergt): op het venster α(12,1]\alpha \in \intoc{\frac12}1 is de convergentie een werkelijk kanstheoretisch verschijnsel — opheffing, geen grootte.
Oplossing

Oplossing van Probleem 22.1.

1. Xn±X_n^{\pm} zijn borelfuncties van XnX_n: ze blijven paarsgewijs onafhankelijk (Oefening 22.3(a)) en identiek verdeeld en integreerbaar, met EX1=EX1+EX1\E X_1 = \E X_1^+ - \E X_1^-. Geldt de stelling voor niet-negatieve veranderlijken, pas haar dan toe op beide helften en trek af: Snn=Sn+nSnnEX1+EX1=m\frac{S_n}n = \frac{S_n^+}n - \frac{S_n^-}n \to \E X_1^+ - \E X_1^- = m b.z.

2. P(XnYn)=P(Xn>n)=P(X1>n)\P(X_n \neq Y_n) = \P(X_n > n) = \P(X_1 > n) (gelijke verdelingen), en nP(X1>n)nP(X1n)EX1<\sum_n\P(X_1 > n) \leq \sum_n\P(X_1 \geq n) \leq \E X_1 < \infty (Oefening 11.3(a)). Borel–Cantelli (1): b.z. is Xn=YnX_n = Y_n voor alle grote nn, dus is SnSnS_n - S_n^* uiteindelijk constant in nn: SnSnn0\frac{S_n - S_n^*}n \to 0 b.z., en de twee genormeerde sommen delen hun asymptotische gedrag.

3. X11X1nX1X_1\mathbf 1_{X_1 \leq n} \nearrow X_1: de monotone convergentie geeft EYnm\E Y_n \to m; en de cesàrogemiddelden van een convergente rij convergeren naar dezelfde limiet: ESnn=1nknEYkm\frac{\E S_n^*}n = \frac1n\sum_{k\leq n}\E Y_k \to m. Het volstaat dus SnESnn0\frac{S^*_n - \E S^*_n}{n} \to 0 b.z. te bewijzen.

4. V(Yn)EYn2=E[X121X1n]\V(Y_n) \leq \E Y_n^2 = \E[X_1^2\mathbf 1_{X_1\leq n}]. Met Tonelli voor reeksen:

nE[X121X1n]n2=E[X12 ⁣ ⁣nmax(X1,1) ⁣1n2]E[X124max(X1,1)]4E[X1]<,\sum_n\frac{\E[X_1^2\mathbf 1_{X_1\leq n}]}{n^2} = \E\Bigl[X_1^2\!\!\sum_{n \geq \max(X_1, 1)}\!\frac1{n^2} \Bigr] \leq \E\Bigl[X_1^2\cdot\frac{4}{\max(X_1,1)}\Bigr] \leq 4\,\E[X_1] < \infty,

met nxn24x\sum_{n\geq x}n^{-2} \leq \frac4x voor x1x \geq 1 (voor x2x \geq 2: 1x12x\leq \frac1{x-1} \leq \frac2x; voor 1x<21 \leq x < 2: π264x\leq \frac{\pi^2}6 \leq \frac4x, want 4x>2\frac4x > 2), en X12/max(X1,1)X1X_1^2/\max(X_1, 1) \leq X_1 in beide gevallen X11X_1 \gtrless 1.

5. De paarsgewijze onafhankelijkheid geeft E[(YiEYi)(YjEYj)]=0\E[(Y_i - \E Y_i)(Y_j - \E Y_j)] = 0 voor iji \neq j (de productformule voor twee veranderlijken), dus tellen de varianties op: V(Sk)=nkV(Yn)\V(S^*_k) = \sum_{n\leq k}\V(Y_n). Chebyshev op elke kjk_j en sommeren:

jP(SkjESkjεkj)1ε2j1kj2nkjV(Yn)=1ε2nV(Yn) ⁣ ⁣j:kjn ⁣1kj2\sum_j\P\Bigl(\abs{S^*_{k_j} - \E S^*_{k_j}} \geq \varepsilon k_j\Bigr) \leq \frac1{\varepsilon^2}\sum_j\frac1{k_j^2}\sum_{n\leq k_j}\V(Y_n) = \frac1{\varepsilon^2}\sum_n\V(Y_n)\!\!\sum_{j : k_j\geq n}\!\frac1{k_j^2}

(Tonelli voor de niet-negatieve dubbele reeks).

6. kj=αjαj2k_j = \lfloor\alpha^j\rfloor \geq \frac{\alpha^j}2 (geldig zodra αj1\alpha^j \geq 1, dus voor alle j0j \geq 0: xx2\lfloor x\rfloor \geq \frac x2 voor x1x \geq 1). Bijgevolg is

j:kjn1kj24j:αjnα2j41α21n2=Cαn2,\sum_{j : k_j \geq n}\frac1{k_j^2} \leq 4\sum_{j : \alpha^j \geq n}\alpha^{-2j} \leq \frac{4}{1 - \alpha^{-2}}\cdot\frac1{n^2} = \frac{C_\alpha}{n^2},

(meetkundige reeks vanaf de eerste jj met αjn\alpha^j \geq n). Samen met de vragen 4 en 5 is de dubbele som eindig; en Borel–Cantelli (1), voor elke rationale ε\varepsilon toegepast en doorsneden, geeft SkjESkjkj0\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \to 0 b.z., en met vraag 3: Skjkjm\frac{S^*_{k_j}}{k_j} \to m b.z.

7. Yn0Y_n \geq 0 maakt nSnn \mapsto S^*_n niet-dalend: voor kjnkj+1k_j \leq n \leq k_{j+1} is

Skjkj+1SnnSkj+1kj,\frac{S^*_{k_j}}{k_{j+1}} \leq \frac{S^*_n}{n} \leq \frac{S^*_{k_{j+1}}}{k_j},

wat na kjkj+1\frac{k_j}{k_{j+1}} en kj+1kj\frac{k_{j+1}}{k_j} in te voegen de weergegeven insluiting is. Omdat kj+1kjα\frac{k_{j+1}}{k_j} \to \alpha, geeft vraag 6 b.z.

mαlim infnSnnlim supnSnnαm.\frac m\alpha \leq \liminf_n\frac{S^*_n}n \leq \limsup_n\frac{S^*_n}n \leq \alpha m .

8. Pas vraag 7 toe voor α=1+1p\alpha = 1 + \frac1p met pNp \in \N^*: aftelbaar veel b.z.-gebeurtenissen; op hun doorsnede geeft pp \to \infty dat limSnn=m\lim\frac{S^*_n}n = m b.z. Met de vragen 1 tot 3 is SnnEX1\frac{S_n}n \to \E X_1 b.z.: de sterke wet van de grote aantallen, onder paarsgewijze onafhankelijkheid.

9. Hypothesen van het onafhankelijkheidstype doken driemaal op: (i) de additiviteit van de varianties (vraag 5) — paarsgewijs volstaat; (ii) de gelijke verdeling, in de sommen bij de afknotting (vraag 2) en de berekening van de verwachtingswaarde (vraag 3) — helemaal geen onafhankelijkheid; (iii) Borel–Cantelli (1) (vragen 2 en 6) — geldig zonder enige onafhankelijkheid. De volledige onderlinge onafhankelijkheid werd nooit ingeroepen: de waarneming van Etemadi.

10. Leg een basis bb en een cijfer rr vast. De cijfers (dk)(d_k) in basis bb van een uniforme ω\omega zijn i.i.d. uniform op {0,,b1}\{0, \dots, b-1\} (elke waarde van de cijfervector beslaat een interval van lengte bmb^{-m}: woordelijk het argument van Stelling 22.6). De sterke wet toegepast op de i.i.d. begrensde veranderlijken 1dk=r\mathbf 1_{d_k = r} geeft: b.z. gaat de frequentie van het cijfer rr naar 1b\frac1b. De aftelbaar vele paren (b,r)(b, r) doorsnijden: bijna elk getal is eenvoudig normaal in elke basis. Een expliciet niet-normaal getal: x=0.1001001002x = 0.100100100\ldots_2 (de frequentie van de enen is 1312\frac13 \neq \frac12). Het contrast is ontnuchterend: bijna alle getallen zijn normaal, en toch is de normaliteit van 2\sqrt2, e\eu of π\pi onbewezen — de maattheorie telt zonder te tonen.

11. Door Oefening 22.10 te itereren levert één uniforme veranderlijke een rij i.i.d. uniforme vectoren UkU_k op [0,1]d\intcc01^d op (splits de cijferverzameling van elke UnU_n uit Stelling 22.6 in dd deelfamilies). Voor gL1([0,1]d)g \in L^1(\intcc01^d) zijn de veranderlijken g(Uk)g(U_k) i.i.d. en integreerbaar met gemiddelde g ⁣dλd\int g\,\dd\lambda_d (overdracht): de sterke wet geeft

1nk=1ng(Uk)nb.z.[0,1]dg ⁣dλd:\frac1n\sum_{k=1}^ng(U_k) \xrightarrow[n\to\infty]{\text{b.z.}} \int_{\intcc01^d}g\,\dd\lambda_d :

de monte-carlo-integratie convergeert bijna zeker, in elke dimensie — de grootte van de fout is de zaak van de centrale limietstelling (Hoofdstuk 23).

12. Zij Ak={Skε}j<k{Sj<ε}A_k = \{\abs{S_k} \geq \varepsilon\} \cap \bigcap_{j<k}\{\abs{S_j} < \varepsilon\}: de AkA_k zijn disjunct met vereniging A={maxknSkε}A = \{\max_{k\leq n}\abs{S_k} \geq \varepsilon\}. Dan is

ESn2k=1nE[Sn21Ak]=k=1nE[(Sk2+2Sk(SnSk)+(SnSk)2)1Ak]k=1nE[Sk21Ak],\E S_n^2 \geq \sum_{k=1}^n\E\bigl[S_n^2\mathbf 1_{A_k}\bigr] = \sum_{k=1}^n\E\Bigl[\bigl(S_k^2 + 2S_k(S_n - S_k) + (S_n - S_k)^2\bigr)\mathbf 1_{A_k}\Bigr] \geq \sum_{k=1}^n\E\bigl[S_k^2\mathbf 1_{A_k}\bigr],

want de kruisterm verdwijnt: Sk1AkS_k\mathbf 1_{A_k} is een borelfunctie van de coalitie (Z1,,Zk)(Z_1, \dots, Z_k), die onafhankelijk is van SnSkS_n - S_k, een functie van (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n) (Stelling 22.5), dus E[Sk1Ak(SnSk)]=E[Sk1Ak]E[SnSk]=0\E[S_k\mathbf 1_{A_k}(S_n - S_k)] = \E[S_k\mathbf 1_{A_k}]\,\E[S_n - S_k] = 0. Op AkA_k is Sk2ε2S_k^2 \geq \varepsilon^2, waaruit ESn2ε2kP(Ak)=ε2P(A)\E S_n^2 \geq \varepsilon^2\sum_k\P(A_k) = \varepsilon^2\P(A); en ESn2=knV(Zk)\E S_n^2 = \sum_{k\leq n}\V(Z_k) (de varianties tellen op). De beslissende stap is de ontbinding: Sk1AkS_k\mathbf 1_{A_k} is een niet-lineaire functie van het hele eerste blok, en haar onafhankelijkheid van het tweede blok is onafhankelijkheid van coalities — de paarsgewijze onafhankelijkheid van de ZiZ_i ontkoppelt enkel paren en zou dat niet rechtvaardigen.

13. Leg NN vast en pas vraag 12 toe op ZN+1,,ZN+mZ_{N+1}, \dots, Z_{N+m}:

P(maxN<kN+mSkSN>ε)1ε2j=N+1N+mV(Zj)rNε2,rN=j>NV(Zj).\P\Bigl(\max_{N < k \leq N+m}\abs{S_k - S_N} > \varepsilon\Bigr) \leq \frac1{\varepsilon^2}\sum_{j=N+1}^{N+m}\V(Z_j) \leq \frac{r_N}{\varepsilon^2}, \qquad r_N = \sum_{j>N}\V(Z_j) .

De gebeurtenissen stijgen met mm; de continuïteit van onderen geeft P(supk>NSkSN>ε)rN/ε2\P(\sup_{k>N}\abs{S_k - S_N} > \varepsilon) \leq r_N/\varepsilon^2, en rN0r_N \to 0 volgens de hypothese. Bijgevolg is voor elke pNp \in \N^* P(N{supk>NSkSN>1p})infNp2rN=0\P\bigl(\bigcap_N\{\sup_{k>N} \abs{S_k - S_N} > \frac1p\}\bigr) \leq \inf_Np^2r_N = 0: bijna zeker is er voor elke pp een NN met supk>NSkSN1p\sup_{k>N}\abs{S_k - S_N} \leq \frac1p (snijd de aftelbaar vele b.z.-gebeurtenissen over pp door), zodat SkSl2p\abs{S_k - S_l} \leq \frac2p voor alle k,l>Nk, l > N: de partiële sommen zijn b.z. cauchy, en dus b.z. convergent.

14. De veranderlijken Zn=xnεnZ_n = x_n\varepsilon_n zijn onafhankelijk (borelfuncties van onafhankelijke veranderlijken, Oefening 22.3(a)) en gecentreerd, met V(Zn)=xn2\V(Z_n) = x_n^2: vraag 13 is van toepassing zodra nxn2<\sum_nx_n^2 < \infty en geeft b.z. convergentie. In het algemeen wordt voor elke NN de convergentie van nxnεn\sum_nx_n\varepsilon_n niet beïnvloed door de waarden van ε1,,εN\varepsilon_1, \dots, \varepsilon_N: de convergentiegebeurtenis ligt in de staart-σ\sigma-algebra van de onafhankelijke rij (εn)(\varepsilon_n), zodat de nul-een-wet van Kolmogorov (Stelling 22.9) haar kans tot 00 of 11 dwingt.

15. (a) Splitsen op het niveau θEZ\theta\E Z en Cauchy–Schwarz op het bovenste stuk geeft

EZ=E[Z1ZθEZ]+E[Z1Z>θEZ]θEZ+EZ2P(Z>θEZ),\E Z = \E\bigl[Z\mathbf 1_{Z \leq \theta\E Z}\bigr] + \E\bigl[Z\mathbf 1_{Z > \theta\E Z}\bigr] \leq \theta\,\E Z + \sqrt{\E Z^2}\, \sqrt{\P(Z > \theta\E Z)} ,

dus (1θ)EZEZ2P(Z>θEZ)(1 - \theta)\E Z \leq \sqrt{\E Z^2\,\P(Z > \theta\E Z)}; kwadrateer. (b) Ontwikkel Tn4=i,j,k,lxixjxkxlE[εiεjεkεl]T_n^4 = \sum_{i,j,k,l}x_ix_jx_kx_l\, \E[\varepsilon_i\varepsilon_j\varepsilon_k\varepsilon_l]: de verwachtingswaarde is 11 wanneer de indices paren vormen (alle vier gelijk, of twee verschillende paren, dat laatste in 33 schikkingen) en 00 anders (een ongepaard teken heeft gemiddelde nul en splitst wegens de onafhankelijkheid af). Bijgevolg is

ETn4=kxk4+3ijxi2xj2=3sn42kxk43sn4.\E T_n^4 = \sum_kx_k^4 + 3\sum_{i\neq j}x_i^2x_j^2 = 3s_n^4 - 2\sum_kx_k^4 \leq 3s_n^4 .

(c) Paley–Zygmund met Z=Tn2Z = T_n^2, EZ=sn2\E Z = s_n^2 en θ=14\theta = \frac14:

P(Tn>sn2)=P(Tn2>sn24)(34)2sn43sn4=316.\P\Bigl(\abs{T_n} > \frac{s_n}2\Bigr) = \P\Bigl(T_n^2 > \frac{s_n^2}4\Bigr) \geq \Bigl(\frac34\Bigr)^2 \frac{s_n^4}{3s_n^4} = \frac3{16} .

Convergeerde de reeks met positieve kans, dan convergeerde ze b.z. (vraag 14), dus was supnTn<\sup_n\abs{T_n} < \infty b.z., en zou een MM voldoen aan P(supnTn>M)<316\P(\sup_n\abs{T_n} > M) < \frac3{16}; maar zodra sn>2Ms_n > 2M is P(Tn>M)P(Tn>sn2)316\P(\abs{T_n} > M) \geq \P(\abs{T_n} > \frac{s_n}2) \geq \frac3{16}: tegenspraak. De divergentie is dus bijna zeker, en met vraag 14 is de dichotomie volledig.

16. Hier is xn=nsx_n = n^{-s} en nn2s<\sum_nn^{-2s} < \infty precies wanneer s>12s > \frac12: volgens de vragen 14 en 15 convergeert nεnns\sum_n\frac{\varepsilon_n}{n^s} b.z. dan en slechts dan als s>12s > \frac12 (voor s12s \leq \frac12 b.z. divergentie). Voor 12<s1\frac12 < s \leq 1 is de convergentie nooit absoluut. De vergelijking is leerzaam: volmaakt alternerende tekens heffen op met sterkte nsn^{-s} voor elke s>0s > 0, terwijl typische toevallige tekens enkel op wortelsterkte opheffen — de toevalswandeling van vraag 21 groeit als n\sqrt n, en de sommatie van Abel zet precies die groei om in de convergentie van εnns\sum\varepsilon_nn^{-s} voor s>12s > \frac12.

17. (a) coshλ=kλ2k(2k)!\cosh\lambda = \sum_k\frac{\lambda^{2k}}{(2k)!} en eλ2/2=kλ2k2kk!\eu^{\lambda^2/2} = \sum_k\frac{\lambda^{2k}}{2^kk!}; en (2k)!2kk!(2k)! \geq 2^kk! geldt term voor term, want (2k)!k!=i=1k(k+i)i=1k(2i)=2kk!\frac{(2k)!}{k!} = \prod_{i=1}^k(k + i) \geq \prod_{i=1}^k(2i) = 2^kk! (elke factor voldoet aan k+i2ik + i \geq 2i voor iki \leq k), zodat zelfs (2k)!2k(k!)22kk!(2k)! \geq 2^k(k!)^2 \geq 2^kk!. (b) Merk op dat a0ba \leq 0 \leq b (ZZ is gecentreerd), en dat wegens de convexiteit van zeλzz \mapsto \eu^{\lambda z} voor z[a,b]z \in \intcc ab

eλzbzbaeλa+zabaeλb,dusEeλZbeλaaeλbba=(1p)ept+pe(1p)t=eφ(t)\eu^{\lambda z} \leq \frac{b - z}{b - a}\,\eu^{\lambda a} + \frac{z - a}{b - a}\,\eu^{\lambda b}, \qquad\text{dus}\qquad \E\,\eu^{\lambda Z} \leq \frac{b\,\eu^{\lambda a} - a\,\eu^{\lambda b}}{b - a} = (1 - p)\eu^{-pt} + p\,\eu^{(1-p)t} = \eu^{\varphi(t)}

met p=aba[0,1]p = \frac{-a}{b-a} \in \intcc01, t=λ(ba)t = \lambda(b - a) en φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t). Dan is φ(0)=0\varphi(0) = 0, verdwijnt φ(t)=p+pet1p+pet\varphi'(t) = -p + \frac{p\eu^t}{1 - p + p\eu^t} in 00, en is φ(t)=ρ(1ρ)14\varphi''(t) = \rho(1 - \rho) \leq \frac14 voor ρ=pet1p+pet[0,1]\rho = \frac{p\eu^t}{1 - p + p\eu^t} \in \intcc01: Taylor tot orde 22 geeft φ(t)t28=λ2(ba)28\varphi(t) \leq \frac{t^2}8 = \frac{\lambda^2(b-a)^2}8.

18. Voor λ>0\lambda > 0 geven Markov toegepast op de positieve veranderlijke eλ(SnESn)\eu^{\lambda(S_n - \E S_n)} (Propositie 22.3) en de productformule voor onafhankelijke veranderlijken

P(SnESnt)eλti=1nEeλ(XiEXi)exp(λt+λ28i(biai)2),\P(S_n - \E S_n \geq t) \leq \eu^{-\lambda t}\prod_{i=1}^n\E\,\eu^{\lambda(X_i - \E X_i)} \leq \exp\Bigl(-\lambda t + \frac{\lambda^2}8\sum_i(b_i - a_i)^2\Bigr),

volgens vraag 17(b) toegepast op elke gecentreerde XiEXi[aiEXi,biEXi]X_i - \E X_i \in \intcc{a_i - \E X_i}{b_i - \E X_i} (dezelfde breedte). De exponent minimaliseren in λ=4tD\lambda = \frac{4t}{D} met D=i(biai)2D = \sum_i(b_i - a_i)^2 levert 2t2D-\frac{2t^2}D. De onderste staart volgt door het resultaat op (Xi)(-X_i) toe te passen.

19. Neem t=nεt = n\varepsilon en D=n(ba)2D = n(b - a)^2:

P(Snnmε)2exp(2n2ε2n(ba)2)=2exp(2nε2(ba)2),\P\Bigl(\Bigl|\frac{S_n}n - m\Bigr| \geq \varepsilon\Bigr) \leq 2\exp\Bigl(\frac{-2n^2\varepsilon^2}{n(b-a)^2}\Bigr) = 2\exp\Bigl(\frac{-2n\varepsilon^2}{(b-a)^2}\Bigr),

wat sommeerbaar is in nn (een reeks van meetkundig type): Borel–Cantelli (Stelling 22.8) geeft dat b.z. uiteindelijk Snnm<ε\abs{\frac{S_n}n - m} < \varepsilon; doorsnijden over ε=1p\varepsilon = \frac1p levert Snnm\frac{S_n}n \to m b.z. Vergelijking: Etemadi vraagt enkel X1L1X_1 \in L^1 en paarsgewijze onafhankelijkheid, en levert geen snelheid; Hoeffding vraagt begrensdheid en volledige onafhankelijkheid, en levert bij elke eindige nn een expliciete exponentiële waarborg — de twee stellingen beantwoorden verschillende vragen over dezelfde limiet.

20. De g(Uk)g(U_k) zijn i.i.d. met waarden in [0,1]\intcc01 en gemiddelde g ⁣dλd\int g\,\dd\lambda_d (overdracht), dus geeft vraag 18 met biai=1b_i - a_i = 1 en t=nεt = n\varepsilon de tweezijdige grens 2e2nε2δ2\eu^{-2n\varepsilon^2} \leq \delta zodra e2nε22δ\eu^{2n\varepsilon^2} \geq \frac2\delta, dat wil zeggen nlog(2/δ)2ε2n \geq \frac{\log(2/\delta)}{2\varepsilon^2}. Voor ε=δ=102\varepsilon = \delta = 10^{-2}:

nlog2002104=5.29830.000226492:n \geq \frac{\log 200}{2\cdot10^{-4}} = \frac{5.2983\ldots}{0.0002} \approx 26\,492 :

ongeveer 2650026\,500 trekkingen waarborgen een nauwkeurigheid van 1%1\% met 99%99\% vertrouwen — in elke dimensie dd, voor elke meetbare integrand met waarden in [0,1]\intcc01. De sterke wet van vraag 11 beloofde convergentie zonder enige waarborg bij eindige nn; een deterministisch rooster met kk punten per as kost kdk^d evaluaties, exponentieel in dd. De concentratie is wat van monte carlo een methode maakt in plaats van een hoop.

21. Onafhankelijkheid en de productformule: EeλSn=(Eeλε1)n=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\E\,\eu^{\lambda\varepsilon_1})^n = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} volgens vraag 17(a). Markov op eλSn\eu^{\lambda S_n}:

P(Snx)eλx+nλ2/2=ex2/(2n)in het optimum λ=xn,\P(S_n \geq x) \leq \eu^{-\lambda x + n\lambda^2/2} = \eu^{-x^2/(2n)} \qquad\text{in het optimum } \lambda = \frac xn,

en de symmetrische grens voor Sn-S_n (dezelfde verdeling) verdubbelt de constante voor Sn\abs{S_n}.

22. Leg η>0\eta > 0 vast en stel xn=(1+η)2nlognx_n = (1 + \eta)\sqrt{2n\log n} voor n2n \geq 2:

P(Snxn)2exp((1+η)2logn)=2n(1+η)2,\P(\abs{S_n} \geq x_n) \leq 2\exp\bigl(-(1 + \eta)^2\log n\bigr) = \frac{2}{n^{(1+\eta)^2}},

sommeerbaar omdat (1+η)2>1(1 + \eta)^2 > 1. Borel–Cantelli: b.z. is Sn<(1+η)2nlogn\abs{S_n} < (1 + \eta)\sqrt{2n\log n} voor alle grote nn, dus lim supnSn2nlogn1+η\limsup_n\frac{\abs{S_n}}{\sqrt{2n\log n}} \leq 1 + \eta b.z.; de b.z.-gebeurtenissen voor η=1p\eta = \frac1p met pNp \in \N^* doorsnijden geeft de bewering. De wandeling van grootte nn heeft een typische amplitude n\sqrt n (haar variantie), en zelfs haar ergste uitstapjes overschrijden die schaal met hoogstens 2logn\sqrt{2\log n}.

23. Met nj=2jn_j = 2^j en x=(1+η)2njloglognjx = (1 + \eta)\sqrt{2n_j\log\log n_j} (gedefinieerd voor j2j \geq 2) geeft vraag 21

P(Snjx)exp((1+η)2loglognj)=(jlog2)(1+η)2,\P\bigl(S_{n_j} \geq x\bigr) \leq \exp\bigl(-(1 + \eta)^2\log\log n_j\bigr) = (j\log 2)^{-(1+\eta)^2},

sommeerbaar in jj omdat (1+η)2>1(1 + \eta)^2 > 1: Borel–Cantelli en η=1p\eta = \frac1p geven lim supjSnj/2njloglognj1\limsup_jS_{n_j}/\sqrt{2n_j \log\log n_j} \leq 1 b.z. Wat voor de volledige bovenste helft ontbreekt, is de brug tussen de ijkpunten: men moet aantonen dat maxnjnnj+1Sn\max_{n_j\leq n\leq n_{j+1}}S_n slechts eindig vaak (1+η)2njloglognj(1+\eta) \sqrt{2n_j\log\log n_j} overschrijdt, wat een maximaalongelijkheid met gaussische staarten vergt (de spiegelingsongelijkheid van Lévy of die van Ottaviani, hier niet bewezen). Vraag 12 is kwantitatief te zwak: ze begrenst de kans door

nj(1+η)22njloglognj=12(1+η)2log(jlog2),\frac{n_j}{(1+\eta)^2\,2n_j\log\log n_j} = \frac{1}{2(1+\eta)^2\log(j\log2)},

wat naar 00 gaat maar niet sommeerbaar is in jj: Borel–Cantelli kan niet besluiten. De onderste helft van de wet van de geïtereerde logaritme past het tweede lemma van Borel–Cantelli toe op de onafhankelijke aangroeiingen Snj+1SnjS_{n_{j+1}} - S_{n_j}, met bijpassende ondergrenzen voor staarten van gaussisch type. Beide verfijningen zijn echte kansrekening van bachelorjaar 3, één cursus verderop; wat deze opgave zonder hulp levert, is de exacte schaal van de geïtereerde logaritme langs meetkundige tijdstippen.

24. Elke p^i\hat p_i is een gemiddelde van nn i.i.d. indicatorveranderlijken met waarden in [0,1]\intcc01 en gemiddelde P(Ai)\P(A_i): Hoeffding geeft P(p^iP(Ai)>ε)2e2nε2\P(\abs{\hat p_i - \P(A_i)} > \varepsilon) \leq 2\eu^{-2n\varepsilon^2}. De somgrens vermenigvuldigt met NN. 2Ne2nε2δ2N\eu^{-2n\varepsilon^2} \leq \delta oplossen: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2}. Numeriek: ln21060.05=ln(4107)17.5\ln\frac{2\cdot10^6}{0.05} = \ln(4\cdot10^7) \approx 17.5, dus n17.5210487600n \geq \frac{17.5}{2\cdot10^{-4}} \approx 87\,600: één kans schatten tot op ±1%\pm1\% vergt ongeveer 1850018\,500 trekkingen (ln(2/δ)/2ε2\ln(2/\delta)/2\varepsilon^2), en een miljoen kansen slechts 4.7\approx 4.7 keer zoveel — de uniformiteit kost lnN\ln N, niet NN: de waarneming die de minimalisering van het empirische risico, en daarmee het machinaal leren, statistisch mogelijk maakt.

25. De veranderlijken Xn=εnnαX_n = \frac{\varepsilon_n} {n^\alpha} zijn onafhankelijk, gecentreerd en begrensd, met nV(Xn)=nn2α\sum_n\V(X_n) = \sum_nn^{-2\alpha}. Is α>12\alpha > \frac12, dan convergeert de variantiereeks, en geeft de eenreeksstelling (Deel VI) de b.z. convergentie van Xn\sum X_n. Is α12\alpha \leq \frac12, dan divergeert de variantiereeks, en geeft de omgekeerde helft (het argument van Paley–Zygmund uit Deel VI, toepasbaar omdat de termen door 11 begrensd zijn) b.z. divergentie. De absolute convergentie vraagt nα<\sum n^{-\alpha} < \infty: α>1\alpha > 1. Op (12,1]\intoc{\frac12}1 convergeert de reeks b.z. hoewel Xn=\sum\abs{X_n} = \infty zeker is: de tekens spannen samen om op te heffen, met kans één — convergentie door opheffing, onzichtbaar voor elke absolute test, en (wegens de nul-een-wet) toch met een deterministisch vonnis.