Mathematics · Livro 5 · Bachelor Year 3

Matemática universitária — Graduação 3

Matemática universitária — Graduação 3 · Bachelor Year 3

22Probabilidade: fundamentos e a lei dos grandes números

O segundo ano construiu a probabilidade em espaços enumeráveis; a teoria da medida agora retira toda restrição. Um espaço de probabilidade é um espaço de medida de massa total 11, as variáveis aleatórias são aplicações mensuráveis, a esperança é a integral de Lebesgue — e, de imediato, todo o arsenal analítico (Capítulos 9, 10 e 11) se aplica ao acaso. Este capítulo instala o dicionário, constrói sequências infinitas de variáveis aleatórias independentes (em [0,1]\intcc01, a partir de dígitos binários: o acaso está escondido dentro da medida de Lebesgue), demonstra os lemas de Borel–Cantelli e a lei zero–um de Kolmogorov, organiza os modos de convergência e demonstra a lei dos grandes números — o teorema que faz as frequências convergirem às probabilidades e torna a estatística possível. O problema de fim de semana dá a demonstração de Etemadi da lei forte em sua forma definitiva L1L^1.

22.1 O dicionário

Definição 22.1

Um espaço de probabilidade é um espaço de medida (Ω,A,P)(\Omega, \mathcal A, \P) com P(Ω)=1\P(\Omega) = 1; os elementos de A\mathcal A são eventos, e uma propriedade vale quase certamente (q.c.) se seu evento tem probabilidade 11. Uma variável aleatória é uma aplicação mensurável X ⁣:ΩRX \colon \Omega \to \R (ou em Rd\R^d: um vetor aleatório); sua lei é a medida de probabilidade imagem PX=XP\P_X = X_*\P em R\R (Exercício 11.9), determinada pela função de distribuição FX(t)=P(Xt)F_X(t) = \P(X \leq t) (Exercício 9.3). XX tem densidade ff se PX=f ⁣dλ\P_X = f\,\dd\lambda; ela é discreta se PX\P_X é uma combinação enumerável de massas de Dirac. A esperança é

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

e o teorema de transferência (Exercício 11.9) a calcula na lei: E[g(X)]=Rg ⁣dPX\E[g(X)] = \int_\R g\,\dd\P_X=g(xk)pk= \sum g(x_k)p_k no caso discreto, =g(x)f(x) ⁣dx= \int g(x)f(x)\dd x no caso com densidade: as fórmulas do segundo ano, agora teoremas de uma única teoria. A variância é V(X)=E[(XEX)2]=E[X2](EX)2\V(X) = \E[(X - \E X)^2] = \E[X^2] - (\E X)^2 para XL2X \in L^2.

Exemplo 22.2

As leis padrão e suas transformadas notáveis: Bernoulli B(p)\mathcal B(p), binomial B(n,p)\mathcal B(n, p), geométrica, Poisson P(λ)\mathcal P(\lambda) (discretas: as tabelas do segundo ano continuam válidas); uniforme em [0,1]\intcc01 (a própria medida de Lebesgue); exponencial E(λ)\mathcal E(\lambda) (densidade λeλx1x>0\lambda\eu^{-\lambda x}\mathbf 1_{x>0}); e a gaussiana N(m,σ2)\mathcal N(m, \sigma^2), de densidade 1σ2πexp((xm)22σ2)\frac1{\sigma\sqrt{2\pi}}\exp\bigl(-\frac{(x - m)^2}{2\sigma^2}\bigr) — uma densidade de probabilidade pelo Problema 10.1, de média mm e variância σ2\sigma^2 (momentos gaussianos, Exercício 11.10).

Proposição 22.3 (Markov e Chebyshev)

Para X0X \geq 0 e a>0a > 0: P(Xa)EXa\P(X \geq a) \leq \frac{\E X}{a}; para XL2X \in L^2: P(XEXa)V(X)a2\P\bigl(\abs{X - \E X} \geq a\bigr) \leq \frac{\V(X)}{a^2}.

Demonstração. Exercício 10.5(a); Chebyshev é Markov aplicado a (XEX)2(X - \E X)^2.

22.2 Independência

Definição 22.4

Sub-σ\sigma-álgebras A1,,AnA\mathcal A_1, \dots, \mathcal A_n \subseteq \mathcal A são independentes se P(A1An)=P(Ai)\P(A_1\cap\dots\cap A_n) = \prod\P(A_i) para todos AiAiA_i \in \mathcal A_i; eventos são independentes se as σ\sigma-álgebras {,Ai,Aic,Ω}\{\varnothing, A_i, A_i^c, \Omega\} o forem; variáveis aleatórias X1,,XnX_1, \dots, X_n, se as σ\sigma-álgebras σ(Xi)=Xi1(B(R))\sigma(X_i) = X_i^{-1}(\mathcal B(\R)) o forem. Uma família infinita é independente se toda subfamília finita o for.

Teorema 22.5

X1,,XnX_1, \dots, X_n são independentes se, e somente se, a lei do vetor (X1,,Xn)(X_1, \dots, X_n) é a medida produto PX1PXn\P_{X_1}\otimes\cdots\otimes\P_{X_n}. Nesse caso, para gi0g_i \geq 0 (ou tais que os produtos sejam integráveis):

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

em particular, E[XY]=EXEY\E[XY] = \E X\,\E Y e V(X1++Xn)=V(Xi)\V(X_1 + \dots + X_n) = \sum\V(X_i) para variáveis L2L^2 independentes.

Demonstração. Se as XiX_i são independentes, as duas medidas de probabilidade P(X1,,Xn)\P_{(X_1,\dots,X_n)} e PXi\bigotimes\P_{X_i} coincidem em todos os produtos B1××BnB_1\times\dots\times B_n de borelianos — um π\pi-sistema que gera B(Rn)\mathcal B(\R^n) (Proposição 11.2(b)) —, logo em toda parte (Teorema 9.7). Reciprocamente, uma lei produto fatoriza todos os eventos iXi1(Bi)\bigcap_iX_i^{-1}(B_i): independência. A fórmula da esperança é então Tonelli/Fubini (Teorema 11.5) através do teorema de transferência; E[XY]=EXEY\E[XY] = \E X\E Y é o caso gi=idg_i = \mathrm{id}, e expandir o quadrado dá a aditividade das variâncias (os termos cruzados E[(XiEXi)(XjEXj)]=0\E[(X_i - \E X_i)(X_j - \E X_j)] = 0).

Teorema 22.6 (Existência de sequências independentes)

Em ([0,1],L,λ)\bigl(\intcc01, \mathcal L, \lambda\bigr) existe uma sequência (Un)n1(U_n)_{n\geq1} de variáveis aleatórias independentes, cada uma uniforme em [0,1]\intcc01. Consequentemente, para leis prescritas (μn)(\mu_n) em R\R quaisquer, existem (Xn)(X_n) independentes com PXn=μn\P_{X_n} = \mu_n.

Demonstração. Dígitos. Para ω[0,1]\omega \in \intcc01, sejam (bk(ω))(b_k(\omega)) seus dígitos binários (ω=bk2k\omega = \sum b_k2^{-k}; escolha a expansão que não termina em uma sequência de 11 — a ambiguidade envolve apenas um conjunto enumerável, logo nulo). Cada bkb_k é variável aleatória ({bk=1}\{b_k = 1\} é uma reunião finita de intervalos diádicos), e o vetor (b1,,bm)(b_1, \dots, b_m) assume cada valor de {0,1}m\{0,1\}^m em um intervalo diádico de comprimento 2m2^{-m}: os bkb_k são Bernoulli(12)(\frac12) independentes.

Reagrupamento. Separe N\N^* em infinitos conjuntos infinitos disjuntos (In)(I_n) (por exemplo, por potências de primos ou por diagonais); seja (kjn)j(k^n_j)_j uma enumeração de InI_n e ponha

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

Cada UnU_n é uniforme: seus dígitos binários são bits independentes e honestos, de modo que P(Un[l2m,(l+1)2m))=2m\P(U_n \in [l2^{-m}, (l+1)2^{-m})) = 2^{-m} para todo intervalo diádico, e os intervalos diádicos determinam a lei (Teorema 9.7). Os UnU_n são independentes: são funções de blocos disjuntos da família independente (bk)(b_k) — formalmente, os eventos {UnDn}\{U_n \in D_n\} para DnD_n diádicos dependem de finitos dígitos de conjuntos disjuntos e se fatorizam; o argumento do π\pi-sistema eleva isso a todos os borelianos.

Leis arbitrárias. Seja Gn(u)=inf{t:Fμn(t)u}G_n(u) = \inf\{t : F_{\mu_n}(t) \geq u\} (a função quantil da função de distribuição FμnF_{\mu_n}); a equivalência-chave Gn(u)t    uFμn(t)G_n(u) \leq t \iff u \leq F_{\mu_n}(t) (continuidade à direita de FF, monotonicidade) mostra que Xn=Gn(Un)X_n = G_n(U_n) é mensurável com 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): lei μn\mu_n; a independência é herdada (funções de variáveis independentes, Exercício 22.3).

Exemplo 22.7 (O problema dos aniversários, honestamente)

Entre nn pessoas com aniversários independentes e uniformes em N=365N = 365 dias, a probabilidade de que todos os aniversários difiram é

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

por condicionamentos sucessivos (ou diretamente: os N(N1)(Nn+1)N(N-1)\cdots(N - n + 1) casos favoráveis sobre o total NnN^n, um argumento de contagem que a fórmula do produto da independência torna rigoroso). Tomando logaritmos e usando ln(1x)=x+O(x2)-\ln(1 - x) = x + O(x^2):

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

O ponto de virada pn=12p_n = \frac12 fica em n2Nln21.18Nn \approx \sqrt{2N\ln2} \approx 1.18\sqrt N: para N=365N = 365, n=23n = 23 (p23=0.4927p_{23} = 0.4927). Duas morais. Primeira: as colisões entre nn itens em NN caixas aparecem na escala nNn \sim \sqrt N, e não nNn \sim N — o escalamento do aniversário, que governa as colisões de hash e o custo N\sqrt N dos ataques de aniversário em criptografia. Segunda: o cálculo é um modelo: os (n2)\binom n2 eventos de colisão de pares não são independentes e, no entanto, a resposta se comporta como se fossem (e(n2)/N\eu^{-\binom n2/N} é exatamente a heurística de pares independentes) — uma primeira instância da aproximação de Poisson, tornada rigorosa no problema de fim de semana do Capítulo 23 (desigualdade de Le Cam).

22.3 Borel–Cantelli e a lei zero–um

Teorema 22.8 (Borel–Cantelli)

Sejam (An)(A_n) eventos e lim supAn=NnNAn\limsup A_n = \bigcap_N \bigcup_{n\geq N}A_n (“AnA_n ocorre infinitas vezes”).

  1. Se P(An)<\sum\P(A_n) < \infty, então P(lim supAn)=0\P(\limsup A_n) = 0.
  2. Se P(An)=\sum\P(A_n) = \infty e os AnA_n são independentes, então P(lim supAn)=1\P(\limsup A_n) = 1.

Demonstração. (1) é o Exercício 9.4. (2): para NMN \leq M, a independência dos complementares (Exercício 22.3) dá

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}; a série diverge). Logo P(nNAn)=1\P\bigl(\bigcup_{n\geq N}A_n\bigr) = 1 para todo NN, e a interseção decrescente em NN continua tendo probabilidade 11 (continuidade por cima, Proposição 9.6).

Teorema 22.9 (Lei zero–um de Kolmogorov)

Sejam (Xn)(X_n) independentes e T=Nσ(XN,XN+1,)\mathcal T = \bigcap_N\sigma(X_N, X_{N+1}, \dots) a σ\sigma-álgebra caudal (eventos insensíveis a qualquer número finito dos XnX_n: convergência de Xn\sum X_n, de Snn\frac{S_n}n, valores de lim sup\limsup, …). Então todo TTT \in \mathcal T tem P(T){0,1}\P(T) \in \{0, 1\}.

Demonstração. Fixe NN. As σ\sigma-álgebras σ(X1,,XN)\sigma(X_1, \dots, X_N) e σ(XN+1,)\sigma(X_{N+1}, \dots) são independentes: eventos que dependem de blocos disjuntos se fatorizam nos π\pi-sistemas geradores (cilindros iN{XiBi}\bigcap_{i\leq N}\{X_i \in B_i\} e condições finitas sobre as variáveis posteriores, respectivamente), e Dynkin (Teorema 9.4, aplicado duas vezes, um lado de cada vez) estende a fatoração. Um evento caudal TT está em σ(XN+1,)\sigma(X_{N+1}, \dots) para todo NN: TT é independente de todo σ(X1,,XN)\sigma(X_1, \dots, X_N), logo da σ\sigma-álgebra que eles geram, σ(X1,X2,)\sigma(X_1, X_2, \dots) (Dynkin mais uma vez: a reunião dos σ(X1,,XN)\sigma(X_1,\dots,X_N) é um π\pi-sistema que a gera). Mas Tσ(X1,X2,)T \in \sigma(X_1, X_2, \dots) também: TT é independente de si mesmo, 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 Modos de convergência

Definição 22.10

XnXX_n \to X quase certamente se P(XnX)=1\P(X_n \to X) = 1; em probabilidade se P(XnXε)0\P(\abs{X_n - X} \geq \varepsilon) \to 0 para todo ε>0\varepsilon > 0; em LpL^p se EXnXp0\E\abs{X_n - X}^p \to 0.

Proposição 22.11

(a) a convergência q.c. implica a convergência em probabilidade; (b) a convergência LpL^p implica a convergência em probabilidade; (c) a convergência em probabilidade implica a convergência q.c. ao longo de uma subsequência; (d) nenhuma outra implicação vale em geral.

Demonstração. (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 sob convergência q.c. (continuidade por cima; o evento limsup exclui a convergência). (b) Markov: P(XnXε)εpEXnXp\P(\abs{X_n - X} \geq \varepsilon) \leq \varepsilon^{-p}\,\E\abs{X_n - X}^p. (c) Tome nkn_k com P(XnkX2k)2k\P(\abs{X_{n_k} - X} \geq 2^{-k}) \leq 2^{-k}; Borel–Cantelli (1) faz XnkX<2k\abs{X_{n_k} - X} < 2^{-k} a partir de certo índice, q.c. (d) A máquina de escrever (Exercício 12.3) em ([0,1],λ)(\intcc01, \lambda) converge em L1L^1 e em probabilidade, mas em ponto algum; n1(0,1/n)0n\mathbf 1_{\intoo0{1/n}} \to 0 q.c. mas não em L1L^1; os detalhes e os demais contraexemplos estão no Exercício 22.6.

22.5 A lei dos grandes números

Ao longo de toda esta seção, (Xn)(X_n) são independentes de mesma lei (i.i.d.), com Sn=X1++XnS_n = X_1 + \dots + X_n.

Teorema 22.12 (Lei fraca dos grandes números)

Se X1L2X_1 \in L^2, com m=EX1m = \E X_1:

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 em probabilidade (e em L2L^2).

Demonstração. ESnn=m\E\frac{S_n}n = m e V(Snn)=nV(X1)n2\V\bigl(\frac{S_n}n\bigr) = \frac{n\V(X_1)}{n^2} (Teorema 22.5); Chebyshev.

Teorema 22.13 (Lei forte dos grandes números)

Se X1L1X_1 \in L^1, então

Snnnq.c.E[X1].\frac{S_n}{n} \xrightarrow[n\to\infty]{\text{q.c.}} \E[X_1].

Demonstramo-la aqui sob a hipótese mais forte X1L4X_1 \in L^4; o caso geral (L1L^1: a demonstração de Etemadi) é o problema de fim de semana.

Demonstração sob EX14<\E X_1^4 < \infty. Centrando (XiXimX_i \mapsto X_i - m), suponha m=0m = 0. Expanda:

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 ,

pois a independência e a centragem matam todo termo que contenha um fator isolado (E[XiXjXkXl]=E[Xi]E[]=0\E[X_iX_jX_kX_l] = \E[X_i]\E[\cdots] = 0, a menos que os índices se emparelhem: os únicos sobreviventes são os nn termos i=j=k=li=j=k=l e os 3n(n1)3n(n-1) termos com dois pares distintos). 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},

somável: Borel–Cantelli (1) dá, para cada racional ε\varepsilon, que Sn/n<ε\abs{S_n/n} < \varepsilon a partir de certo índice, q.c.; intersectando em εQ+\varepsilon \in \Q_+^* (uma quantidade enumerável de eventos de probabilidade 11): Sn/n0S_n/n \to 0 q.c.

Exemplo 22.14 (O que a lei forte compra)

(a) Frequências: para lançamentos de moeda i.i.d., a frequência observada de caras converge q.c. a pp — a justificativa empírica da própria probabilidade. (b) Monte Carlo: para gL1([0,1])g \in L^1(\intcc01) e (Un)(U_n) i.i.d. uniformes (Teorema 22.6), 1nkng(Uk)01g\frac1n\sum_{k\leq n}g(U_k) \to \int_0^1g q.c.: integrais por amostragem, em qualquer dimensão, à taxa n1/2\sim n^{-1/2} independente da dimensão, tornada precisa no Capítulo 23. (c) Números normais: quase todo número real tem, em sua expansão binária, frequência assintótica 12\frac12 de uns (aplique a lei forte às variáveis dígito do Teorema 22.6) — o teorema de Borel, uma afirmação sobre os números de todo dia demonstrada por medida; o Problema 22.1 a completa em todas as bases.

Método 22.15

A ordem de trabalho para enunciados assintóticos sobre sequências aleatórias: (1) o evento é caudal? Então sua probabilidade é 00 ou 11 (Teorema 22.9), e resta apenas decidir qual. (2) Para demonstrar enunciados q.c.: Borel–Cantelli — somabilidade das probabilidades dos eventos “ruins”, via cotas do tipo Markov/Chebyshev sobre os momentos que existirem; a independência só é necessária na direção recíproca. (3) Subsequência + sanduíche: demonstre a convergência ao longo de uma subsequência tratável e controle a oscilação intermediária por monotonicidade ou desigualdades maximais — o esqueleto da demonstração de Etemadi. (4) Para limites em distribuição, espere o Capítulo 23.

22.6 Exercícios

Exercício 22.1

(a) Seja XX com função de distribuição FF contínua e estritamente crescente. Mostre que F(X)F(X) é uniforme em [0,1]\intcc01 e que G(U)FG(U) \sim F para UU uniforme, G=F1G = F^{-1}: simulação por inversão. (b) Calcule a função de distribuição e a densidade de X2X^2 para XX uniforme em [1,1]\intcc{-1}1, e de 1λlnU-\frac1\lambda\ln U para UU uniforme em (0,1)\intoo01.

Solução

Solução de Exercício 22.1.

(a) Para 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 (a continuidade e a monotonicidade estrita fazem de FF uma bijeção sobre (0,1)\intoo01 com {F(X)u}={XF1(u)}\{F(X) \leq u\} = \{X \leq F^{-1}(u)\}): F(X)F(X) é uniforme. Reciprocamente, P(G(U)t)=P(UF(t))=F(t)\P(G(U) \leq t) = \P(U \leq F(t)) = F(t): para simular uma lei, aplique a inversa da função de distribuição a uma amostra uniforme.

(b) Y=X2Y = X^2, XX uniforme em [1,1]\intcc{-1}1: para t[0,1]t \in \intcc01, FY(t)=P(tXt)=tF_Y(t) = \P(-\sqrt t \leq X \leq \sqrt t) = \sqrt t: densidade 12t1(0,1)\frac1{2\sqrt t}\mathbf 1_{\intoo01}. E 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}: a exponencial E(λ)\mathcal E(\lambda) — a inversão em ação.

Exercício 22.2

(a) Calcule a média e a variância das leis de Poisson P(λ)\mathcal P(\lambda) e geométrica via o teorema de transferência. (b) Mostre que uma variável aleatória positiva TT com P(T>t)>0\P(T > t) > 0 para todo tt satisfaz a propriedade de falta de memória P(T>t+sT>t)=P(T>s)\P(T > t + s \mid T > t) = \P(T > s) para todos s,t0s, t \geq 0 se, e somente se, TT é exponencial. (A função de sobrevivência satisfaz a equação funcional de Cauchy; a monotonicidade substitui a continuidade.)

Solução

Solução de Exercício 22.2.

(a) Poisson: EX=k0keλλkk!=λ\E X = \sum_{k\geq0}k\,\eu^{-\lambda} \frac{\lambda^k}{k!} = \lambda, E[X(X1)]=λ2\E[X(X-1)] = \lambda^2, de modo que V=λ2+λλ2=λ\V = \lambda^2 + \lambda - \lambda^2 = \lambda. Geométrica (P(X=k)=p(1p)k1\P(X = k) = p(1-p)^{k-1}): EX=1p\E X = \frac1p, V=1pp2\V = \frac{1-p}{p^2} (derive duas vezes a série geométrica).

(b) G(t)=P(T>t)G(t) = \P(T > t) é não crescente com G(0+)G(0^+)\dots G ⁣:[0,)(0,1]G \colon \intco0\infty \to \intoc01; a ausência de memória se lê G(t+s)=G(t)G(s)G(t + s) = G(t)G(s). Então G(nt)=G(t)nG(n t) = G(t)^n e G(t/n)=G(t)1/nG(t/n) = G(t)^{1/n}: G(q)=G(1)qG(q) = G(1)^q para q0q \geq 0 racionais; escrevendo G(1)=eλG(1) = \eu^{-\lambda} ((0,1)\in \intoo01: G(1)=1G(1) = 1 forçaria G1G \equiv 1, impossível para uma variável aleatória finita; G(1)=0G(1) = 0 está excluído por hipótese) e espremendo um tt arbitrário entre racionais (monotonicidade): G(t)=eλtG(t) = \eu^{-\lambda t} — a lei exponencial. A recíproca é um cálculo.

Exercício 22.3 ★★

(a) Mostre que, se X1,,XnX_1, \dots, X_n são independentes e as fif_i são funções borelianas, então as fi(Xi)f_i(X_i) são independentes. (b) Mostre que eventos A1,,AnA_1, \dots, A_n são independentes se, e somente se, seus complementares o são, e se, e somente se, as indicadoras 1Ai\mathbf 1_{A_i} são variáveis aleatórias independentes. (c) (Duas a duas é mais fraco) Duas moedas honestas: A=A = a primeira dá cara, B=B = a segunda dá cara, C=C = as duas coincidem. Mostre que A,B,CA, B, C são independentes duas a duas, mas não são independentes.

Solução

Solução de Exercício 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 borelianas), e sub-σ\sigma-álgebras de σ\sigma-álgebras independentes são independentes (a identidade que as define vale 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}): as três afirmações asseguram a independência das mesmas σ\sigma-álgebras. (Que a fatoração sobre os AiA_i se propague aos complementares é o argumento de λ\lambda-sistemas dentro da equivalência da Definição 22.4 — ou inclusão-exclusão direta.)

(c) P(A)=P(B)=P(C)=12\P(A) = \P(B) = \P(C) = \frac12; AB=AC=BCA\cap B = A\cap C = B\cap C nos pares: cada interseção é “duas caras” ou análoga, de probabilidade 14\frac14: independentes dois a dois. Mas P(ABC)=P(HH)=1418\P(A\cap B\cap C) = \P(\text{HH}) = \frac14 \neq \frac18: não independentesCC é determinado por AA e BB.

Exercício 22.4 ★★

(a) (O macaco infinito) Uma sequência i.i.d. de teclas uniformes em um alfabeto finito contém q.c. todo texto finito infinitas vezes: demonstre-o com Borel–Cantelli (2) em blocos disjuntos. (b) (Corridas) Para bits honestos i.i.d., seja RnR_n o comprimento da corrida de uns que começa na posição nn. Mostre que, q.c., Rn(1+ε)log2nR_n \geq (1+\varepsilon)\log_2n apenas um número finito de vezes, e Rnlog2nR_n \geq \log_2 n infinitas vezes (as duas metades de Borel–Cantelli; para a segunda, passe a blocos disjuntos para ganhar independência): a corrida mais longa nos primeiros nn dígitos cresce como log2n\log_2n.

Solução

Solução de Exercício 22.4.

(a) Seja o texto TT de comprimento LL e q=aLq = a^{-L} (aa o tamanho do alfabeto). Os eventos Ek={E_k = \{as posições kL+1,,(k+1)LkL+1, \dots, (k+1)L soletram T}T\} são independentes (blocos disjuntos de letras i.i.d.), cada um de probabilidade q>0q > 0: P(Ek)=\sum\P(E_k) = \infty, e Borel–Cantelli (2) dá uma infinidade de ocorrências q.c.

(b) Superior: 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)}, somável: por Borel–Cantelli (1), q.c. só finitos nn desses. Inferior: empacote blocos disjuntos — o jj-ésimo de comprimento j=log2sj\ell_j = \lceil\log_2s_j\rceil começando em sj=i<jis_j = \sum_{i<j}\ell_i; os eventos “o bloco jj é só de uns” são independentes, de probabilidade 2j1sj1jlog2j2^{-\ell_j} \asymp \frac1{s_j} \asymp \frac1{j\log_2 j}, cuja soma diverge: Borel–Cantelli (2) dá uma infinidade de blocos só de uns, isto é, Rsjlog2sjR_{s_j} \geq \log_2 s_j uma infinidade de vezes. Juntos: o comprimento máximo de uma sequência de uns nos nn primeiros dígitos é (1+o(1))log2n(1 + o(1))\log_2n q.c.

Exercício 22.5 ★★

Sejam (Xn)(X_n) independentes. (a) Mostre que o raio de convergência de Xnzn\sum X_n z^n é uma constante q.c. (possivelmente 00 ou \infty). (b) Mostre que P(Xn converge){0,1}\P(\sum X_n \text{ converge}) \in \{0, 1\} e P(Sn/nm){0,1}\P(S_n/n \to m) \in \{0,1\}. (c) Dê um evento sobre (Xn)(X_n) que não seja caudal e verifique que a lei zero–um pode falhar para ele.

Solução

Solução de Exercício 22.5.

(a) R=(lim supXn1/n)1R = \bigl(\limsup\abs{X_n}^{1/n}\bigr)^{-1} não muda se finitos XnX_n forem modificados: para todo NN, RR é σ(XN,XN+1,)\sigma(X_N, X_{N+1}, \dots)-mensurável, isto é, mensurável na cauda. Então cada evento {Rc}\{R \leq c\} tem probabilidade 00 ou 11 (Teorema 22.9), de modo que a função de distribuição de RR assume apenas os valores 0,10, 1: ela salta num único ponto c0[0,+]c_0 \in \intcc0{+\infty}, e R=c0R = c_0 q.c.

(b) A convergência de Xn\sum X_n e a de Snn\frac{S_n}n são insensíveis à alteração de finitos termos (para a segunda: os termos modificados contribuem com O(1/n)0O(1/n) \to 0): eventos de cauda; lei zero–um.

(c) {X1>0}\{X_1 > 0\} depende de X1X_1: para sinais i.i.d. (P(X1=±1)=12\P(X_1 = \pm1) = \frac12), sua probabilidade vale 12{0,1}\frac12 \notin \{0,1\} — nenhuma contradição: não é evento de cauda.

Exercício 22.6 ★★

Em ([0,1],λ)(\intcc01, \lambda), exiba — com demonstrações — variáveis aleatórias tais que: (a) Xn0X_n \to 0 em probabilidade e em todo LpL^p, mas em ponto algum q.c.; (b) Xn0X_n \to 0 q.c. mas em nenhum LpL^p; (c) Xn0X_n \to 0 em L1L^1 mas não em L2L^2; (d) e mostre: se XnXX_n \to X em probabilidade e XnYL1\abs{X_n} \leq Y \in L^1, então XnXX_n \to X em L1L^1 (subsequências + convergência dominada + o truque da subsubsequência).

Solução

Solução de Exercício 22.6.

Trabalhe em ([0,1],λ)(\intcc01, \lambda). (a) A máquina de escrever 1In\mathbf 1_{I_n} (Exercício 12.3): Xnpp=λ(In)0\norm{X_n}_p^p = \lambda(I_n) \to 0 (todo p<p < \infty), logo também em probabilidade; em todo ω\omega os valores 00 e 11 recorrem ambos: nenhuma convergência pontual em parte alguma. (b) Xn=n1(0,1/n)0X_n = n\mathbf 1_{\intoo0{1/n}} \to 0 fora de 00, mas 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, EXn2=1\E X_n^2 = 1. (d) De qualquer subsequência extraia (convergência em probabilidade) uma subsequência ulterior convergindo q.c. (Proposição 22.11(c)); a convergência dominada dá a convergência L1L^1 ao longo dela, com o mesmo limite XX. Assim, toda subsequência da sequência numérica EXnX\E\abs{X_n - X} tem uma sub-subsequência tendendo a 00: a sequência inteira tende a 00.

Exercício 22.7 ★★

Uma pesquisa de opinião estima uma proporção desconhecida pp pela frequência empírica p^n\hat p_n de nn sorteios independentes. (a) Chebyshev: mostre que P(p^npε)14nε2\P(\abs{\hat p_n - p} \geq \varepsilon) \leq \frac1{4n\varepsilon^2} (use p(1p)14p(1-p) \leq \frac14). (b) Quantos sorteios garantem um erro 3%\leq 3\% com probabilidade 95%\geq 95\% por essa cota? (A resposta verdadeira, via o Capítulo 23, é cerca de 10701070: Chebyshev é honesto, mas grosseiro.)

Solução

Solução de Exercício 22.7.

(a) p^n=Snn\hat p_n = \frac{S_n}n com SnS_n binomial: V(p^n)=p(1p)n14n\V(\hat p_n) = \frac{p(1-p)}n \leq \frac1{4n}, e Chebyshev (Proposição 22.3) dá a cota. (b) Resolva 14n(0.03)20.05\frac1{4n(0.03)^2} \leq 0.05: n140.00090.055556n \geq \frac{1}{4\cdot0.0009\cdot0.05} \approx 5556. O teorema central do limite justificará n1070n \approx 1070 para a mesma garantia: Chebyshev paga sua generalidade com um fator 5\approx 5.

Exercício 22.8 ★★★

(Bernstein) Para fC([0,1])f \in \mathcal C(\intcc01), defina o polinômio de Bernstein 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) Reconheça Bnf(x)=E[f(Snn)]B_nf(x) = \E\bigl[f\bigl(\frac {S_n}n\bigr)\bigr] para SnS_n binomial B(n,x)\mathcal B(n, x). (b) Demonstre BnffB_nf \to f uniformemente em [0,1]\intcc01: separe em {Snnxδ}\{\abs{\frac{S_n}n - x} \leq \delta\} e seu complementar, usando a continuidade uniforme e Chebyshev com a cota uniforme V(Snn)14n\V(\frac{S_n}n) \leq \frac1{4n}. (c) Conclua: uma segunda demonstração, probabilística, do teorema de aproximação de Weierstrass (Corolário 7.16), com a taxa explícita Bnff32ωf(n1/2)\norm{B_nf - f}_\infty \leq \frac32\,\omega_f(n^{-1/2}) para o módulo de continuidade ωf\omega_f — demonstre ao menos a forma O(ωf(n1/2))O(\omega_f(n^{-1/2})).

Solução

Solução de Exercício 22.8.

(a) Se SnB(n,x)S_n \sim \mathcal B(n, x), o teorema de transferência dá 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) Seja ω=ωf\omega = \omega_f o módulo de continuidade (f(u)f(v)ω(uv)\abs{f(u) - f(v)} \leq \omega(\abs{u - v}), e ω(cδ)(1+c)ω(δ)\omega(c \delta) \leq (1 + c)\,\omega(\delta) encadeando passos). Então, para qualquer δ>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)

(se uxδ\abs{u - x} \leq \delta, claro; caso contrário, ω(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)). Tome esperanças em 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) ;

com δ=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 (continuidade uniforme no compacto): um teorema de Weierstrass probabilístico, com taxa explícita e uniforme.

Exercício 22.9 ★★★

(Colecionador de figurinhas) Figurinhas de nn tipos são sorteadas uniformemente com reposição; seja TnT_n o número de sorteios até que todos os tipos apareçam. (a) Escreva Tn=k=1nτkT_n = \sum_{k=1}^{n}\tau_k com τk\tau_k geométricas de parâmetro nk+1n\frac{n - k + 1}n, os τk\tau_k independentes, e deduza ETn=nHnnlnn\E T_n = n\,H_n \sim n\ln n (com HnH_n o número harmônico) e V(Tn)π26n2\V(T_n) \leq \frac{\pi^2}6n^2. (b) Chebyshev: Tnnlnn1\frac{T_n}{n\ln n} \to 1 em probabilidade. (c) Refine com Borel–Cantelli: mostre diretamente que P(Tn>βnlnn)n1β\P(T_n > \beta n\ln n) \leq n^{1 - \beta} para β>1\beta > 1 (cota da união sobre o evento de algum tipo faltar após βnlnn\beta n\ln n sorteios, usando 1xex1 - x \leq \eu^{-x}) e deduza que, ao longo de n=2mn = 2^m, q.c. TnβnlnnT_n \leq \beta n\ln n a partir de certo índice, para todo β>2\beta > 2.

Solução

Solução de Exercício 22.9.

(a) Depois de k1k - 1 tipos coletados, cada sorteio é novo com probabilidade pk=nk+1np_k = \frac{n-k+1}n: τk\tau_k é geométrica (pk)(p_k), e os τk\tau_k são independentes (os sorteios o são). Somas: ETn=knnk+1=nHnnlnn\E T_n = \sum_k\frac n{n-k+1} = nH_n \sim n\ln n; 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, e nHnnlnn1\frac{nH_n}{n\ln n} \to 1: Tnnlnn1\frac{T_n}{n\ln n} \to 1 em probabilidade.

(c) Cota da união: Tn>tT_n > t significa que algum tipo não foi visto após t\lceil t\rceil sorteios, de modo que P(Tn>t)n(11n)tnet/n\P(T_n > t) \leq n(1 - \frac1n)^{t} \leq n\,\eu^{-t/n}; em t=βnlnnt = \beta n\ln n: n1β\leq n^{1 - \beta}. Para β>1\beta > 1, m2m(1β)<\sum_m 2^{m(1-\beta)} < \infty: Borel–Cantelli dá, ao longo de n=2mn = 2^m, q.c. TnβnlnnT_n \leq \beta n\ln n a partir de certa ordem — em particular, para todo β>2\beta > 2 como enunciado (qualquer β>1\beta > 1 serve ao longo da subsequência).

Exercício 22.10 ★★

Usando a construção por dígitos (Teorema 22.6): (a) verifique por cálculo direto que U=b2k2kU = \sum b_{2k}2^{-k} (os dígitos de índice par de um ω\omega uniforme) é uniforme e independente de V=b2k12kV = \sum b_{2k-1}2^{-k}; (b) deduza uma bijeção mensurável a menos de conjuntos nulos entre [0,1]\intcc01 e [0,1]2\intcc01^2 que preserva a medida, e comente: um único número aleatório uniforme contém dois (e enumeráveis) independentes — compare com a curva de Peano (Problema 6.1), que obteve sobrejetividade, mas não preservação de medida nem injetividade.

Solução

Solução de Exercício 22.10.

(a) Os dígitos de índice par (b2k)k(b_{2k})_k são bits honestos i.i.d. (uma subfamília da família independente de dígitos), de modo que U=kb2k2kU = \sum_kb_{2k}2^{-k} dá a todo intervalo diádico sua probabilidade correta (como no Teorema 22.6): uniforme; do mesmo modo VV; e (U,V)(U, V) dependem de blocos de dígitos disjuntos: independentes (fatoração em retângulos diádicos, depois Dynkin).

(b) Φ(ω)=(U(ω),V(ω))\Phi(\omega) = (U(\omega), V(\omega)) é mensurável com Φλ=λλ=λ2\Phi_*\lambda = \lambda\otimes\lambda = \lambda_2 (coincidência em retângulos diádicos + unicidade). Intercalar dígitos define uma inversa definida fora do conjunto (nulo) dos racionais diádicos em qualquer dos fatores: uma bijeção que preserva a medida entre subconjuntos de medida total de [0,1]\intcc01 e [0,1]2\intcc01^2. Contraste com Peano (Problema 6.1): a continuidade forçava a sobrejetividade sem a injetividade; trocar a continuidade pela mera mensurabilidade compra um isomorfismo de medida — a dimensão é invisível à teoria da medida e visível à topologia.

Exercício 22.11 ★★

(Recordes) Sejam (Xn)n1(X_n)_{n\geq1} i.i.d. com função de distribuição contínua, e diga que ocorre um recorde no instante nn se Xn>max(X1,,Xn1)X_n > \max(X_1, \dots, X_{n-1}) (o instante 11 é recorde). Seja RnR_n a indicadora de recorde. (a) Mostre que P(Rn=1)=1n\P(R_n = 1) = \frac1n (por simetria, cada uma das n!n! ordenações de X1,,XnX_1, \dots, X_n é igualmente provável, e os empates têm probabilidade 00). (b) Mostre que os RnR_n são independentes (conte as ordenações compatíveis com posições de recorde prescritas, ou argumente que a ordem relativa de X1,,Xn1X_1, \dots, X_{n-1} é independente do posto de XnX_n entre eles). (c) Deduza de Borel–Cantelli (Teorema 22.8, as duas metades) que ocorrem infinitos recordes q.c., mas que recordes em instantes consecutivos n,n+1n, n+1 ocorrem infinitas vezes com probabilidade — decida qual! — e calcule nP(Rn=1,Rn+1=1)\sum_n\P(R_n = 1, R_{n+1} = 1).

Solução

Solução de Exercício 22.11.

(a) A continuidade da distribuição torna nulos os eventos de empate (como nos argumentos de estatísticas de ordem do capítulo), e as n!n! ordenações relativas de (X1,,Xn)(X_1, \dots, X_n) são permutáveis, logo equiprováveis. Rn=1R_n = 1 significa que o máximo ocupa a última posição: probabilidade (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n.

(b) Fixe nn e condicione à ordenação relativa de X1,,Xn1X_1, \dots, X_{n-1}: inserir XnX_n nas nn posições de posto possíveis é uniforme e independente dessa ordenação (permutabilidade da nn-upla). Logo RnR_n (o evento “XnX_n ocupa a primeira posição”) é independente de todo o histórico de recordes (R1,,Rn1)(R_1, \dots, R_{n-1}), que é função da ordenação relativa das n1n - 1 primeiras variáveis. A indução dá a independência plena, com P(Rn=1)=1n\P(R_n = 1) = \frac1n.

(c) P(Rn=1)=1n=\sum\P(R_n = 1) = \sum\frac1n = \infty com independência: a segunda metade de Borel–Cantelli dá recordes uma infinidade de vezes q.c. (os recordes nunca cessam — mas rareiam logaritmicamente: E[#recordesn]=Hnlnn\E[\#\text{recordes} \leq n] = H_n \approx \ln n). Recordes consecutivos: P(Rn=Rn+1=1)=1n(n+1)\P(R_n = R_{n+1} = 1) = \frac1{n(n+1)} (independência), e

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

a primeira metade de Borel–Cantelli se aplica — só finitos pares de recordes consecutivos ocorrem, q.c.

Exercício 22.12 ★★

(A corrida de caras mais longa) Lance uma moeda honesta infinitas vezes e seja LnL_n o comprimento da maior corrida de caras consecutivas entre os primeiros nn lançamentos. (a) Mostre que, para todo ε>0\varepsilon > 0, q.c. Ln(1+ε)log2nL_n \leq (1 + \varepsilon)\log_2n a partir de certo índice (a probabilidade de alguma corrida de comprimento \ell começar entre os primeiros nn lançamentos é no máximo n2n2^{-\ell}; Borel–Cantelli ao longo de n=2kn = 2^k). (b) Mostre que q.c. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n a partir de certo índice (reparta os primeiros nn lançamentos em n/\lfloor n/\ell\rfloor blocos disjuntos de comprimento =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil; os blocos são independentes, cada um todo de caras com probabilidade 22^{-\ell}, e a probabilidade de nenhum ser todo de caras é no máximo exp(n2/)\exp(-n2^{-\ell}/\ell); some ao longo de n=2kn = 2^k, de novo). (c) Conclua que Lnlog2n1\frac{L_n}{\log_2n} \to 1 q.c.: em um milhão de lançamentos honestos deve-se esperar uma corrida de cerca de 2020 caras — e um conjunto de dados sem ela provavelmente foi fabricado.

Solução

Solução de Exercício 22.12.

(a) Uma sequência de comprimento \ell começando na posição ini \leq n tem probabilidade 22^{-\ell}; cota da união: P(Ln)n2\P(L_n \geq \ell) \leq n2^{-\ell}. Com n=(1+ε)log2n\ell_n = (1 + \varepsilon)\log_2n: P(Lnn)nε\P(L_n \geq \ell_n) \leq n^{-\varepsilon}. Ao longo de n=2kn = 2^k: k2kε<\sum_k2^{-k\varepsilon} < \infty, de modo que q.c. L2k<(1+ε)kL_{2^k} < (1+\varepsilon)k a partir de certa ordem (Borel–Cantelli); para nn geral, tome 2k1<n2k2^{k-1} < n \leq 2^k e use a monotonicidade de LnL_n mais 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, e o fator extra é absorvido aumentando ligeiramente ε\varepsilon.

(b) Com =(1ε)log2n\ell = \lceil(1 - \varepsilon)\log_2n\rceil e m=n/m = \lfloor n/\ell\rfloor blocos disjuntos: os blocos são independentes, cada um só de caras com probabilidade 2n(1ε)/22^{-\ell} \geq n^{-(1-\varepsilon)}/2, de modo que

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)

para uma constante c>0c > 0 e nn grande. Essas probabilidades são somáveis ao longo de n=2kn = 2^k (de fato, ao longo de todos os nn): Borel–Cantelli dá q.c. Ln(1ε)log2nL_n \geq (1 - \varepsilon)\log_2n a partir de certa ordem (a monotonicidade preenche entre os 2k2^k como em (a), inofensivamente).

(c) Ambas as cotas ao longo de uma sequência ε=1j\varepsilon = \frac1j, intersectando uma infinidade enumerável de eventos de medida total: Lnlog2n1\frac{L_n}{\log_2n} \to 1 q.c. Para n=106n = 10^6: log2n19.9\log_2n \approx 19.9 — uma sequência de 20\approx 20 caras não é uma anomalia suspeita, e sim uma certeza matemática, e sua ausência é indício de um humano forjando “aleatoriedade” (humanos raramente ousam escrever mais de 55 ou 66 caras seguidas).

22.7 Problema: a demonstração de Etemadi da lei forte

Problema 22.1

Problema de fim de semana — a lei forte dos grandes números para variáveis i.i.d. integráveis

A lei forte de Kolmogorov — SnnEX1\frac{S_n}n \to \E X_1 q.c. para XnL1X_n \in L^1 i.i.d. — teve por muito tempo apenas demonstrações intrincadas; em 1981, N. Etemadi encontrou uma de economia notável, que não usa nada além deste capítulo (e chega a enfraquecer a independência para independência duas a duas). Vamos segui-la. Sejam (Xn)(X_n) independentes duas a duas, de mesma lei, integráveis; m=EX1m = \E X_1, Sn=X1++XnS_n = X_1 + \dots + X_n.

Parte I — Reduções.

  1. Mostre que basta tratar Xn0X_n \geq 0 (separe Xn=Xn+XnX_n = X_n^+ - X_n^-: verifique que as duas metades são de novo i.i.d. integráveis e independentes duas a duas). Suponha daqui em diante Xn0X_n \geq 0.
  2. (Truncamento) Sejam Yn=Xn1XnnY_n = X_n\,\mathbf 1_{X_n \leq n} e Sn=Y1++YnS_n^* = Y_1 + \dots + Y_n. Mostre que

    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

    (Exercício 11.3) e deduza, via Borel–Cantelli, que SnSnn0\frac{S_n - S_n^*}{n} \to 0 q.c.: basta demonstrar Snnm\frac{S^*_n}n \to m q.c.

  3. Mostre que EYn=E[X11X1n]m\E Y_n = \E\bigl[X_1\mathbf 1_{X_1\leq n}\bigr] \to m (convergência monótona), logo 1nknEYkm\frac1n\sum_{k\leq n}\E Y_k \to m (Cesàro): basta demonstrar SnESnn0\frac{S_n^* - \E S_n^*}{n} \to 0 q.c.

Parte II — A estimativa de variância.

  1. Mostre que

    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]

    e, usando o bolo de camadas (Proposição 11.8), a cota-chave

    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

    (troque a soma com a esperança — Tonelli para séries — e limite nx1n22max(x,1)\sum_{n \geq x}\frac1{n^2} \leq \frac2{\max(x,1)} para a estimativa interna x2nxn22xx^2\sum_{n\geq x}n^{-2} \leq 2x).

Parte III — Convergência ao longo de subsequências geométricas. Fixe α>1\alpha > 1 e ponha kj=αjk_j = \lfloor\alpha^j\rfloor.

  1. Usando a independência duas a duas (as variâncias se somam, Teorema 22.5 — verifique que a aditividade das variâncias só precisa da independência duas a duas) e Chebyshev, mostre que, para todo ε>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. Mostre que j:kjnkj2Cαn2\sum_{j : k_j \geq n}k_j^{-2} \leq \frac{C_\alpha}{n^2} (série geométrica; cuidado com a parte inteira: kjαj2k_j \geq \frac{\alpha^j}2 para o cuidado do tipo αj2\alpha^j \geq 2) e conclua, com a questão 4 e Borel–Cantelli:

    SkjESkjkjjq.c.0,logoSkjkjm q.c.\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \xrightarrow[j\to\infty]{\text{q.c.}} 0, \qquad\text{logo}\qquad \frac{S^*_{k_j}}{k_j} \to m \ \text{q.c.}

Parte IV — Sanduíche e conclusão.

  1. Para kjnkj+1k_j \leq n \leq k_{j+1}, use a monotonicidade de SnS^*_n (as parcelas são não negativas!) para mostrar

    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}},

    e deduza, q.c.:

    mαlim infSnnlim supSnnαm.\frac m\alpha \leq \liminf\frac{S^*_n}n \leq \limsup\frac{S^*_n}n \leq \alpha\,m .
  2. Faça α1\alpha \downarrow 1 ao longo de uma sequência e conclua Snnm\frac{S_n^*}n \to m q.c., logo (Parte I) a lei forte dos grandes números:

     Snnnq.c.E[X1]. \boxed{\ \frac{S_n}{n} \xrightarrow[n\to\infty]{\text{q.c.}} \E[X_1].\ }
  3. Onde exatamente a independência duas a duas (em vez da independência plena) foi suficiente? Liste os três pontos em que hipóteses do tipo independência foram invocadas.

Parte V — Dividendos.

  1. (Números normais de Borel) Mostre que λ\lambda-quase todo x[0,1]x \in \intcc01 é normal em toda base b2b \geq 2: cada dígito 0,,b10, \dots, b-1 aparece com frequência assintótica 1b\frac1b (fixe bb e um dígito, aplique a lei forte às variáveis indicadoras — justifique que os dígitos na base bb de uma variável uniforme são i.i.d. uniformes em {0,,b1}\{0,\dots,b-1\}, como no Teorema 22.6 — e intersecte os enumeráveis eventos de probabilidade um). Exiba um número explicitamente não normal e reflita: o teorema afirma a normalidade de quase todos os números e, no entanto, demonstrar a normalidade de 2\sqrt2 ou de π\pi continua em aberto.
  2. (Monte Carlo, com garantia) Justifique por completo o método do Exemplo 22.14(b) para gL1([0,1]d)g \in L^1(\intcc01^d): construa a amostra uniforme i.i.d. em [0,1]d\intcc01^d a partir do Teorema 22.6 e do Exercício 22.10, e enuncie o que a lei forte entrega.

Parte VI — O que a independência plena compra: desigualdades maximais e séries aleatórias. Etemadi gasta apenas a independência duas a duas; as partes restantes exploram a versão plena (mútua). Sejam (Zn)(Z_n) variáveis independentes centradas de L2L^2 e Sk=Z1++ZkS_k = Z_1 + \dots + Z_k (notação nova, sem relação com os XnX_n acima).

  1. (Desigualdade maximal de Kolmogorov) Para ε>0\varepsilon > 0, demonstre

    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) :

    o preço de Chebyshev compra o máximo (reparta o evento segundo o primeiro índice kk com Skε\abs{S_k} \geq \varepsilon; nessa peça, escreva Sn2Sk2+2Sk(SnSk)S_n^2 \geq S_k^2 + 2S_k(S_n - S_k) e use a independência das coalizões (Z1,,Zk)(Z_1, \dots, Z_k) e (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n), Teorema 22.5). Aponte o passo em que a independência duas a duas já não bastaria.

  2. (Teorema de uma série de Khinchin–Kolmogorov) Deduza: se nV(Zn)<\sum_n\V(Z_n) < \infty, então nZn\sum_nZ_n converge quase certamente (mostre que, q.c., as somas parciais formam uma sequência de Cauchy: faça mm \to \infty na desigualdade maximal aplicada a ZN+1,,ZN+mZ_{N+1}, \dots, Z_{N+m} e, depois, faça NN \to \infty).
  3. (Séries de Rademacher) Sejam (εn)(\varepsilon_n) sinais i.i.d., P(εn=±1)=12\P(\varepsilon_n = \pm1) = \frac12 (Teorema 22.6), e sejam (xn)(x_n) números reais. Mostre que nxnεn\sum_nx_n\varepsilon_n converge q.c. assim que nxn2<\sum_nx_n^2 < \infty; mostre também que, quaisquer que sejam os (xn)(x_n), a probabilidade de nxnεn\sum_nx_n\varepsilon_n convergir é 00 ou 11 (Teorema 22.9).
  4. A recíproca, de maneira elementar. Ponha Tn=knxkεkT_n = \sum_{k\leq n}x_k\varepsilon_k e sn2=knxk2s_n^2 = \sum_{k\leq n}x_k^2, e suponha sns_n \to \infty. (a) Demonstre a desigualdade de Paley–Zygmund: para Z0Z \geq 0 com EZ2<\E Z^2 < \infty e 0<θ<10 < \theta < 1,

    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}

    (separe EZ\E Z no nível θEZ\theta\E Z e aplique Cauchy–Schwarz à peça superior). (b) Mostre que ETn43sn4\E T_n^4 \leq 3s_n^4. (c) Deduza P(Tn>sn2)316\P\bigl(\abs{T_n} > \frac{s_n}2\bigr) \geq \frac3{16} e conclua que nxnεn\sum_nx_n\varepsilon_n diverge q.c.; donde a dicotomia

    nxnεn converge q.c.    nxn2<.\sum_nx_n\varepsilon_n\ \text{converge q.c.} \iff \sum_nx_n^2 < \infty .
  5. (Série harmônica aleatória) Conclua que nεnns\sum_n\frac{\varepsilon_n}{n^s} converge q.c. se, e somente se, s>12s > \frac12. Para 12<s1\frac12 < s \leq 1, a série converge q.c. enquanto nns=\sum_nn^{-s} = \infty: sinais aleatórios produzem cancelamento de intensidade raiz quadrada — compare com a série alternada n(1)nns\sum_n\frac{(-1)^n}{n^s}, que converge para todo s>0s > 0.

Parte VII — Concentração: a desigualdade de Hoeffding. A lei forte diz que Snnm\frac{S_n}n \to m; as desigualdades de concentração dizem quão improvável é um desvio em cada nn fixado.

  1. (Lema de Hoeffding) (a) Mostre que coshλeλ2/2\cosh\lambda \leq \eu^{\lambda^2/2} para todo λR\lambda \in \R, comparando as duas séries termo a termo. (b) Seja ZZ centrada com aZba \leq Z \leq b, a<ba < b. Mostre que

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

    (limite eλz\eu^{\lambda z} em [a,b]\intcc ab por sua corda, tome esperanças e estude φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t) com p=abap = \frac{-a}{b-a} e t=λ(ba)t = \lambda(b - a): mostre que φ(0)=φ(0)=0\varphi(0) = \varphi'(0) = 0 e φ14\varphi'' \leq \frac14).

  2. (Desigualdade de Hoeffding) Sejam X1,,XnX_1, \dots, X_n independentes com aiXibia_i \leq X_i \leq b_i e Sn=X1++XnS_n = X_1 + \dots + X_n. Demonstre, para 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),

    e a mesma cota para a cauda inferior (Chebyshev exponencial: limite Eeλ(SnESn)\E\,\eu^{\lambda(S_n - \E S_n)} usando a independência e a questão 17 e, depois, otimize em λ>0\lambda > 0).

  3. (A lei forte, caso limitado, com taxa) Sejam as XiX_i i.i.d. com valores em [a,b]\intcc ab e m=EX1m = \E X_1. Mostre que

    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)

    e recupere Snnm\frac{S_n}n \to m q.c. por Borel–Cantelli: uma segunda demonstração da lei forte para variáveis limitadas — sem truncamento, com taxa exponencial em cada nn finito, mas com parcelas limitadas e independência plena. Compare as hipóteses com as de Etemadi.

  4. (Monte Carlo, com garantia em nn fixado) Sejam g ⁣:[0,1]d[0,1]g \colon \intcc01^d \to \intcc01 mensurável e (Uk)(U_k) a amostra uniforme i.i.d. da questão 11. Dado ε,δ>0\varepsilon, \delta > 0, mostre que

    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,

    e avalie o limiar para ε=δ=102\varepsilon = \delta = 10^{-2}. A cota não envolve dd: compare com a questão 11 e com as grades determinísticas.

Parte VIII — Qual é o tamanho de um passeio aleatório? Rumo ao logaritmo iterado. Seja Sn=ε1++εnS_n = \varepsilon_1 + \dots + \varepsilon_n o passeio aleatório simples construído a partir de sinais honestos i.i.d.

  1. (Caudas sub-gaussianas) Mostre que EeλSn=(coshλ)nenλ2/2\E\,\eu^{\lambda S_n} = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} e deduza, para x>0x > 0,

    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. Deduza, via Borel–Cantelli,

    lim supnSn2nlogn1q.c.\limsup_{n\to\infty}\frac{\abs{S_n}} {\sqrt{2n\log n}} \leq 1 \quad\text{q.c.}

    (para η>0\eta > 0, some as cotas de cauda em x=(1+η)2nlognx = (1 + \eta)\sqrt{2n\log n} e, depois, intersecte em η=1p\eta = \frac1p). Em particular, o passeio vive na escala do teorema central do limite, n\sqrt n, a menos de um fator logarítmico — muito abaixo da cota grosseira Snn\abs{S_n} \leq n.

  3. Ao longo da subsequência de dobramento nj=2jn_j = 2^j, mostre que

    lim supjSnj2njloglognj1q.c.,\limsup_{j\to\infty}\frac{S_{n_j}} {\sqrt{2n_j\log\log n_j}} \leq 1 \quad\text{q.c.},

    e reflita: a lei do logaritmo iterado (Khinchin; Hartman–Wintner para parcelas L2L^2 centradas gerais) afirma que

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

    Explique com precisão o que separa a estimativa em subsequência que acabamos de demonstrar da metade superior desse enunciado (é preciso controlar maxnjnnj+1Sn\max_{n_j \leq n \leq n_{j+1}}S_n dentro de cada bloco, o que exige uma desigualdade maximal em escala exponencial) e verifique quantitativamente que a desigualdade da questão 12 é fraca demais para esse fim. A metade inferior repousa sobre o segundo lema de Borel–Cantelli aplicado a blocos independentes; as duas metades são material honesto de terceiro ano para um curso dedicado de probabilidade.

  4. (Desvio uniforme sobre uma classe finita) Sejam A1,,ANA_1, \dots, A_N eventos de um experimento repetível, e estime cada probabilidade por sua frequência empírica p^i\hat p_i em nn repetições i.i.d. Combinando a desigualdade de Hoeffding com uma cota da união, mostre que

    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},

    e deduza a regra de tamanho de amostra: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2} garante que as NN estimativas sejam simultaneamente ε\varepsilon-precisas com probabilidade 1δ\geq 1 - \delta. Calcule nn para N=106N = 10^6, ε=0.01\varepsilon = 0.01, δ=0.05\delta = 0.05: o preço logarítmico da uniformidade.

  5. (A janela harmônica aleatória) Combinando as duas metades da teoria de séries aleatórias, mostre que, para sinais i.i.d. (εn)(\varepsilon_n), a série nεnnα\sum_n\frac{\varepsilon_n}{n^\alpha} converge q.c. se α>12\alpha > \frac12 e diverge q.c. se α12\alpha \leq \frac12; contraste com a convergência absoluta (que exige α>1\alpha > 1): na janela α(12,1]\alpha \in \intoc{\frac12}1, a convergência é um fenômeno genuinamente probabilístico — cancelamento, não tamanho.
Solução

Solução de Problema 22.1.

1. Xn±X_n^{\pm} são funções borelianas de XnX_n: elas permanecem independentes dois a dois (Exercício 22.3(a)) e identicamente distribuídas, integráveis, com EX1=EX1+EX1\E X_1 = \E X_1^+ - \E X_1^-. Se o teorema vale para variáveis não negativas, aplique-o às duas metades e subtraia: 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 q.c.

2. P(XnYn)=P(Xn>n)=P(X1>n)\P(X_n \neq Y_n) = \P(X_n > n) = \P(X_1 > n) (leis idênticas), e 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 (Exercício 11.3(a)). Borel–Cantelli (1): q.c. Xn=YnX_n = Y_n para todo nn grande, de modo que SnSnS_n - S_n^* é constante em nn a partir de certa ordem: SnSnn0\frac{S_n - S_n^*}n \to 0 q.c., e as duas somas normalizadas partilham o comportamento assintótico.

3. X11X1nX1X_1\mathbf 1_{X_1 \leq n} \nearrow X_1: o TCM dá EYnm\E Y_n \to m; as médias de Cesàro de uma sequência convergente convergem ao mesmo limite: ESnn=1nknEYkm\frac{\E S_n^*}n = \frac1n\sum_{k\leq n}\E Y_k \to m. Logo basta demonstrar SnESnn0\frac{S^*_n - \E S^*_n}{n} \to 0 q.c.

4. V(Yn)EYn2=E[X121X1n]\V(Y_n) \leq \E Y_n^2 = \E[X_1^2\mathbf 1_{X_1\leq n}]. Por Tonelli para séries,

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,

usando nxn24x\sum_{n\geq x}n^{-2} \leq \frac4x para x1x \geq 1 (para x2x \geq 2: 1x12x\leq \frac1{x-1} \leq \frac2x; para 1x<21 \leq x < 2: π264x\leq \frac{\pi^2}6 \leq \frac4x, pois 4x>2\frac4x > 2), e X12/max(X1,1)X1X_1^2/\max(X_1, 1) \leq X_1 nos dois casos X11X_1 \gtrless 1.

5. A independência dois a dois dá E[(YiEYi)(YjEYj)]=0\E[(Y_i - \E Y_i)(Y_j - \E Y_j)] = 0 para iji \neq j (a fórmula do produto para duas variáveis), de modo que as variâncias se somam: V(Sk)=nkV(Yn)\V(S^*_k) = \sum_{n\leq k}\V(Y_n). Chebyshev em cada kjk_j e somando:

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 para a série dupla não negativa).

6. kj=αjαj2k_j = \lfloor\alpha^j\rfloor \geq \frac{\alpha^j}2 (válido assim que αj1\alpha^j \geq 1, isto é, para todo j0j \geq 0: xx2\lfloor x\rfloor \geq \frac x2 para x1x \geq 1). Logo

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},

(série geométrica a partir do primeiro jj com αjn\alpha^j \geq n). Combinando com as questões 4–5, a soma dupla é finita; Borel–Cantelli (1), aplicado para cada ε\varepsilon racional e intersectado, dá SkjESkjkj0\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \to 0 q.c. e, com a questão 3: Skjkjm\frac{S^*_{k_j}}{k_j} \to m q.c.

7. Yn0Y_n \geq 0 torna nSnn \mapsto S^*_n não decrescente: para kjnkj+1k_j \leq n \leq k_{j+1},

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

que é o sanduíche em destaque após inserir kjkj+1\frac{k_j}{k_{j+1}} e kj+1kj\frac{k_{j+1}}{k_j}. Como kj+1kjα\frac{k_{j+1}}{k_j} \to \alpha, a questão 6 dá q.c.

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. Aplique a questão 7 para α=1+1p\alpha = 1 + \frac1p, pNp \in \N^*: uma infinidade enumerável de eventos q.c.; em sua interseção, fazendo pp \to \infty: limSnn=m\lim\frac{S^*_n}n = m q.c. Com as questões 1–3, SnnEX1\frac{S_n}n \to \E X_1 q.c.: a lei forte dos grandes números, sob independência dois a dois.

9. As hipóteses do tipo independência apareceram três vezes: (i) a aditividade das variâncias (questão 5) — a independência dois a dois basta; (ii) a distribuição idêntica, nas somas truncadas (questão 2) e no cálculo da média (questão 3) — nenhuma independência; (iii) Borel–Cantelli (1) (questões 2 e 6) — válido sem independência alguma. A independência mútua plena nunca foi invocada: a observação de Etemadi.

10. Fixe uma base bb e um dígito rr. Os dígitos na base bb (dk)(d_k) de um ω\omega uniforme são i.i.d. uniformes em {0,,b1}\{0, \dots, b-1\} (cada valor do vetor de dígitos ocupa um intervalo de comprimento bmb^{-m}: o argumento do Teorema 22.6 literalmente). A lei forte aplicada às variáveis i.i.d. limitadas 1dk=r\mathbf 1_{d_k = r} dá: q.c., a frequência do dígito rr tende a 1b\frac1b. Intersectando sobre a infinidade enumerável de pares (b,r)(b, r): quase todo número é simplesmente normal em toda base. Um número não normal explícito: x=0.1001001002x = 0.100100100\ldots_2 (frequência de uns 1312\frac13 \neq \frac12). O contraste é humilhante: quase todos os números são normais e, no entanto, para 2\sqrt2, e\eu ou π\pi a normalidade permanece indemonstrada — a teoria da medida conta sem exibir.

11. Pelo Exercício 22.10 iterado, uma única variável uniforme fornece uma sequência de vetores uniformes i.i.d. UkU_k em [0,1]d\intcc01^d (reparta o conjunto de dígitos de cada UnU_n do Teorema 22.6 em dd subfamílias). Para gL1([0,1]d)g \in L^1(\intcc01^d), as variáveis g(Uk)g(U_k) são i.i.d. integráveis de média g ⁣dλd\int g\,\dd\lambda_d (transferência): a lei forte dá

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

A integração de Monte Carlo converge quase certamente, em toda dimensão — o tamanho do erro é assunto do teorema central do limite (Capítulo 23).

12. Seja Ak={Skε}j<k{Sj<ε}A_k = \{\abs{S_k} \geq \varepsilon\} \cap \bigcap_{j<k}\{\abs{S_j} < \varepsilon\}: os AkA_k são disjuntos, de união A={maxknSkε}A = \{\max_{k\leq n}\abs{S_k} \geq \varepsilon\}. Então

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],

porque o termo cruzado se anula: Sk1AkS_k\mathbf 1_{A_k} é função boreliana da coalizão (Z1,,Zk)(Z_1, \dots, Z_k), que é independente de SnSkS_n - S_k, função de (Zk+1,,Zn)(Z_{k+1}, \dots, Z_n) (Teorema 22.5), de modo que 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. Em AkA_k, Sk2ε2S_k^2 \geq \varepsilon^2, donde ESn2ε2kP(Ak)=ε2P(A)\E S_n^2 \geq \varepsilon^2\sum_k\P(A_k) = \varepsilon^2\P(A); e ESn2=knV(Zk)\E S_n^2 = \sum_{k\leq n}\V(Z_k) (as variâncias se somam). O passo decisivo é a fatoração: Sk1AkS_k\mathbf 1_{A_k} é função não linear de todo o primeiro bloco, e sua independência do segundo bloco é a independência por coalizões — a independência dois a dois dos ZiZ_i apenas descorrelaciona pares e não a justificaria.

13. Fixe NN e aplique a questão 12 a 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) .

Os eventos crescem com mm; a continuidade por baixo dá P(supk>NSkSN>ε)rN/ε2\P(\sup_{k>N}\abs{S_k - S_N} > \varepsilon) \leq r_N/\varepsilon^2, e rN0r_N \to 0 por hipótese. Logo, para cada 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: quase certamente, para todo pp existe NN com supk>NSkSN1p\sup_{k>N}\abs{S_k - S_N} \leq \frac1p (intersecte a infinidade enumerável de eventos q.c. sobre pp), de modo que SkSl2p\abs{S_k - S_l} \leq \frac2p para todo k,l>Nk, l > N: as somas parciais são q.c. de Cauchy, logo q.c. convergentes.

14. As variáveis Zn=xnεnZ_n = x_n\varepsilon_n são independentes (funções borelianas de variáveis independentes, Exercício 22.3(a)), centradas, com V(Zn)=xn2\V(Z_n) = x_n^2: a questão 13 se aplica quando nxn2<\sum_nx_n^2 < \infty e dá a convergência q.c. Em geral, para cada NN a convergência de nxnεn\sum_nx_n\varepsilon_n não é afetada pelos valores de ε1,,εN\varepsilon_1, \dots, \varepsilon_N: o evento de convergência está na σ\sigma-álgebra de cauda da sequência independente (εn)(\varepsilon_n), de modo que a lei zero–um de Kolmogorov (Teorema 22.9) força sua probabilidade a ser 00 ou 11.

15. (a) Separando no nível θEZ\theta\E Z e usando Cauchy–Schwarz na peça superior,

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)} ,

de modo que (1θ)EZEZ2P(Z>θEZ)(1 - \theta)\E Z \leq \sqrt{\E Z^2\,\P(Z > \theta\E Z)}; eleve ao quadrado. (b) Expanda 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]: a esperança vale 11 quando os índices se emparelham (todos os quatro iguais, ou dois pares distintos, este último em 33 arranjos) e 00 nos demais casos (um sinal desemparelhado tem média nula e se fatora por independência). Logo

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 com Z=Tn2Z = T_n^2, EZ=sn2\E Z = s_n^2, θ=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} .

Se a série convergisse com probabilidade positiva, ela convergiria q.c. (questão 14), de modo que supnTn<\sup_n\abs{T_n} < \infty q.c., e algum MM satisfaria P(supnTn>M)<316\P(\sup_n\abs{T_n} > M) < \frac3{16}; mas, assim que sn>2Ms_n > 2M, P(Tn>M)P(Tn>sn2)316\P(\abs{T_n} > M) \geq \P(\abs{T_n} > \frac{s_n}2) \geq \frac3{16}: contradição. Logo a divergência é quase certa e, com a questão 14, a dicotomia está completa.

16. Aqui xn=nsx_n = n^{-s} e nn2s<\sum_nn^{-2s} < \infty exatamente quando s>12s > \frac12: pelas questões 14–15, nεnns\sum_n\frac{\varepsilon_n}{n^s} converge q.c. se, e somente se, s>12s > \frac12 (para s12s \leq \frac12, divergência q.c.). Para 12<s1\frac12 < s \leq 1 a convergência nunca é absoluta. A comparação é instrutiva: sinais perfeitamente alternados se cancelam na força nsn^{-s} para todo s>0s > 0, ao passo que sinais aleatórios típicos se cancelam apenas na força da raiz quadrada — o passeio aleatório da questão 21 cresce como n\sqrt n, e a soma de Abel converte exatamente esse crescimento em convergência de εnns\sum\varepsilon_nn^{-s} para s>12s > \frac12.

17. (a) coshλ=kλ2k(2k)!\cosh\lambda = \sum_k\frac{\lambda^{2k}}{(2k)!} e eλ2/2=kλ2k2kk!\eu^{\lambda^2/2} = \sum_k\frac{\lambda^{2k}}{2^kk!}; e (2k)!2kk!(2k)! \geq 2^kk! vale termo a termo, pois (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! (cada fator satisfaz k+i2ik + i \geq 2i para iki \leq k), de modo que, de fato, (2k)!2k(k!)22kk!(2k)! \geq 2^k(k!)^2 \geq 2^kk!. (b) Note a0ba \leq 0 \leq b (ZZ é centrada) e, pela convexidade de zeλzz \mapsto \eu^{\lambda z}, para z[a,b]z \in \intcc ab:

eλzbzbaeλa+zabaeλb,logoEeλ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{logo}\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)}

com p=aba[0,1]p = \frac{-a}{b-a} \in \intcc01, t=λ(ba)t = \lambda(b - a), φ(t)=pt+log(1p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t). Então φ(0)=0\varphi(0) = 0, φ(t)=p+pet1p+pet\varphi'(t) = -p + \frac{p\eu^t}{1 - p + p\eu^t} se anula em 00 e φ(t)=ρ(1ρ)14\varphi''(t) = \rho(1 - \rho) \leq \frac14 para ρ=pet1p+pet[0,1]\rho = \frac{p\eu^t}{1 - p + p\eu^t} \in \intcc01: Taylor de ordem 22φ(t)t28=λ2(ba)28\varphi(t) \leq \frac{t^2}8 = \frac{\lambda^2(b-a)^2}8.

18. Para λ>0\lambda > 0, Markov aplicado à variável positiva eλ(SnESn)\eu^{\lambda(S_n - \E S_n)} (Proposição 22.3) e a fórmula do produto para variáveis independentes dão

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),

pela questão 17(b) aplicada a cada XiEXi[aiEXi,biEXi]X_i - \E X_i \in \intcc{a_i - \E X_i}{b_i - \E X_i} centrada (mesma largura). Minimizando o expoente em λ=4tD\lambda = \frac{4t}{D}, D=i(biai)2D = \sum_i(b_i - a_i)^2, obtém-se 2t2D-\frac{2t^2}D. A cauda inferior segue aplicando o resultado a (Xi)(-X_i).

19. Tome t=nεt = n\varepsilon e 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),

que é somável em nn (uma série do tipo geométrico): Borel–Cantelli (Teorema 22.8) dá que q.c. Snnm<ε\abs{\frac{S_n}n - m} < \varepsilon a partir de certa ordem; intersectando sobre ε=1p\varepsilon = \frac1p, obtém-se Snnm\frac{S_n}n \to m q.c. Comparação: Etemadi pede apenas X1L1X_1 \in L^1 e independência dois a dois, e não fornece taxa; Hoeffding pede limitação e independência plena, e fornece uma garantia exponencial explícita em todo nn finito — os dois teoremas respondem a perguntas diferentes sobre o mesmo limite.

20. As g(Uk)g(U_k) são i.i.d. com valores em [0,1]\intcc01 e média g ⁣dλd\int g\,\dd\lambda_d (transferência), de modo que a questão 18 com biai=1b_i - a_i = 1, t=nεt = n\varepsilon dá a cota bilateral 2e2nε2δ2\eu^{-2n\varepsilon^2} \leq \delta assim que e2nε22δ\eu^{2n\varepsilon^2} \geq \frac2\delta, isto é, nlog(2/δ)2ε2n \geq \frac{\log(2/\delta)}{2\varepsilon^2}. Para ε=δ=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 :

cerca de 2650026\,500 amostras garantem 1%1\% de precisão com 99%99\% de confiança — em toda dimensão dd, para todo integrando mensurável com valores em [0,1]\intcc01. A lei forte da questão 11 prometia convergência sem garantia alguma em nn finito; uma malha determinística com kk pontos por eixo custa kdk^d avaliações, exponencial em dd. A concentração é o que faz de Monte Carlo um método e não uma esperança.

21. Independência e a fórmula do produto: 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} pela questão 17(a). Markov em eλSn\eu^{\lambda S_n}:

P(Snx)eλx+nλ2/2=ex2/(2n)no oˊtimo λ=xn,\P(S_n \geq x) \leq \eu^{-\lambda x + n\lambda^2/2} = \eu^{-x^2/(2n)} \qquad\text{no ótimo } \lambda = \frac xn,

e a cota simétrica para Sn-S_n (mesma lei) duplica a constante para Sn\abs{S_n}.

22. Fixe η>0\eta > 0 e ponha xn=(1+η)2nlognx_n = (1 + \eta)\sqrt{2n\log n} para 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}},

somável, já que (1+η)2>1(1 + \eta)^2 > 1. Borel–Cantelli: q.c. Sn<(1+η)2nlogn\abs{S_n} < (1 + \eta)\sqrt{2n\log n} para todo nn grande, de modo que lim supnSn2nlogn1+η\limsup_n\frac{\abs{S_n}}{\sqrt{2n\log n}} \leq 1 + \eta q.c.; intersectar os eventos q.c. para η=1p\eta = \frac1p, pNp \in \N^*, dá o afirmado. O passeio de tamanho nn tem amplitude típica n\sqrt n (sua variância), e mesmo suas piores excursões excedem essa escala em, no máximo, 2logn\sqrt{2\log n}.

23. Com nj=2jn_j = 2^j e x=(1+η)2njloglognjx = (1 + \eta)\sqrt{2n_j\log\log n_j} (definido para j2j \geq 2), a questão 21 dá

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},

somável em jj, pois (1+η)2>1(1 + \eta)^2 > 1: Borel–Cantelli e η=1p\eta = \frac1p dão lim supjSnj/2njloglognj1\limsup_jS_{n_j}/\sqrt{2n_j \log\log n_j} \leq 1 q.c. O que falta para a metade superior completa é a ponte entre os pontos de controle: é preciso mostrar que maxnjnnj+1Sn\max_{n_j\leq n\leq n_{j+1}}S_n excede (1+η)2njloglognj(1+\eta)\sqrt{2n_j\log\log n_j} apenas finitas vezes, o que exige uma desigualdade maximal com caudas gaussianas (a desigualdade de reflexão de Lévy ou a desigualdade de Ottaviani, não demonstradas aqui). A questão 12 é quantitativamente fraca demais: ela limita a probabilidade por

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)},

que tende a 00, mas não é somável em jj: Borel–Cantelli não consegue concluir. A metade inferior da lei do logaritmo iterado aplica o segundo lema de Borel–Cantelli aos incrementos independentes Snj+1SnjS_{n_{j+1}} - S_{n_j}, usando cotas inferiores correspondentes para caudas de tipo gaussiano. Ambos os refinamentos são probabilidade genuína do terceiro ano, um curso adiante; o que este problema entrega sem auxílio é a escala exata do logaritmo iterado ao longo de tempos geométricos.

24. Cada p^i\hat p_i é uma média de nn variáveis indicadoras i.i.d. com valores em [0,1]\intcc01 e média P(Ai)\P(A_i): Hoeffding dá P(p^iP(Ai)>ε)2e2nε2\P(\abs{\hat p_i - \P(A_i)} > \varepsilon) \leq 2\eu^{-2n\varepsilon^2}. A cota da união multiplica por NN. Resolvendo 2Ne2nε2δ2N\eu^{-2n\varepsilon^2} \leq \delta: nln(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2}. Numericamente: ln21060.05=ln(4107)17.5\ln\frac{2\cdot10^6}{0.05} = \ln(4\cdot10^7) \approx 17.5, de modo que n17.5210487600n \geq \frac{17.5}{2\cdot10^{-4}} \approx 87\,600: estimar uma probabilidade com ±1%\pm1\% exige cerca de 1850018\,500 amostras (ln(2/δ)/2ε2\ln(2/\delta)/2\varepsilon^2), e um milhão de probabilidades apenas 4.7\approx 4.7 vezes mais — a uniformidade custa lnN\ln N, não NN: a observação que torna a minimização do risco empírico, e com ela o aprendizado de máquina, estatisticamente possível.

25. As variáveis Xn=εnnαX_n = \frac{\varepsilon_n} {n^\alpha} são independentes, centradas, limitadas, com nV(Xn)=nn2α\sum_n\V(X_n) = \sum_nn^{-2\alpha}. Se α>12\alpha > \frac12: a série das variâncias converge, e o teorema de uma série (Parte VI) dá a convergência q.c. de Xn\sum X_n. Se α12\alpha \leq \frac12: a série das variâncias diverge, e a metade recíproca (o argumento de Paley–Zygmund da Parte VI, aplicável já que as parcelas são limitadas por 11) dá a divergência q.c. A convergência absoluta pede nα<\sum n^{-\alpha} < \infty: α>1\alpha > 1. Em (12,1]\intoc{\frac12}1, a série converge q.c. embora Xn=\sum\abs{X_n} = \infty com certeza: os sinais conspiram para se cancelar, com probabilidade um — convergência por cancelamento, invisível a qualquer teste absoluto e, pela lei zero–um, com um veredicto determinístico ainda assim.