Resumo de Distribuições Minitab
Resumo de Distribuições Minitab
Resumo de Distribuições Minitab
Frases
Roteiro
1. Distribuições de Probabilidade
2. Aproximação da Binomial
3. Geração de Números Aleatórios
4. Ajuste de Distribuições
5. Referências
1
Distribuições de Probabilidade
Funcionalidades
2
Distribuições Discretas
Minitab: Outras:
• Binomial • Binomial Negativa
• Hipergeométrica • Geométrica
• Discreta • Bernoulli
• Inteira
• Poisson
P{X ≤ a} = p
Quantil referente a
uma probabilidade
Valores de entrada
Distribuição Binomial
• Variância: np (1 − p )
3
Distribuição Binomial – Exemplo
Montar tabela com os valores de distribuição
binomial com parâmetros n = 30 e p = 0,4
• Para gerar os valores dos pontos (k):
Calc > Make Patterned Data > Simple Set of Numbersè
Valores de k
Parâmetros
Valores de entrada
Valores de saída
Parâmetros
Valores de entrada
Valores de saída
4
Binomial (30; 0,4) – Tabela
k P(X=k) P(X<=k)
0 0,000 0,000
1 0,000 0,000
2 0,000 0,000
3 0,000 0,000
4 0,001 0,002
5 0,004 0,006
6 0,012 0,017
7 0,026 0,044
8 0,050 0,094
9 0,082 0,176
10 0,115 0,291
11 0,140 0,431
12 0,147 0,578
13 0,136 0,715
14 0,110 0,825
15 0,078 0,903
16 0,049 0,952
17 0,027 0,979
18 0,013 0,992
19 0,005 0,997
20 0,002 0,999
21 0,001 1,000
22 0,000 1,000
23 0,000 1,000
0,14
0,12
0,10
P(X=k)
0,08
0,06
0,04
0,02
0,00
0 5 10 15 20 25 30
k
5
Quantis
Quantis da binomial(30; 0,4) para 0,056;
0,058; 0,060; ... ; 0,0672
• Para gerar os valores das probabilidades
Calc > Make Patterned Data > Simple Set of Numbersè
Valores de
probabilidade
Parâmetros
Valores de entrada
Valores de saída
P{X ≤ a} = p
0,058 8
0,060 8
0,062 8
0,064 8
0,066 8
0,068 8
0,070 8
0,072 8
0,074 8
0,076 8
0,078 8
0,080 8
0,082 8
0,084 8
0,086 8
0,088 8
0,090 8
0,092 8
0,094 8
0,096 9
0,098 9
0,100 9
0,102 9
6
Distribuições Contínuas
• Qui-quadrado • Beta
• Normal • Exponencial
• F • Gama
• t de Student • Valor extremo
• Uniforme • Lognormal
• Logística
• Weibull
• outras
P{X ≤ a} = p
Quantil referente a
uma probabilidade
Valores de entrada
Distribuição Normal
• Notação: X ~ N(µ, s 2)
1 1 x − µ 2
• Função de densidade: f ( x) = exp−
2π σ 2 σ
• Esperança: µ
• Variância: s 2
7
Cálculo de Probabilidades
• Para valores únicos, são mais fáceis de serem
obtidos através da janela Session:
√ Editor > Enable Commands
• Para cálculo de densidade:
PDF valor do quantil;
nome distribuição parâmetro 1 parâmetro 2.
• Para cálculo da função de distribuição acumulada
CDF valor do quantil;
nome distribuição parâmetro 1 parâmetro 2.
• Para cálculo de quantil:
InvCDF valor da probabilidade;
nome distribuição parâmetro 1 parâmetro 2.
Normal Padrão
MTB > PDF -1;
SUBC> Normal 0,0 1,0.
x f( x )
-1 0,241971
P{X ≤ −1}
SUBC> Normal 0,0 1,0.
x P( X <= x )
-1 0,158655
P{X ≤ a} = 0,05
SUBC> Normal 0,0 1,0.
P( X <= x ) x
0,05 -1,64485
X ~ N (100,100)
MTB > cdf 95;
SUBC> normal 100 10.
x P( X <= x )
95 0,308538
= 1 – 0,158655
P{X > 95} = 0,841345
8
Gráficos Densidade Normal
• Sobrepor gráficos das funções de densidade
N (100, 100), N (100, 25) e N(90, 25).
• Criar coluna das abcissas X:
Calc > Make Patterned Data > Simple Set of Numbersè
Valores de X
Parâmetros
Valores de entrada
Valores de saída
9
Funções de Densidade Normal
0,09 Variable
N(100, 100)
0,08 N(100, 25)
N(90, 25)
Densidade de Probabilidade
0,07
0,06
0,05
0,04
0,03
0,02
0,01
0,00
Parâmetro de locação: µ
Parâmetro de escala: s 2
Parâmetros
Valores de entrada
Valores de saída
10
• Repetir o procedimento de cálculo de
probabilidades para Ac N(100, 25) e Ac N(90, 25).
• Fazer gráfico sobreposto:
Graph > Scatterplot > Simple è
0,6
0,4
0,2
0,0
Distribuição Exponencial
11
Distribuição Exponencial
• Variância: β 2
Cálculos de Probabilidade
P{X > 2}
SUBC> exponencial 0,5.
P( X <= x ) x
0,5 0,346574
12
Variáveis Discretas – Probabilidade
P{ X = 7 }
Probability Density Function
x P( X = x )
7 0,0212030
P{ X < 8 }
SUBC> binomial 10 0,35.
x P( X <= x )
7 0,995179
= P{ 3 < X ≤ 6 }
Binomial with n = 10 and p = 0,35
x P( X <= x )
3 0,513827
= 0,513827 –0,973976 6 0,973976
= 0,460149
Função Gama
∞
• Dada por: Γ(r ) = ∫ x r −1e − x dx
0
13
Exemplo Função Gama
Data Display
K1 1,16671
Distribuição de Weibull
• Notação: X ~ Weibull (ß, d)
β x β −1 − x
β
δ
x≥0
• Função de densidade: f ( x) = δ δ e
0 c.c.
• Esperança: δ Γ(1+ 1
β )
• Variância: δ 2 {Γ(1 + 2 β ) − [Γ(1 + 1
β )]2 }
• Parâmetro de forma: ß
• Parâmetro de escala: d
First value: 0
Last value: 15
In steps of: 0,1
14
• Calcular as probabilidades das colunas
W(1, 1), W(3,4; 2) e W(4,5; 6,2):
Calc > Probability Distributions > Weibull è
Parâmetros
Valores de entrada
Valores de saída
0,7
W(ß, d)
0,6
0,5
0,4
0,3
0,2
0,1
0,0
0 2 4 6 8 10 12 14 16
X_w
Parâmetro de forma: ß
Parâmetro de escala: d
15
Distribuição Gama
• Notação: X ~ Gama (r, ?)
λr x r −1e − λx
x>0
• Função de densidade: f ( x) = Γ ( r )
0 c.c.
• Esperança: r
λ
r
• Variância:
λ2
• Parâmetro de forma: r
• Parâmetro de escala: ?
• Variância: rβ 2
• Parâmetro de forma: r
• Parâmetro de escala: ß
First value: 0
Last value: 6
In steps of: 0,1
16
• Calcular as probabilidades das colunas
G(1, 1), G(2; 0,5) e G(2; 1):
Calc > Probability Distributions > Weibull è
Parâmetros
Valores de entrada
Valores de saída
Densidades Gama
Variable
0,9 G(1, 1)
G(2; 0,5)
0,8 G(2, 1)
Densidade de Probabilidade
0,7
G(r, ß)
0,6
0,5
0,4
0,3
0,2
0,1
0,0
0 1 2 3 4 5 6
X_g
Parâmetro de forma: r
Parâmetro de escala: ß
17
Aproximação da Binomial
• Se n ? 8 e p ? 0, com np ? ? então a
variável aleatória XB ~ Binomial (n, p) pode
ser aproximada por XP ~ Poisson (?), neste
caso ? ˜ np.
Exemplo
18
• Gerar coluna das abcissas X_b:
Calc > Make Patterned Data > Simple Set of Numbersè
First value: 0
Last value: 70
In steps of: 1
• Calcular as probabilidade binomiais
19
Distribuição Acumulada Exata e Aproximada Função de Probabilidade Exata e Aproximada
Variable
1,0 0,5
P(X_b=x)
P(X_p=x)
0,9 0,4
0,8 0,3
Y-Data
Y-Data
0,7 0,2
0,6 0,1
0,5 0,0
0 2 4 6 8 10 0 1 2 3 4 5
x_b x_b
Exemplo
• Seja XB ~ binomial (70; 0,98)
• Calcular as seguintes probabilidades exatas,
comparando-as com sua aproximação pela
Poisson
• P{XB = 68}
• P{XB = 65}
• P{65 < XB = 68}
20
XB ~ binomial (70; 0,98)
MTB > pdf 68;
SUBC> binomial 70 0,98.
P{ X B = 68 }
Probability Density Function
x P( X = x )
68 0,244540
P{ X B ≤ 65 }
Cumulative Distribution Function
x P( X <= x )
65 0,0132298
P{ 65 < X B ≤ 68 }
SUBC> binomial 70 0,98.
x P( X = x )
2 0,241665
P{ X B ≤ 65 }
Poisson with mean = 1,4
≈ P{ YP ≥ 5 } = 1 − P{YP ≤ 4}
x P( X <= x )
4 0,985747
1 0,591833
= 1 – 0,985747
= 0,014253
P{ 65 < X B ≤ 68 }
≈ P{ 2 ≤ YP < 5 } = P{1 < YP ≤ 4} = 0,985747 – 0,591833
= 0,393914
Comparação
Aproximação
Probabilidades Exata
Poisson
21
Aproximação pela Normal
Correção de Continuidade
22
Comparação
Aproximação
Probabilidades Exata
Poisson
First value: 0
Last value: 100
In steps of: 1
• Calcular na coluna
XB(100; 0,5) as
probabilidades acumuladas
binomiais (exatas)
• Calcular na coluna
XN(50; 25) as
probabilidades acumuladas
aproximadas pela normal
(com correção de
continuidade)
23
• Construir gráfico das distribuições acumuladas
sobrepostas XB(100, 0,5) e XN(50, 25)
1,0
0,8
Probabilidade A cumulada
0,6
0,4
Variable
XB(100,0,5) * X_b
XN(50,25) * X_co
0,2
0,0
40 45 50 55 60
Valores de X
Abcissas entre 40 e 60
1,0
0,8
0,6
Y-Data
0,4
Variable
0,2 P(X_b<=k) * x_b
P(X_p<=x) * x_b
XN(7;0,693) * x_b
0,0
0 2 4 6 8 10
X-Data
Abcissas entre 0 e 10
24
Comparação entre Aproximações (2)
1,0
0,8
0,6
Y-Data
0,4 Variable
XB(100; 0,5) * x_b
XP(50) * x_bn
XN(50,25) * x_cor
0,2
0,0
30 40 50 60 70
X-Data
Abcissas entre 30 e 70
Aproximações – Comentários
Mantidas as condições de qualidade de cada
aproximação, temos:
• A aproximação de Poisson é melhor que a
aproximação normal para valores de p nas
proximidades de 0 e 1.
• A aproximação Normal é melhor que a
aproximação de Poisson para valores
intermediários de p.
25
Geração de Números Aleatórios
26
Exemplo 1 – Geração de Normais
Valores exatos:
Média: 86,6
Desvio padrão: 4,6
27
Construção Histogramas
Histogramas de N1 e N2
N1
0,08
0,06
Densidade de Probabilidade
0,04
0,02
0,00
N2
0,08
0,06
0,04
0,02
0,00
78 81 84 87 90 93 96 99
28
Histogramas de N1 e N2
Normal
N1
0,08
Densidade de Probabilidade
0,06
0,04
0,02
0,00
N2
0,08
0,06
0,04
0,02
0,00
78 81 84 87 90 93 96 99
29
Média e Desvio Padrão Amostrais
MTB > describe 'E1'-'E3';
SUBC> mean;
SUBC> StDev.
Valores exatos:
Média: 128
Desvio padrão: 128
Histogramas de E1, E2 e E3
0 80 160 240 320 400 480 560
E1 E2
0,0048
0,0036
0,0024
0,0012
Density
0,0000
E3
0,0048
0,0036
0,0024
0,0012
0,0000
0 80 160 240 320 400 480 560
30
Histogramas de E1, E2 e E3
Exponential
0 80 160 240 320 400 480 560
E1 E2
0,008
0,006
0,004
0,002
Density
0,000
E3
0,008
0,006
0,004
0,002
0,000
0 80 160 240 320 400 480 560
31
Média e Desvio Padrão Amostrais
MTB > describe 'W1' 'W2';
SUBC> mean;
SUBC> StDev.
K3 131,115
Data Display
K5 60,4552
Histogramas de W1 e W2
W1
0,008
0,006
0,004
0,002
Density
0,000
W2
0,008
0,006
0,004
0,002
0,000
0 60 120 180 240
32
Histogramas de W1 e W2
Weibull
W1
0,008
0,006
0,004
0,002
Density
0,000
W2
0,008
0,006
0,004
0,002
0,000
0 60 120 180 240
Ajuste de Distribuições
Verificação de Normalidade
33
• Teste de normalidade da coluna N1:
Graph > Probability Plot > Single è
Probability Plot of N1
Normal - 95% CI Parâmetros estimados
99
Mean 86,87
StDev 4,235
95 N 29 Estatística
AD 0,149
90
P-Value 0,959 de teste
80
70
Percent
60
50
40 P-valor p/ o teste:
30
20 H0: dados normais
10
5
H1: Dados não-normais
1
75 80 85 90 95 100
N1
60
50
40 P-valor p/ o teste:
30
20
H0: dados normais
10
34
Exemplo 4 – Ajuste de Modelo
• Repetir o teste de normalidade para as colunas
E2, E3, W1 e W2
Graph > Probability Plot > Single è
1 1 P-Value <0,005
-200 0 200 400 -200 0 200 400 600 W1
Mean 134,9
W1 W2
99 99 StDev 63,81
N 43
90 90 AD 1,281
P-Value <0,005
50 50 W2
Mean 144,5
10 10 StDev 56,17
N 43
1 1 AD 0,543
Sem evidência pra
-100 0 100 200 300 0 100 200 300 P-Value 0,154
rejeitar H0
35
Probability Plot of E1; E2; E3 H0: modelo exponencial
Exponential - 95% CI
H1: outro modelo
E1 E2 E1
Mean 116,3
90 90
N 34
50 50 AD 0,621
P-Value 0,341
E2
10 10 Mean 119,4
N 34 Não há evidências
AD 1,253
Percent
1
1 10 100 1000
1
0,01 0,10 1,00 10,00 100,00 1000,00
P-Value 0,057 para rejeitar H0
E3
E3
Mean 124,1
90 N 34
AD 0,935
50 P-Value 0,136
10
1
1 10 100 1000
10 100 1000
H1: outro modelo
W1 W2 W1
99 99
Shape 2,263
Scale 151,8
90
80
90
80
N 43 Há evidências
AD 1,521
70
60
70
60 P-Value <0,010 para rejeitar H0
50 50 W2
40 40 Shape 2,786
Percent
30 30 Scale 161,1
20 20 N 43
AD 0,678 Não há evidências
10 10 P-Value 0,074
para rejeitar H0
5 5
3 3
2 2
1 1
10 100 1000
36
Exemplo 5 – Amostra Binomial
Padronização
• No caso da binomial a padronização exata é
dada por: xi − np
zi =
np(1 − p)
Data Display
K1 28,0000
K2 4,09878
37
• No Minitab, a padronização é obtida por:
Calc > Standardize è
1 1 StDev 1,065
-4 -2 0 2 4 -2 0 2 4 N 25
Bp3
AD
P-Value
0,492
0,199
Não há evidências
99
90
Bp3 para rejeitar H0
Mean 0,02928
StDev 1,237
50 N 25
AD 0,235
10 P-Value 0,768
1
-4 -2 0 2 4
38
Referências
Bibliografia Recomendada
39