Inferência Estatística

Semanas 1–2 · 2026 · Modelos probabilísticos, distribuições e convergência

O que é a inferência estatística?

A inferência estatística é o processo de utilizar dados de uma pequena amostra para tirar conclusões ou fazer generalizações sobre uma população mais vasta.

A amostragem parte de uma população com distribuição desconhecida e produz uma amostra observada. A inferência estatística parte da amostra para obter conclusões sobre a população: sob um modelo probabilístico, estimam-se propriedades da distribuição e preveem-se novas observações.

Habitualmente, supomos que os nossos dados seguem uma distribuição de probabilidade e procuramos estimar propriedades da distribuição subjacente (média, variância, covariâncias, outros parâmetros...), para fazer previsões sobre a população a partir da amostra.

A inferência estatística é a base de todas as disciplinas empíricas.

«Uma vez que todos os modelos estão errados, o cientista deve estar atento aos erros que realmente importam.»

Retrato de George E. P. Box.

— George E. P. Box
Science and Statistics (1976), p. 792 · tradução

Fotografia: DavidMCEddy · CC BY-SA 3.0

Uma torre instável de Jenga assente num único bloco com a inscrição Inferência estatística. Os blocos estão identificados com Física, Aprendizagem automática, Economia, Matemática financeira, Biologia, Psicologia e Empregos tecnológicos excessivamente remunerados.

Introdução: trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2020, com uma reta ajustada. Dimensão da amostra 104, declive 0.099 e correlação amostral 0.79. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2020
\(n=104\) observações sintéticas
Declive \(\widehat\beta\)0.099
Correlação \(r\)0.79

Correlação elevada.

As notas dos trabalhos permitem prever as notas dos exames.

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2021, com uma reta ajustada. Dimensão da amostra 98, declive 0.097 e correlação amostral 0.82. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2021
\(n=98\) observações sintéticas
Declive \(\widehat\beta\)0.097
Correlação \(r\)0.82

Observa-se novamente um padrão semelhante.

Existe uma relação causal entre as notas dos trabalhos e as notas dos exames?

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2022, com uma reta ajustada. Dimensão da amostra 112, declive 0.107 e correlação amostral 0.84. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2022
\(n=112\) observações sintéticas
Declive \(\widehat\beta\)0.107
Correlação \(r\)0.84

O padrão mantém-se.

Haverá alguma lei da natureza segundo a qual bons resultados nos trabalhos implicam bons resultados nos exames?

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2023, com uma reta ajustada. Dimensão da amostra 96, declive 0.092 e correlação amostral 0.80. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2023
\(n=96\) observações sintéticas
Declive \(\widehat\beta\)0.092
Correlação \(r\)0.80

A correlação positiva continua visível.

Como poderia isto alguma vez mudar?

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2024, com uma reta ajustada. Dimensão da amostra 108, declive 0.049 e correlação amostral 0.33. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2024
\(n=108\) observações sintéticas
Declive \(\widehat\beta\)0.049
Correlação \(r\)0.33

O que aconteceu aqui?

O declive diminuiu. A correlação é mais fraca. As notas dos trabalhos aumentaram.

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Trabalhos e exames

Classificações sintéticas de trabalhos e exames em 2025, com uma reta ajustada. Dimensão da amostra 101, declive 0.015 e correlação amostral 0.08. O eixo horizontal vai de 0 a 100 e o vertical de 0 a 10.

2025
\(n=101\) observações sintéticas
Declive \(\widehat\beta\)0.015
Correlação \(r\)0.08

Quase não resta correlação linear.

O trabalho continua a medir aquilo que o exame mede?

Dados sintéticos. Inspirado em T. Osborne, QuSoft (2026).

Programa e avaliação

A

Teoria da probabilidade

Variáveis aleatórias, distribuições, convergência, teoremas limite

B

Modelos não paramétricos

Estimação, ajustamento, ordens, associação

C

Modelos paramétricos

Suficiência, verosimilhança, estimação, testes

D

Estatística bayesiana

Distribuições a priori e a posteriori, conjugação, estimadores de Bayes

Avaliação · 2026/27

Exame: 100%

Antigo · 2025/2026 · Já não se aplica

Exame 60% · Projeto/apresentação 25% · Relatório escrito 15%

Semanas 1–2: as primeiras quatro aulas

Distribuições unidimensionais

Espaços de probabilidade, elementos aleatórios, funções de distribuição, quantis, densidades

Vetores aleatórios

Distribuições conjuntas, condicionamento, σ-álgebras geradas, independência

Estatística descritiva

Valores esperados, momentos, covariância, famílias usuais, transformadas

Convergência

Quase certa, em \(L^\infty\), em probabilidade, em \(L^p\), em distribuição

Espaços de probabilidade e distribuições unidimensionais

01 · SEMANA 1 · AULA 1

Espaços de probabilidade, elementos aleatórios, funções de distribuição e quantis.

Breve história da teoria da probabilidade

Séc. IV a.C.

Aristóteles: a probabilidade não faz parte da matemática. A matemática descreve o que é logicamente necessário, enquanto a probabilidade trata do que é plausível. Visão predominante até ao século XVII.

1654

Blaise Pascal e Pierre de Fermat estabelecem fundamentos da probabilidade matemática na sua correspondência sobre problemas de jogos de azar.

1898–1902

Émile Borel e Henri Lebesgue introduzem as \(\sigma\)-álgebras e a teoria da medida para formalizar a integração e o cálculo de volumes.

1900

David Hilbert apresenta uma lista de problemas em aberto no Congresso Internacional de Matemáticos, em Paris. Um deles: encontrar uma formalização axiomática rigorosa da teoria da probabilidade que permita demonstrar os teoremas limite da física estatística.

1924

Stefan Banach e Alfred Tarski mostram que uma bola em \(\mathbb R^3\) pode ser decomposta num número finito de partes e reconstituída em duas cópias de igual volume, por rotações e translações. Se se atribuir um volume a subconjuntos arbitrários de \(\mathbb{R}^3\), a noção de volume torna-se paradoxal. Daí o recurso a \(\sigma\)-álgebras.

1933

Andrey Kolmogorov axiomatiza a probabilidade usando a teoria da medida. Formulação moderna: uma \(\sigma\)-álgebra e três axiomas da probabilidade.

\(\sigma\)-álgebras

Seja \(\Omega\) um conjunto não vazio e seja \(\mathcal P(\Omega)\) o seu conjunto das partes.

Definição 1.1

Uma família \(\mathcal F\subseteq\mathcal P(\Omega)\) chama-se \(\sigma\)-álgebra se

\[ \begin{aligned} \text{(i)}\quad&\Omega\in\mathcal F,\\ \text{(ii)}\quad&A\in\mathcal F\ \Longrightarrow\ A^c\in\mathcal F,\\ \text{(iii)}\quad&A_1,A_2,\ldots\in\mathcal F\ \Longrightarrow\ \bigcup_{n=1}^{\infty}A_n\in\mathcal F. \end{aligned} \]

O par \((\Omega,\mathcal F)\) chama-se espaço mensurável. Os elementos de \(\mathcal F\) chamam-se acontecimentos.

Espaços de probabilidade e axiomas de Kolmogorov

Definição 1.2

Uma aplicação \(\mathbb P:\mathcal F\to\mathbb R\) chama-se medida de probabilidade em \((\Omega,\mathcal F)\) se satisfaz os axiomas de Kolmogorov, isto é,

\[ \begin{aligned} \text{(i)}\quad &\mathbb P(A)\in [0, 1] &&\text{para todo o }A\in\mathcal F,\\[2mm] \text{(ii)}\quad &\mathbb P(\Omega)=1,\\[2mm] \text{(iii)}\quad &\mathbb P\!\left(\bigcup_{n=1}^{\infty}A_n\right) =\sum_{n=1}^{\infty}\mathbb P(A_n) &&\text{para }A_n\in\mathcal F\text{ disjuntos dois a dois}. \end{aligned} \]

O triplo \((\Omega,\mathcal F,\mathbb P)\) chama-se espaço de probabilidade.

Consequências dos axiomas de Kolmogorov

Proposição 1.3

Para acontecimentos \(A,B\) num espaço de probabilidade \((\Omega,\mathcal F,\mathbb P)\),

\[ \begin{aligned} \mathbb P(A^c)&=1-\mathbb P(A),\\ \mathbb P(A\cup B)&=\mathbb P(A)+\mathbb P(B)-\mathbb P(A\cap B). \end{aligned} \]

Se \(A\subseteq B\), então \(\mathbb P(B\setminus A)=\mathbb P(B)-\mathbb P(A)\) e \(\mathbb P(A)\le \mathbb P(B)\).

Demonstração. A aditividade finita decorre da aditividade numerável e de \(\mathbb P(\varnothing)=0\). Aplique-se às decomposições disjuntas

\[ \Omega=A\mathbin{\dot\cup}A^c,\qquad A\cup B=A\mathbin{\dot\cup}(B\setminus A),\qquad B=(A\cap B)\mathbin{\dot\cup}(B\setminus A). \]

Se \(A\subseteq B\), a última decomposição é \(B=A\mathbin{\dot\cup}(B\setminus A)\). A não negatividade dá \(\mathbb P(B)\ge \mathbb P(A)\).

Sucessões monótonas de acontecimentos

Definição 1.4

Sejam \((\Omega,\mathcal F,\mathbb P)\) um espaço de probabilidade e \(A,B,A_1,A_2,\ldots\in\mathcal F\).

  • Escrevemos \(A_n\nearrow A\) se \(A_n\subseteq A_{n+1}\) para todo o \(n\ge1\) e \(A=\bigcup_{n\ge1}A_n\).
  • Escrevemos \(A_n\searrow B\) se \(A_n\supseteq A_{n+1}\) para todo o \(n\ge1\) e \(B=\bigcap_{n\ge1}A_n\).

Continuidade da probabilidade

Proposição 1.5

Sejam \(A,B,A_1,A_2,\ldots\) acontecimentos num espaço de probabilidade. Então

\[ A_n\nearrow A\ \Longrightarrow\ \mathbb P(A_n)\nearrow \mathbb P(A), \qquad A_n\searrow B\ \Longrightarrow\ \mathbb P(A_n)\searrow \mathbb P(B). \]

Demonstração. No caso crescente, sejam \(A_0=\varnothing\) e \(D_n=A_n\setminus A_{n-1}\). Estes conjuntos são disjuntos, pelo que

\[ \mathbb P(A)=\sum_{n=1}^{\infty}\mathbb P(D_n) =\lim_{N\to\infty}\sum_{n=1}^{N}\mathbb P(D_n) =\lim_{N\to\infty}\mathbb P(A_N). \]

Caso decrescente. De \(A_n^c\nearrow B^c\) resulta \(\mathbb P(A_n)=1-\mathbb P(A_n^c)\searrow 1-\mathbb P(B^c)=\mathbb P(B)\).

Definição 1.6 · σ-álgebra de Borel

Num espaço métrico \(S\), a menor \(\sigma\)-álgebra que contém todos os abertos de \(S\) chama-se \(\sigma\)-álgebra de Borel e denota-se por \(\mathcal B(S)\).

Aplicações mensuráveis

Sejam \((S,\mathcal S)\) e \((T,\mathcal T)\) espaços mensuráveis.

Definição 1.7

Uma aplicação \(g:S\to T\) chama-se mensurável se

\[ g^{-1}(B)=\{s\in S:g(s)\in B\}\in\mathcal S \qquad\text{para todo o }B\in\mathcal T. \]

Exemplo 1.8

(i) Aplicações constantes. Fixado \(t_0\in T\), \(g(s)=t_0\) é mensurável: \(g^{-1}(B)=S\) se \(t_0\in B\), e \(g^{-1}(B)=\varnothing\) caso contrário.

(ii) Aplicações contínuas. Se os espaços métricos \(S\) e \(T\) estão munidos das suas \(\sigma\)-álgebras de Borel, toda a aplicação contínua \(g:S\to T\) é mensurável.

As imagens inversas preservam complementares e uniões numeráveis. Assim, \(\{B\subseteq T:g^{-1}(B)\in\mathcal B(S)\}\) é uma \(\sigma\)-álgebra. Pela continuidade, contém todos os abertos de \(T\) e, portanto, \(\mathcal B(T)\).

Elementos aleatórios

Sejam \((\Omega,\mathcal F,\mathbb P)\) um espaço de probabilidade e \((S,\mathcal S)\) um espaço mensurável.

Definição 1.9

Uma aplicação mensurável \(X:(\Omega,\mathcal F)\to(S,\mathcal S)\) chama-se elemento aleatório em \((S,\mathcal S)\). Equivalentemente,

\[ X^{-1}(B)=\{\omega\in\Omega:X(\omega)\in B\}\in\mathcal F \qquad\text{para todo o }B\in\mathcal S. \]

Tipos de elementos aleatórios

Definição 1.10
Espaço de estados \(S\) Terminologia Valor no resultado \(\omega\)
\(\mathbb R\) Variável aleatória Um número real \(X(\omega)\)
\(\mathbb R^d\), \(d\ge2\) Vetor aleatório \((X_1(\omega),\ldots,X_d(\omega))\)
\(C([0,1])\) Função aleatória Uma função contínua \(t\mapsto X(\omega)(t)\)

Munimos \(C([0,1])\) da norma do supremo e da sua \(\sigma\)-álgebra de Borel.

Mensurabilidade de vetores aleatórios

Proposição 1.11

Sejam \((\Omega,\mathcal F)\) um espaço mensurável e \(d\ge1\) um inteiro. Munimos \(\mathbb R\) e \(\mathbb R^d\) das suas \(\sigma\)-álgebras de Borel.

Uma aplicação \(\mathbf X=(X_1,\ldots,X_d):\Omega\to\mathbb R^d\) é mensurável se e só se cada coordenada \(X_j\) é mensurável.

Demonstração \((\Rightarrow)\). Fixe-se \(j\in\{1,\ldots,d\}\). A projeção \(\pi_j:\mathbb R^d\to\mathbb R\), \(\pi_j(x)=x_j\), é contínua.

Para as \(\sigma\)-álgebras de Borel, \(\pi_j\) é mensurável pelo Exemplo 1.8 (ii).

Se \(\mathbf X\) é mensurável, então, para todo o \(B\in\mathcal B(\mathbb R)\),

\[ X_j^{-1}(B)=\mathbf X^{-1}\!\left(\pi_j^{-1}(B)\right)\in\mathcal F. \]

Logo, \(X_j\) é mensurável para todo o \(j\).

Mensurabilidade de vetores aleatórios II

Proposição 1.11 · recíproca

Se cada coordenada \(X_j\) é mensurável, então \(\mathbf X=(X_1,\ldots,X_d)\) é mensurável.

Demonstração \((\Leftarrow)\). Defina-se

\[ \mathcal C=\{A\subseteq\mathbb R^d:\mathbf X^{-1}(A)\in\mathcal F\}. \]

As imagens inversas preservam complementares e uniões numeráveis, e \(\mathbf X^{-1}(\mathbb R^d)=\Omega\). Assim, \(\mathcal C\) é uma \(\sigma\)-álgebra.

Para intervalos abertos \(I_1,\ldots,I_d\), seja \(R=\prod_{j=1}^d I_j\). Então

\[ \mathbf X^{-1}(R)=\bigcap_{j=1}^dX_j^{-1}(I_j)\in\mathcal F. \]

Os retângulos abertos \(I_1 \times \ldots \times I_d\) com extremos racionais formam uma base numerável da topologia de \(\mathbb R^d\), pelo que geram \(\mathcal B(\mathbb R^d)\). Logo, \(\mathcal B(\mathbb R^d)\subseteq\mathcal C\), o que prova a mensurabilidade.

Distribuição de um elemento aleatório

Seja \(X:(\Omega,\mathcal F,\mathbb P)\to(S,\mathcal S)\) um elemento aleatório.

Definição 1.12

A distribuição de \(X\) define-se como a medida de probabilidade \(\mathbb P_X\) em \((S,\mathcal S)\) dada por

\[ \mathbb P_X(B)=\mathbb P(X\in B)=\mathbb P\!\left(X^{-1}(B)\right), \qquad B\in\mathcal S. \]

Demonstração de que \(\mathbb P_X\) é uma medida de probabilidade. A mensurabilidade de \(X\) dá \(X^{-1}(B)\in\mathcal F\) para todo o \(B\in\mathcal S\). Logo, \(\mathbb P_X\) está bem definida, \(\mathbb P_X(B)\ge0\) e \(\mathbb P_X(S)=\mathbb P(\Omega)=1\).

Para \(B_1,B_2,\ldots\in\mathcal S\) disjuntos dois a dois, seja \(A_n=X^{-1}(B_n)\). Estes acontecimentos são disjuntos dois a dois, pelo que a aditividade numerável de \(\mathbb P\) dá

\[ \mathbb P_X\!\left(\bigcup_{n\ge1}B_n\right) =\mathbb P\!\left(\bigcup_{n\ge1}A_n\right) =\sum_{n\ge1}\mathbb P(A_n) =\sum_{n\ge1}\mathbb P_X(B_n). \]

Transformações mensuráveis

Proposição 1.13

Sejam \(X:(\Omega,\mathcal F,\mathbb P)\to(S,\mathcal S)\) um elemento aleatório e \(g:(S,\mathcal S)\to(T,\mathcal T)\) uma aplicação mensurável.

Então \(Y=g\circ X\) é um elemento aleatório em \((T,\mathcal T)\) com distribuição

\[ \mathbb P_Y(C)=\mathbb P_X\!\left(g^{-1}(C)\right) \qquad\text{para todo o }C\in\mathcal T. \]

Demonstração. Seja \(C\in\mathcal T\). Como \(g\) é mensurável, \(g^{-1}(C)\in\mathcal S\). Da mensurabilidade de \(X\) resulta

\[ Y^{-1}(C)=(g\circ X)^{-1}(C)=X^{-1}\!\left(g^{-1}(C)\right)\in\mathcal F. \]

Logo, \(Y\) é um elemento aleatório. Pelas definições de \(\mathbb P_Y\) e \(\mathbb P_X\),

\[ \mathbb P_Y(C)=\mathbb P\!\left(Y^{-1}(C)\right) =\mathbb P\!\left(X^{-1}\!\left(g^{-1}(C)\right)\right) =\mathbb P_X\!\left(g^{-1}(C)\right). \]

Exemplo: lançamento de moedas

Exemplo 1.14

Lançam-se duas moedas equilibradas, independentemente. Denotamos cara por \(H\) e coroa por \(T\).

\[ \Omega=\{TT,HT,TH,HH\},\qquad \mathcal F=\mathcal P(\Omega),\qquad \mathbb P(A)=\frac{|A|}{4}\quad(A\in\mathcal F). \]

Seja \(X:\Omega\to\mathbb R\) o número de caras.

\(\omega\) \(TT\) \(HT\) \(TH\) \(HH\)
\(X(\omega)\) 0 1 1 2

Por exemplo, \(X^{-1}(\{1\})=\{HT,TH\}\). Da tabela, obtemos a distribuição \(\mathbb P_X\):

\[ \mathbb P_X(\{0\})=\frac14,\qquad \mathbb P_X(\{1\})=\frac12,\qquad \mathbb P_X(\{2\})=\frac14. \]

Por outro lado, para a variável aleatória \(Y=\mathbf1_{[1,\infty)}(X)\), temos \(\mathbb P_Y(\{1\})=\mathbb P_X(\{1,2\})=3/4\).

Funções de distribuição

Definição 1.15

Para uma variável aleatória real \(X\), a sua função de distribuição (FD) define-se por \(F_X(x)=\mathbb P(X\le x)\).

Dizemos que uma função \(F\) é contínua à direita, se \(\displaystyle\lim_{h\searrow0}F(x+h)=F(x)\) para todo o \(x\).

Teorema 1.16

Uma função \(F:\mathbb R\to[0,1]\) é uma FD se e só se:

  1. \(F\) é não decrescente.
  2. \(F\) é contínua à direita.
  3. \(\displaystyle\lim_{x\to-\infty}F(x)=0\) e \(\displaystyle\lim_{x\to+\infty}F(x)=1\).

Toda a função \(F\) nestas condições determina uma única medida de probabilidade de Borel em \(\mathbb R\).

Continuidade à direita e limites nos extremos

Proposição 1.17

Sejam \(X\) uma variável aleatória real e \(F_X(x)=\mathbb P(X\le x)\). Então \(F_X\) é contínua à direita e

\[ \lim_{x\to-\infty}F_X(x)=0, \qquad \lim_{x\to+\infty}F_X(x)=1. \]

Demonstração. Fixe-se \(x\in\mathbb R\) e seja \(h_n\searrow0\) com \(h_n>0\). Pela Proposição 1.5,

\[ X^{-1}((-\infty,x+h_n])\searrow X^{-1}((-\infty,x]) \quad\Longrightarrow\quad F_X(x+h_n)\searrow F_X(x). \]

Para os limites nos extremos, o mesmo resultado dá

\[ \begin{aligned} X^{-1}((-\infty,-n])\searrow\varnothing &\quad\Longrightarrow\quad F_X(-n)\searrow0,\\ X^{-1}((-\infty,n])\nearrow\Omega &\quad\Longrightarrow\quad F_X(n)\nearrow1. \end{aligned} \]

A inclusão \(X^{-1}((-\infty,x])\subseteq X^{-1}((-\infty,y])\) para \(x\le y\) implica que \(F_X\) é não decrescente. Assim, os limites ao longo dos inteiros estendem-se a \(x\to\pm\infty\).

Limites à esquerda e probabilidades pontuais

Proposição 1.18

Seja \(X\) uma variável aleatória real com FD \(F_X\). Para todo o \(x\in\mathbb R\), existe o limite à esquerda \(F_X(x^-)=\lim_{t\nearrow x}F_X(t)\) e

\[ \begin{aligned} F_X(x^-)&=\mathbb P(X<x),\\ F_X(x)-F_X(x^-)&=\mathbb P(X=x). \end{aligned} \]

Demonstração. A função \(F_X\), limitada e não decrescente, tem limite à esquerda em cada \(x\). Como \(X^{-1}((-\infty,x-1/n])\nearrow X^{-1}((-\infty,x))\), a continuidade da probabilidade dá

\[ F_X(x^-)=\lim_{n\to\infty}\mathbb P(X\le x-1/n)=\mathbb P(X<x). \]

A segunda identidade decorre da aditividade aplicada à decomposição disjunta

\[ X^{-1}((-\infty,x])=X^{-1}((-\infty,x))\mathbin{\dot\cup}X^{-1}(\{x\}). \]

Átomos e saltos da função de distribuição

Definição 1.19

Um ponto \(x\in\mathbb R\) chama-se átomo da distribuição \(\mathbb P_X\) se \(\mathbb P(X=x)>0\).

Pela Proposição 1.18, os átomos são precisamente as descontinuidades por salto de \(F_X\). O salto em \(x\) tem amplitude \(\mathbb P(X=x)\).

Exemplo 1.20

Se \(F_X(2^-)=0.55\) e \(F_X(2)=0.70\), então

\[ \mathbb P(X=2)=0.15. \]

Uma função de distribuição com coordenadas verticais à escala. Um salto de 0.55 para 0.70 em x igual a 2 tem amplitude 0.15.

Probabilidades de intervalos a partir da FD

Proposição 1.21

Seja \(X\) uma variável aleatória real com FD \(F_X\). Para quaisquer \(a,b\in\mathbb R\) com \(a<b\),

\[ \begin{aligned} \mathbb P(a<X\le b)&=F_X(b)-F_X(a),\\ \mathbb P(a\le X\le b)&=F_X(b)-F_X(a^-). \end{aligned} \]

Demonstração. As seguintes decomposições são disjuntas:

\[ \begin{aligned} X^{-1}((-\infty,b])&=X^{-1}((-\infty,a])\mathbin{\dot\cup}X^{-1}((a,b]),\\ X^{-1}((-\infty,b])&=X^{-1}((-\infty,a))\mathbin{\dot\cup}X^{-1}([a,b]). \end{aligned} \]

Aplique-se a aditividade, usem-se \(\mathbb P(X\le a)=F_X(a)\) e \(\mathbb P(X<a)=F_X(a^-)\) e rearranjem-se os termos.

Número de descontinuidades da FD

Proposição 1.22

Se \(X\) é uma variável aleatória real, a sua FD \(F_X\) tem, no máximo, um conjunto numerável de descontinuidades.

Demonstração. Para cada inteiro \(m\ge1\), defina-se

\[ D_m=\{x\in\mathbb R:\mathbb P(X=x)\ge1/m\}. \]

O conjunto \(D_m\) tem no máximo \(m\) pontos: quaisquer \(m+1\) acontecimentos singulares deste tipo, disjuntos, teriam probabilidade total superior a \(1\).

Todo o átomo pertence a algum \(D_m\), pelo que

\[ \{x\in\mathbb R:\mathbb P(X=x)>0\}=\bigcup_{m=1}^{\infty}D_m \]

é no máximo numerável. Como \(F_X\) é não decrescente e contínua à direita, toda a descontinuidade é um salto.

A inversa generalizada

Definição 1.23

Seja \(F\) uma função com as propriedades do Teorema 1.16. A sua inversa generalizada define-se por

\[ F^{-1}(u)=\inf\{x\in\mathbb R:F(x)\ge u\},\qquad 0<u<1. \]

Proposição 1.24

Seja \(F\) nas condições do Teorema 1.16. Para todo o \(u\in(0,1)\), \(F^{-1}(u)\) é finito e

\[ F\!\left(F^{-1}(u)\right)\ge u. \]

Demonstração. Pelos limites nos extremos, \(\{x\in\mathbb R:F(x)\ge u\}\) é não vazio e minorado. Tome-se \(q=F^{-1}(u)\) e \(x_n\searrow q\) com \(F(x_n)\ge u\). Pela continuidade à direita,

\[ F(q)=\lim_{n\to\infty}F(x_n)\ge u. \]

Uma função de distribuição prescrita

Proposição 1.25

Seja \(F\) nas condições do Teorema 1.16 e seja \(Q=F^{-1}\).

Então \(Q:(0,1)\to\mathbb R\) é mensurável. Se \(U\) é uniforme em \((0,1)\), \(Q(U)\) tem FD \(F\).

Demonstração. Para \(u\in(0,1)\) e \(x\in\mathbb R\),

\[ Q(u)\le x\quad\Longleftrightarrow\quad u\le F(x). \]

A implicação direta decorre de \(F(x)\ge F(Q(u))\ge u\). A recíproca decorre da definição de ínfimo.

Assim, \(Q^{-1}(({-\infty},x])=(0,1)\cap(0,F(x)]\) é boreliano. Estas semirretas geram \(\mathcal B(\mathbb R)\), pelo que \(Q\) é mensurável. Além disso,

\[ \mathbb P(Q(U)\le x)=\mathbb P(U\le F(x))=F(x). \]

Passando aos complementares, obtém-se também \(F(x)<u\Longleftrightarrow x<Q(u)\).

Continuidade dos quantis

Proposição 1.26

Para uma FD \(F\), a inversa generalizada \(Q=F^{-1}:(0,1)\to\mathbb R\) é não decrescente e contínua à esquerda.

Demonstração. Para \(0<u<v<1\),

\[ \{x\in\mathbb R:F(x)\ge v\}\subseteq\{x\in\mathbb R:F(x)\ge u\} \quad\Longrightarrow\quad Q(u)\le Q(v). \]

Para provar a continuidade à esquerda, seja \(u_n\nearrow u\in(0,1)\). Então \(Q(u_n)\nearrow\ell\le Q(u)\) e

\[ F(\ell)\ge F(Q(u_n))\ge u_n \quad\Longrightarrow\quad F(\ell)\ge u. \]

A relação inversa dá \(Q(u)\le\ell\). Portanto,

\[ \lim_{v\nearrow u}Q(v)=Q(u). \]

Saltos da FD e patamares dos quantis

Exemplo 1.27

Para a variável das duas moedas, \(\mathbb P(X=0)=1/4\), \(\mathbb P(X=1)=1/2\) e \(\mathbb P(X=2)=1/4\).

A função de distribuição do número de caras em dois lançamentos de uma moeda equilibrada e a sua inversa generalizada, contínua à esquerda. Um salto de um meio em x=1 corresponde a um patamar para um quarto menor que u menor ou igual a três quartos. Um patamar da FD de 0 a 1 corresponde a um salto à direita do quantil em u=um quarto.

\[ F(x^-)<u\le F(x)\quad\Longrightarrow\quad F^{-1}(u)=x. \]

O patamar \(F(x)=1/4\) em \([0,1)\) dá um salto à direita de \(F^{-1}\) em \(u=1/4\).

Distribuições discretas

Definição 1.28

A medida de Dirac em \(a\) define-se por \(\delta_a(A)=\mathbf1_A(a)\).

Uma distribuição chama-se discreta se, para pontos distintos \(\xi_j\) num conjunto finito ou numerável,

\[ \mathbb P_X=\sum_j p_j\delta_{\xi_j}, \qquad p_j\ge0,\qquad\sum_jp_j=1. \]

A função de probabilidade define-se por \(p_X(\xi_j)=p_j\) e \(p_X(x)=0\) nos outros pontos.

O suporte, conjunto dos átomos, define-se por \(\operatorname{supp}(X):=\{x\in\mathbb R:p_X(x)>0\}\). Assim,

\[ \mathbb P(X\in A)=\sum_{\xi_j\in A}p_j, \qquad F_X(x)=\sum_{\xi_j\le x}p_j. \]

Para suporte finito, a FD é uma função em escada com salto \(p_X(x)\) em cada \(x\in\operatorname{supp}(X)\).

Distribuições absolutamente contínuas

Definição 1.29

Uma distribuição chama-se absolutamente contínua se existe uma função mensurável \(f_X:\mathbb R\to[0,\infty)\) tal que

\[ \int_{\mathbb R}f_X(x)\,dx=1, \qquad \mathbb P(X\in A)=\int_A f_X(x)\,dx \quad(A\in\mathcal B(\mathbb R)). \]

A função \(f_X\) chama-se função de densidade de probabilidade e

\[ F_X(x)=\int_{-\infty}^{x}f_X(t)\,dt. \]

A densidade chama-se derivada de Radon–Nikodym \(d\mathbb P_X/dx\) e é definida a menos de igualdade quase em toda a parte.

Derivadas de Radon–Nikodym

Proposição 1.30

Seja \(X\) de distribuição absolutamente contínua, com densidade \(f_X\). A sua FD \(F_X\) é contínua e \(F_X'(x)=f_X(x)\) para quase todo o \(x\) (Lebesgue). A igualdade vale também em cada ponto de continuidade de \(f_X\).

Demonstração nos pontos de continuidade de \(f_X\). Para \(h\ne0\),

\[ \frac{F_X(x+h)-F_X(x)}{h} =\frac1h\int_x^{x+h}f_X(t)\,dt \ \longrightarrow\ f_X(x). \]

O teorema fundamental para integrais de Lebesgue dá a continuidade de \(F_X\) e a afirmação geral quase em toda a parte. Em particular,

\[ \mathbb P(X=x)=\int_{\{x\}}f_X(t)\,dt=0 \qquad\text{para todo o }x. \]

Densidade não é massa pontual: \(U\sim U(0,1/2)\) tem \(f_U(x)=2\) em \((0,1/2)\).

Ficha 1: distribuição uniforme em \([-2,2]\)

Exercício 1.31 · Ficha 1, questão 11

Seja \(X\) uma variável com função de densidade de probabilidade

\[ f_X(x)=c\,\mathbf1_{[-2,2]}(x), \]

  1. Determine \(c\).
  2. Determine e esboce a função de distribuição \(F_X\).
  3. Determine o quantil de ordem \(0.95\).

Uma distribuição mista

Exemplo 1.32

Para \(0<\pi<1\), seja \(\mathbb P_X=\pi\delta_0+(1-\pi)U(0,1)\). Então

\[ F_X(x)= \begin{cases} 0,&x<0,\\ \pi+(1-\pi)x,&0\le x<1,\\ 1,&x\ge1. \end{cases} \]

O átomo em zero tem massa \(\pi\). As probabilidades em \((0,1)\) obtêm-se integrando a constante \(1-\pi\).

Esta distribuição não é discreta nem absolutamente contínua.

Vetores aleatórios e independência

02 · SEMANA 1 · AULA 2

Distribuições conjuntas, independência e vetores aleatórios discretos e absolutamente contínuos.

Vetores aleatórios e funções de distribuição conjuntas

Definição 2.1

Para um vetor aleatório \(\mathbf X=(X_1,\ldots,X_d)\), a distribuição conjunta define-se por

\[ \mathbb P_{\mathbf X}(A)=\mathbb P(\mathbf X\in A), \qquad A\in\mathcal B(\mathbb R^d). \]

Para \(x=(x_1,\ldots,x_d)\), defina-se \(S_x=\prod_{j=1}^d(-\infty,x_j]\). A função de distribuição conjunta define-se por

\[ F_{\mathbf X}(x)=\mathbb P(X_1\le x_1,\ldots,X_d\le x_d)=\mathbb P_{\mathbf X}(S_x). \]

É não decrescente em cada coordenada e determina univocamente a distribuição conjunta.

Função de distribuição conjunta dos resultados de dois dados equilibrados com 36 pares ordenados equiprováveis. A superfície em escada tem altura ij/36 em cada retângulo [i,i+1) × [j,j+1), para i e j de 1 a 6. Os eixos horizontais são x e y e o eixo vertical é F_X(x,y). O ponto assinalado em (3,4) tem altura 1/3.

FD conjuntas: limites laterais

Proposição 2.2

Sejam \(\mathbf X\) um vetor aleatório em \(\mathbb R^d\), \(x\in\mathbb R^d\) e \(h_n>0\) com \(h_n\searrow0\). Escreva-se \(\mathbf1=(1,\ldots,1)\) e \(S_x=\prod_{j=1}^d(-\infty,x_j]\). Então

\[ \begin{aligned} F_{\mathbf X}(x+h_n\mathbf1)&\searrow F_{\mathbf X}(x),\\ F_{\mathbf X}(x-h_n\mathbf1)&\nearrow \mathbb P_{\mathbf X}(S_x^\circ), \end{aligned} \] onde \(S_x^\circ=\prod_j(-\infty,x_j)\) denota o interior de \(S_x\).

Demonstração. Os retângulos inferiores satisfazem

\[ S_{x+h_n\mathbf1}\searrow S_x, \qquad S_{x-h_n\mathbf1}\nearrow S_x^\circ. \]

Aplique-se a continuidade da probabilidade às suas imagens inversas por \(\mathbf X\).

Continuidade de uma FD conjunta

Proposição 2.3

Sejam \(\mathbf X\) um vetor aleatório em \(\mathbb R^d\) e \(x\in\mathbb R^d\). Para \(S_x=\prod_{j=1}^d(-\infty,x_j]\), com fronteira \(\partial S_x=S_x\setminus S_x^\circ\),

\[ F_{\mathbf X}\text{ é contínua em }x \quad\Longleftrightarrow\quad \mathbb P_{\mathbf X}(\partial S_x)=0. \]

Demonstração. Se a fronteira tem probabilidade nula, os limites por cima e por baixo do slide anterior são ambos iguais a \(F_{\mathbf X}(x)\). Para \(\|y-x\|_\infty<h\),

\[ F_{\mathbf X}(x-h\mathbf1) \le F_{\mathbf X}(y) \le F_{\mathbf X}(x+h\mathbf1). \]

Para qualquer \(\varepsilon>0\), escolha-se \(h>0\) de modo que ambos os extremos difiram de \(F_{\mathbf X}(x)\) por menos de \(\varepsilon\). Então \(|F_{\mathbf X}(y)-F_{\mathbf X}(x)|<\varepsilon\) sempre que \(\|y-x\|_\infty<h\), o que prova a continuidade em \(x\).

Reciprocamente, a continuidade dá \(\mathbb P_{\mathbf X}(S_x^\circ)=\mathbb P_{\mathbf X}(S_x)\), pelo que a fronteira tem probabilidade nula.

Distribuições marginais

Definição 2.4

Para um vetor aleatório \(\mathbf X=(X_1,\ldots,X_d)\), as distribuições \(\mathbb P_{X_j}\) das suas coordenadas chamam-se distribuições marginais.

Denotando a projeção na coordenada por \(\pi_j(x)=x_j\),

\[ \mathbb P_{X_j}(B)=\mathbb P_{\mathbf X}\!\left(\pi_j^{-1}(B)\right), \qquad B\in\mathcal B(\mathbb R). \]

Esta identidade decorre de \(X_j=\pi_j\circ\mathbf X\).

Para um par \((X,Y)\), temos \((X,Y)^{-1}((-\infty,x]\times(-\infty,n])\nearrow X^{-1}((-\infty,x])\), pelo que

\[ F_X(x)=\lim_{y\to+\infty}F_{X,Y}(x,y). \]

Analogamente, \(F_Y(y)=\lim_{x\to+\infty}F_{X,Y}(x,y)\).

Projeção para os dois dados equilibrados. Trinta e seis pontos equiprováveis formam uma grelha de seis por seis. A projeção pi_1 envia cada coluna para a primeira coordenada, somando seis massas de 1/36 numa massa marginal de 1/6. A coluna de primeira coordenada 3 e a sua imagem estão realçadas.

Probabilidade condicionada

Definição 2.5

Seja \((\Omega,\mathcal F,\mathbb P)\) um espaço de probabilidade. Para \(B\in\mathcal F\) com \(\mathbb P(B)>0\), a probabilidade condicionada por \(B\) define-se por

\[ \mathbb P(A\mid B):=\frac{\mathbb P(A\cap B)}{\mathbb P(B)}, \qquad A\in\mathcal F. \]

Verificação. Trata-se de uma medida de probabilidade: é não negativa, \(\mathbb P(\Omega\mid B)=1\) e, para \(A_n\) disjuntos,

\[ \mathbb P\!\left(\bigcup_n A_n\,\middle|\,B\right) =\frac{\sum_n\mathbb P(A_n\cap B)}{\mathbb P(B)} =\sum_n\mathbb P(A_n\mid B). \]

Em particular, \(\mathbb P(B\mid B)=1\) e \(\mathbb P(A\cap B)=\mathbb P(A\mid B)\mathbb P(B)\).

Acontecimentos independentes

Definição 2.6

Os acontecimentos \(A,B\in\mathcal F\) chamam-se independentes se

\[ \mathbb P(A\cap B)=\mathbb P(A)\mathbb P(B). \]

Para \(\mathbb P(B)>0\), isto equivale a \(\mathbb P(A\mid B)=\mathbb P(A)\). No entanto, a definição aplica-se também quando \(\mathbb P(B)=0\) ou \(1\).

Dezasseis resultados equiprováveis formam uma grelha de quatro linhas e quatro colunas. O acontecimento A contém as oito células azuis das duas linhas superiores. O acontecimento B contém as oito células das duas colunas da esquerda. Ao condicionar por B, restam oito resultados, quatro dos quais pertencem a A. Assim, ℙ(A)=8/16=1/2 e ℙ(A condicionado por B)=4/8=1/2, ilustrando a independência.

Acontecimentos disjuntos de probabilidade positiva não são independentes:

\[ A\cap B=\varnothing,\quad \mathbb P(A)\mathbb P(B)>0 \quad\Longrightarrow\quad \mathbb P(A\cap B)=0\ne \mathbb P(A)\mathbb P(B). \]

Um critério de independência

Proposição 2.7

Para acontecimentos \(A,B\) com \(0<\mathbb P(B)<1\), \(A\) e \(B\) são independentes se e só se

\[ \mathbb P(A\mid B)=\mathbb P(A\mid B^c). \]

Demonstração. Sejam \(\alpha=\mathbb P(A\mid B)\) e \(\beta=\mathbb P(A\mid B^c)\).

Como \(A=(A\cap B)\mathbin{\dot\cup}(A\cap B^c)\),

\[ \mathbb P(A)=\alpha \mathbb P(B)+\beta \mathbb P(B^c). \]

Se \(\alpha=\beta\), obtemos \(\mathbb P(A)=\alpha\) e, portanto, \(\mathbb P(A\cap B)=\mathbb P(A)\mathbb P(B)\).

Reciprocamente, se \(A\) e \(B\) são independentes, então

\[ \mathbb P(A\cap B^c)=\mathbb P(A)-\mathbb P(A\cap B)=\mathbb P(A)\mathbb P(B^c), \]

pelo que \(\alpha=\beta=\mathbb P(A)\). O mesmo cálculo mostra que a independência se preserva por passagem ao complementar.

Um dado: independência e dependência

Exemplo 2.8

Em \(\Omega=\{1,\ldots,6\}\), com resultados equiprováveis, sejam

\[ A=\{2,4,6\},\qquad B=\{1,2\},\qquad C=\{1,2,3\}. \]

Par Probabilidade da interseção Produto das probabilidades Independentes?
\(A,B\) \(1/6\) \((1/2)(1/3)=1/6\) Sim
\(A,C\) \(1/6\) \((1/2)(1/2)=1/4\) Não

Duas representações dos seis resultados equiprováveis de um dado. Os resultados pares de A estão a azul. Contornos dourados selecionam B={1,2} à esquerda e C={1,2,3} à direita. Ao condicionar por B, um dos dois resultados restantes é par, pelo que ℙ(A condicionado por B)=1/2=ℙ(A). Ao condicionar por C, um dos três resultados restantes é par, pelo que ℙ(A condicionado por C)=1/3, diferente de ℙ(A)=1/2.

A dependência na segunda linha exprime-se também por

\[ \mathbb P(A\mid C)=\frac13\ne\frac23=\mathbb P(A\mid C^c). \]

Independência de intervalos

Exemplo 2.9

Em \([0,1]\), com probabilidade uniforme, sejam \(I=[0,1/3]\) e \(J=[a,b]\), onde \(0\le a<b\le1\).

Se \(J\subseteq I\) ou \(a\ge1/3\), a igualdade \(\mathbb P(I\cap J)=\mathbb P(I)\mathbb P(J)\) é impossível. Assim, a independência exige \(a<1/3<b\) e

\[ \mathbb P(I\cap J)=\frac13-a,\qquad \mathbb P(I)\mathbb P(J)=\frac13(b-a). \]

\(I\) e \(J\) são independentes se e só se

\[ 0\le a<\frac13,\qquad b=1-2a. \]

Se \(a=b\), \(J\) tem probabilidade nula e é independente de qualquer acontecimento.

A \(\sigma\)-álgebra gerada por variáveis aleatórias

Definição 2.10

Para uma classe \(\mathcal C\) de subconjuntos de \(\Omega\), \(\sigma(\mathcal C)\) denota a menor \(\sigma\)-álgebra que contém \(\mathcal C\).

Para variáveis aleatórias reais \(X_1,X_2,\ldots\), define-se

\[ \sigma(X_1,X_2,\ldots) =\sigma\!\left(\{X_i^{-1}(B):B\in\mathcal B(\mathbb R),\ i\ge1\}\right). \]

É a menor \(\sigma\)-álgebra que torna cada \(X_i\) mensurável.

Para uma só variável, \(\sigma(X)=\{X^{-1}(B):B\in\mathcal B(\mathbb R)\}\), pois as imagens inversas já formam uma \(\sigma\)-álgebra.

Boa definição. Interseções arbitrárias de \(\sigma\)-álgebras em \(\Omega\) são \(\sigma\)-álgebras. A interseção de todas as que contêm \(\mathcal C\) é, portanto, a menor. Esta família é não vazia, pois contém \(\mathcal P(\Omega)\).

A \(\sigma\)-álgebra gerada por uma família finita

Proposição 2.11

Seja \(\mathcal C=\{A_1,\ldots,A_m\}\subseteq\mathcal F\), com \(m\ge1\). Sejam \(A_i^{(1)}=A_i\) e \(A_i^{(0)}=A_i^c\).

Para \(\varepsilon\in\{0,1\}^m\), defina-se \(C_\varepsilon=\bigcap_{i=1}^m A_i^{(\varepsilon_i)}\). Então

\[ \sigma(\mathcal C) =\left\{\bigcup_{\varepsilon\in J}C_\varepsilon:J\subseteq\{0,1\}^m\right\}, \]

onde a união para \(J=\varnothing\) é \(\varnothing\).

Demonstração. Os conjuntos não vazios \(C_\varepsilon\) formam uma partição de \(\Omega\). As suas uniões formam uma \(\sigma\)-álgebra \(\mathcal G\): os complementares selecionam os restantes conjuntos, e as uniões numeráveis são ainda uniões destes conjuntos.

Como \(A_i=\bigcup_{\varepsilon:\varepsilon_i=1}C_\varepsilon\), temos \(\sigma(\mathcal C)\subseteq\mathcal G\). Reciprocamente, cada \(C_\varepsilon\in\sigma(\mathcal C)\), pelo que qualquer união destes conjuntos, em número finito, pertence a \(\sigma(\mathcal C)\).

Exemplo: acontecimentos determinados por uma observação

Exemplo 2.12

Em \(\Omega=[0,1]\), sejam \(X_1(\omega)=\omega\) e

\[ X_2(\omega)= \begin{cases} 1,&0\le\omega<1/2,\\ 2,&1/2\le\omega\le1. \end{cases} \]

Então

\[ \sigma(X_1)=\mathcal B([0,1]), \]

\[ \sigma(X_2)= \{\varnothing,\ [0,1/2),\ [1/2,1],\ [0,1]\}. \]

Toda a imagem inversa por \(X_2\) é uma união dos seus dois conjuntos de nível. Ocorrem as quatro uniões possíveis.

Independência de classes e sucessões

Sejam \(\mathcal C_i\subseteq\mathcal F\) classes de acontecimentos, para \(i=1,2,\ldots\).

Definição 2.13

As classes chamam-se independentes se, para todo o conjunto finito não vazio \(I\) de índices distintos e toda a escolha \(A_i\in\mathcal C_i\),

\[ \mathbb P\!\left(\bigcap_{i\in I}A_i\right)=\prod_{i\in I}\mathbb P(A_i). \]

As variáveis aleatórias \(X_1,X_2,\ldots\) chamam-se independentes se as \(\sigma\)-álgebras \(\sigma(X_1),\sigma(X_2),\ldots\) são independentes.

Independência e FD conjuntas

Teorema 2.14

Para um vetor aleatório real \(\mathbf X=(X_1,\ldots,X_n)\), são equivalentes as condições seguintes.

(i) \(X_1,\ldots,X_n\) são independentes.

(ii) Para todo o \(\mathbf x=(x_1,\ldots,x_n)\in\mathbb R^n\),

\[ F_{\mathbf X}(x_1,\ldots,x_n)=F_{X_1}(x_1)\cdots F_{X_n}(x_n). \]

(iii) Para todos os borelianos \(B_1,\ldots,B_n\subseteq\mathbb R\),

\[ \mathbb P_{\mathbf X}(B_1\times\cdots\times B_n)=\mathbb P_{X_1}(B_1)\cdots \mathbb P_{X_n}(B_n). \]

Demonstração. (i) ⇒ (ii). Tomem-se os acontecimentos \(X_i^{-1}(({-\infty},x_i])\).

(ii) ⇒ (iii). A igualdade vale para \(B_i=(-\infty,x_i]\). Fixados os outros conjuntos, ambos os membros são medidas finitas em \(B_i\). A unicidade estende a igualdade a todos os borelianos, uma coordenada de cada vez.

Dos retângulos à independência

Demonstração do Teorema 2.14, (iii) ⇒ (i).

Seja \(I=\{i_1,\ldots,i_k\}\subseteq\{1,\ldots,n\}\) não vazio, e tomem-se \(A_i\in\sigma(X_i)\) para \(i\in I\).

Pela Definição 2.10, escreva-se \(A_i=X_i^{-1}(B_i)\) com \(B_i\in\mathcal B(\mathbb R)\). Tome-se \(B_i=\mathbb R\) para \(i\notin I\).

A condição (iii) e \(\mathbb P_{X_i}(\mathbb R)=1\) dão

\[ \begin{aligned} \mathbb P\!\left(\bigcap_{i\in I}A_i\right) &=\mathbb P\!\left(\mathbf X^{-1}(B_1\times\cdots\times B_n)\right)\\ &=\mathbb P_{\mathbf X}(B_1\times\cdots\times B_n)\\ &=\mathbb P_{X_1}(B_1)\cdots \mathbb P_{X_n}(B_n)\\ &=\mathbb P(A_{i_1})\cdots \mathbb P(A_{i_k}). \end{aligned} \]

Isto vale para todo o \(I\) e qualquer escolha dos acontecimentos \(A_i\). Logo, as \(\sigma\)-álgebras \(\sigma(X_1),\ldots,\sigma(X_n)\) são independentes, como se pretendia.

Independência dois a dois e mútua

Exemplo 2.15

Lancem-se duas moedas equilibradas e independentes. Sejam \(X,Y\in\{H,T\}\) os resultados, com \(H\) para cara e \(T\) para coroa. Defina-se \(Z=1\) se \(X=Y\), e \(Z=0\) caso contrário.

Para \(a\in\{H,T\}\) e \(b\in\{0,1\}\), exatamente um dos quatro resultados equiprováveis de \((X,Y)\) dá \((X,Z)=(a,b)\), pelo que

\[ \mathbb P(X=a,Z=b)=\frac14=\mathbb P(X=a)\mathbb P(Z=b). \]

O mesmo cálculo vale para \((Y,Z)\), pelo que todos os pares são independentes. No entanto,

\[ \mathbb P(X=H,Y=H,Z=1)=\frac14\ne\frac18 =\mathbb P(X=H)\mathbb P(Y=H)\mathbb P(Z=1). \]

Funções de probabilidade conjuntas

Definição 2.16

Sejam \(X\) e \(Y\) concentradas em conjuntos no máximo numeráveis \(D_X\) e \(D_Y\). Então \((X,Y)\) está concentrado em \(D=D_X\times D_Y\).

A sua função de probabilidade conjunta define-se por

\[ p_{X,Y}(x,y)=\mathbb P(X=x,Y=y). \]

Proposição 2.17

Se \((X,Y)\) está concentrado num conjunto no máximo numerável \(D\subseteq\mathbb R^2\), então, para todo o boreliano \(A\subseteq\mathbb R^2\),

\[ \mathbb P((X,Y)\in A)=\sum_{(x,y)\in A\cap D}p_{X,Y}(x,y), \qquad \sum_{(x,y)\in D}p_{X,Y}(x,y)=1. \]

Demonstração. Aplique-se a aditividade numerável aos acontecimentos disjuntos \((X,Y)^{-1}(\{(x,y)\})\), \((x,y)\in D\), cuja união tem probabilidade um.

Funções de probabilidade marginais

Proposição 2.18

Sejam \(X,Y\) variáveis concentradas em conjuntos no máximo numeráveis \(D_X,D_Y\). Para todo o \(x\in D_X\) e todo o \(y\in D_Y\),

\[ p_X(x)=\sum_{v\in D_Y}p_{X,Y}(x,v), \qquad p_Y(y)=\sum_{u\in D_X}p_{X,Y}(u,y). \]

Demonstração. Fixado \(x\), os acontecimentos \(X^{-1}(\{x\})\cap Y^{-1}(\{v\})\), \(v\in D_Y\), são disjuntos. A sua união é \(X^{-1}(\{x\})\cap Y^{-1}(D_Y)\). Como \(\mathbb P(Y^{-1}(D_Y))=1\), a aditividade numerável dá

\[ p_X(x)=\mathbb P\bigl(X^{-1}(\{x\})\bigr) =\sum_{v\in D_Y}\mathbb P\bigl(X^{-1}(\{x\})\cap Y^{-1}(\{v\})\bigr). \]

Trocando \(X\) e \(Y\), obtém-se a segunda identidade. As marginais são as somas das linhas e colunas da tabela da função de probabilidade conjunta.

Funções de probabilidade condicionadas

Definição 2.19

Seja \((X,Y)\) discreto com função de probabilidade conjunta \(p_{X,Y}\). Para \(p_X(x)>0\), a função de probabilidade de \(Y\) condicionada por \(X=x\) define-se por

\[ p_{Y\mid X=x}(y) :=\mathbb P(Y=y\mid X=x) :=\frac{p_{X,Y}(x,y)}{p_X(x)}. \]

A identidade das marginais prova a normalização:

\[ \sum_{y\in\operatorname{supp}(Y)} p_{Y\mid X=x}(y) =\frac{\sum_{y\in\operatorname{supp}(Y)} p_{X,Y}(x,y)}{p_X(x)}=1. \]

Para \(p_Y(y)>0\), a função de probabilidade condicionada no sentido inverso é

\[ p_{X\mid Y=y}(x):=\frac{p_{X,Y}(x,y)}{p_Y(y)}. \]

Independência e fatorização da função de probabilidade

Proposição 2.20

Sejam \(X,Y\) variáveis aleatórias concentradas em conjuntos no máximo numeráveis \(D_X,D_Y\). São independentes se e só se

\[ p_{X,Y}(x,y)=p_X(x)p_Y(y)\quad\text{para quaisquer }x,y. \]

Demonstração. A independência dá a identidade para os acontecimentos \(X^{-1}(\{x\})\) e \(Y^{-1}(\{y\})\).

Reciprocamente, somando uma função de probabilidade fatorizada num retângulo boreliano \(A\times B\), obtém-se

\[ \begin{aligned} \mathbb P(X\in A,Y\in B) &=\sum_{x\in A\cap D_X}\sum_{y\in B\cap D_Y}p_X(x)p_Y(y)\\ &=\mathbb P(X\in A)\mathbb P(Y\in B). \end{aligned} \]

Função de probabilidade: marginais, condicionamento e independência

Exemplo 2.21
\(Y=0\) \(Y=1\) \(p_X(x)\)
\(X=0\) 0.10 0.20 0.30
\(X=1\) 0.30 0.40 0.70
\(p_Y(y)\) 0.40 0.60 1

\[ \mathbb P(Y=1\mid X=0)=\frac{0.20}{0.30}=\frac23. \]

As variáveis não são independentes porque

\[ p_{X,Y}(0,0)=0.10\ne(0.30)(0.40)=p_X(0)p_Y(0). \]

Funções de densidade conjuntas

Definição 2.22

Um par \((X,Y)\) chama-se absolutamente contínuo se existe uma função mensurável \(f_{X,Y}:\mathbb R^2\to[0,\infty)\) tal que, para todo o boreliano \(A\subseteq\mathbb R^2\),

\[ \mathbb P((X,Y)\in A)=\iint_A f_{X,Y}(x,y)\,dx\,dy. \]

A função \(f_{X,Y}\) chama-se função de densidade de probabilidade conjunta.

Tomando \(A=\mathbb R^2\), obtém-se \(\iint_{\mathbb R^2}f_{X,Y}(x,y)\,dx\,dy=1\).

Exemplo 2.23 · Marginais absolutamente contínuas sem densidade conjunta

Sejam \(X\sim U(0,1)\) e \(Y=X\). Então \(X\) e \(Y\) têm distribuições absolutamente contínuas, mas \((X,Y)\) não admite densidade conjunta.

De facto, a diagonal \(D=\{(x,y)\in\mathbb R^2:x=y\}\) tem medida de Lebesgue bidimensional nula, mas \(\mathbb P((X,Y)\in D)=1\).

Funções de densidade marginais

Proposição 2.24

Se \((X,Y)\) tem densidade conjunta \(f_{X,Y}\), então \(X\) e \(Y\) têm densidades marginais dadas, quase em toda a parte, por

\[ f_X(x)=\int_{\mathbb R}f_{X,Y}(x,y)\,dy, \qquad f_Y(y)=\int_{\mathbb R}f_{X,Y}(x,y)\,dx. \]

Demonstração. Para todo o boreliano \(A\subseteq\mathbb R\), o teorema de Tonelli dá

\[ \mathbb P(X\in A) =\int_A\left[\int_{\mathbb R}f_{X,Y}(x,y)\,dy\right]dx. \]

O integral interior é mensurável e não negativo, com integral total igual a um. É, portanto, finito quase em toda a parte e define uma densidade de \(X\). Trocando as coordenadas, prova-se a afirmação para \(Y\).

Funções de densidade condicionadas

Definição 2.25

Seja \((X,Y)\) um par com densidade conjunta \(f_{X,Y}\) e funções marginais \(f_X,f_Y\) dadas pela Proposição 2.24.

Para \(0<f_X(x)<\infty\), a função de densidade de probabilidade de \(Y\) condicionada por \(X=x\) define-se por

\[ f_{Y\mid X=x}(y):=\frac{f_{X,Y}(x,y)}{f_X(x)}. \]

Analogamente, \(f_{X\mid Y=y}(x):=f_{X,Y}(x,y)/f_Y(y)\) para \(0<f_Y(y)<\infty\).

Verificação. Para estes valores de \(x\), \(f_{Y\mid X=x}\) é mensurável e não negativa, e

\[ \int_{\mathbb R}f_{Y\mid X=x}(y)\,dy =\frac{\int_{\mathbb R}f_{X,Y}(x,y)\,dy}{f_X(x)}=1. \]

Distribuições condicionadas

Proposição 2.26

Seja \((X,Y)\) um par com densidade conjunta \(f_{X,Y}\). Defina-se \(f_{Y\mid X=x}\) como na Definição 2.25, usando uma densidade de probabilidade fixa quando \(f_X(x)\in\{0,\infty\}\).

Para quaisquer conjuntos borelianos \(A,B\subseteq\mathbb R\),

\[ \mathbb P(X\in A,Y\in B) =\int_A\left[\int_Bf_{Y\mid X=x}(y)\,dy\right]d\mathbb P_X(x). \]

Demonstração. Seja \(D=\{x\in\mathbb R:0<f_X(x)<\infty\}\). Como \(\int_{\mathbb R}f_X(x)\,dx=1\), \(\mathbb P_X(D)=1\).

Em \(D\), \(f_{Y\mid X=x}(y)f_X(x)=f_{X,Y}(x,y)\). Assim, por Tonelli,

\[ \begin{aligned} \int_A\left[\int_Bf_{Y\mid X=x}(y)\,dy\right]d\mathbb P_X(x) &=\int_{A\cap D}\int_Bf_{X,Y}(x,y)\,dy\,dx\\ &=\mathbb P(X\in A,Y\in B). \end{aligned} \]

Independência e fatorização da densidade

Proposição 2.27

Seja \((X,Y)\) um par com densidade conjunta \(f_{X,Y}\) e densidades marginais \(f_X,f_Y\). Então \(X\) e \(Y\) são independentes se e só se

\[ f_{X,Y}(x,y)=f_X(x)f_Y(y) \quad\text{para quase todo o }(x,y)\text{ (Lebesgue)}. \]

Demonstração. Seja \(g(x,y)=f_X(x)f_Y(y)\). Para conjuntos borelianos \(A,B\subseteq\mathbb R\), Tonelli dá

\[ \iint_{A\times B}g(x,y)\,dx\,dy=\mathbb P_X(A)\mathbb P_Y(B). \]

Tomando \(A=B=\mathbb R\), conclui-se que \(g\) é uma densidade de probabilidade.

(⇐) Se \(f_{X,Y}=g\) quase em toda a parte, então \(\mathbb P(X\in A,Y\in B)=\mathbb P_X(A)\mathbb P_Y(B)\).

(⇒) Se \(X\) e \(Y\) são independentes, as medidas com densidades \(f_{X,Y}\) e \(g\) coincidem nos retângulos borelianos. Esta classe é fechada para interseções finitas e gera \(\mathcal B(\mathbb R^2)\). Pela unicidade das medidas de probabilidade e das densidades, \(f_{X,Y}=g\) quase em toda a parte.

Exemplo: uma densidade num triângulo

Exemplo 2.28

Sejam \(D=\{(x,y)\in\mathbb R^2:0<y<x<1\}\) e \(f_{X,Y}(x,y)=2\cdot\mathbf1_D(x,y)\). Então

\[ \int_0^1\int_0^x2\,dy\,dx=\int_0^12x\,dx=1. \]

As densidades marginais são

\[ f_X(x)=2x\cdot\mathbf1_{(0,1)}(x), \qquad f_Y(y)=2(1-y)\cdot\mathbf1_{(0,1)}(y). \]

Para \(0<x<1\),

\[ f_{Y\mid X=x}(y)=\frac1x\cdot\mathbf1_{(0,x)}(y). \]

Assim, \(Y\mid X=x\) tem distribuição uniforme em \((0,x)\). As variáveis não são independentes.

Fórmula de Bayes

Preparação facultativa para a Ficha 1

Proposição 2.29

Para acontecimentos \(A,B\) num espaço de probabilidade, com \(\mathbb P(A)>0\) e \(\mathbb P(B)>0\),

\[ \mathbb P(B\mid A)=\frac{\mathbb P(A\mid B)\mathbb P(B)}{\mathbb P(A)}. \]

Demonstração. Pela definição de probabilidade condicionada,

\[ \mathbb P(A\cap B)=\mathbb P(A\mid B)\mathbb P(B)=\mathbb P(B\mid A)\mathbb P(A). \]

Divida-se por \(\mathbb P(A)\).

Ficha 1:

Exercício 2.30 · Ficha 1, questão 5

Exercício facultativo

Sejam \(H\) o acontecimento de ter hipertensão e \(D\) o de consumir álcool regularmente. Suponha que

\[ \mathbb P(H)=0.05,\qquad \mathbb P(D\mid H)=0.75,\qquad \mathbb P(D\mid H^c)=0.50. \]

  1. Determine a proporção de pessoas que consomem álcool, \(\mathbb P(D)\).
  2. Determine a proporção de consumidores de álcool que têm hipertensão, \(\mathbb P(H\mid D)\).

Valor esperado, famílias usuais e transformadas

03 · SEMANA 2 · AULA 3

Valor esperado, variância, covariância. Distribuições binomial, de Poisson, multinomial, normal, exponencial e gama. Funções características.

Valor esperado

Definição 3.1

Uma variável aleatória real \(X\) chama-se integrável se \(\int_\Omega|X|\,d\mathbb P<\infty\). Se \(X\) for integrável, o seu valor esperado define-se por

\[ \mathbb E(X)=\int_\Omega X(\omega)\,d\mathbb P(\omega) =\int_{\mathbb R}x\,d\mathbb P_X(x). \]

Observação: A condição \(\int_\Omega|X|\,d\mathbb P<\infty\) é necessária por razões de teoria da medida. Se \(X\ge0\) (isto é, se \(X\) tomar apenas valores não negativos), o integral também está definido sem esta hipótese, mas pode tomar um valor em \([0,\infty]\).

Valor esperado de funções de uma variável aleatória

Teorema 3.2 · "Transferência" (LOTUS, do inglês “Law of the unconscious statistician”)

Sejam \(X\) uma variável aleatória real e \(g:\mathbb R\to\mathbb R\) uma aplicação mensurável. Se \(\mathbb E(|g(X)|)<\infty\), então

\[ \mathbb E(g(X))=\int_{\mathbb R}g(x)\,d\mathbb P_X(x). \]

Se \(X\) for discreta ou absolutamente contínua, o integral pode escrever-se, respetivamente, usando a função de probabilidade ou a função de densidade de probabilidade:

\[ \mathbb E(g(X)) =\sum_{x\in\operatorname{supp}(X)} g(x)p_X(x) \quad\text{ou}\quad \mathbb E(g(X)) =\int_{\mathbb R}g(x)f_X(x)\,dx. \]

Propriedades do valor esperado

Proposição 3.3

Sejam \(X,Y,X_1,\ldots,X_m\) variáveis aleatórias integráveis em \((\Omega,\mathcal F,\mathbb P)\), \(a_1,\ldots,a_m\in\mathbb R\) e \(A\in\mathcal F\). Então o valor esperado \(\mathbb E\) tem as seguintes propriedades.

\[ \begin{aligned} &\text{Linearidade:}&& \mathbb E\!\left(\sum_{j=1}^m a_jX_j\right)=\sum_{j=1}^m a_j\mathbb E(X_j),\\[3mm] &\text{Monotonia:}&& X\le Y\ \text{q.c.}\ \Longrightarrow\ \mathbb E(X)\le \mathbb E(Y),\\[3mm] &\text{Indicadoras:}&& \mathbb E(\mathbf1_A)=\mathbb P(A). \end{aligned} \]

Dizemos que uma propriedade se verifica quase certamente (q.c.), se se verifica com probabilidade um. Precisamente, dizemos que \(X\le Y\) se verifica quase certamente, se \(\mathbb P (\{\omega\in \Omega \mid X(\omega) \le Y(\omega) \}) = 1\). Analogamente, dizemos que \(X = Y\) se verifica quase certamente, se \(\mathbb P (\{\omega\in \Omega \mid X(\omega) = Y(\omega) \}) = 1\). Denota-se por \([X]\) a classe de equivalência de todas as variáveis aleatórias \(Y\) iguais a \(X\) quase certamente.

Momentos e espaços \(L^p\)

Definição 3.4

Para \(1\le p<\infty\), sejam

\[ L^p(\mathbb P)=\{[X]:\mathbb E(|X|^p)<\infty\}, \qquad \|[X]\|_p=(\mathbb E(|X|^p))^{1/p}. \]

Aqui, \(X\) percorre as variáveis aleatórias reais em \((\Omega,\mathcal F,\mathbb P)\). Escreve-se \([X]\) para denotar a classe de equivalência definida no slide anterior, identificando duas variáveis aleatórias quando são iguais quase certamente. Esta identificação é importante porque \(\|[X]\|_p\) é uma norma nas classes de variáveis aleatórias, mas apenas uma seminorma nas variáveis aleatórias.

Para um inteiro \(k\ge1\) e \(X \in L^k(\mathbb P)\), o momento de ordem \(k\) e o momento centrado de ordem \(k\) definem-se, respetivamente, por

\[ \mathbb E(X^k) \qquad \text{e}\qquad \mathbb E((X-\mathbb E(X))^k). \]

Comparação de normas \(L^p\)

Proposição 3.5

Num espaço de probabilidade, se \(p>q\ge1\) e \(X\in L^p(\mathbb P)\), então \(X\in L^q(\mathbb P)\) e \(\|X\|_q\le\|X\|_p\).

Demonstração. A função \(t\mapsto t^{q/p}\) é côncava em \([0,\infty)\). A desigualdade de Jensen aplicada a \(|X|^p\) dá \(\mathbb E(|X|^q)\le(\mathbb E(|X|^p))^{q/p}<\infty\). Tomando raízes de ordem \(q\), obtém-se a conclusão.

Proposição 3.6 · Desigualdade de Jensen

Seja \(Z\) integrável com valores num intervalo \(I\), e seja \(g:I\to\mathbb R\) convexa. Então

\[ g(\mathbb E(Z))\le\mathbb E(g(Z)). \]

O segundo membro pode ser \(+\infty\). Para \(g\) côncava, a desigualdade inverte-se, admitindo \(-\infty\) no segundo membro.

Variância

Definição 3.7

Para \(X\in L^2(\mathbb P)\), a variância define-se por

\[ \mathbb V(X)=\mathbb E((X-\mathbb E(X))^2). \]

Proposição 3.8

Para \(X\in L^2(\mathbb P)\) e \(a,b\in\mathbb R\),

\[ \mathbb V(X)=\mathbb E(X^2)-(\mathbb E(X))^2,\qquad \mathbb V(aX+b)=a^2\mathbb V(X). \]

Além disso, \(\mathbb V(X)=0\) se e só se \(X=\mathbb E(X)\) quase certamente.

Demonstração. Expanda-se o quadrado e aplique-se a linearidade para obter a primeira identidade. A segunda decorre de \(aX+b-\mathbb E(aX+b)=a(X-\mathbb E(X))\).

Por fim, a variável não negativa \((X-\mathbb E(X))^2\) tem valor esperado nulo se e só se é nula quase certamente.

Ficha 1: média e variância

Exercício 3.9 · Ficha 1, questão 8

Seja \(X\) uma variável com a seguinte distribuição:

\[ \mathbb P(X=-2)=\frac13,\qquad \mathbb P(X=3)=\frac12,\qquad \mathbb P(X=1)=\frac16. \]

  1. Calcule \(\mathbb E(X)\) e \(\mathbb E(X^2)\).
  2. Calcule \(\mathbb V(X)\).
  3. Determine \(\mathbb E(2X+5)\), \(\mathbb V(2X)\) e \(\mathbb V(2X+5)\).

Ficha 1: variância da distribuição uniforme

Exercício 3.10 · Ficha 1, questão 11, continuação

Seja \(X\sim U(-2,2)\).

  1. Calcule \(\mathbb E(X)\), \(\mathbb E(X^2)\) e \(\mathbb V(X)\).
  2. Determine o desvio padrão de \(X\).
  3. Determine \(\mathbb E(2X+5)\) e \(\mathbb V(3X+4)\).

Covariância

Definição 3.11

Para \(X,Y\in L^2(\mathbb P)\), a sua covariância define-se por

\[ \operatorname{Cov}(X,Y) =\mathbb E((X-\mathbb E(X))(Y-\mathbb E(Y))). \]

As variáveis aleatórias \(X, Y\) chamam-se não correlacionadas se \(\operatorname{Cov}(X,Y)=0\).

Proposição 3.12

Para \(X,Y\in L^2(\mathbb P)\) e \(a,b,c,d\in\mathbb R\),

\[ \begin{aligned} \operatorname{Cov}(X,Y)&=\mathbb E(XY)-\mathbb E(X)\mathbb E(Y),\\ \operatorname{Cov}(aX+b,cY+d)&=ac\,\operatorname{Cov}(X,Y),\\ |\operatorname{Cov}(X,Y)|&\le\sqrt{\mathbb V(X)\mathbb V(Y)}. \end{aligned} \]

Demonstração. Expanda-se o produto, use-se a linearidade e aplique-se Cauchy–Schwarz às variáveis centradas.

Covariância e produtos internos

Proposição 3.13

No espaço vetorial real \(L^2(\mathbb P)\), define-se um produto interno por

\[ \langle X,Y\rangle_2:=\int_\Omega X(\omega)Y(\omega)\,d\mathbb P(\omega). \]

No subespaço \(L^2_0(\mathbb P):=\{X\in L^2(\mathbb P):\mathbb E(X)=0\}\), a covariância coincide com este produto interno:

\[ \operatorname{Cov}(X,Y)=\langle X,Y\rangle_2\qquad(X,Y\in L^2_0(\mathbb P)). \]

Demonstração. Por Cauchy–Schwarz, \(\mathbb E(|XY|)<\infty\). A simetria e a bilinearidade decorrem do integral. Além disso, \(\langle X,X\rangle_2=\mathbb E(X^2)\ge0\), com igualdade exatamente quando \(X=0\) quase certamente, isto é, \(X=0\) em \(L^2(\mathbb P)\).

Para variáveis centradas, a fórmula da covariância é \(\mathbb E(XY)\). Em todo o espaço \(L^2(\mathbb P)\), a covariância não é um produto interno: \(\operatorname{Cov}(1,1)=0\), embora \(1\ne0\).

Identidade de Hoeffding para a covariância

Teorema 3.14

Sejam \(X,Y\in L^2(\mathbb P)\), com FD conjunta \(F_{X,Y}\) e FD marginais \(F_X,F_Y\). Então

\[ \operatorname{Cov}(X,Y) =\int_{\mathbb R^2} \bigl(F_{X,Y}(x,y)-F_X(x)F_Y(y)\bigr)\,dx\,dy. \]

Para variáveis independentes, a função integranda é nula em todos os pontos. A recíproca não se conclui: um integral nulo pode resultar da compensação entre valores positivos e negativos da função integranda.

Variância de uma soma

Proposição 3.15

Para \(X_1,\ldots,X_n\in L^2(\mathbb P)\),

\[ \mathbb V\!\left(\sum_{j=1}^nX_j\right) =\sum_{j=1}^n\mathbb V(X_j) +2\sum_{1\le i<j\le n}\operatorname{Cov}(X_i,X_j). \] Se as variáveis são não correlacionadas duas a duas, então

\[ \mathbb V\!\left(\sum_{j=1}^nX_j\right)=\sum_{j=1}^n\mathbb V(X_j). \]

Demonstração. Mostrámos em Proposição 3.12 que a covariância é bilinear. Esta propriedade permite expandir \(\mathbb E\!\left(\left(\sum_j(X_j-\mathbb E(X_j))\right)^2\right)\). Os termos cruzados são o dobro das covariâncias entre pares.

Valor esperado de um produto

Proposição 3.16

Se \(X\) e \(Y\) são independentes e integráveis, então \(XY\) é integrável e

\[ \mathbb E(XY)=\mathbb E(X)\mathbb E(Y). \]

Demonstração. Por Teorema 2.14, a distribuição conjunta é o produto das marginais. O teorema de Tonelli dá

\[ \begin{aligned} \mathbb E(|XY|) &=\int_{\mathbb R}\int_{\mathbb R}|xy|\,d\mathbb P_Y(y)\,d\mathbb P_X(x)\\ &=\mathbb E(|X|)\mathbb E(|Y|)<\infty. \end{aligned} \]

O teorema de Fubini aplica-se, portanto, ao produto com sinal:

\[ \mathbb E(XY) =\int_{\mathbb R}x\left(\int_{\mathbb R}y\,d\mathbb P_Y(y)\right)d\mathbb P_X(x) =\mathbb E(X)\mathbb E(Y). \]

Não correlação e independência

Proposição 3.17

Sejam \(X,Y\in L^2(\mathbb P)\).

(i) \(X\) e \(Y\) são não correlacionadas se e só se \(\operatorname{Cov}(f(X),g(Y))=0\) para quaisquer funções afins \(f,g:\mathbb R\to\mathbb R\).

(ii) \(X\) e \(Y\) são independentes se e só se \(\operatorname{Cov}(f(X),g(Y))=0\) para quaisquer funções Borel-mensuráveis \(f,g:\mathbb R\to\mathbb R\) tais que \(f(X),g(Y)\in L^2(\mathbb P)\).

Em cada condição, o integral que se exige ser nulo é

\[ \int_\Omega \bigl(f(X)-\mathbb E(f(X))\bigr) \bigl(g(Y)-\mathbb E(g(Y))\bigr)\,d\mathbb P. \]

A restrição de integrabilidade em (ii) pode ser substituída pela condição de \(f,g\) serem limitadas.

Demonstração da comparação

Demonstração de Proposição 3.17.

(i) Se \(f(x)=ax+b\) e \(g(y)=cy+d\), por Proposição 3.12,

\[ \operatorname{Cov}(f(X),g(Y))=ac\,\operatorname{Cov}(X,Y). \]

Esta expressão é nula para quaisquer \(a,b,c,d\) exatamente quando \(\operatorname{Cov}(X,Y)=0\), tomando \(a=c=1\) e \(b=d=0\).

(ii), implicação direta. A mensurabilidade dá \(\sigma(f(X))\subseteq\sigma(X)\) e \(\sigma(g(Y))\subseteq\sigma(Y)\), pelo que \(f(X)\) e \(g(Y)\) são independentes. Por Proposição 3.16, a sua covariância é nula.

Implicação recíproca. Para quaisquer conjuntos de Borel \(A,B\subseteq\mathbb R\), escolham-se as funções limitadas e mensuráveis \(f=\mathbf1_A\) e \(g=\mathbf1_B\). Então

\[ 0=\operatorname{Cov}(\mathbf1_A(X),\mathbf1_B(Y)) =\mathbb P(X\in A,Y\in B)-\mathbb P(X\in A)\mathbb P(Y\in B). \]

Por Teorema 2.14, \(X\) e \(Y\) são independentes.

Variáveis não correlacionadas, mas dependentes

Exemplo 3.18

Sejam \(\mathbb P(X=-1)=\mathbb P(X=0)=\mathbb P(X=1)=1/3\) e \(Y=X^2\). Então

\[ \mathbb E(X)=\mathbb E(XY)=\mathbb E(X^3)=0, \]

pelo que \(\operatorname{Cov}(X,Y)=0\). No entanto, \(Y\) é uma função de \(X\) e

\[ \mathbb P(X=0,Y=0)=\frac13\ne\frac19=\mathbb P(X=0)\mathbb P(Y=0). \]

Momentos por integrais de cauda

Proposição 3.19

Para uma variável aleatória não negativa \(X\) e um real \(r>0\), com igualdade em \([0,\infty]\),

\[ \mathbb E(X^r) =r\int_0^\infty t^{r-1}\mathbb P(X>t)\,dt. \]

Demonstração. Tonelli permite a troca abaixo porque a função integranda é não negativa.

\[ \begin{aligned} \mathbb E(X^r) &=\int_\Omega X(\omega)^r\,d\mathbb P(\omega) =\int_\Omega\int_0^{X(\omega)}r t^{r-1}\,dt\,d\mathbb P(\omega)\\ &=\int_0^\infty\int_{\{\omega\in\Omega\mid X(\omega)>t\}}r t^{r-1}\,d\mathbb P(\omega)\,dt\\ &=r\int_0^\infty\mathbb P(\{\omega\in\Omega\mid X(\omega)>t\})t^{r-1}\,dt =r\int_0^\infty t^{r-1}\mathbb P(X>t)\,dt. \end{aligned} \]

Fórmulas de cauda para a média

Proposição 3.20

Para uma variável aleatória não negativa \(X\),

\[ \mathbb E(X)=\int_0^\infty \mathbb P(X>t)\,dt =\int_0^\infty(1-F_X(t))\,dt. \]

Proposição 3.21 · Fórmula da soma das probabilidades de cauda

Para uma variável aleatória \(X\) com valores em \(\{0,1,2,\ldots\}\),

\[ \mathbb E(X)=\sum_{k=1}^\infty \mathbb P(X\ge k) =\sum_{k=0}^\infty \mathbb P(X>k). \]

Ambas as igualdades admitem o valor \(+\infty\). O integral decorre do caso \(r=1\) em Proposição 3.19. A soma obtém-se tomando valores esperados em \(X=\sum_{k\ge1}\mathbf1_{[k,\infty)}(X)\).

Somas de cauda: momentos superiores

Proposição 3.22

Seja \(X\) uma variável com valores em \(\{0,1,2,\ldots\}\) e seja \(r\ge1\) um inteiro. Então, com igualdade em \([0,\infty]\),

\[ \mathbb E(X^r) =\sum_{k=0}^{\infty} \bigl((k+1)^r-k^r\bigr)\mathbb P(X>k). \]

Demonstração. Para \(t\in[k,k+1)\), \(\mathbb P(X>t)=\mathbb P(X>k)\). Decompondo o integral não negativo em Proposição 3.19, obtém-se

\[ \begin{aligned} \mathbb E(X^r) &=\sum_{k=0}^\infty\mathbb P(X>k)\int_k^{k+1}r t^{r-1}\,dt\\ &=\sum_{k=0}^\infty\bigl((k+1)^r-k^r\bigr)\mathbb P(X>k). \end{aligned} \]

Valor esperado como integral de quantis

Proposição 3.23

Sejam \(X\) uma variável aleatória real integrável e \(F_X^{-1}\) a sua inversa generalizada. Então

\[ \mathbb E(X)=\int_0^1F_X^{-1}(u)\,du. \]

Demonstração. Seja \(Q=F_X^{-1}\) e seja \(U\) uniforme em \((0,1)\). Por Proposição 1.25, \(Q\) é mensurável e \(Q(u)\le x\) se e só se \(u\le F_X(x)\). Logo, para todo o \(x\in\mathbb R\),

\[ \mathbb P(Q(U)\le x)=\mathbb P(U\le F_X(x))=F_X(x). \]

Assim, \(Q(U)\) e \(X\) têm a mesma distribuição. Em particular, \(\mathbb E(|Q(U)|)=\mathbb E(|X|)<\infty\). Aplicando Teorema 3.2 a \(U\) e \(Q\), obtém-se

\[ \mathbb E(X)=\mathbb E(Q(U))=\int_0^1Q(u)\,du. \]

A distribuição binomial

Definição 3.24

Sejam \(n\ge1\) um inteiro e \(0\le p\le1\). Diz-se que uma variável \(X\) tem distribuição binomial, e escreve-se \(X\sim B(n,p)\), se

\[ \mathbb P(X=k)=\binom nkp^k(1-p)^{n-k}, \qquad k=0,\ldots,n. \]

Usamos a convenção \(0^0=1\) para incluir \(p=0\) e \(p=1\). A distribuição \(B(1,p)\) chama-se Bernoulli\((p)\).

Proposição 3.25

Sejam \(X_1,\ldots,X_n\) variáveis independentes de Bernoulli\((p)\) e seja \(S_n=X_1+\cdots+X_n\). Então \(S_n\sim B(n,p)\) e

\[ \mathbb E(S_n)=np,\qquad \mathbb V(S_n)=np(1-p). \]

Reciprocamente, qualquer variável \(B(n,p)\) tem a mesma distribuição que \(S_n\).

Contagem dos resultados de Bernoulli

Demonstração de Proposição 3.25. Suponha-se primeiro que \(0<p<1\).

Cada sequência \((\varepsilon_1,\ldots,\varepsilon_n)\in\{0,1\}^n\) com exatamente \(k\) uns tem probabilidade \(p^k(1-p)^{n-k}\), por independência.

Há \(\binom nk\) escolhas para as posições dos uns. Os acontecimentos correspondentes são disjuntos, pelo que

\[ \mathbb P(S_n=k)=\binom nkp^k(1-p)^{n-k}. \]

Para \(p=0\) ou \(p=1\), \(S_n=np\) quase certamente.

Quatro passos independentes de Bernoulli formam uma rede de Galton. Os passos para a esquerda têm probabilidade um menos p, e os passos para a direita, probabilidade p. Os números de caminhos nos cinco pontos finais são 1, 4, 6, 4, 1. Destaca-se um dos seis caminhos com dois passos para a direita.

Os números contam caminhos. Os seis caminhos até \(k=2\) têm probabilidade \(p^2(1-p)^2\).

Como \(X_j^2=X_j\), \(\mathbb E(X_j)=p\) e \(\mathbb V(X_j)=p(1-p)\). Pela linearidade e por Proposição 3.15,

\[ \mathbb E(S_n)=np, \qquad \mathbb V(S_n)=np(1-p). \]

Ficha 2: uma probabilidade binomial

Exercício 3.26 · Ficha 2, questão 1

Cada uma de 25 pessoas selecionadas independentemente tem probabilidade \(0.10\) de ter problemas renais. Seja \(X\) o número de casos na amostra.

  1. Identifique a distribuição de \(X\) e os seus parâmetros.
  2. Determine a probabilidade de haver, no máximo, cinco casos.
  3. Determine a probabilidade de haver entre seis e nove casos, inclusive.

A distribuição trinomial

Exemplo 3.27

Suponha-se que lançamos um dado viciado de 3 faces 10 vezes independentemente, ou seja, \(n=3\), \(d=10\). As faces 1,2,3 têm probabilidades \(\mathbf p=(0.2,0.3,0.5)\). O vetor \(\mathbf X=(X_1,X_2,X_3)\) regista as contagens de cada face, pelo que \(X_1+X_2+X_3=10\).

Uma grelha triangular com as 66 contagens possíveis em dez lançamentos de um dado de três faces, com probabilidades 0.2, 0.3 e 0.5. A área dos círculos representa a probabilidade. Destaca-se a contagem 2,3,5.

Cada ponto representa uma contagem possível das três faces. A área do círculo é proporcional à sua probabilidade.

Para o ponto destacado,

\[ \begin{gathered} \mathbb P(\mathbf X=(2,3,5))\\ =\frac{10!}{2!\,3!\,5!}\,0.2^2\,0.3^3\,0.5^5\\ =0.08505. \end{gathered} \]

A distribuição multinomial

Definição 3.28

Sejam \(n,d\ge1\) inteiros e \(\mathbf p=(p_1,\ldots,p_n)\), com \(p_j\ge0\) e \(p_1+\cdots+p_n=1\).

Um vetor aleatório \(\mathbf X=(X_1,\ldots,X_n)\) tem distribuição multinomial, e escreve-se \(\mathbf X\sim\operatorname{Multinomial}(d,\mathbf p)\), se

\[ \mathbb P(\mathbf X=\mathbf x) =\frac{d!}{x_1!\cdots x_n!}\,p_1^{x_1}\cdots p_n^{x_n} \]

para \(\mathbf x\in\mathbb N_0^n\) com \(x_1+\cdots+x_n=d\), sendo a probabilidade nula nos restantes casos.

Interpretação por contagem. Lance-se um dado de \(n\) faces \(d\) vezes, independentemente. Denote-se por \(C_\ell\) o resultado do lançamento \(\ell\), com \(\mathbb P(C_\ell=j)=p_j\) para cada face \(j\). Seja \(X_j\) o número de ocorrências da face \(j\).

Há exatamente \(d!/(x_1!\cdots x_n!)\) sequências com contagens \(\mathbf x\). Cada uma tem probabilidade \(p_1^{x_1}\cdots p_n^{x_n}\), obtendo-se a fórmula apresentada.

Momentos da distribuição multinomial

Proposição 3.29

Se \(\mathbf X=(X_1,\ldots,X_n)\sim\operatorname{Multinomial}(d,\mathbf p)\), então, para cada \(j\),

\[ X_j\sim B(d,p_j), \qquad \mathbb E(X_j)=dp_j, \qquad \mathbb V(X_j)=dp_j(1-p_j). \]

Para \(i\ne j\), \(\operatorname{Cov}(X_i,X_j)=-dp_ip_j\).

Demonstração. Em cada um dos \(d\) lançamentos independentes, a face \(j\) tem probabilidade \(p_j\). Logo, \(X_j\sim B(d,p_j)\) por Proposição 3.25.

Para \(i\ne j\), \(X_iX_j\) conta pares ordenados de lançamentos distintos com faces \(i\) e \(j\). Pela linearidade e independência,

\[ \begin{aligned} \mathbb E(X_iX_j) &=\sum_{\substack{1\le\ell,m\le d\\\ell\ne m}}\mathbb P(C_\ell=i,C_m=j) =d(d-1)p_ip_j,\\ \operatorname{Cov}(X_i,X_j) &=d(d-1)p_ip_j-(dp_i)(dp_j)=-dp_ip_j. \end{aligned} \]

A distribuição de Poisson

Definição 3.30

Para \(\lambda>0\), definimos a distribuição de Poisson com taxa \(\lambda\), denotada por \(\operatorname{Poisson}(\lambda)\), da seguinte forma.

\[ X\sim\operatorname{Poisson}(\lambda) \quad\Longleftrightarrow\quad \mathbb P(X=k)=e^{-\lambda}\frac{\lambda^k}{k!}, \quad k=0,1,\ldots \]

Uma variável aleatória com distribuição de Poisson toma valores em \(\mathbb N_0\).

A média e a variância são ambas iguais a \(\lambda\):

\[ \mathbb E(X)=\mathbb V(X)=\lambda. \]

Se \(X_j\sim\operatorname{Poisson}(\lambda_j)\) são independentes, então

\[ \sum_{j=1}^nX_j\sim\operatorname{Poisson}\!\left(\sum_{j=1}^n\lambda_j\right). \]

A distribuição exponencial

Definição 3.31

Para \(\lambda>0\), definimos a distribuição exponencial com taxa \(\lambda\), denotada por \(\operatorname{Exp}(\lambda)\), da seguinte forma.

\[ X\sim\operatorname{Exp}(\lambda) \quad\Longleftrightarrow\quad f_X(x)=\lambda e^{-\lambda x}\mathbf1_{[0,\infty)}(x). \]

Uma variável aleatória com distribuição exponencial toma valores em \([0,\infty)\).

Por integração, obtém-se

\[ \mathbb P(X>x)=e^{-\lambda x}\quad(x\ge0), \qquad \mathbb E(X)=\frac1\lambda, \qquad \mathbb V(X)=\frac1{\lambda^2}. \]

No exemplo de chegadas de chamadas que se segue, \(1/\lambda\) é o tempo médio de espera.

Chegadas de chamadas: contagem

Exemplo 3.32 · Modelo de um centro de atendimento

Sejam \(\lambda>0\) e \(n>\lambda\) inteiro. Nos instantes \(j/n\) minutos (\(j=1,2,\ldots\)), há uma chamada com probabilidade \(\lambda/n\), independentemente. Seja \(N_n(t)\) o número de chamadas até \(t\). Então \(N_n(1)\sim B(n,\lambda/n)\).

Probabilidades do número de chamadas para 6, 30 e 1000 intervalos curtos por minuto, com uma taxa média de três chamadas por minuto. As barras binomiais em azul-petróleo aproximam-se das probabilidades fixas da Poisson(3), assinaladas por pontos em coral. Os painéis usam os mesmos eixos.

Para \(k\in\mathbb N_0\) fixo e \(n\ge k\), quando \(n\to\infty\),

\[ \mathbb P(N_n(1)=k) =\frac{\lambda^k}{k!} \underbrace{\frac{n!}{(n-k)!\,n^k}}_{\longrightarrow 1} \underbrace{\left(1-\frac{\lambda}{n}\right)^{n-k}}_{\longrightarrow e^{-\lambda}} \longrightarrow e^{-\lambda}\frac{\lambda^k}{k!}. \]

Gráficos: \(\lambda=3\) chamadas/minuto. Barras: probabilidades binomiais. Pontos: probabilidades da \(\operatorname{Poisson}(3)\).

Chegadas de chamadas: tempo de espera

Exemplo 3.33 · Modelo de um centro de atendimento, continuação

No Exemplo 3.32, seja \(X_n\) o tempo de espera, em minutos, pela primeira chamada. Até ao instante \(t\ge0\), realizam-se \(\lfloor nt\rfloor\) provas.

Não há chamadas entre o instante zero e t. A primeira chamada chega no instante X_n, à direita de t. Assim, X_n maior do que t e N_n(t) igual a zero descrevem o mesmo acontecimento.

Esperar mais do que \(t\) significa que nenhuma destas provas produziu uma chamada:

\[ \begin{aligned} \mathbb P(X_n>t) &=\mathbb P(N_n(t)=0)=\left(1-\frac{\lambda}{n}\right)^{\lfloor nt\rfloor}\\ &=\left[\left(1-\frac{\lambda}{n}\right)^n\right]^{\lfloor nt\rfloor/n} \xrightarrow[n\to\infty]{} e^{-\lambda t}. \end{aligned} \]

Esta é a probabilidade de sobrevivência da \(\operatorname{Exp}(\lambda)\). Para \(\lambda=3\) chamadas/minuto, a espera média é de \(1/3\) minuto, ou \(20\) segundos.

Ausência de memória

Proposição 3.34

Se \(X\sim\operatorname{Exp}(\lambda)\) com \(\lambda>0\), então, para quaisquer \(s,t\ge0\),

\[ \mathbb P(X>s+t\mid X>s) =\mathbb P(X>t). \]

Equivalentemente, a distribuição de \(X-s\) condicionada por \(X>s\) é \(\operatorname{Exp}(\lambda)\).

Demonstração. Como \(\mathbb P(X>s)=e^{-\lambda s}>0\),

\[ \mathbb P(X>s+t\mid X>s) =\frac{\mathbb P(X>s+t)}{\mathbb P(X>s)} =\frac{e^{-\lambda(s+t)}}{e^{-\lambda s}} =e^{-\lambda t}. \]

Esta é a probabilidade de sobrevivência \(\mathbb P(X>t)\) para todo o \(t\ge0\), pelo que determina a mesma distribuição.

Ficha 2: uma probabilidade de Poisson

Exercício 3.35 · Ficha 2, questão 3

O número anual de novos casos segue uma distribuição de Poisson com média \(13\).

Determine a probabilidade de, num dado ano, haver:

  1. Exatamente dez novos casos.
  2. Pelo menos oito novos casos.
  3. Entre nove e quinze novos casos, inclusive.

A distribuição normal

Definição 3.36

Para \(\mu\in\mathbb R\) e \(\sigma>0\), \(X\) tem distribuição normal, também chamada distribuição gaussiana, e escreve-se \(X\sim N(\mu,\sigma^2)\), se a sua densidade de probabilidade é

\[ f_X(x)=\frac{1}{\sigma\sqrt{2\pi}} \exp\!\left(-\frac12\left(\frac{x-\mu}{\sigma}\right)^2\right). \]

Densidades normais N(0,1), N(2,1) e N(0,4). Alterar a média desloca a curva. Uma variância maior dá uma curva mais larga e mais baixa.

O parâmetro \(\mu\) é a média e \(\sigma\) é o desvio padrão:

\[ \mathbb E(X)=\mu,\qquad \mathbb V(X)=\sigma^2,\qquad \frac{X-\mu}{\sigma}\sim N(0,1). \]

Da binomial à normal

Distribuições binomiais estandardizadas para 4, 16 e 64 provas, com p igual a um meio. Os histogramas aproximam-se da densidade normal reduzida, representada por uma curva suave em coral. Os eixos têm as mesmas escalas nos três painéis.
Teorema 3.37 · De Moivre-Laplace

Fixe-se \(0<p<1\). Se \(X_d\sim B(d,p)\), então, para todo o \(z\in\mathbb R\),

\[ \lim_{d\to\infty}\mathbb P\!\left(\frac{X_d-dp}{\sqrt{dp(1-p)}}\le z\right) =\int_{-\infty}^z\frac{e^{-t^2/2}}{\sqrt{2\pi}}\,dt. \]

Barras: \(B(d,1/2)\) estandardizada, com altura \(\mathbb P(X_d=k)/(2/\sqrt d)\). A área de cada barra é a sua probabilidade. Curva: a densidade \(N(0,1)\). Este é um caso particular do teorema do limite central, que será enunciado e demonstrado na semana 3.

Ficha 2: estandardização

Exercício 3.38 · Ficha 2, questão 4

Seja \(X\) o peso de uma pessoa, em quilogramas. Suponha que \(X\sim N(60,12^2)\).

  1. Determine \(\mathbb P(X=60)\) e \(\mathbb P(X>70)\).
  2. Determine \(\mathbb P(55<X<65)\).
  3. Determine o quantil de ordem \(0.75\) de \(X\).

A distribuição gama

Definição 3.39

Para forma \(\kappa>0\) e escala \(\theta>0\),

\[ X\sim\Gamma(\kappa,\theta) \quad\Longleftrightarrow\quad f_X(x)= \frac{x^{\kappa-1}e^{-x/\theta}} {\Gamma(\kappa)\theta^\kappa}\mathbf1_{(0,\infty)}(x). \]

\[ \mathbb E(X)=\kappa\theta,\qquad \mathbb V(X)=\kappa\theta^2. \]

Se \(\kappa=m\in\mathbb N\), então \(\Gamma(m,\theta)\) é a distribuição de uma soma de \(m\) variáveis \(\operatorname{Exp}(1/\theta)\) independentes. Na parametrização forma–taxa, \(\lambda=1/\theta\).

Funções geradoras de momentos

Definição 3.40

Para uma variável aleatória real \(X\), a sua função geradora de momentos (FGM) define-se por

\[ M_X(t)=\mathbb E(e^{tX}) \]

para os valores de \(t\) em que o valor esperado é finito.

Teorema 3.41

Seja \(X\) uma variável real e suponha-se \(M_X\) finita num intervalo aberto \(I\) que contém \(0\). Então todos os momentos absolutos de \(X\) são finitos e

\[ M_X^{(k)}(0)=\mathbb E(X^k),\qquad k=1,2,\ldots \]

Se uma variável aleatória real \(Y\) satisfaz \(M_Y(t)=M_X(t)\) para todo o \(t\in I\), então \(\mathbb P_Y=\mathbb P_X\).

A distribuição de Cauchy, por exemplo, não tem FGM finita para nenhum \(t\ne0\).

Funções características

Definição 3.42

A função característica de uma variável aleatória real \(X\) define-se por

\[ \varphi_X(t)=\mathbb E(e^{itX}) =\int_{\mathbb R}e^{itx}\,d\mathbb P_X(x), \qquad t\in\mathbb R. \]

Aqui, \(i^2=-1\). Como \(|e^{itX}|=1\), o valor esperado existe para todo o real \(t\). Toda a função característica satisfaz

\[ \varphi_X(0)=1,\qquad |\varphi_X(t)|\le1,\qquad \varphi_X(-t)=\overline{\varphi_X(t)}, \]

e é uniformemente contínua em \(\mathbb R\).

Aplicações afins e somas independentes

Proposição 3.43

Para uma variável aleatória real \(X\) e \(a,b\in\mathbb R\),

\[ \varphi_{aX+b}(t)=e^{ibt}\varphi_X(at)\qquad(t\in\mathbb R). \]

Para variáveis aleatórias reais independentes \(X_1,\ldots,X_n\),

\[ \varphi_{\sum_{j=1}^nX_j}(t) =\prod_{j=1}^n\varphi_{X_j}(t)\qquad(t\in\mathbb R). \]

Demonstração. Substitua-se \(aX+b\) na definição e coloque-se \(e^{ibt}\) em evidência. Para a soma, a independência dá

\[ \mathbb E\!\left(e^{it\sum_jX_j}\right) =\mathbb E\!\left(\prod_je^{itX_j}\right) =\prod_j\mathbb E(e^{itX_j}). \]

As identidades para as FGM valem onde todas as FGM envolvidas são finitas.

Derivadas de funções características

Proposição 3.44

Sejam \(k\ge1\) um inteiro e \(X\) uma variável aleatória real com \(\mathbb E(|X|^k)<\infty\). Então \(\varphi_X\) é \(k\) vezes continuamente diferenciável e, para \(1\le j\le k\),

\[ \varphi_X^{(j)}(t) =\mathbb E\!\left((iX)^je^{itX}\right), \qquad t\in\mathbb R. \]

Em particular,

\[ \varphi_X^{(j)}(0)=i^j\mathbb E(X^j). \]

Demonstração. Como \(|(ix)^je^{itx}|=|x|^j\le1+|x|^k\), a diferenciação dominada sob o sinal de integral é válida até à ordem \(k\). O teorema da convergência dominada de Lebesgue mostra também que as derivadas são contínuas. Tome-se \(t=0\) para obter a última identidade.

Funções características de distribuições usuais

Distribuição \(\varphi_X(t)\)
\(B(n,p)\) \((1-p+pe^{it})^n\)
\(\operatorname{Poisson}(\lambda)\) \(\exp(\lambda(e^{it}-1))\)
\(U(a,b)\) \(\dfrac{e^{itb}-e^{ita}}{it(b-a)}\) para \(t\ne0\). Valor \(1\) em \(t=0\)
\(N(\mu,\sigma^2)\) \(\exp(i\mu t-\tfrac12\sigma^2t^2)\)
\(\operatorname{Exp}(\lambda)\) \(\dfrac{\lambda}{\lambda-it}\)
\(\Gamma(\kappa,\theta)\) \((1-i\theta t)^{-\kappa}\)

Unicidade das funções características

Definição 3.45

Diz-se que as variáveis aleatórias reais \(X\) e \(Y\) têm a mesma distribuição, e escreve-se \(X\overset{d}=Y\), se \(\mathbb P_X=\mathbb P_Y\).

Teorema 3.46

Para variáveis aleatórias reais \(X\) e \(Y\),

\[ X\overset{d}=Y \quad\Longleftrightarrow\quad \varphi_X(t)=\varphi_Y(t) \quad\text{para todo o }t\in\mathbb R. \]

Em conjunto com a fórmula do produto, a unicidade identifica distribuições de somas sem calcular integrais de convolução.

Tipos de convergência

04 · SEMANA 2 · AULA 4

Tipos de convergência: visão geral

Existem vários sentidos distintos em que se pode dizer que uma sucessão de variáveis aleatórias \(X_n\) converge para uma variável aleatória \(X\). Nesta aula, vamos estudar os tipos mais importantes de convergência de variáveis aleatórias.

Para variáveis aleatórias reais e \(1\le p<\infty\):

Tipo Notação Quando \(n\to\infty\) Mesmo espaço de probabilidade?
Quase certa \(X_n\xrightarrow{\mathrm{q.c.}}X\) \(X_n(\omega)\to X(\omega)\) fora de um conjunto de probabilidade nula Sim
\(L^\infty\) \(X_n\xrightarrow{L^\infty}X\) \(\|X_n-X\|_\infty\to0\) Sim
Em probabilidade \(X_n\overset{\mathbb P}{\rightharpoonup}X\) \(\mathbb P(|X_n-X|>\varepsilon)\to0\) para todo o \(\varepsilon>0\) Sim
\(L^p\) \(X_n\xrightarrow{L^p}X\) \(\mathbb E(|X_n-X|^p)\to0\) Sim
Em distribuição \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\) \(F_{X_n}(x)\to F_X(x)\) em todo o ponto de continuidade \(x\) de \(F_X\) Não

Por isso, não devemos dizer apenas que \(X_n\) converge para \(X\), pois isso seria ambíguo. É importante especificar o tipo de convergência.

Convergência quase certa

Definição 4.1

Para variáveis aleatórias reais \(X,X_1,X_2,\ldots\) num mesmo espaço de probabilidade, dizemos que \(X_n\) convergem quase certamente para \(X\), e escrevemos \(X_n\xrightarrow{\mathrm{q.c.}}X\), se

\[ \mathbb P\!\left(\{\omega\in\Omega:X_n(\omega)\to X(\omega)\}\right)=1. \]

Definição 4.2

Para acontecimentos \(A_1,A_2,\ldots\), o seu limite superior define-se por

\[ \limsup_{n\to\infty} A_n=\bigcap_{m=1}^{\infty}\bigcup_{n\ge m}A_n. \]

É o acontecimento de ocorrerem infinitos acontecimentos da sucessão \(A_n\).

Convergência quase certa: um critério

Proposição 4.3

Sejam \(X,X_1,X_2,\ldots\) variáveis aleatórias reais num mesmo espaço de probabilidade. Então \(X_n\xrightarrow{\mathrm{q.c.}}X\) se e só se

\[ \mathbb P\!\left(\limsup_{n\to\infty}\bigl(|X_n-X|\bigr)^{-1}((\varepsilon,\infty))\right)=0 \quad\text{para todo o }\varepsilon>0. \]

Demonstração de “\(\implies\)”. Suponha-se que \(X_n\xrightarrow{\mathrm{q.c.}}X\). Fora de um conjunto \(A\) de probabilidade zero, \(X_n(\omega)\) converge (pontualmente) para \(X(\omega)\). Para cada \(\varepsilon>0\) e \(\omega \in \Omega\setminus A\), só pode haver um número finito de \(n\ge1\) tais que \(|X_n(\omega)-X(\omega)| \ge \varepsilon\). Assim, \(\omega \in \Omega\setminus A\) não pode pertencer ao limite superior dos conjuntos \(\bigl(|X_n-X|\bigr)^{-1}((\varepsilon,\infty))\). Portanto, o limite superior está contido em \(A\) e tem necessariamente probabilidade zero.

Convergência quase certa: recíproca

Proposição 4.3 · Recíproca

Suponha-se que, para todo o \(\varepsilon>0\),

\[ \mathbb P\!\left(\limsup_{n\to\infty}\bigl(|X_n-X|\bigr)^{-1}((\varepsilon,\infty))\right)=0. \]

Então \(X_n\xrightarrow{\mathrm{q.c.}}X\).

Demonstração. Pela condição assumida e pela subaditividade numerável,

\[ A:=\bigcup_{m=1}^{\infty}\limsup_{n\to\infty}\bigl(|X_n-X|\bigr)^{-1}((1/m,\infty)) \qquad\text{satisfaz }\mathbb P(A)=0. \]

Para \(\omega\notin A\) e cada \(m\), vale \(|X_n(\omega)-X(\omega)|\le1/m\) para todos os \(n\) suficientemente grandes. Logo, \(X_n(\omega)\) converge pontualmente para \(X(\omega)\) para \(\omega\notin A\), o que prova \(X_n\xrightarrow{\mathrm{q.c.}}X\).

O espaço \(L^\infty(\mathbb P)\)

Definição 4.4

Para uma variável aleatória real \(X\) em \((\Omega,\mathcal F,\mathbb P)\), define-se

\[ \|X\|_\infty :=\operatorname*{ess\,sup}|X| :=\inf\{a\in[0,\infty):\mathbb P(|X|>a)=0\}, \]

com \(\inf\varnothing:=\infty\). A variável \(X\) chama-se essencialmente limitada se \(\|X\|_\infty<\infty\).

O espaço das variáveis aleatórias essencialmente limitadas denota-se por

\[ L^\infty(\mathbb P):=\{X:\|X\|_\infty<\infty\}. \]

Como em \(L^p(\mathbb P)\), identificam-se variáveis iguais quase certamente. Nestas classes de equivalência, \(\|\cdot\|_\infty\) chama-se norma do supremo essencial.

Convergência em \(L^\infty\)

Definição 4.5

Para \(X,X_1,X_2,\ldots\in L^\infty(\mathbb P)\), dizemos que \(X_n\) converge em \(L^\infty\) para \(X\), e escrevemos \(X_n\xrightarrow{L^\infty}X\), se

\[ \|X_n-X\|_\infty\to0. \]

Proposição 4.6

Para \(X,X_1,X_2,\ldots\in L^\infty(\mathbb P)\),

\[ X_n\xrightarrow{L^\infty}X \quad\Longrightarrow\quad X_n\xrightarrow{\mathrm{q.c.}}X. \]

Demonstração. Para cada \(n\), escolha-se um conjunto \(A_n\) de probabilidade zero fora do qual \(|X_n-X|\le\|X_n-X\|_\infty\). A união \(A=\bigcup_{n\ge1}A_n\) tem probabilidade zero.

Para \(\omega\notin A\), tem-se \(|X_n(\omega)-X(\omega)|\le\|X_n-X\|_\infty\to0\), o que prova \(X_n\xrightarrow{\mathrm{q.c.}}X\).

Convergência em probabilidade

Definição 4.7

Para variáveis aleatórias reais \(X,X_1,X_2,\ldots\) num mesmo espaço de probabilidade, dizemos que \(X_n\) converge em probabilidade para \(X\), e escrevemos \(X_n\overset{\mathbb P}{\rightharpoonup}X\), se

\[ \mathbb P(|X_n-X|>\varepsilon)\to0 \quad\text{para todo o }\varepsilon>0. \]

Proposição 4.8

Se as variáveis aleatórias reais \(X_n\) convergem para \(X\) quase certamente num mesmo espaço de probabilidade, então \(X_n\overset{\mathbb P}{\rightharpoonup}X\).

Demonstração. Fixado \(\varepsilon>0\), as indicadoras \(\mathbf1_{(\varepsilon,\infty)}(|X_n-X|)\) convergem para zero quase certamente e são limitadas por \(1\). Aplicando o teorema da convergência dominada de Lebesgue, obtemos

\[ \mathbb P(|X_n-X|>\varepsilon) =\mathbb E\!\left(\mathbf1_{(\varepsilon,\infty)}(|X_n-X|)\right)\to0. \]

Convergência em \(L^p\)

Definição 4.9

Sejam \(1\le p<\infty\) e \(X,X_1,X_2,\ldots\in L^p(\mathbb P)\). Dizemos que \(X_n\) converge em \(L^p\) para \(X\), e escrevemos \(X_n\xrightarrow{L^p}X\), se

\[ \|X_n-X\|_p =(\mathbb E(|X_n-X|^p))^{1/p}\to0. \]

Se \(q>p\ge1\), então, num espaço de probabilidade,

\[ \|Y\|_p\le\|Y\|_q. \]

Consequentemente,

\[ \text{convergência em }L^\infty\Rightarrow\text{convergência em }L^q \Rightarrow\text{convergência em }L^p. \]

Desigualdade de Markov

Proposição 4.10 · Desigualdade de Markov

Seja \(Z\) uma variável aleatória não negativa e seja \(f:[0,\infty)\to[0,\infty)\) não decrescente, com \(\mathbb E(f(Z))<\infty\). Para todo o \(\varepsilon>0\) tal que \(f(\varepsilon)>0\),

\[ \mathbb P(Z\ge\varepsilon)\le\frac{\mathbb E(f(Z))}{f(\varepsilon)}. \]

Em particular, isto aplica-se a \(f(t) = t\) e \(f(t) = t^2\).

Demonstração. A desigualdade pontual \(f(\varepsilon)\mathbf1_{[\varepsilon,\infty)}(Z)\le f(Z)\) e a monotonia do valor esperado dão

\[ f(\varepsilon)\,\mathbb P(Z\ge\varepsilon) =\mathbb E(f(\varepsilon)\mathbf1_{[\varepsilon,\infty)}(Z)) \le\mathbb E(f(Z)). \]

Agora, divide-se por \(f(\varepsilon)>0\).

Para \(f(t)=t\), obtém-se \(\mathbb P(Z\ge\varepsilon)\le\mathbb E(Z)/\varepsilon\).

Desigualdade de Chebyshev

Proposição 4.11 · Desigualdade de Chebyshev

Seja \(X\in L^2(\mathbb P)\). Para todo o \(\varepsilon>0\),

\[ \mathbb P\!\left(|X-\mathbb E(X)|\ge\varepsilon\right) \le\frac{\mathbb V(X)}{\varepsilon^2}. \]

Demonstração. Aplique-se a Proposição 4.10 a \(Z=|X-\mathbb E(X)|\), com \(f(t)=t^2\):

\[ \mathbb P\!\left(|X-\mathbb E(X)|\ge\varepsilon\right) \le\frac{\mathbb E(|X-\mathbb E(X)|^2)}{\varepsilon^2} =\frac{\mathbb V(X)}{\varepsilon^2}. \]

Convergência em \(L^p\) e em probabilidade

Proposição 4.12

Sejam \(1\le p<\infty\) e \(X,X_1,X_2,\ldots\in L^p(\mathbb P)\). Então

\[ X_n\xrightarrow{L^p}X \quad\Longrightarrow\quad X_n\overset{\mathbb P}{\rightharpoonup}X. \]

Demonstração. Para todo o \(\varepsilon>0\), aplique-se a Proposição 4.10 com \(Z=|X_n-X|\) e \(f(t)=t^p\):

\[ \begin{aligned} \mathbb P(|X_n-X|>\varepsilon) &\le\frac{\mathbb E(f(|X_n-X|))}{f(\varepsilon)}\\ &=\frac{\mathbb E(|X_n-X|^p)}{\varepsilon^p}\longrightarrow0. \end{aligned} \]

Sob convergência em \(L^p\), o numerador tende para zero, o que prova a convergência em probabilidade.

A sucessão da máquina de escrever

Exemplo 4.13

Considere-se \(\Omega=[0,1)\), \(\mathcal F=\mathcal B([0,1))\) e seja \(\mathbb P\) a medida de Lebesgue. Para inteiros \(m\ge0\) e \(k=0,\ldots,2^m-1\), defina-se

\[ I_{m,k}:=[k2^{-m},(k+1)2^{-m}), \qquad X_{m,k}:=\mathbf1_{I_{m,k}}. \]

Enumerem-se todas as variáveis linha a linha, da esquerda para a direita, para definir a sucessão

\[ (Y_n)_{n\ge1}:=(X_{0,0},X_{1,0},X_{1,1},X_{2,0},\ldots). \]

As linhas m=0,1,2,3 dividem o intervalo de zero a um num intervalo, duas metades, quatro quartos e oito oitavos. Leia-se cada linha da esquerda para a direita. Cada ponto pertence a exatamente um intervalo em cada linha.

Máquina de escrever: convergência em \(L_p\) sem convergência quase certa

Proposição 4.14

A sucessão \((Y_n)_{n\ge1}\) do Exemplo 4.13 converge para \(0\) em \(L^p\) para todo o \(1\le p<\infty\) e em probabilidade, mas não converge quase certamente.

Demonstração. Se \(Y_n=X_{m,k}\), então \(\|Y_n\|_p^p=\int_\Omega\mathbf1_{I_{m,k}}\,d\mathbb P=2^{-m}\).

Cada linha é finita, pelo que \(m\to\infty\) quando \(n\to\infty\). Logo, \(\|Y_n\|_p\to0\). A convergência em probabilidade decorre da Proposição 4.12.

Fixe-se \(\omega\in\Omega\). Em cada linha \(m\ge1\), exatamente uma indicadora vale \(1\) em \(\omega\), e as restantes valem \(0\). Sejam \(n_m(\omega)\) e \(r_m(\omega)\) os índices em \((Y_n)\) deste \(1\) e do primeiro \(0\) nessa linha. Ambas as sucessões de índices são estritamente crescentes, e

\[ Y_{n_m(\omega)}(\omega)=1,\qquad Y_{r_m(\omega)}(\omega)=0 \qquad(m\ge1). \]

Logo, \((Y_n(\omega))\) não converge para nenhum \(\omega\in\Omega\).

Convergência quase certa sem convergência em \(L^1\)

Exemplo 4.15

Seja \(U\sim U(0,1)\) e

\[ X_n=n\,\mathbf1_{(0,1/n]}(U). \]

Para todo o \(\omega\) com \(U(\omega)>0\), tem-se \(X_n(\omega)=0\) para todos os \(n\) suficientemente grandes. Como \(\mathbb P(U>0)=1\),

\[ X_n\xrightarrow{\mathrm{q.c.}}0. \]

Mas

\[ \mathbb E(|X_n|)=n\mathbb P(0<U\le1/n)=1, \]

pelo que \(X_n\) não converge para \(0\) em \(L^1\).

Gráficos de u para n vezes a indicadora de (0,1/n], para n=2,4,8, com os mesmos eixos. As alturas aumentam e as larguras diminuem. Cada retângulo sombreado tem área um.

Convergência dominada em \(L^p\)

Teorema 4.16

Sejam \(1\le p<\infty\) e \(X,X_1,X_2,\ldots,Y\) variáveis aleatórias reais num mesmo espaço de probabilidade, com \(Y\ge0\). Suponha-se que \(X_n\overset{\mathbb P}{\rightharpoonup}X\) e

\[ |X_n|\le Y\quad\text{q.c. para todo o }n, \qquad Y\in L^p(\mathbb P). \]

Então \(X\in L^p(\mathbb P)\) e \(X_n\xrightarrow{L^p}X\).

Demonstração. Não exigida no exame Passo 1) Para todo o \(\varepsilon>0\), a desigualdade triangular e \(|X_n|\le Y\) quase certamente dão

\[ \mathbb P(|X|>Y+\varepsilon) \le\mathbb P(|X-X_n|>\varepsilon)\longrightarrow0. \]

A probabilidade à esquerda não depende de \(n\), pelo que é zero. Tomando \(\varepsilon=1/m\) para \(m=1,2,\ldots\) e usando a subaditividade numerável, obtemos \(|X|\le Y\) quase certamente.

Assim, \(\mathbb E(|X|^p)\le\mathbb E(Y^p)<\infty\), pelo que \(X\in L^p(\mathbb P)\).

Convergência dominada: o limite em \(L^p\)

Demonstração do Teorema 4.16. Não exigida no exame Passo 2) Defina-se \(Z_n:=|X_n-X|^p\) e \(Z:=(2Y)^p\). Pelo passo 1, \(0\le Z_n\le Z\) quase certamente, com \(\mathbb E(Z)<\infty\).

Fixem-se \(\delta,M>0\). Em \(Z\le M\), usa-se \(Z_n\le\delta+M\mathbf1_{(\delta,\infty)}(Z_n)\). Em \(Z>M\), usa-se \(Z_n\le Z\). Tomando valores esperados, obtemos

\[ \mathbb E(Z_n)\le\delta+M\,\mathbb P(Z_n>\delta) +\mathbb E\!\left(Z\mathbf1_{(M,\infty)}(Z)\right). \]

Pela convergência em probabilidade,

\[ \mathbb P(Z_n>\delta) =\mathbb P(|X_n-X|>\delta^{1/p})\longrightarrow0. \]

Primeiro, toma-se o limite superior quando \(n\to\infty\). Depois, faz-se \(M\to\infty\): o último valor esperado tende para zero pelo teorema da convergência dominada de Lebesgue. Obtemos assim \(\limsup_{n\to\infty}\mathbb E(Z_n)\le\delta\).

Como \(\delta>0\) é arbitrário, \(\mathbb E(|X_n-X|^p)\to0\), como se pretendia.

Convergência fraca

Seja \(S\) um espaço métrico e denote-se por \(C_b(S)\) o espaço das funções contínuas e limitadas \(f\colon S\to \mathbb R\).

Definição 4.17

Dizemos que as medidas de probabilidade \(\mathbb P_n\) em \((S,\mathcal B(S))\) convergem fracamente para \(\mathbb P\), e escrevemos \(\mathbb P_n\rightharpoonup \mathbb P,\) se

\[ \int_Sg\,d\mathbb P_n\longrightarrow\int_Sg\,d\mathbb P \qquad\text{para todo o }g\in C_b(S). \]

Para elementos aleatórios, dizemos que \(X_n\) converge em distribuição para \(X\), e escrevemos \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X,\) se \(\mathbb P_{X_n}\rightharpoonup \mathbb P_X\).

As variáveis aleatórias podem estar definidas em espaços de probabilidade diferentes.

Continuidade uniforme e conjuntos de continuidade

Definição 4.18

Uma função \(g:S\to\mathbb R\) num espaço métrico \((S,d)\) chama-se uniformemente contínua se, para todo o \(\varepsilon>0\), existe \(\delta>0\) tal que

\[ d(x,y)<\delta\quad\Longrightarrow\quad |g(x)-g(y)|<\varepsilon \qquad(x,y\in S). \]

Definição 4.19

Denota-se por \(A^\circ\) o interior de \(A\), por \(\overline A\) o seu fecho e por \(\partial A:=\overline A\setminus A^\circ\) a sua fronteira.

Para uma medida de probabilidade de Borel \(\mathbb P\) em \(S\), um boreliano \(A\subseteq S\) chama-se conjunto de continuidade de \(\mathbb P\) se \(\mathbb P(\partial A)=0\).

Teorema de Portmanteau

Teorema 4.20

Para medidas de probabilidade de Borel \(\mathbb P,\mathbb P_1,\mathbb P_2,\ldots\) num espaço métrico \(S\), as condições seguintes são equivalentes:

  1. \(\mathbb P_n\rightharpoonup\mathbb P\).
  2. \(\int_Sg\,d\mathbb P_n\to\int_Sg\,d\mathbb P\) para toda a função limitada e uniformemente contínua \(g:S\to\mathbb R\).
  3. \(\limsup_{n\to\infty}\mathbb P_n(F)\le\mathbb P(F)\) para todo o conjunto fechado \(F\).
  4. \(\liminf_{n\to\infty}\mathbb P_n(G)\ge\mathbb P(G)\) para todo o conjunto aberto \(G\).
  5. \(\mathbb P_n(A)\to\mathbb P(A)\) para todo o conjunto de continuidade \(A\) de \(\mathbb P\).

Portmanteau: conjuntos fechados

Demonstração do Teorema 4.20, (1) implica (2) implica (3). Não exigida no exame A primeira implicação decorre de a continuidade uniforme ser uma condição mais forte do que a mera continuidade.

Para um conjunto fechado não vazio \(F\subseteq S\) e inteiros \(m\ge1\), defina-se

\[ d(x,F):=\inf_{y\in F}d(x,y),\qquad g_m(x):=\max\{1-m\,d(x,F),0\}. \]

Então \(|g_m(x)-g_m(y)|\le m\,d(x,y)\), pelo que \(g_m\) é uniformemente contínua. Além disso, \(0\le g_m\le1\) e \(g_m\searrow\mathbf1_F\), pois \(F\) é fechado. Por (2),

\[ \limsup_{n\to\infty}\mathbb P_n(F) \le\lim_{n\to\infty}\int_Sg_m\,d\mathbb P_n =\int_Sg_m\,d\mathbb P. \]

Faça-se \(m\to\infty\). Pelo teorema da convergência dominada de Lebesgue, o último integral tende para \(\mathbb P(F)\). Para \(F=\varnothing\), a desigualdade é imediata.

Portmanteau: abertos e conjuntos de continuidade

Demonstração do Teorema 4.20, (3) e (4) são equivalentes. Não exigida no exame Se \(G\) é aberto, então \(G^c\) é fechado, pelo que (3) dá

\[ \liminf_{n\to\infty}\mathbb P_n(G) =1-\limsup_{n\to\infty}\mathbb P_n(G^c) \ge1-\mathbb P(G^c)=\mathbb P(G). \]

O mesmo cálculo com complementares prova a implicação recíproca.

Estas condições implicam (5). Para um boreliano \(A\), usa-se \(A^\circ\subseteq A\subseteq\overline A\) para obter

\[ \mathbb P(A^\circ) \le\liminf_{n\to\infty}\mathbb P_n(A) \le\limsup_{n\to\infty}\mathbb P_n(A) \le\mathbb P(\overline A). \]

Se \(\mathbb P(\partial A)=0\), os dois extremos são iguais a \(\mathbb P(A)\). Logo, \(\mathbb P_n(A)\to\mathbb P(A)\).

Portmanteau: conjuntos de nível

Demonstração do Teorema 4.20, (5) implica (1). Não exigida no exame Suponha-se (5). Fixe-se uma função contínua \(h:S\to[0,1]\) e, para \(0<t<1\), defina-se

\[ A_t:=h^{-1}((t,\infty)). \]

Por continuidade, um ponto com \(h(x)>t\) tem uma vizinhança contida em \(A_t\), enquanto um ponto com \(h(x)<t\) tem uma vizinhança disjunta de \(A_t\). Assim,

\[ \partial A_t\subseteq h^{-1}(\{t\}). \]

Pela Proposição 1.22, a distribuição de \(h\) sob \(\mathbb P\) tem, no máximo, um conjunto numerável de átomos. Para todos os outros \(t\), tem-se \(\mathbb P(\partial A_t)=0\), pelo que (5) dá

\[ \mathbb P_n(A_t)\longrightarrow\mathbb P(A_t). \]

Portmanteau: integrais de funções contínuas

Demonstração do Teorema 4.20, conclusão. Não exigida no exame Com \(h\) e \(A_t\) como acima, a fórmula do integral de cauda Proposição 3.20 dá

\[ \int_Sh\,d\mathbb P_n=\int_0^1\mathbb P_n(A_t)\,dt \longrightarrow\int_0^1\mathbb P(A_t)\,dt=\int_Sh\,d\mathbb P. \]

O limite decorre do teorema da convergência dominada de Lebesgue: os integrandos pertencem a \([0,1]\) e convergem fora de um conjunto numerável, que tem medida de Lebesgue zero.

Para uma função geral \(g\in C_b(S)\), escolha-se \(M>0\) com \(|g|\le M\) e defina-se \(h=(g+M)/(2M)\). Então \(h:S\to[0,1]\) é contínua e

\[ \int_Sg\,d\mathbb P_n=2M\int_Sh\,d\mathbb P_n-M \longrightarrow2M\int_Sh\,d\mathbb P-M=\int_Sg\,d\mathbb P. \]

Isto é convergência fraca, o que completa todas as equivalências.

Convergência em distribuição: critério das FD

Teorema 4.21

Sejam \(X,X_1,X_2,\ldots\) variáveis aleatórias reais com FD \(F,F_1,F_2,\ldots\). Então

\[ X_n\xrightarrow{\mathmakebox[1.5em]{d}}X \quad\Longleftrightarrow\quad F_n(x)\to F(x) \]

para todo o ponto de continuidade \(x\) de \(F\).

Num ponto de salto de \(F\), não se exige a convergência de \(F_n(x)\).

Demonstração, implicação direta. Suponha-se que \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\) e que \(F\) é contínua em \(x\). Pela Proposição 1.18,

\[ \mathbb P_X(\partial(-\infty,x]) =\mathbb P_X(\{x\})=F(x)-F(x^-)=0. \]

A condição (5) do Teorema 4.20 dá \(F_n(x)=\mathbb P_{X_n}((-\infty,x])\to\mathbb P_X((-\infty,x])=F(x)\).

O critério das FD: intervalos

Demonstração do Teorema 4.21, recíproca. Suponha-se que \(F_n(x)\to F(x)\) em todo o ponto de continuidade de \(F\).

Pela Proposição 1.22, as descontinuidades de \(F\) formam um conjunto numerável. Escolha-se um ponto de continuidade em cada intervalo aberto não vazio com extremos racionais. O conjunto \(D\) destes pontos é numerável e denso em \(\mathbb R\).

Para \(a,b\in D\) com \(a<b\),

\[ \begin{aligned} \mathbb P_{X_n}((a,b]) &=F_n(b)-F_n(a)\\ &\longrightarrow F(b)-F(a)=\mathbb P_X((a,b]). \end{aligned} \]

A mesma convergência vale para qualquer união finita destes intervalos. Para o verificar, escreve-se essa união como uma união disjunta de um número finito de intervalos semiabertos com extremos em \(D\) e usa-se a aditividade finita.

Observação. Todo o aberto \(G\subseteq\mathbb R\) é uma união numerável de intervalos deste tipo. Aproximando por uniões finitas crescentes, obtém-se \(\liminf_{n\to\infty}\mathbb P_{X_n}(G)\ge\mathbb P_X(G)\), pelo que o Teorema 4.20 dá \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\).

Um conjunto com fronteira de probabilidade positiva

Exemplo 4.22

Sejam \(X_n\equiv1/n\) e \(X\equiv0\). Para toda a função contínua e limitada \(g\),

\[ \mathbb E(g(X_n))=g(1/n)\to g(0)=\mathbb E(g(X)), \]

pelo que \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\).

No entanto, para \(A=\{0\}\),

\[ \mathbb P(X_n\in A)=0 \quad\text{mas}\quad \mathbb P(X\in A)=1. \]

Não há contradição: \(0\in\partial A\), pelo que \(\mathbb P(X\in\partial A)=1\).

Convergência em probabilidade e em distribuição

Proposição 4.23

Sejam \(X,X_1,X_2,\ldots\) variáveis aleatórias reais num mesmo espaço de probabilidade. Então

\[ X_n\overset{\mathbb P}{\rightharpoonup}X \quad\Longrightarrow\quad X_n\xrightarrow{\mathmakebox[1.5em]{d}}X. \] Para uma constante \(c\in\mathbb R\), vale também a recíproca:

\[ X_n\xrightarrow{\mathmakebox[1.5em]{d}}c \quad\Longleftrightarrow\quad X_n\overset{\mathbb P}{\rightharpoonup}c. \]

Distribuições iguais sem convergência em probabilidade

Exemplo 4.24

Sejam \(X,X_1,X_2,\ldots\) variáveis de Rademacher independentes:

\[ \mathbb P(X_j=-1)=\mathbb P(X_j=1)=\frac12. \]

\[ \begin{aligned} X_n&\overset{d}=X &&\text{para todo o }n,\\ \mathbb P(|X_n-X|>1)&=\mathbb P(X_n\ne X)=\frac12 &&\text{para todo o }n. \end{aligned} \]

Assim, \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\), mas \(X_n\overset{\mathbb P}{\not\rightharpoonup}X\).

Teorema de continuidade de Lévy

Teorema 4.25

Sejam \(X_1,X_2,\ldots\) variáveis aleatórias reais com funções características \(\varphi_1,\varphi_2,\ldots\). Se

\[ \varphi_n(t)\to\varphi(t)\qquad\text{para todo o }t\in\mathbb R \]

e \(\varphi\) é contínua em \(0\), então \(\varphi\) é a função característica de alguma distribuição \(\mathbb P\) e

\[ \mathbb P_{X_n}\rightharpoonup \mathbb P. \]

Reciprocamente, \(X_n\xrightarrow{\mathmakebox[1.5em]{d}}X\) implica \(\varphi_n(t)\to\varphi_X(t)\) para todo o \(t\).

Teorema da aplicação contínua

Teorema 4.26

Sejam \(S,T\) espaços métricos, \(X,X_1,X_2,\ldots\) elementos aleatórios em \(S\) e \(g:S\to T\) uma aplicação mensurável. Seja \(D_g\) o seu conjunto de pontos de descontinuidade. Se

\[ X_n\xrightarrow{\mathmakebox[1.5em]{d}}X \qquad\text{e}\qquad \mathbb P(X\in D_g)=0, \]

então

\[ g(X_n)\xrightarrow{\mathmakebox[1.5em]{d}}g(X). \]

Para a convergência quase certa ou em probabilidade, vale a conclusão análoga quando \(g\) é contínua em \(X\) quase certamente.

Relações de implicação

Implicações válidas entre convergência quase certa, em L infinito, em L q, em L p, em probabilidade e em distribuição.

As setas \(L^q\Rightarrow L^p\) pressupõem \(q>p\ge1\) e \(\mathbb P(\Omega)=1\).

Recíprocas: contraexemplos e condições adicionais

Recíproca inválida Contraexemplo Condição adicional suficiente
\(L^p\not\Rightarrow\) quase certamente Sucessão da máquina de escrever. Proposição 4.14 —
Quase certamente \(\not\Rightarrow L^1\) Sucessão de picos. Exemplo 4.15 Uma variável dominante integrável. Teorema 4.16 com \(p=1\)
Em distribuição \(\not\Rightarrow\) em probabilidade Cópias independentes. Exemplo 4.24 Um limite constante. Proposição 4.23

Uma versão preliminar destes slides foi preparada com o apoio de agentes de programação baseados em IA.