3.3 Lei fraca dos grandes números

Nessa seção iremos mostrar um dos resultados mais importantes da Teoria da Probabilidade. O que nossa intuição tem a nos dizer sobre a probabilidade de obtermos um resultado em um dado é 1/61/6? Uma possível explicação seria por simetria, mas e o que podemos dizer no caso de um dado viciado?

Se dizemos a alguém que a probabilidade de obter 66 em um certo dado é 1/101/10, naturalmente a pessoa pode se perguntar como descobrimos isso. Um bom jeito de obter tal medida seria jogar o dado várias vezes independentemente e calcular em qual proporção dos ensaios ele retornou um seis.

O objetivo desta seção é confirmar a validade desse experimento de maneira quantitativa.

{theorem}

Se X1,X2,X_{1},X_{2},\dots são i.i.d.s em 2\mathcal{L}^{2} e definimos

Sn=i=1nXi,S_{n}=\sum_{i=1}^{n}X_{i}, (3.32)

então para todo ε>0\varepsilon>0

limnP[|SnnE(X1)|>ε]=0.\lim_{n\to\infty}P\Big{[}\Big{|}\frac{S_{n}}{n}-E(X_{1})\Big{|}>\varepsilon% \Big{]}=0. (3.33)

Ou seja, SnnE(X1)\tfrac{S_{n}}{n}\to E(X_{1}) em medida (que também chamamos de “em probabilidade”).

Demonstração.

Sabemos que

P[|SnnE(X1)|>ε]\Var(Snn)ε2,P\Big{[}\Big{|}\frac{S_{n}}{n}-E(X_{1})\Big{|}>\varepsilon\Big{]}\leq\frac{% \Var(\tfrac{S_{n}}{n})}{\varepsilon^{2}}, (3.34)

pois E(Sn/n)=1/nE(X1++Xn)=E(X1)E(S_{n}/n)=1/nE(X_{1}+\dots+X_{n})=E(X_{1}).

Mas como \Var(Sn/n)=1/n2\Var(X1++Xn)=(n/n2)\Var(X1)\Var(S_{n}/n)=1/n^{2}\Var(X_{1}+\dots+X_{n})=(n/n^{2})\Var(X_{1}), temos o resultado. ∎

Observe que nós apenas utilizamos que as variáveis XiX_{i} eram independentes duas a duas.

Além disso, obtivemos o seguinte resultado quantitativo que vale mesmo para valores finitos de nn:

{scholia}

Se X1,X2,X_{1},X_{2},\dots são i.i.d.s em 2\mathcal{L}^{2} e definimos Sn=i=1nXiS_{n}=\sum_{i=1}^{n}X_{i} como acima, então, para todo ε>0\varepsilon>0 e n1n\geq 1, temos

P[|SnnE(X1)|>ε]\Var(X1)ε2n.P\Big{[}\Big{|}\frac{S_{n}}{n}-E(X_{1})\Big{|}>\varepsilon\Big{]}\leq\frac{% \Var(X_{1})}{\varepsilon^{2}n}. (3.35)
{corollary}

Se A1,A2,A_{1},A_{2},\dots são eventos independentes dois a dois com P(Ai)=p[0,1]P(A_{i})=p\in[0,1] para todo ii, então

limnP[|#{in;ωAi}np|>ε]=0,\lim_{n\to\infty}P\Big{[}\Big{|}\frac{\#\{i\leq n;\omega\in A_{i}\}}{n}-p\Big{% |}>\varepsilon\Big{]}=0, (3.36)

ou em outras palavras a proporção de ensaios onde o evento AiA_{i} ocorre converge em probabilidade para pp.

Demonstração.

Basta tomar Xi=\1AiX_{i}=\1_{A_{i}} no Teorema 3.3. ∎

{exercise}

Sejam (Xi)i1(X_{i})_{i\geq 1} variáveis \iidcom distribuição Ber(p)(p), p[0,1]p\in[0,1]. Mostre que

limN1Ni=1NXiXi+1=p2, em probabilidade.\lim_{N\to\infty}\frac{1}{N}\sum_{i=1}^{N}X_{i}X_{i+1}=p^{2},\text{ em % probabilidade.} (3.37)
{exercise}

Sejam X1,,XnX_{1},\dots,X_{n} e Y1,,YnY_{1},\dots,Y_{n} variáveis independentes com distribuição \Ber(p)\Ber(p). Defina agora Zi,j=XiYjZ_{i,j}=X_{i}Y_{j}, para i,j{1,,n}i,j\in\{1,\dots,n\} e

  1.  a)

    calcule a esperança de Sn=1n2i=1nj=1nZi,jS_{n}=\tfrac{1}{n^{2}}\sum_{i=1}^{n}\sum_{j=1}^{n}Z_{i,j} e

  2.  b)

    estime P[|SnE(Sn)|>a]P[|S_{n}-E(S_{n})|>a] usando o método do segundo momento. Como esse resultado se compara com o caso em que os Zi,jZ_{i,j} são i.i.d.?

{exercise}

Considere uma rua infinita com casas ii\in\mathbb{Z}. Para todo ii\in\mathbb{Z}, existia uma rua entre as casas ii e i+1i+1, mas após uma grande tempestade essas ruas foram danificadas. Mais precisamente, para cada ii\in\mathbb{Z}, temos variáveis aleatórias XiX_{i} que são i.i.d. com distribuição Ber(p)\text{Ber}(p), onde Xi=1X_{i}=1 indica que o trecho da rua entre as casas ii e i+1i+1 foi danificado e não pode ser utilizado. Defina, para ii\in\mathbb{Z}, RiR_{i} como sendo o número de casas que continuaram acessíveis à casa ii após a tempestade. Por exemplo, se X2X_{-2} e X0=1X_{0}=1 e X1=0X_{-1}=0, temos que a casa 0 somente pode acessar a casa 1-1, logo R0=1R_{0}=1. Nesse contexto,

  1.  a)

    Calcule a distribuição e a esperança de R0R_{0},

  2.  b)

    Use o método do segundo momento para estimar a probabilidade

    P[|1ni=1nRiE(R0)|>a].P\Big{[}\Big{|}\frac{1}{n}\sum_{i=1}^{n}R_{i}-E(R_{0})\Big{|}>a\Big{]}. (3.38)