Aula 5 Distribui esdeProbabilidade II 10-06-2010

35
Distribuições Contínuas Normal Gama Valores Extremos Exponencial Discretas Binomial Poisson Geométrica

Transcript of Aula 5 Distribui esdeProbabilidade II 10-06-2010

Page 1: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Distribuições

Contínuas

• Normal• Gama• Valores Extremos• Exponencial

Discretas

• Binomial• Poisson• Geométrica

Page 2: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Funções

• A diferença entre as distribuições contínuas e discretas é que as distribuições discretas envolvem somas sobre uma função de probabilidade descontínua, enquanto a probabilidade para variáveis aleatórias contínuas envolve integração sobre funções contínuas denominadas FUNÇÃO DENSIDADE DE PROBABILIDADE, ou (PDFs).

• Convencionalmente, a PDF para uma variável aleatória X é denominada f(x).

Page 3: Aula 5 Distribui esdeProbabilidade II 10-06-2010

PDF ou FDP

• Uma função densidade de probabilidade, PDF, tem que satisfazer a seguinte condição.

• Além disso, f(x) não pode ser negativo para todos os valores de x.

∫ =x

xf 1)(

Page 4: Aula 5 Distribui esdeProbabilidade II 10-06-2010

A f(x) tem significado quando pensamos em calcular probabilidades para valores de uma variável aleatória em uma vizinhança não infinitesimal em torno de um ponto, por exemplo X=1.

Page 5: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• Uma idéia relacionada com a PDF éaquela de uma função de distribuição cumulativa (CDF). A CDF é uma função da variável aleatória X, dada pela integral da PDF até um valor particular de x. Convencionalmente, CDFs são denominadas F(x):

∫ ≤=≤=

xXdxxfxXxF )(}Pr{)(

Page 6: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Distribuição Normal

Uma variável aleatória contínua érepresentada por uma distribuição normal se a sua distribuição for:

• Simétrica• A forma gráfica é similar a um sino

Page 7: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Distribuição Normal

( )( )

2

2

2

X

e2

1Xf σ

µ−−

πσ=

+∞<<∞− xpara

µ média da populaçãoσ desvio-padrão da população

µµµµ

σ

Page 8: Aula 5 Distribui esdeProbabilidade II 10-06-2010

A função densidade pode ser compreendia como uma extensão natural do histograma.

Page 9: Aula 5 Distribui esdeProbabilidade II 10-06-2010
Page 10: Aula 5 Distribui esdeProbabilidade II 10-06-2010
Page 11: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• Cada para de parâmetros (µ; σ) define uma distribuição normal distinta;

Page 12: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• Convencionalmente, a variável aleatória descrita por uma distribuição Gaussiana padrão é denominada de Z.

• Sua densidade de probabilidade se simplifica a:

Onde Z é

( )

−=

2

2exp

2

1)(

Zxf

π

σµ−= X

Z

Page 13: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• Quando estimamos os coeficientes da Gassiana pelos dados, então utilizamos a seguinte notação (notem que agora a variável transformada é denotada como “z”):

s

xxz

−=

Page 14: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Exemplo 1:• Suponha que uma distribuição Gaussiana para o mês de

janeiro em uma certa localidade seja caracterizada por µ=22.2º C e σ=4.4º C.

• Suponha que você esteja interessado em avaliar a probabilidade de que um certo mês de janeiro tenha temperaturas menores ou iguais a 21.4º C.

• O primeiro passo para a solução desse problema écalcular o valor padronizado z.

z = (21.4º C – 22.2º C)/4.4º C = -0.18.

Assim, a probabilidade de uma temperatura igual ou mais fria que 21.4º C é a mesma que a probabilidade de um valor de Z igual ou menor que -0.18:

Pr{X≤21.4º C} = Pr{Z≤ -0.18}.

Page 15: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Como encontrar a probabilidade

• Identifique na tabela que contenha probabilidades acumuladas para a distribuição Gaussiana o valor z.

• Encontre a linha que indica -0.1 e a coluna que mostra 0.08

• Você encontra a probabilidade de 0.4286.

Page 16: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Como funciona a tabela

z = z(linha) + z(coluna)

Linh

aColuna

Page 17: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Exemplo 2:

Considere uma variável aleatória X com µ=15 e σ=25.

Qual a probabilidade de que X assuma valores entre 16 ≤ X ≤ 20 ?

Page 18: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• A probabilidade desejada pode ser obtida, a partir da variável transformada (Z).

Para X = 16 ⇒

Para X = 20 ⇒

P(0,04 ≤ Z ≤ 0,20) = P(Z ≤ 0,2)-P(Z ≤ 0,04)

( )04,0

25

1516Z =−=

( )20,0

25

1520Z =−=

Page 19: Aula 5 Distribui esdeProbabilidade II 10-06-2010

P(Z ≤ 0,2) =

P(Z ≤ 0,04) =

Page 20: Aula 5 Distribui esdeProbabilidade II 10-06-2010

P(Z ≤ 0,2) = 0,5793

P(Z ≤ 0,04) = 0,5160

P(0,04 ≤ Z ≤ 0,20) = P(Z ≤ 0,2)-P(Z ≤ 0,04)

P(0,04 P(0,04 ≤≤ Z Z ≤≤ 0,20) = 0,5793 0,20) = 0,5793 –– 0,5160 = 0,0633 ou 6,33%0,5160 = 0,0633 ou 6,33%

Page 21: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Exemplo 3

Queremos saber agora qual a temperatura que corresponde ao primeiro decil, sendo que esta distribuicao de temperatura tem µ=22.2º C e σ=4.4º C.

Para isso, utilizamos a tabela para achar o valor de probabilidade 0.10. z = -1.28

Page 22: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Mas

X = (-1.28)x (4.4) + 22.2 = 16.6º C

σµ−= X

Z

µσ += ZXLogo

Page 23: Aula 5 Distribui esdeProbabilidade II 10-06-2010

DISTRIBUIÇÃO GAMA• Muitas variáveis atmosféricas são distintamente

assimétricas, e possuem uma assimetria para a direita. • Um exemplo clássico é a precipitação. Sabemos que

não existem precipitações negativas, certo?• Então suponha que uma certa localidade tenha uma

média de precipitação diária de 1.96 in e desvio-padrão de 1.12 in.

•• Utilizando a tabela de distribuiUtilizando a tabela de distribuiçções acumuladas ões acumuladas Gaussiana podemos calcular a probabilidade de Gaussiana podemos calcular a probabilidade de precipitaprecipitaçções negativas como ões negativas como

Pr { Z Pr { Z ≤≤ (0.00 (0.00 –– 1.96)/1.12} = Pr{Z1.96)/1.12} = Pr{Z≤≤ --1.75} = 0.040.1.75} = 0.040.• Esta probabilidade calculada não é especialmente

grande, mas por outro lado não pode ser considerada zero.

• AGORA, SABEMOS PELA NATUREZA QUE PRECIPITAÇÕES NEGATIVAS SÃO IMPOSSÍVEIS!

Page 24: Aula 5 Distribui esdeProbabilidade II 10-06-2010

DISTRIBUIÇÃO GAMA

)(

)/exp()/()(

1

αβββ α

Γ−=

− xxxf onde x, α, β > 0

α PARÂMETRO DE FORMA , β PARÂMETRO DE ESCALA. Γ(α) FUNÇÃO GAMMA, definida pela integral

∫∞

−−=Γ0

1)( dttetαα

Page 25: Aula 5 Distribui esdeProbabilidade II 10-06-2010

∫∞

−−=Γ0

1)( dtet tααΓ(α) FUNÇÃO GAMMA

)()1( ααα Γ=+Γ

Exemplo 1

Γ(3.5) = 2.5 x Γ(2.5)

= 2.5 x 1.5 x Γ(1.5)

= 2.5 x 1.5 x 0.8862 =3.323.

!)1( αα =+ΓOu

Exemplo 2

Γ(4) = Γ(3 + 1)

= 3 !

= 3 x 2 = 6

Page 26: Aula 5 Distribui esdeProbabilidade II 10-06-2010

α PARÂMETRO DE FORMA e β PARÂMETRO DE ESCALA

Page 27: Aula 5 Distribui esdeProbabilidade II 10-06-2010

α=k PARÂMETRO DE FORMA e β= θ PARÂMETRO DE ESCALA

ββ== θ θ θ θ θ θ θ θ

Page 28: Aula 5 Distribui esdeProbabilidade II 10-06-2010

A PDF da distribuição Gamma pode apresentar uma grande variedade de formas, dependendo, portanto, do parâmetro de forma α.

Para valores de α muito altos, a distribuição gamma tende à Gaussiana

O parâmetro de escala β, tem a função de ESTICAR OU ENCOLHER (isto é escalonar) a função de densidade gamma para a direita ou esquerda, dependendo das magnitudes gerais dos valores dos dados representados.

Page 29: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Calculando os fatores da função gama

• Existem duas maneiras:1) Estimadores de Thom – Thom (1958)) 2) Aproximação polinomial

Page 30: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Estimadores de Thom – Thom

D

D

4

3/411ˆ

++=α

αβ

ˆˆ x=

∑=

−=n

iix

nxD

1)ln(

1)ln(

parâmetro de forma α

parâmetro de escala β

Page 31: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Aproximação polinomial

D

DD 20544274.01648852.05000876.0ˆ

−+=α

32968477.1179728.17

29775373.0059950.9898919.8ˆ

DDD

DD

++

++=α

Para 0≤ D ≤0.5772

Para 0.5772 ≤ D ≤ 17

αβ

ˆˆ x= ∑

=−=

n

iix

nxD

1)ln(

1)ln(

parâmetro de forma α

parâmetro de escala β

Page 32: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Para calcular a Probabilidade

é sempre necessário fazer uma transformação para re-escalonar a variável X de interesse (caracterizada por uma gama com parâmetro de escala arbitrário β) para a variável

βξ x=

Page 33: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Exemplo 1Suponha que a a média aritmética da ppt é igual a 1.96 polegadas.

Sera Sera queque a a precipitaprecipita çção de janeiro de 1987, 3.15 polegadas ão de janeiro de 1987, 3.15 polegadas ééusual?usual?

1) Calculando o valor da média dos logaritmos dos totais mensais = 0.5346

2) Obtemos o valor de D = 0.1393) De acordo com o método de Thom α=3.76 e β=0.52 polegadas.

ζ = (X/β) = 3.15 /0.52 = 6.06.

O passo seguinte é encontrar no corpo da tabela onde se encontra a probabilidade 6.06 para α=3.76. Tabela 2 - Wilks

Page 34: Aula 5 Distribui esdeProbabilidade II 10-06-2010

• Este valor se encontra entre os valores tabulados F(5.214)=0.90 e F(6.354)=0.90.

• Logo uma interpolação nos dá F(6.06)=0.874. • Portanto a probabilidade de chover em janeiro

menos ou igual a 3.15 polegadas em Ithaca éde 0.874.

• O complementar (1- 0.874) = 0.126 é a probabilidade de chover mais do que esse valor (a qual é equivalente a aproximadamente 1 chance em 8 (1/8).

Page 35: Aula 5 Distribui esdeProbabilidade II 10-06-2010

Lista de Exercício 3Entrega: 17 de Junho

• Utilizando como referência a temperatura média do ar e o respectivo desvio padrão calculado na lista de exercício 1 e que a distribuição de temperatura do ar pode ser representada por uma distribuição normal, calcule.

• 1) Qual a probabilidade de termos temperaturas acima de 25,30 e 35oC ?

• 2) Qual a probabilidade de a temperatura do ar estiver entre o intervalo de 25 a 30oC ?

• 3) Qual a probabilidade de termos temperaturas abaixo de 10,15 e 20oC ?

• 4) Qual é a temperatura que representa um decil?