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=X∗P\P_X = X_*\P em R\R (Exercício 11.9), determinada pela função de distribuição FX(t)=P(X≤t)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(X≥0 ou X∈L1(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[(X−EX)2]=E[X2]−(EX)2\V(X) = \E[(X - \E X)^2] = \E[X^2] - (\E X)^2 para X∈L2X \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⁡(−(x−m)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 X≥0X \geq 0 e a>0a > 0: P(X≥a)≤EXa\P(X \geq a) \leq \frac{\E X}{a}; para X∈L2X \in L^2: P(∣X−EX∣≥a)≤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 (X−EX)2(X - \E X)^2. ∎

22.2 Independência

Definição 22.4

Sub-σ\sigma-álgebras A1,…,An⊆A\mathcal A_1, \dots, \mathcal A_n \subseteq \mathcal A são independentes se P(A1∩⋯∩An)=∏P(Ai)\P(A_1\cap\dots\cap A_n) = \prod\P(A_i) para todos Ai∈AiA_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)=Xi−1(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 PX1⊗⋯⊗PXn\P_{X_1}\otimes\cdots\otimes\P_{X_n}. Nesse caso, para gi≥0g_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]=EX EY\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 ⋂iXi−1(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[(Xi−EXi)(Xj−EXj)]=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)n≥1(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 (ω=∑bk2−k\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 2−m2^{-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=∑j≥1bkjn 2−j.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∈[l2−m,(l+1)2−m))=2−m\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 {Un∈Dn}\{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  ⟺  u≤Fμ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(Xn≤t)=P(Un≤Fμ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=1n−1(1−kN),p_n = \prod_{k=1}^{n-1}\Bigl(1 - \frac kN\Bigr),

por condicionamentos sucessivos (ou diretamente: os N(N−1)⋯(N−n+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⁡(1−x)=x+O(x2)-\ln(1 - x) = x + O(x^2):

ln⁡pn=−n(n−1)2N+O(n3N2),logopn≈e−n2/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 n≈2Nln⁡2≈1.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 n∼Nn \sim \sqrt N, e não n∼Nn \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 sup⁡An=⋂N⋃n≥NAn\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 sup⁡An)=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 sup⁡An)=1\P(\limsup A_n) = 1.

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

P(⋂n=NMAnc)=∏n=NM(1−P(An))≤exp⁡(−∑n=NMP(An))→M→∞0\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

(1−x≤e−x1 - x \leq \eu^{-x}; a série diverge). Logo P(⋃n≥NAn)=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 T∈TT \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 ⋂i≤N{Xi∈Bi}\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(T∩T)=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

Xn→XX_n \to X quase certamente se P(Xn→X)=1\P(X_n \to X) = 1; em probabilidade se P(∣Xn−X∣≥ε)→0\P(\abs{X_n - X} \geq \varepsilon) \to 0 para todo ε>0\varepsilon > 0; em LpL^p se E∣Xn−X∣p→0\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(∣Xn−X∣≥ε)≤P(sup⁡m≥n∣Xm−X∣≥ε)↓P(lim sup⁡{∣Xm−X∣≥ε})=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(∣Xn−X∣≥ε)≤ε−p E∣Xn−X∣p\P(\abs{X_n - X} \geq \varepsilon) \leq \varepsilon^{-p}\,\E\abs{X_n - X}^p. (c) Tome nkn_k com P(∣Xnk−X∣≥2−k)≤2−k\P(\abs{X_{n_k} - X} \geq 2^{-k}) \leq 2^{-k}; Borel–Cantelli (1) faz ∣Xnk−X∣<2−k\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 X1∈L2X_1 \in L^2, com m=EX1m = \E X_1:

P(∣Snn−m∣≥ε)≤V(X1)n ε2→n→∞0:\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 :

Snn→m\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 X1∈L1X_1 \in L^1, então

Snn→n→∞q.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 X1∈L4X_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 (Xi↦Xi−mX_i \mapsto X_i - m), suponha m=0m = 0. Expanda:

E[Sn4]=∑i,j,k,lE[XiXjXkXl]=n E[X14]+3n(n−1) (E[X12])2≤C n2,\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(n−1)3n(n-1) termos com dois pares distintos). Markov:

P(∣Snn∣≥ε)=P(Sn4≥n4ε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/n→0S_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 g∈L1([0, 1])g \in L^1(\intcc01) e (Un)(U_n) i.i.d. uniformes (Teorema 22.6), 1n∑k≤ng(Uk)→∫01g\frac1n\sum_{k\leq n}g(U_k) \to \int_0^1g q.c.: integrais por amostragem, em qualquer dimensão, à taxa ∼n−1/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=F−1G = 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λln⁡U-\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(X≤F−1(u))=F(F−1(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}={X≤F−1(u)}\{F(X) \leq u\} = \{X \leq F^{-1}(u)\}): F(X)F(X) é uniforme. Reciprocamente, P(G(U)≤t)=P(U≤F(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(−t≤X≤t)=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λln⁡U≤t)=P(U≥e−λt)=1−e−λ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+s∣T>t)=P(T>s)\P(T > t + s \mid T > t) = \P(T > s) para todos s,t≥0s, 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=∑k≥0k e−λλkk!=λ\E X = \sum_{k\geq0}k\,\eu^{-\lambda} \frac{\lambda^k}{k!} = \lambda, E[X(X−1)]=λ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(1−p)k−1\P(X = k) = p(1-p)^{k-1}): EX=1p\E X = \frac1p, V=1−pp2\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 q≥0q \geq 0 racionais; escrevendo G(1)=e−λG(1) = \eu^{-\lambda} (∈(0, 1)\in \intoo01: G(1)=1G(1) = 1 forçaria G≡1G \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)⊆Xi−1(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; A∩B=A∩C=B∩CA\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(A∩B∩C)=P(HH)=14≠18\P(A\cap B\cap C) = \P(\text{HH}) = \frac14 \neq \frac18: não independentes — CC é 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+ε)log⁡2nR_n \geq (1+\varepsilon)\log_2n apenas um número finito de vezes, e Rn≥log⁡2nR_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 log⁡2n\log_2n.

Solução

Solução de Exercício 22.4.

(a) Seja o texto TT de comprimento LL e q=a−Lq = 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+ε)log⁡2n)≤2−(1+ε)log⁡2n=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=⌈log⁡2sj⌉\ell_j = \lceil\log_2s_j\rceil começando em sj=∑i<jℓis_j = \sum_{i<j}\ell_i; os eventos “o bloco jj é só de uns” são independentes, de probabilidade 2−ℓj≍1sj≍1jlog⁡2j2^{-\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 é, Rsj≥log⁡2sjR_{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))log⁡2n(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/n→m)∈{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 sup⁡∣Xn∣1/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 {R≤c}\{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) Xn→0X_n \to 0 em probabilidade e em todo LpL^p, mas em ponto algum q.c.; (b) Xn→0X_n \to 0 q.c. mas em nenhum LpL^p; (c) Xn→0X_n \to 0 em L1L^1 mas não em L2L^2; (d) e mostre: se Xn→XX_n \to X em probabilidade e ∣Xn∣≤Y∈L1\abs{X_n} \leq Y \in L^1, então Xn→XX_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): ∥Xn∥pp=λ(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 ∥Xn∥p≥n1−1/p≥1\norm{X_n}_p \geq n^{1 - 1/p} \geq 1. (c) Xn=n 1(0, 1/n)X_n = \sqrt n\,\mathbf 1_{\intoo0{1/n}}: E∣Xn∣=n−1/2→0\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 E∣Xn−X∣\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^n−p∣≥ε)≤14nε2\P(\abs{\hat p_n - p} \geq \varepsilon) \leq \frac1{4n\varepsilon^2} (use p(1−p)≤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(1−p)n≤14n\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)2≤0.05\frac1{4n(0.03)^2} \leq 0.05: n≥14⋅0.0009⋅0.05≈5556n \geq \frac{1}{4\cdot0.0009\cdot0.05} \approx 5556. O teorema central do limite justificará n≈1070n \approx 1070 para a mesma garantia: Chebyshev paga sua generalidade com um fator ≈5\approx 5.

Exercício 22.8 ★★★

(Bernstein) Para f∈C([0, 1])f \in \mathcal C(\intcc01), defina o polinômio de Bernstein Bnf(x)=∑k=0n(nk)xk(1−x)n−kf(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 Bnf→fB_nf \to f uniformemente em [0, 1]\intcc01: separe em {∣Snn−x∣≤δ}\{\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 ∥Bnf−f∥∞≤32 ωf(n−1/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(n−1/2))O(\omega_f(n^{-1/2})).

Solução

Solução de Exercício 22.8.

(a) Se Sn∼B(n,x)S_n \sim \mathcal B(n, x), o teorema de transferência dá E[f(Snn)]=∑k(nk)xk(1−x)n−kf(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)∣≤ω(∣u−v∣)\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+(u−x)2δ2)ω(δ)\abs{f(u) - f(x)} \leq \Bigl(1 + \frac{(u - x)^2}{\delta^2}\Bigr)\omega(\delta)

(se ∣u−x∣≤δ\abs{u - x} \leq \delta, claro; caso contrário, ω(∣u−x∣)≤(1+∣u−x∣δ)ω(δ)≤(1+(u−x)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 δ=n−1/2\delta = n^{-1/2}: ∥Bnf−f∥∞≤54 ω(n−1/2)≤32 ω(n−1/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 n−k+1n\frac{n - k + 1}n, os τk\tau_k independentes, e deduza ETn=n Hn∼nln⁡n\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: Tnnln⁡n→1\frac{T_n}{n\ln n} \to 1 em probabilidade. (c) Refine com Borel–Cantelli: mostre diretamente que P(Tn>βnln⁡n)≤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 βnln⁡n\beta n\ln n sorteios, usando 1−x≤e−x1 - x \leq \eu^{-x}) e deduza que, ao longo de n=2mn = 2^m, q.c. Tn≤βnln⁡nT_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 k−1k - 1 tipos coletados, cada sorteio é novo com probabilidade pk=n−k+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=∑knn−k+1=nHn∼nln⁡n\E T_n = \sum_k\frac n{n-k+1} = nH_n \sim n\ln n; V(Tn)=∑1−pkpk2≤n2∑j=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(∣Tn−nHn∣≥εnln⁡n)≤π2n2/6ε2n2ln⁡2n→0\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 nHnnln⁡n→1\frac{nH_n}{n\ln n} \to 1: Tnnln⁡n→1\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(1−1n)t≤n e−t/n\P(T_n > t) \leq n(1 - \frac1n)^{t} \leq n\,\eu^{-t/n}; em t=βnln⁡nt = \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≤βnln⁡nT_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=∑b2k2−kU = \sum b_{2k}2^{-k} (os dígitos de índice par de um ω\omega uniforme) é uniforme e independente de V=∑b2k−12−kV = \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=∑kb2k2−kU = \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)n≥1(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,…,Xn−1)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,…,Xn−1X_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 (n−1)!n!=1n\frac{(n-1)!}{n!} = \frac1n.

(b) Fixe nn e condicione à ordenação relativa de X1,…,Xn−1X_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,…,Rn−1)(R_1, \dots, R_{n-1}), que é função da ordenação relativa das n−1n - 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[#recordes≤n]=Hn≈ln⁡n\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(1n−1n+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+ε)log⁡2nL_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 n2−ℓn2^{-\ell}; Borel–Cantelli ao longo de n=2kn = 2^k). (b) Mostre que q.c. Ln≥(1−ε)log⁡2nL_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−ε)log⁡2n⌉\ell = \lceil(1 - \varepsilon)\log_2n\rceil; os blocos são independentes, cada um todo de caras com probabilidade 2−ℓ2^{-\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 Lnlog⁡2n→1\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 i≤ni \leq n tem probabilidade 2−ℓ2^{-\ell}; cota da união: P(Ln≥ℓ)≤n2−ℓ\P(L_n \geq \ell) \leq n2^{-\ell}. Com ℓn=(1+ε)log⁡2n\ell_n = (1 + \varepsilon)\log_2n: P(Ln≥ℓn)≤n−ε\P(L_n \geq \ell_n) \leq n^{-\varepsilon}. Ao longo de n=2kn = 2^k: ∑k2−kε<∞\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 2k−1<n≤2k2^{k-1} < n \leq 2^k e use a monotonicidade de LnL_n mais log⁡22k−1≤log⁡2n\log_22^{k-1} \leq \log_2n: Ln≤L2k<(1+ε)k≤(1+ε)kk−1log⁡2nL_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−ε)log⁡2n⌉\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 2−ℓ≥n−(1−ε)/22^{-\ell} \geq n^{-(1-\varepsilon)}/2, de modo que

P(Ln<ℓ)≤(1−2−ℓ)m≤exp⁡(−m2−ℓ)≤exp⁡(−c nεlog⁡2n)\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−ε)log⁡2nL_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: Lnlog⁡2n→1\frac{L_n}{\log_2n} \to 1 q.c. Para n=106n = 10^6: log⁡2n≈19.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 — Snn→EX1\frac{S_n}n \to \E X_1 q.c. para Xn∈L1X_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 Xn≥0X_n \geq 0 (separe Xn=Xn+−Xn−X_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 Xn≥0X_n \geq 0.
  2. (Truncamento) Sejam Yn=Xn 1Xn≤nY_n = X_n\,\mathbf 1_{X_n \leq n} e Sn∗=Y1+⋯+YnS_n^* = Y_1 + \dots + Y_n. Mostre que

    ∑n≥1P(Xn≠Yn)=∑n≥1P(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 Sn−Sn∗n→0\frac{S_n - S_n^*}{n} \to 0 q.c.: basta demonstrar Sn∗n→m\frac{S^*_n}n \to m q.c.

  3. Mostre que EYn=E[X11X1≤n]→m\E Y_n = \E\bigl[X_1\mathbf 1_{X_1\leq n}\bigr] \to m (convergência monótona), logo 1n∑k≤nEYk→m\frac1n\sum_{k\leq n}\E Y_k \to m (Cesàro): basta demonstrar Sn∗−ESn∗n→0\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[X12 1X1≤n]\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

    ∑n≥1V(Yn)n2≤∑n≥11n2 E[X121X1≤n]≤C E[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 ∑n≥x1n2≤2max⁡(x,1)\sum_{n \geq x}\frac1{n^2} \leq \frac2{\max(x,1)} para a estimativa interna x2∑n≥xn−2≤2xx^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=⌊αj⌋k_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:

    ∑j≥1P(∣Skj∗−ESkj∗kj∣≥ε)≤1ε2∑j≥11kj2∑n≤kjV(Yn)=1ε2∑n≥1V(Yn)∑j : kj≥n1kj2.\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:kj≥nkj−2≤Cα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 αj≥2\alpha^j \geq 2) e conclua, com a questão 4 e Borel–Cantelli:

    Skj∗−ESkj∗kj→j→∞q.c.0,logoSkj∗kj→m 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 kj≤n≤kj+1k_j \leq n \leq k_{j+1}, use a monotonicidade de Sn∗S^*_n (as parcelas são não negativas!) para mostrar

    kjkj+1 Skj∗kj  ≤  Sn∗n  ≤  kj+1kj Skj+1∗kj+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 inf⁡Sn∗n≤lim sup⁡Sn∗n≤α 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 Sn∗n→m\frac{S_n^*}n \to m q.c., logo (Parte I) a lei forte dos grandes números:

     Snn→n→∞q.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 b≥2b \geq 2: cada dígito 0,…,b−10, \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,…,b−1}\{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 g∈L1([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(max⁡1≤k≤n∣Sk∣≥ε)  ≤  1ε2∑k=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 Sn2≥Sk2+2Sk(Sn−Sk)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 m→∞m \to \infty na desigualdade maximal aplicada a ZN+1,…,ZN+mZ_{N+1}, \dots, Z_{N+m} e, depois, faça N→∞N \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=∑k≤nxkεkT_n = \sum_{k\leq n}x_k\varepsilon_k e sn2=∑k≤nxk2s_n^2 = \sum_{k\leq n}x_k^2, e suponha sn→∞s_n \to \infty. (a) Demonstre a desigualdade de Paley–Zygmund: para Z≥0Z \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 ETn4≤3sn4\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<s≤1\frac12 < s \leq 1, a série converge q.c. enquanto ∑nn−s=∞\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 Snn→m\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 a≤Z≤ba \leq Z \leq b, a<ba < b. Mostre que

    E eλZ≤exp⁡(λ2(b−a)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⁡(1−p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t) com p=−ab−ap = \frac{-a}{b-a} e t=λ(b−a)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 ai≤Xi≤bia_i \leq X_i \leq b_i e Sn=X1+⋯+XnS_n = X_1 + \dots + X_n. Demonstre, para t>0t > 0,

    P(Sn−ESn≥t)≤exp⁡(−2t2∑i=1n(bi−ai)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 E eλ(Sn−ESn)\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(∣Snn−m∣≥ε)≤2exp⁡(−2nε2(b−a)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 Snn→m\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

    n≥log⁡(2/δ)2ε2  ⟹  P(∣1n∑k=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 ε=δ=10−2\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 E eλSn=(cosh⁡λ)n≤enλ2/2\E\,\eu^{\lambda S_n} = (\cosh\lambda)^n \leq \eu^{n\lambda^2/2} e deduza, para x>0x > 0,

    P(Sn≥x)≤e−x2/(2n),P(∣Sn∣≥x)≤2 e−x2/(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 sup⁡n→∞∣Sn∣2nlog⁡n≤1q.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+η)2nlog⁡nx = (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 ∣Sn∣≤n\abs{S_n} \leq n.

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

    lim sup⁡j→∞Snj2njlog⁡log⁡nj≤1q.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 sup⁡n→∞Sn2nlog⁡log⁡n=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 max⁡nj≤n≤nj+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(max⁡i≤N ∣p^i−P(Ai)∣>ε)  ≤  2N e−2nε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: n≥ln⁡(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+n−Sn−n→EX1+−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(Xn≠Yn)=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(X1≥n)≤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 Sn−Sn∗S_n - S_n^* é constante em nn a partir de certa ordem: Sn−Sn∗n→0\frac{S_n - S_n^*}n \to 0 q.c., e as duas somas normalizadas partilham o comportamento assintótico.

3. X11X1≤n↗X1X_1\mathbf 1_{X_1 \leq n} \nearrow X_1: o TCM dá EYn→m\E Y_n \to m; as médias de Cesàro de uma sequência convergente convergem ao mesmo limite: ESn∗n=1n∑k≤nEYk→m\frac{\E S_n^*}n = \frac1n\sum_{k\leq n}\E Y_k \to m. Logo basta demonstrar Sn∗−ESn∗n→0\frac{S^*_n - \E S^*_n}{n} \to 0 q.c.

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

∑nE[X121X1≤n]n2=E[X12 ⁣ ⁣∑n≥max⁡(X1,1) ⁣1n2]≤E[X12⋅4max⁡(X1,1)]≤4 E[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 ∑n≥xn−2≤4x\sum_{n\geq x}n^{-2} \leq \frac4x para x≥1x \geq 1 (para x≥2x \geq 2: ≤1x−1≤2x\leq \frac1{x-1} \leq \frac2x; para 1≤x<21 \leq x < 2: ≤π26≤4x\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 X1≷1X_1 \gtrless 1.

5. A independência dois a dois dá E[(Yi−EYi)(Yj−EYj)]=0\E[(Y_i - \E Y_i)(Y_j - \E Y_j)] = 0 para i≠ji \neq j (a fórmula do produto para duas variáveis), de modo que as variâncias se somam: V(Sk∗)=∑n≤kV(Yn)\V(S^*_k) = \sum_{n\leq k}\V(Y_n). Chebyshev em cada kjk_j e somando:

∑jP(∣Skj∗−ESkj∗∣≥εkj)≤1ε2∑j1kj2∑n≤kjV(Yn)=1ε2∑nV(Yn) ⁣ ⁣∑j:kj≥n ⁣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 αj≥1\alpha^j \geq 1, isto é, para todo j≥0j \geq 0: ⌊x⌋≥x2\lfloor x\rfloor \geq \frac x2 para x≥1x \geq 1). Logo

∑j:kj≥n1kj2≤4∑j:αj≥nα−2j≤41−α−2⋅1n2=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 αj≥n\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á Skj∗−ESkj∗kj→0\frac{S^*_{k_j} - \E S^*_{k_j}}{k_j} \to 0 q.c. e, com a questão 3: Skj∗kj→m\frac{S^*_{k_j}}{k_j} \to m q.c.

7. Yn≥0Y_n \geq 0 torna n↦Sn∗n \mapsto S^*_n não decrescente: para kj≤n≤kj+1k_j \leq n \leq k_{j+1},

Skj∗kj+1≤Sn∗n≤Skj+1∗kj,\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 inf⁡nSn∗n≤lim sup⁡nSn∗n≤α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, p∈N∗p \in \N^*: uma infinidade enumerável de eventos q.c.; em sua interseção, fazendo p→∞p \to \infty: lim⁡Sn∗n=m\lim\frac{S^*_n}n = m q.c. Com as questões 1–3, Snn→EX1\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,…,b−1}\{0, \dots, b-1\} (cada valor do vetor de dígitos ocupa um intervalo de comprimento b−mb^{-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.100100100…2x = 0.100100100\ldots_2 (frequência de uns 13≠12\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 g∈L1([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á

1n∑k=1ng(Uk)→n→∞q.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={max⁡k≤n∣Sk∣≥ε}A = \{\max_{k\leq n}\abs{S_k} \geq \varepsilon\}. Então

ESn2≥∑k=1nE[Sn21Ak]=∑k=1nE[(Sk2+2Sk(Sn−Sk)+(Sn−Sk)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 Sn−SkS_n - S_k, função de (Zk+1,…,Zn)(Z_{k+1}, \dots, Z_n) (Teorema 22.5), de modo que E[Sk1Ak(Sn−Sk)]=E[Sk1Ak] E[Sn−Sk]=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≥ε2∑kP(Ak)=ε2P(A)\E S_n^2 \geq \varepsilon^2\sum_k\P(A_k) = \varepsilon^2\P(A); e ESn2=∑k≤nV(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(max⁡N<k≤N+m∣Sk−SN∣>ε)≤1ε2∑j=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(sup⁡k>N∣Sk−SN∣>ε)≤rN/ε2\P(\sup_{k>N}\abs{S_k - S_N} > \varepsilon) \leq r_N/\varepsilon^2, e rN→0r_N \to 0 por hipótese. Logo, para cada p∈N∗p \in \N^*, P(⋂N{sup⁡k>N∣Sk−SN∣>1p})≤inf⁡Np2rN=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 sup⁡k>N∣Sk−SN∣≤1p\sup_{k>N}\abs{S_k - S_N} \leq \frac1p (intersecte a infinidade enumerável de eventos q.c. sobre pp), de modo que ∣Sk−Sl∣≤2p\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+EZ2 P(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−θ)EZ≤EZ2 P(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,lxixjxkxl E[ε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+3∑i≠jxi2xj2=3sn4−2∑kxk4≤3sn4.\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 sup⁡n∣Tn∣<∞\sup_n\abs{T_n} < \infty q.c., e algum MM satisfaria P(sup⁡n∣Tn∣>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=n−sx_n = n^{-s} e ∑nn−2s<∞\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 s≤12s \leq \frac12, divergência q.c.). Para 12<s≤1\frac12 < s \leq 1 a convergência nunca é absoluta. A comparação é instrutiva: sinais perfeitamente alternados se cancelam na força n−sn^{-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 ∑εnn−s\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+i≥2ik + i \geq 2i para i≤ki \leq k), de modo que, de fato, (2k)!≥2k(k!)2≥2kk!(2k)! \geq 2^k(k!)^2 \geq 2^kk!. (b) Note a≤0≤ba \leq 0 \leq b (ZZ é centrada) e, pela convexidade de z↦eλzz \mapsto \eu^{\lambda z}, para z∈[a, b]z \in \intcc ab:

eλz≤b−zb−a eλa+z−ab−a eλb,logoE eλZ≤b eλa−a eλbb−a=(1−p)e−pt+p e(1−p)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=−ab−a∈[0, 1]p = \frac{-a}{b-a} \in \intcc01, t=λ(b−a)t = \lambda(b - a), φ(t)=−pt+log⁡(1−p+pet)\varphi(t) = -pt + \log(1 - p + p\eu^t). Então φ(0)=0\varphi(0) = 0, φ′(t)=−p+pet1−p+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 ρ=pet1−p+pet∈[0, 1]\rho = \frac{p\eu^t}{1 - p + p\eu^t} \in \intcc01: Taylor de ordem 22 dá φ(t)≤t28=λ2(b−a)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λ(Sn−ESn)\eu^{\lambda(S_n - \E S_n)} (Proposição 22.3) e a fórmula do produto para variáveis independentes dão

P(Sn−ESn≥t)≤e−λt∏i=1nE eλ(Xi−EXi)≤exp⁡(−λt+λ28∑i(bi−ai)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 Xi−EXi∈[ai−EXi, bi−EXi]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(bi−ai)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(b−a)2D = n(b - a)^2:

P(∣Snn−m∣≥ε)≤2exp⁡(−2n2ε2n(b−a)2)=2exp⁡(−2nε2(b−a)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. ∣Snn−m∣<ε\abs{\frac{S_n}n - m} < \varepsilon a partir de certa ordem; intersectando sobre ε=1p\varepsilon = \frac1p, obtém-se Snn→m\frac{S_n}n \to m q.c. Comparação: Etemadi pede apenas X1∈L1X_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 bi−ai=1b_i - a_i = 1, t=nεt = n\varepsilon dá a cota bilateral 2e−2nε2≤δ2\eu^{-2n\varepsilon^2} \leq \delta assim que e2nε2≥2δ\eu^{2n\varepsilon^2} \geq \frac2\delta, isto é, n≥log⁡(2/δ)2ε2n \geq \frac{\log(2/\delta)}{2\varepsilon^2}. Para ε=δ=10−2\varepsilon = \delta = 10^{-2}:

n≥log⁡2002⋅10−4=5.2983…0.0002≈26 492:n \geq \frac{\log 200}{2\cdot10^{-4}} = \frac{5.2983\ldots}{0.0002} \approx 26\,492 :

cerca de 26 50026\,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: E eλSn=(E eλε1)n=(cosh⁡λ)n≤enλ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(Sn≥x)≤e−λx+nλ2/2=e−x2/(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+η)2nlog⁡nx_n = (1 + \eta)\sqrt{2n\log n} para n≥2n \geq 2:

P(∣Sn∣≥xn)≤2exp⁡(−(1+η)2log⁡n)=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+η)2nlog⁡n\abs{S_n} < (1 + \eta)\sqrt{2n\log n} para todo nn grande, de modo que lim sup⁡n∣Sn∣2nlog⁡n≤1+η\limsup_n\frac{\abs{S_n}}{\sqrt{2n\log n}} \leq 1 + \eta q.c.; intersectar os eventos q.c. para η=1p\eta = \frac1p, p∈N∗p \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, 2log⁡n\sqrt{2\log n}.

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

P(Snj≥x)≤exp⁡(−(1+η)2log⁡log⁡nj)=(jlog⁡2)−(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 sup⁡jSnj/2njlog⁡log⁡nj≤1\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 max⁡nj≤n≤nj+1Sn\max_{n_j\leq n\leq n_{j+1}}S_n excede (1+η)2njlog⁡log⁡nj(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+η)2 2njlog⁡log⁡nj=12(1+η)2log⁡(jlog⁡2),\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+1−SnjS_{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^i−P(Ai)∣>ε)≤2e−2nε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 2Ne−2nε2≤δ2N\eu^{-2n\varepsilon^2} \leq \delta: n≥ln⁡(2N/δ)2ε2n \geq \frac{\ln(2N/\delta)}{2\varepsilon^2}. Numericamente: ln⁡2⋅1060.05=ln⁡(4⋅107)≈17.5\ln\frac{2\cdot10^6}{0.05} = \ln(4\cdot10^7) \approx 17.5, de modo que n≥17.52⋅10−4≈87 600n \geq \frac{17.5}{2\cdot10^{-4}} \approx 87\,600: estimar uma probabilidade com ±1%\pm1\% exige cerca de 18 50018\,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 ln⁡N\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)=∑nn−2α\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.

Termos definidos neste capítulo

Ver todos os 395 termos do glossário