An engineer is interested in the effect of cutting speed
, metal hardness
, and cutting angle
on the life of a cutting tool. Two levels of each factor are chosen, and two replicates of a
factorial design are run. The tool life data (in hours) are shown in the following table.
Treatment
Combination
Replicate
I
II
221
311
325
435
354
348
552
472
440
453
406
377
605
500
392
419
(a) Analyze the data from this experiment.
(b) Find an appropriate regression model that explains tool life in terms of the variables used in the experiment.
(c) Analyze the residuals from this experiment.
|
Replicate | |||
| I | II | |||
|
|
221 | 311 | ||
|
|
325 | 435 | ||
|
|
354 | 348 | ||
|
|
552 | 472 | ||
|
|
440 | 453 | ||
|
|
406 | 377 | ||
|
|
605 | 500 | ||
|
|
392 | 419 | ||
(b) Find an appropriate regression model that explains tool life in terms of the variables used in the experiment.
(c) Analyze the residuals from this experiment.
Passo 1
Vou traduzir claramente o que o enunciado pede e listar os dados fornecidos (tudo retirado da tabela que você mandou). Temos três fatores binários: velocidade de corte
, dureza do metal
e ângulo de corte
. Foram feitas duas réplicas para cada uma das
combinações de níveis. O que se pede:
- (a) Analisar os dados (ou seja, identificar efeitos importantes, fazer ANOVA).
- (b) Encontrar um modelo de regressão apropriado para a vida da ferramenta em função dos fatores.
- (c) Analisar os resíduos do experimento.
Dados (cada par é réplica I e II; eu peguei os valores da tabela do enunciado):
Observações importantes extraídas da tabela:
- Cada combinação tem duas observações (duas réplicas).
- Notação: por exemplo,
= nível alto de (os outros fatores no nível baixo), = e no nível alto, = no nível alto, etc.
A partir daqui a resolução começa com o cálculo das médias, efeitos, ANOVA e análise de resíduos.
Passo 2
Vou transformar os dados em médias por combinação (útil porque há nº igual de réplicas). Calculei a média de cada tratamento (peguei os dois valores da tabela e fiz a média):
Agora o valor médio geral (grand mean), que é a média de todas as 16 observações (equivalente à média das 8 médias pois há mesmas réplicas):
Passo 3
Explico como medir efeitos num projeto
e calculo os efeitos principais e de interação. Convenção: codificamos níveis baixos como
e níveis altos como
. A estimativa do efeito (tamanho da mudança de
para
) para um fator num
com médias por tratamento é dada por
onde
então o fator é
vezes a soma com sinais
conforme presença do fator. Vou listar os sinais para cada tratamento (ordem:
):
- Para
: . - Para
: . - Para
: . - Para
: . - Para
: . - Para
: . - Para
: .
Faço a soma ponderada (
) e divido por
. Cálculos (mostro as somas e o efeito final):
Interações:
Interpretação didática rápida:
- O efeito de
é grande e positivo: quando passa de baixo para alto, a vida média aumenta cerca de horas. - O efeito de
também é grande ( ). - A interação
é grande em magnitude e negativa ( ), indicando que o efeito de depende fortemente do nível de (ou vice-versa). - Os outros efeitos são menores.
Passo 4
Ligando esses efeitos ao modelo de regressão em variáveis codificadas
:
- No modelo linear com todas as interações possíveis,
Relação entre efeitos (definidos antes) e os coeficientes
: com codificação
, a diferença entre
e
é
. Assim
Calculo os coeficientes (cada
é metade do efeito):
Passo 5
Agora construo a ANOVA a partir dos contrastes (forma padrão para
). Primeiro obtenho o somatório dos quadrados de erro “pure error” pela variação entre réplicas (cada tratamento tem 2 réplicas). Para cada tratamento com medidas
,
Calculei a soma destes valores para os 8 tratamentos (peguei as diferenças das réplicas da tabela):
Logo
Para cada efeito podemos calcular a soma de quadrados a partir do contraste sobre as 16 observações. Fórmula prática:
onde
(soma sobre as 16 observações) e
. Como eu já tinha
para as 8 médias, o contraste sobre as 16 observações é o dobro daquele valor (pois há 2 réplicas):
. Usei isso para calcular as somas de quadrados. Resultados:
Total regressão (soma desses SS)
. Soma total
regressão
erro
(com
graus de liberdade).
Faço o teste F para cada efeito (cada um tem 1 grau de liberdade);
:
Tomando nível de significância
(teste usual), o valor crítico
. Portanto os efeitos significativos são:
-
(pois ), -
(pois ), -
(porque ).
Os demais efeitos não são estatisticamente significativos ao nível
.
Passo 6
A partir da ANOVA, construo um modelo reduzido contendo os termos significativos (método usual: manter intercepto e os termos que deram significância). O modelo reduzido em codificação
é:
Substituindo os coeficientes calculados:
Explicação do que significa: por exemplo,
quando
está no nível alto;
é o produto dos codificados e representa a interação entre
e
. Se você quiser a mudança de
para
em
(mantendo os outros fixos), o ganho esperado seria
horas, que é o efeito que calculamos.
Passo 7
Agora analiso os resíduos do modelo reduzido. Primeiro prevejo
para cada combinação de níveis (mesma predição para as duas réplicas da mesma combinação) e calculo os resíduos
. Usei o modelo reduzido acima. Resultados (mostro tratamento,
e resíduos nas duas réplicas; peguei os
da sua tabela):
Calculei a soma dos quadrados dos resíduos do modelo reduzido:
Grau de liberdade do erro no modelo reduzido:
(quatro parâmetros: intercepto,
). Portanto
Passo 8
Verificação de adequação do modelo reduzido: faço um teste F parcial para ver se os efeitos que retirei do modelo completo (os quatro termos:
) são coletivamente significativos. A diferença de SSE entre modelo reduzido e completo é:
com
. O F de comparação é
Esse
é muito menor que o crítico
(que é bem maior que 1), portanto não há evidência de que os termos removidos tragam melhoria significativa ao ajuste. Conclusão: o modelo reduzido é adequado e parcimonioso.
Passo 9
Análise dos resíduos em termos práticos (o que você precisa checar na prova):
- Soma dos resíduos é zero (verificação básica de ajuste com intercepto) — confirmei que a soma dos resíduos é exatamente zero.
- Verifique outliers: calculei os resíduos padronizados dividindo cada resíduo pela raiz de
( ). Os resíduos padronizados foram aproximadamente:
Nenhum desses excede em magnitude
; portanto não há outliers evidentes.
- Homocedasticidade (variância constante): os resíduos são de magnitudes parecidas (não há tendência clara de aumento da variância com o valor previsto). Com 16 pontos não se pode ser muito exigente, mas inspeção numérica não mostra heterocedasticidade forte.
- Normalidade: com 16 observações não posso provar normalidade, mas os resíduos padronizados estão razoavelmente centrados e sem outliers extremos; um gráfico QQ (se você fizesse) provavelmente não mostraria grandes desvios. Logo, nenhuma evidência forte contra supostos normais para
. - Conclusão prática sobre resíduos: não há motivos fortes para transformação; o modelo reduzido é aceitável.
Passo 10
Interpretação final e recomendações para a prova (resumo didático):
- Os fatores que mais influenciam a vida da ferramenta são
(dureza) e (ângulo), e existe uma interação forte entre (velocidade) e . Isso significa que o efeito da velocidade sobre a vida depende do ângulo. - O modelo em variáveis codificadas
que explica bem os dados é
onde
.
- Os resíduos do modelo não mostram problemas evidentes (sem outliers extremos, variância razoavelmente constante, soma dos resíduos zero). Um teste parcial confirmado que os termos removidos do modelo completo não são necessários.
Resposta
(a) Os efeitos estatisticamente significativos ao nível de
são
,
e a interação
.
(b) Modelo adequado (em codificação ):
(b) Modelo adequado (em codificação
(c) Análise de resíduos:
com
(
); resíduos padronizados dentro de
aproximadamente, sem outliers extremos e sem evidência forte de heterocedasticidade ou não-normalidade; o modelo reduzido é adequado (teste parcial para os termos removidos deu
, não significativo).