21 Análise de Variância
21 Análise de Variância
21 Análise de Variância
2 Tipos de Experimentos
Fatores controláveis a níveis fixos
Ex.: 5 valores de temperatura
Fatores Controláveis a níveis aleatórios
Ex.: 3 fábricas escolhidas ao acaso
EXEMPLO
Um pesquisador deseja investigar o efeito da temperatura do forno
sobre o número de bactérias contadas após o processo de
esterilização. Os dados revelaram o seguinte Estudo experimental em
solados de borracha.
Exemplo (níveis aleatórios)
Um pesquisador
U i d d deseja
j iinvestigar
i se a permeabilidade
bilid d ddas llentes d
de
uso flexível fabricadas em sua indústria permanece uniforme ou não.
Escolhe-se aleatoriamente três lotes de produção e realizam-se
ensaios:
Disposição dos Dados
A Análise de Variância (ANOVA)
• Formulação matemática
do problema:
Modelo Estatístico:
xij = μ + αj + εij
• Hipóteses:
– H0: não há diferenças significativas entre os grupos (α1=α2=α3=αk=0);
– H1: há diferenças significativas entre os grupos.
• D
Decomposição
i ã ddos resíduos:
íd
Σ(x
( ij - x..))2 = Σ(x.
( j - x..))2 + Σ(x
( ij - x.j)2
• Graus de Liberdade:
(N - 1) = (K - 1) + (N - K)
A Análise de Variância (ANOVA)
• Médias Quadradas: Variabilidade Global
=
– MQG = SQG / (K - 1) Variabilidade entre Grupos + Variabilidade dentro do Grupo
– MQR = SQR / (N - K)
Di t ib i ã F
Distribuição
Se não há diferenças significativas entre
os grupos, então: α=0.05
E [MQG] < E [MQR].
[MQR]
Fcrítico (α,GDL1,GDL2)
• Estatística F – Cálculo de Fcal
GDL1-Numerador
GDL1-Denominador
F = MQG / MQR, segue Distribuição F
Fórmulas
Fó l para os cálculos
ál l :
TC = T..2/N
SQT = Σ(xij2) - TC
SQG = Σ(T.j2/nj) - TC
SQR = SQT – SQG
TABELA ANOVA:
EXEMPLO
Um pesquisador
U i d d deseja
j iinvestigar
i o efeito
f i d da temperatura d
do fforno
sobre o número de bactérias contadas em equipamentos cirúrgicos
após o processo de esterilização. Os dados revelaram o seguinte:
EXEMPLO
Cálculos
Cál l Iniciais:
I i i i
TC = T..2 / N = (818,7)2 / 60 = 11.171,1
SQT = Σ(xij2) - TC = 11.608,2
11 608 2 - 11.171,1
11 171 1 = 437
437,11
SQG = Σ(T.j2/nj) - TC = [(206,1)2/12] + ... + [(146,2)2/12] - 11.171,1 =
= 222,3
SQR = SQT - SQG = 437,1 - 222,3 = 214,8
Fonte SQ GDL MQ Teste F
Entre Grupos (Lotes) 222,3 4 55,6 14,2
Dentro dos Grupos (Residual) 214,8 55 3,9
Total 437 1
437,1 59
Æ Há diferenças significativas
entre grupos.
EXEMPLO
4. Usar barras contínuas sobre as médias que não diferem entre sí.
Um pesquisador
U i d d deseja
j iinvestigar
i se a permeabilidade
bilid d ddas llentes d
de
uso flexível fabricadas em sua indústria permanece uniforme ou não.
Escolhe-se aleatoriamente três lotes de produção e realizam-se
ensaios:
Exemplo (níveis aleatórios)
Cálculos
Cál l Iniciais:
I i i i
TC = T..² / N = (1099)² / 18 = 67.100,06
SQT = Σ(xij²)) - TC = 67.163,00
67 163 00 - 67.100,06
67 100 06 = 6262,64
64
SQG = Σ(T.j²/nj) - TC = [(375)2/6] + ... + [(308)2/5] - 67.100,06 = 32,53
SQR = SQT - SQG = 62,94 - 32,53 = 30,41
Æ Há diferenças significativas
entre grupos.
Exemplo (níveis aleatórios)
Pó i
Próximo passo: Estimar
E i componentes de
d variação
i ã
E [MQR] = σ2
E [MQG] = σ2 + ncσα2 ____ σα2 = (E [MQG] - σ2)/nc
2.1)
2 1) Q
Quatro catalisadores
li d que podem
d afetar
f a concentração
ã dde um
componente em uma mistura química estão sendo investigados. As
seguintes concentrações foram obtidas:
Exercícios
P d
Pede-se :
Cálculos
Cál l IIniciais:
i i i
TC = T..2 / N = Σ(xij2)
SQT = Σ(xij2)
( j ) - TC =
SQG = Σ(T.j2/nj) - TC =
SQR = SQT - SQG =
TABELA ANOVA:
Fonte SQ GDL MQ Teste F
p
Entre Grupos
Total
>
Exercícios
Comparação
C ã múltipla
últi l de
d médias
édi
1. Calcular o desvio padrão das médias
sx = √MQR / √nc =
onde nc = (n1 +n2 + ... + nk) / k
2. Calcular o limite de decisão
Ld = 3 x sx =
3. Escrever as médias em ordem crescente ou decrescente e compará-las
duas a duas. A diferença será significativa se for maior que o Ld.
X(1) - X(2) =
X(1) - X(3) =
X(1) - X(4) =
X(2) - X(3) =
X(2) - X(4) =
X(3) - X(4) =
4. Usar barras contínuas sobre as médias que não diferem entre si.
Exercícios
2.2
2 2RResultados
l d d de corpos d
de prova d
de concreto com adição
di ã dde
Microssílica indicaram os seguintes resultados de resistência à
compressão:
Adição Resistência (Mpa)
0% 28,1 26,5 24,3 23,8 28,5
5% 35,3 34,3 37,5 38,0 33,9
10% 39,8 44,1 42,3 39,2 44,8
15% 39,1 40,8 43,0 40,1 43,5
2.4
2 4UUma metalúrgica
lú i tem um grande d número
ú d
de ffornos usados
d para
fundição de metais. A temperatura desses fornos deveria ser a mesma.
Para testar essa hipótese foram feitas medições em 4 fornos escolhidos
aleatoriamentes. Analise os resultados e conclua a respeito de
possíveis diferenças entre os fornos.
2.5
2 5UUm engenheiro
h i iindustrial
d ti ld desenvolveu
l um modelo
d l estocástico
t á ti d de
simulação que prevê a produtividade mensal em função do intervalo de
tempo entre manutenções preventivas. Se esse intervalo for muito curto,
as máquinas estarão constantemente em manutenção e a produtividade
será baixa. Se o intervalo for muito longo, haverá quebras, exigindo
manutenção corretiva, mais demorada, novamente prejudicando a
produtividade Os resultados da simulação aparecem a seguir
produtividade. seguir.
Intervalo Temperatura (°C)
4 136 137 135 140 136
6 145 146 147 147 147
8 146 144 148 145 145
10 134 131 136 134 133
12 117 119 117 115 116
2.6
2 6EEm uma iindústria
dú i química
í i um catalisador
li d é utilizado
ili d para acelerar
l
um processo de deposição metálica. Foi feito um experimento
variando-se a concentração desse catalisador e anotando-se o tempo
necessário para complementar o processo. Analise os dados usando a
Tabela Anova. Depois faça uma comparação múltipla de médias, plote
um ggráfico de linhas e conclua a respeito
p da concentração
ç ideal.
Concentração Tempos
25 11,2 10,4 10,1 10,6
20 10,6 11,1 10,8 11,7
15 12 5
12,5 12 0
12,0 13 2
13,2 12 6
12,6
10 18,8 19,0 18,4 19,6
Projetos Fatoriais com
D i Fatores
Dois F t
Os Experimentos Fatoriais de Dois Fatores
(Two-way Anova)
•Os
O experimentos
i t fatoriais
f t i i maisi simples
i l envolveml d
dois
i ffatores;
t
•Fator A com “a” níveis e Fator B com “b” níveis.
•Cada repetição
p completa
p do experimento
p envolve “ab” ensaios.
1 2 ... b
: : : :
Resposta
60
50
40
30
20
10
0
A1 A2
B1 B2
Projetos Fatoriais com Dois Fatores
Com Interação
Em alguns
E l experimentos
i a dif
diferença na resposta observada
b d quando d se
modifica os níveis de um dos fatores irá depender do nível do outro
fator. Por exemplo:
Resposta
p
60
50
40
30
20
10
0
A1 A2
B1 B2
Modelo Estatístico
M d l Estatístico:
Modelo E í i
• S
Suposições
i õ
εijk → N(0,σ)
• Hipoteses:
– Para o Fator A:
• H0: τi = 0;
• H1: τi ≠ 0 p
para algum
g i.
– Para o Fator B:
• H0: βj = 0;
• H1: βj ≠ 0 para algum j.j
– Para a interação AB:
• H0: βτij = 0;
• H1 βτij
H1: β ij ≠ 0 para algum
l ij
ij.
Formulário para os Cálculos da Significância
de A, B, AB
• TC = (T...)²
abn
Verificação:
• SQA = Σ(Ti..)
Σ(Ti )² - TC
SQT = SQA + SQB + SQAB +SQR
bn
• SQB = Σ(T.j.)² - TC
an
• SQAB = Σ(Tij.)² - TC - SQA - SQB
n
• SQR = Σ(yijk)²
(y j ) - Σ(Ti..)²
( )
n
• SQT = Σ(yijk)² - TC
Análise de Variância para Projetos
Cruzados de 2 Fatores
Suspeita-se
S i que a máxima
á i voltagem
l d
de saída
íd dde um tipo
i d de b
bateria
i é
afetada pelo material usado nas placas e pela temperatura. Quatro
repetições completas de um experimento fatorial completo foram
rodadas em laboratório e os seguintes dados foram obtidos:
Material Temperatura (B) Ti..
(A) 50 65 80
Erro
o 18644
86 27 691
69
Total 71611 35
200
•O efeito do Material é significativo; 150
Material1
oltagem
•O efeito da Temperatura é significativo; 100 Material 2
Material 3
Vo
•O efeito da interação é significativo. 50
0
50 65 80
T
Temperatura
t
Comparação Múltipla de Médias (CMM)
• Se há efeitos
S f i significativos,
i ifi i em gerall procede-se
d a uma CMM.
CMM
• Como a interação é significativa, as comparações devem ser feitas
fixando-se o nível de um dos fatores.
• Por exemplo, podemos investigar se há DS entre as médias
obtidas com os três tipos de materiais para a Temperatura de 65°C.
• Li i
Limites de
d decisão
d i ã
– Ld = 3 Sx = 39,3
• Comparação duas a duas:
εijk → N(0,σ)
Comparação Múltipla de Médias (CMM)
•Experimento
E i sem Repetição
R i ã
–Lembrando,
Lembrando, o número de GDL do termo de erro vem dado por: ab(n – 1)
–Se não há repetições do experimento, isto é, se n = 1, não sobram GDL
para calcular de modo independente a MQR.
•Tabela Anova
Um pesquisador
U i d acredita
dit que a resistência
i tê i à ttração
ã dde certos
t corpos dde prova d
de
argamassa depende da % de microssílica utilizada na sua fabricação e do operador
que confecciona os CPs. Os dados revelaram:
% Microssílica
Operador 0 5 10 15 20 Totais
1 4 5 6 5 3 23
2 1 3 4 3 2 13
3 1 1 3 2 1 8
Totais 6 9 13 10 6 44
% de Madeira
0 5 10 15 Totais
Escura
5% 14 13 23 27 32 33 34 38
10% 31 28 38 39 43 40 41 39
Totais 514
Pergunta-se:
a) Qual a variável de resposta?
b) Quais os fatores controláveis e Qual o número de níveis dos fatores controláveis?
c) Faça a análise de variância e indique quais os efeitos significativos?
d) Faça um gráfico de dois fatores
e) O que fazer oara assegurar a qualidade? (Resistência maior-é-melhor)
maior é melhor) O que fazer para
assegurar economia? (Supor que um aumento na % de madeira escura ou na
temperatura dos rolos implica maior custo)
Generalização dos Projetos
Fatoriais
(Multi-factor Anova)
Generalização dos Projetos Fatoriais
• Os resultados
O l d d do P
Projeto
j F
Fatorial
i ldde d
dois
i ffatores podem
d ser
estendidos para o caso onde há vários fatores.
-Fator A, a níveis
-Fator B, b níveis
-Fator C, c níveis
:
-nn observações por parcela
• Por exemplo,
P l sejaj um projeto
j t ffatorial
t i l com ttrês
ê ffatores:
t
• Modelo Estatístico:
• S
Suposições
i õ
εijkl → N(0,σ)
• Hipoteses:
p
• Para o Fator A:
– H0: τi = 0;
– H1: τi ≠ 0 para algum i.
:
• Para o Fator AB:
H0: τβij
– H0 βij = 0
0;
– H1: τβij ≠ 0 para algum ij.
:
• Para a interação ABC:
– H0: βτγijk = 0;
– H1: βτiγjk ≠ 0 para algum ij.
Formulário para os Cálculos da Significância
de A, B, AB
Æ TC = (T...)2
(T )2 Æ SQAC = Σ(Ti.k.)² - TC - SQA - SQC
abcn bn
Æ SQA = Σ(Ti...)²
( ) - TC Æ SQBC = Σ(T.jk.)²
( j ) - TC - SQB - SQC
bcn an
Æ SQB = Σ(T.j..)² - TC Æ SQR = Σ(yijkl)² - Σ(Tijk.)²
acn n
Æ SQC = Σ(T..k.)² - TC Æ SQT = Σ(yijkl)2 - TC
abn
Æ SQAB = Σ(Tij.)² - TC - SQA - SQB
cn
Æ SQABC = Σ(Tijk.)
Σ(Tijk.)² - TC - SQA - SQB - SQC - SQAB - SQAC - SQBC
n
Verificação:
SQT = SQA + SQB + SQAB + ... + SQR
Análise de Variância para Projetos
Cruzados de 2 Fatores
F t
Fonte SQ GDL MQ T t F
Teste
A SQA a-1 MQA MQA / MQR
Um ffabricante
U bi t dde refrigerantes
fi t está
tá estudando
t d d o efeito
f it da
d %d de carbonatação
b t ã (A)(A),
pressão de enchimento (B) e velocidade da linha (C) sobre o volume do refrigerante. Os
dados revelaram:
Pressão de Enchimento (C)
25 psi 30 psi
T b l Anova
Tabela A
C 22 04
22,04 1 22 04
22,04 31 1*
31,1*
AC 0,58
, 2 0,29
, 0,4
,
Total 336,63 23
Exercícios
4.2
4 2SSupõe-se
õ que a tensão
ã dde cisalhamento
i lh suportada
d por peças
coadas depende do fornecedor de adesivo e da pressão e temperatura
usadas no processo de colagem. Analise os dados a seguir,
respondendo as mesmas questões enunciadas no exercício 4.1
Fornec.: 1 2
Temp: 250 260 270 250 260 270
Pressões:
4.3
4 3SSabe-se
b que a tensão
ã dde cisalhamento
i lh suportada
d por peças
coladas depende do fornecedor de adesivos da temperatura usadas no
processo de colagem. Além disso, diferentes fornecedores sugerem
diferentes temperaturas ótimas de colagem. Indique qual o modelo
estatístico desse experimento e, depois, analise os dados a seguir,
respondendo
p as mesmas q questões enunciadas no exercício 4.1
Fornecedor 1 2
Temperatura 250 260 270 250 260 270