4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a...

36
4. PRINCIPAIS MODELOS 4. PRINCIPAIS MODELOS DISCRETOS DISCRETOS 2012 2012

Transcript of 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a...

Page 1: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

4. PRINCIPAIS MODELOS 4. PRINCIPAIS MODELOS

DISCRETOSDISCRETOS

20122012

Page 2: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Principais modelos probabilísticos discretos

4.1. Modelo Bernoulli

Muitos experimentos admitem apenas dois resultados.

Exemplos:

1. Uma peça é classificada como defeituosa ou não defeituosa;

2. O resultado de um exame médico para detecção de uma doença é positivo ounegativo.

3. Um entrevistado concorda ou não com uma afirmação feita;

2

3. Um entrevistado concorda ou não com uma afirmação feita;

4. No lançamento de um dado ocorre ou não face 6;

5. Um item produzido é classificado como conforme ou não conforme.

Situações com alternativas dicotômicas podem ser representadas genericamentepor resposta do tipo sucesso ou insucesso (fracasso ou falha).

Esses experimentos recebem o nome de ensaios de Bernoulli e originam uma v.a.com distribuição de Bernoulli.

Page 3: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Distribuição de Bernoulli

X é uma v.a. que assume apenas dois valores: 1 se ocorrer sucesso (S) e 0 seocorrer fracasso (F). Sendo p a probabilidade de sucesso, 0 < p <1.

X(S) = 1 e X(F) = 0. A distribuição de probabilidade é dada por

x

P(X=x)

0 1

1 – p p

Notação: X ~ Bernoulli (p) indica que a v.a. X tem distribuição de Bernoulli. O

=−

===−

c.c.,0

1. ,0 se,)1()()(

1xpp

xXPxf

xx

3

Notação: X ~ Bernoulli (p) indica que a v.a. X tem distribuição de Bernoulli. Oparâmetro da distribuição é p.

Se X ~ Bernoulli(p), então

E(X) = p

e Var(X) = p (1 – p).

Repetições independentes de um ensaio de Bernoulli dão origem aomodelo binomial.

Page 4: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

4.2. Modelo binomial

Exemplo. Uma moeda é lançada 3 vezes e a probabilidade de cara é p em cadalançamento. Determinar a distribuição de probabilidade da variável número decaras nos 3 lançamentos (X).

Denotemos S: sucesso, ocorre cara (c) e F: fracasso, ocorre coroa (k).

O espaço amostral para este experimento é

Ω = FFF, FFS, FSF, SFF, FSS, SFS, SSF,SSS.

Fazemos Xi ~ Bernoulli(p), i = 1,2,3. Logo, X = X1 + X2 + X3 representa o númerode caras nos 3 lançamentos.

4

Ω Probabilidade X1 X2 X3 X=X1+X2+X3

FFF (1-p)3 0 0 0 0

FFS (1-p)2p 0 0 1 1

FSF (1-p)2p 0 1 0 1

SFF (1-p)2p 1 0 0 1

FSS (1-p)p2 0 1 1 2

SFS (1-p)p2 1 0 1 2

SSF (1-p)p2 1 1 0 2

SSS P3 1 1 1 3

Page 5: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

.)(P)3(P

e )1(3),,(P)2(P

,)1(3),,(P)1(P

,)1()(P)0(P

3

2

2

3

pSSSX

ppSSFSFSFSSX

ppSFFFSFFFSX

pFFFX

===

−===

−===

−===

Calculamos

A distribuição de probabilidade da v.a. X é dada por

3210x

5

3223 )1(3)1(3)1()()( ppppppxXPxf −−−==

f(x) pode ser escrita como

.)!3(!

!33 que em

c.c.,0

,3,2,1,0 se,)1(3

)(3

xxx

xppxxf

xx

−=

=−

=

Page 6: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Distribuição binomial

Repetição de n ensaios de Bernoulli independentes, todos com a mesmaprobabilidade de sucesso p. A variável aleatória que conta o número de sucessosnos n ensaios de Bernoulli é denominada de variável aleatória binomial comparâmetros n e p.

binomial. ecoeficient o representa !

que em

,c.c.,0

,,,1,0 se,)1()(

nn

nxppx

n

xfxnx

=

=−

=

−L

6

binomial. ecoeficient o representa )!(!

! que em

xnx

n

x

n

−=

Notação: X ~ B(n,p) para indicar que a v.a. X tem distribuição Binomial comparâmetros n e p.

Se X ~ B(n, p), então

E(X) = np e

Var(X) = np(1 – p).

Page 7: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0 2 4 6 80.

00.

20.

4

x

P(X

=x)

p=0,1

0 2 4 6 8

0.00

0.20

x

P(X

=x)

p=0,3

Distribuição B(n =10, p)

7

0 2 4 6 8

0.00

0.15

x

P(X

=x)

p=0,5

0 2 4 6 80.

000.

20

x

P(X

=x)

p=0,8

Page 8: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0 5 10 15 200.

000.

20

x

P(X

=x)

p=0,1

0 5 10 15 20

0.00

0.15

x

P(X

=x)

p=0,3

Distribuição B(n = 20, p)

8

0 5 10 15 20

0.00

0.10

x

P(X

=x)

p=0,5

0 5 10 15 200.

000.

15

x

P(X

=x)

p=0,8

Page 9: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0 10 20 300.

000.

15

x

P(X

=x)

p=0,1

0 10 20 30

0.00

0.10

x

P(X

=x)

p=0,3

Distribuição B(n = 30, p)

9

0 10 20 30

0.00

0.10

x

P(X

=x)

p=0,5

0 10 20 300.

000.

15

x

P(X

=x)

p=0,8

Page 10: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

O professor da disciplina de Estatística elaborou uma prova de múltipla escolha,composta de 10 questões cada uma com 5 alternativas. Aprovação na disciplinarequer pelo menos 6 questões corretas. Se um aluno responde a todas asquestões baseado em palpite (“chute”), qual a probabilidade de ser aprovado?

Solução. X é a v.a. número de questões respondidas corretamente nas 10questões. Eventos: S: “questão respondida corretamente” e F: “questãorespondida incorretamente”.

P(S) = 1 / 5 e P(F) = 4 / 5. Logo, X ~ B(10, p).

Exemplo

10

=

=

c.c.,0

,10,,1,0 se,5

4

5

110

)(

10

Lxxxf

xx

A probabilidade de aprovação é

.00637,00,99363061

)5(1)5(P1)6(P1)6(P

=−=

−=≤−=<−=≥ FXXX

Page 11: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

x f(x) F(x) 0 0,107374 0,10737 1 0,268435 0,37581 2 0,301990 0,67780 3 0,201327 0,87913 4 0,088080 0,96721 5 0,026424 0,99363

Exemplo

11

5 0,026424 0,99363 6 0,005505 0,99914 7 0,000786 0,99992 8 0,000074 1,00000 9 0,000004 1,00000 10 0,000000 1,00000

Em R: dbinom(0:10,10,1/5) e pbinom(0:10, 10, 1/5).

Page 12: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0.20

0.30

P(X

=x)

B(10,p=0,20)

Exemplo

12

0 2 4 6 8 10

0.00

0.10

x

P(X

=x)

Page 13: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Exemplo

Um fabricante adquire certo tipo de componente de um fornecedor. Segundoeste fornecedor, a proporção de componentes defeituosos é 2%.

(a) O fabricante seleciona 15 componentes de um lote para inspeção. Qual aprobabilidade de que seja encontrado pelo menos um componente defeituosoneste lote?

(b) O fabricante adquire 10 lotes por mês e de cada lote são selecionados 15componentes para inspeção, como no item (a). Qual a probabilidade de quesejam encontrados três lotes com pelo menos um componente defeituoso?

13

Solução. (a) Definimos o evento sucesso (S) como “o componente selecionado édefeituoso”. Pelo enunciado, P(S) = p = 0,02. A v.a. X é definida como sendo onúmero de componentes defeituosos (sucessos) em n = 15 componentes.Supondo independência, X ~ B(n = 15, p = 0,02).

Devemos calcular P(X ≥ 1), que é dada por

.261,098,01

)02,01(02,00

151

)0(P11P1)1(P

15

0150

=−=

−××

−=

=−=<−=≥

X)(XX

Em Excel:

= 1 - DISTRBINOM(0; 15; 0,02; FALSO)

Page 14: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Exemplo

Solução. (b) Definimos o evento sucesso (S) como “o lote contém pelo menos umcomponente defeituoso”. De acordo com o item (a), P(S) = p = 0,261. A v.a. Y édefinida como sendo o número de lotes com pelo um componente defeituoso(sucessos) em n = 10 lotes. Supondo independência, Y ~ B(n = 10, p = 0,261).

Devemos calcular P(Y = 3), que é dada por

.257,0)261,01(261,03

10)3(P 3103 =−××

== −

Y

Em Excel:

= DISTRBINOM(3; 10; 0,261; FALSO)

( a ) ( b )

14

0 5 1 0 1 5

0.0

0.2

0.4

0.6

( a )

x

P(X

= x

)

0 2 4 6 8 1 0

0.00

0.05

0.10

0.15

0.20

0.25

( b )

y

P(Y

= y

)

Page 15: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

4.3. Modelo hipergeométrico

Um conjunto de N elementos é dividido em duas classes. Uma classe com M (M < N)elementos (sucessos) e a outra com N – M elementos (fracassos).

Por exemplo, no caso de N itens produzidos, podem ser considerados M itens defeituosose N – M itens não defeituosos.

Uma amostra de tamanho n (n < N) é sorteada sem reposição. A v.a. X é definida como onúmero de elementos com a característica de interesse (sucesso) na amostra de tamanho n.

(1) n elementos são selecionados de um conjunto de N elementos. (2) xsucessos são escolhidos de uma classe com M sucessos. (3) Finalmente, n – xfracassos são escolhidos de uma classe com N – M fracassos.

15

−−=

=

c.c.,0

,,min, ),(,0max se,)(

MnMNnx

n

N

xn

MN

x

M

xfL

Notação: X ~ H(N, M, n) indica que a v.a. X tem distribuição hipergeométricacom parâmetros N, M e n.

Se X ~ H(N, M, n), então .1

1)( e )(

=

=

N

nN

N

M

N

MnXVar

N

MnXE

A função de probabilidade da v.a. X é

Page 16: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Exemplo (Hines et al., 2006, p. 105)

Em um departamento de inspeção de recebimento, lotes de eixo de bomba sãorecebidos periodicamente. Os lotes contêm 100 unidades e o seguinte plano deamostragem de aceitação é usado. Seleciona-se uma amostra de 10 unidadessem reposição. O lote é aceito se a amostra tiver, no máximo, um eixodefeituoso. Suponha que um lote seja recebido e que 5% dos itens sejamdefeituosos. Qual a probabilidade de que o lote seja aceito ?

X: número de defeituosos na amostra ⇒ X ~ H(N = 100, M = 5, n = 10).

)1(P)0(P)1(P)loteoaceitar(P =+==≤= XXX

16

.923,0

10

100

9

95

1

5

10

100

10

95

0

5

=

+

=

Em R: dhyper(0,5,95,10) + dhyper(1,5,95,10) ou phyper(1,5,95,10).

Em Excel: =DIST.HIPERGEOM(0;10;5;100) + DIST.HIPERGEOM(1;10;5;100).

Page 17: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0.3

0.4

0.5

0.6

P(X

=x)

Exemplo (Hines et al., 2006, p. 105)

17

0 1 2 3 4 5

0.0

0.1

0.2

x

P(X

=x)

Page 18: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Muitos experimentos consistem em observar a ocorrência de eventos emdeterminada unidade (de tempo, volume, comprimento, área, ...)

4.4. Modelo de Poisson

Exemplos

1. Número de consultas a uma base de dados em um minuto.

2. Número de acidentes de trabalho por semana em uma fábrica.

3. Número de pequenas manchas por m2 no esmaltado de uma geladeira.

18

4. Número de chamadas que chegam a uma central telefônica de uma empresa a cada 10 min.

5. Número de carros que chegam ao campus entre 7:00 e 8:00h.

6. Número de microorganismos por cm3 de água contaminada.

7. Número de defeitos em cada teclado produzido por uma fábrica.

Page 19: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

O fenômeno estudado ocorre em intervalos (de tempo, por exemplo).

O intervalo pode ser dividido em subintervalos com comprimentossuficientemente pequenos tais que

• a probabilidade de ocorrência de mais um evento em umsubintervalo é pequena,

• a probabilidade de ocorrência de um evento em um subintervaloseja a mesma para todos os subintervalos e proporcional ao

Suposições básicas

19

seja a mesma para todos os subintervalos e proporcional aocomprimento do subintervalo e

• a contagem em cada subintervalo seja independente de outrossubintervalos.

Pode ser provado que a distribuição do número de ocorrências é Poisson.

Page 20: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Uma v. a. discreta X tem distribuição de Poisson com parâmetro µ se suafunção de probabilidade é dada por

==

,c.c.,0

,,2,1,0 se,!)(

Lxx

exf

xµµ

Distribuição de Poisson

20

em que x é número de eventos em t unidades de medida,

λ é o número médio de eventos (taxa) em uma unidade de medida (t = 1) e

µ = λ t é o número médio de eventos em t unidades de medida.

Notação: X ~ Po(µ) indica que a v.a. X tem distribuição de Poisson comparâmetro µ.

Propriedades: E(X) = µ e Var(X) = µ.

Page 21: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Distribuição Po(µ)

0 5 1 0 1 5 2 0 2 5

0.0

0.2

0.4

P o (0 ,8 )

x

P(X

= x

)

0 5 1 0 1 5 2 0 2 5

0.00

0.10

0.20

P o (2 )

x

P(X

= x

)

21

x x

0 5 1 0 1 5 2 0 2 5

0.00

0.10

P o (5 ,7 )

x

P(X

= x

)

0 5 1 0 1 5 2 0 2 5

0.00

0.04

0.08

0.12

P o (1 0 )

x

P(X

= x

)

Page 22: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

As chegadas a um posto de atendimento ocorrem de formaindependente seguindo a distribuição de Poisson. Suponha que amédia de chegadas é 3 a cada 4 minutos. Qual é a probabilidade deque este posto receba no máximo 2 solicitações em um intervalo de 2minutos?Solução. Se X é número de chegadas a este posto a cada 2 minutos,

então X ~ Po(µ). Aqui, t = 2 min e λ = ¾ = 0,75. Logo, µ = 0,75 × 2 = 1,5. Ou seja,

X ~ Po(1,5) e

Exemplo

P o (1 ,5 )

22

X ~ Po(1,5) e

....3,2,1,0 se,!

5,1)(

5,1

==−

xx

exf

x

.809,0)2

5,15,11(

)2(P)1(P)0(P)2()2(P

Calculamos

25,1 =++=

=+=+===≤

−e

XXXFX

Em R: ppois(2,1.5); em Excel: =POISSON(2;1,5; VERDADEIRO).

0 2 4 6 8 1 0

0.00

0.05

0.10

0.15

0.20

0.25

0.30

x

P(X

= x

)

Page 23: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

O número de pedidos de empréstimos que um banco recebe por dia é umavariável aleatória, sendo que em média são recebidos 7,5 empréstimos pordia. Determine as probabilidades de que, em um dia qualquer, o bancoreceba

(a) Exatamente 2 pedidos de empréstimo;

(b) No máximo 2 pedidos de empréstimo;

(c) No mínimo 8 pedidos de empréstimo.

Exemplo

Solução. Supomos que X (número de

0.15

Po(7,5)

23

Solução. Supomos que X (número depedidos de empréstimos que o bancorecebe por dia) tem distribuição Poissoncom média µ = 7,5. Logo,

.,2,1,0 se,!

5,7)(

5,7

L==−

xx

exf

x

0 5 10 15 20

0.00

0.05

0.10

x

P(X

= x

)

Page 24: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

x f(x)=P(X=x)

0 0,000553

1 0,004148

2 0,015555

3 0,038889

4 0,072916

5 0,109375

6 0,136718

7 0,146484

8 0,137329

9 0,114440

10 0,085830

11 0,058521

12 0,036575

13 0,021101

Exemplo

Em R: dpois(0:27, 7.5); em Excel:

24

13 0,021101

14 0,011304

15 0,005652

16 0,002649

17 0,001169

18 0,000487

19 0,000192

20 0,000072

21 0,000026

22 0,000009

23 0,000003

24 0,000001

25 0,000000

26 0,000000

27 0,000000

Page 25: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

0,0156,2

)5,7()2(P)(

25,7

===−

eXa

e 0,0203 0,0155550,0041480,000553

)2(P)1(P)0(P)2()2(P)(

=++=

=+=+===≤ XXXFXb

Exemplo

Calculamos

25

e 0,0203 0,0155550,0041480,000553 =++=

0,4754.0,5246385 1

0,146484)0,000553(1

)(P1)7(1)8(P1)8(P)(7

0

=−=

++−=

=−=−=<−=≥ ∑=

L

x

xXFXXc

Page 26: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Contaminação é um problema de fabricação de discos ópticos. O número departículas de contaminação que ocorrem em um disco óptico tem uma distribuiçãode Poisson e o número médio de partículas por cm2 de superfície é 0,1. A área dasuperfície do disco em estudo é 100 cm2. Encontre a probabilidade de que 12partículas sejam encontradas em um disco.

Solução. Se X é o número de partículas nasuperfície do disco, então X ~ Po(µ).Temos t = 100 cm2 e λ = 0,1 por cm2.Logo, µ = t × λ = 100 × 0,1 = 10. Ou seja,

Exemplo

0.12

Po(10)

26

Logo, µ = t × λ = 100 × 0,1 = 10. Ou seja,X ~ Po(10) e

L,2,1,0,!

10)(

10

==−

xx

exf

x

0,095.!12

10)12(P

1210

===−

eX

Calculamos

Em R: dpois(12,10);

em Excel: =POISSON(12;10;FALSO).0 5 10 15 20 25

0.00

0.02

0.04

0.06

0.08

0.10

x

P(X

= x

)

Page 27: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Exemplo. Em uma fábrica, dados históricos mostram que em trêssemanas típicas os números médios de acidentes são 2,5 na primeirasemana, 2 na segunda semana e 1,5 na terceira semana. Suponha que onúmero de acidentes por semana segue uma distribuição de Poisson.Qual a probabilidade de que ocorram 4 acidentes em três semanas

Resultado. Se X1,..., Xn são variáveis aleatórias independentes com distribuição de Poisson com parâmetros µ1,..., µn respectivamente, então a variável aleatória Y = X1 + ... + Xn tem distribuição Poisson com parâmetro µ = µ1 + ... + µn.

27

Qual a probabilidade de que ocorram 4 acidentes em três semanastípicas?

.1339,0!4

6)4(P

64

===−

eY

Solução. Xi representa o número de acidentes na i-ésima semana, i =1,2,3, com Xi ~ Po(µi). Supomos que X1, X2 e X3 são independentes.Portanto, Y = X1 + X2 + X3 tem distribuição Poisson com parâmetro µ =2,5 + 2 + 1,5 = 6. Calculamos

Page 28: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

4.5. Modelo geométrico

Ensaios de Bernoulli são realizados de forma independente e cada umcom probabilidade de sucesso igual a p.

Estamos interessados no número de ensaios que antecedem aocorrência do 1o sucesso.

A v.a. X que conta este número tem distribuição geométrica com parâmetro p,notando que X ∈ 0, 1 , 2, ....

Se “S” e “F” representam os eventos sucesso e fracasso e X = x, temos asequência

28

Sendo assim,

sequência

.SF...FF

fracassos

4434421x

.)1(...)1()1()(P

fracassos

ppppxX

x

×−××−×−==444444 3444444 21

Page 29: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Distribuição geométrica

Se ensaios de Bernoulli independentes e com probabilidade de sucessoigual a p são realizados, o número de ensaios que antecedem o primeirosucesso tem uma distribuição geométrica com parâmetro p e sua funçãode probabilidade é dada por

Notação: X ~ Geo(p).

.10 e ,...2,1,0 se ,)1()(P)( <<=−=== pxppxXxf x

29

Se X ~ Geo(p), então

E(X) = (1 – p) / p e

Var(X) = (1 – p) / p2.

Propriedade: Se X ~ Geo(p), então P(X > k + m | X > m) = P(X > k).

É a única distribuição discreta com esta propriedade (“falta de memória”).

Page 30: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

G e o ( 0 .0 1 )

x

P(X

=x)

0 1 2 3 4 5 6 7 8 9 10

0.00

920.

0096

0.01

00

G e o ( 0 .1 )

x

P(X

=x)

0 1 2 3 4 5 6 7 8 9 10

0.04

0.06

0.08

0.10

Distribuição Geo(p)

30

G e o ( 0 .5 )

x

P(X

=x)

0 1 2 3 4 5 6 7 8 9

0.0

0.2

0.4

G e o ( 0 .9 )

x

P(X

=x)

0 1 2 3 4 5

0.0

0.2

0.4

0.6

0.8

Page 31: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Outra definição de distribuição geométrica

Se ensaios de Bernoulli independentes e com probabilidade de sucessoigual a p são realizados, o número de ensaios Y até que ocorra o primeirosucesso tem uma distribuição geométrica com parâmetro p e sua funçãode probabilidade é dada por

Relação entre as duas definições:

.10 e ,...2,1y se ,)1()(P)( 1 <<=−=== − pppyYyf y

31

Y = X + 1,

E(Y) = E(X) + 1 = (1 – p ) / p + 1 = 1 / p e

Var(Y) = Var(X) = (1 – p) / p2.

Obs. Qual a relação entre a distribuição geométrica e os álbuns defigurinhas?

Page 32: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Certo experimento deve ser realizado até que seja obtido um resultado bemsucedido. As realizações são independentes e o custo de cada experimento é$25.000, sendo que se o resultado for um insucesso, há um custo adicionalde $5.000 para o preparo da próxima realização.

(a) Obtenha o custo esperado do experimento.

(b) Se o orçamento não pode ultrapassar $500.000, qual a probabilidade de queeste valor seja ultrapassado.

Solução. Definimos Y como sendo o número de realizações até que ocorra oprimeiro resultado bem sucedido, notando que Y ∈ 1, 2, ... e tem distribuição

Exemplo (Hines et al., 2006, p. 101)

32

primeiro resultado bem sucedido, notando que Y ∈ 1, 2, ... e tem distribuiçãogeométrica com parâmetro p e f(y) = (1 – p)y–1p (veja lâmina 29). Pelo enunciado,o custo é uma v.a., função de Y, dada por

C(Y) = 25000 Y + 5000 (Y – 1) = 30000 Y – 5000.

Usando propriedades do valor esperado obtemos

E[C(Y)] = 30000 E(Y) – 5000 = 30000 / p – 5000.

Se p = 0,25, o custo esperado vale $115.000.

Page 33: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Na letra (b) devemos calcular P(C(Y) > 500000). Usando a expressão de C(Y),

Exemplo (Hines et al., 2006, p. 101)

P(C(Y) > 500000) = P(30000 Y – 5000) > 500000)

= P(Y > 505000 / 30000) = P(Y > 16,8)

= 1 – P(Y ≤ 16,8) = 1 – P(Y ≤ 16)

.10 ,)1(116

1

1 <<−−= ∑=

−ppp

k

k

1.0

33

0.0 0 .1 0 .2 0 .3 0 .4

0.0

0.2

0.4

0.6

0.8

1.0

p

P(C

(Y)

> 5

0000

)

Se p = 0,25,

P(C(Y) > 500000) = 0,010.

Em R: 1 – pgeom(15, 0.25).

Page 34: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

4.6. Modelo binomial negativa

Ensaios de Bernoulli são realizados de forma independente e cada umcom probabilidade de sucesso igual a p.

Interesse no número de ensaios que até que ocorram r sucessos, r ≥ 1.

A v.a. X que conta este número tem distribuição binomial negativa comparâmetros r e p, notando que X ∈ r, r + 1 , r + 2, ....

Se “S” e “F” representam os eventos sucesso e fracasso e X = x, temossequências do tipo

34

sequências do tipo

,SF...SFFSF

ensaios 1 em sucessos 1

44 344 21−− xr

cada uma com probabilidade = pr (1 – p)x – r.

.)!()!1(

)!1(

1

1sequencias de Número

rxr

x

r

x

−−

−=

−=

Page 35: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

Distribuição binomial negativa

Se ensaios de Bernoulli independentes e com probabilidade de sucessoigual a p são realizados, o número de ensaios até que ocorram rsucessos tem uma distribuição binomial negativa com parâmetros r e p.Sua função de probabilidade é dada por

Notação: X ~ BN(r, p).

.10 e ,...2,1, se ,)1(1

1)()( <<++=−

−=== −

prrrxppr

xxXPxf

rxr

35

Notação: X ~ BN(r, p).

Se X ~ BN(r, p), então

E(X) = r / p e

Var(X) = r (1 – p) / p2.

Obs. (a) r = 1: distribuição geométrica na lâmina 29.

(b) Em Excel: função DIST.BIN.NEG.

Page 36: 4. PRINCIPAIS MODELOS DISCRETOS - USPwiki.icmc.usp.br/images/3/3d/ModelosDiscretos.pdfDeterminar a distribuição de probabilidade da variável número de caras nos 3 lançamentos

5 1 0 1 5 2 0 2 5

0.00

0.06

0.12

B N (r = 2 , p = 0 .3 )

x

P(X

= x

)

1 0 2 0 3 0 4 0

0.00

0.04

0.08

B N (r = 4 , p = 0 .3 )

x

P(X

= x

)

Distribuição BN(r, p)

36

x

5 1 0 1 5 2 0 2 5

0.00

0.15

0.30

B N (r = 2 , p = 0 .6 )

x

P(X

= x

)

x

1 0 2 0 3 0 4 0

0.00

0.05

0.10

0.15

B N (r = 4 , p = 0 .5 )

x

P(X

= x

)