Transcript Introdução aos experimentos fatoriais
Introdução aos experimentos fatoriais
Blocagem – cap. 5 (5.6)
Blocagem em um experimento fatorial: modelo • Considere um experimento fatorial a dois fatores (A e B) com é:
n
replicações. O modelo estatístico linear para esse experimento
y ijk
i j
( )
ij
ijk
i k j
1, 2,..., 1, 2,..., 1, 2,...,
a b n
Suponha que para uma realização um material particular é exigido e que esse material está disponível em lotes que não são suficientemente grandes para permitir que todas as
abn
combinações de tratamentos sejam realizadas com o mesmo lote.
Porém, se um lote contém material suficiente para
ab
realizações, então um plano alternativo é rodar cada uma das
n
replicações usando um lote diferente de material.
Blocagem em um experimento fatorial: modelo • Consequentemente, os lotes de material representam uma restrição em aleatorização ou um
bloco,
e uma replicação de um experimento fatorial completo é realizada para cada bloco.
• O modelo de efeitos para esse novo plano é:
y ijk
i
j
( )
ij
k
ijk
,
i k j
1 , 2 ,...,
a
1 , 2 ,..., 1 , 2 ,...,
b n
com
k
representando o efeito do
k
ésimo bloco.
Dentro de cada bloco, é claro, as realizações são feitas de modo completamente aleatorizado.
Tabela ANOVA
FV Blocos A B AB Erro Total SQ 1
ab
k y
2 ..
k
y
2 ...
abn
1
bn
i y i
2 ..
y
2 ...
abn
1
an
j y
.
2
j
.
y
2 ...
abn
1
n i
j
2
y ij
.
y
2 ...
abn
SS A
SS B
Diferença
y ijk
y
2 ...
abn gl n-1 a-1 b-1 (a-1)(b-1) (ab-1)(n-1) abn-1
QM QMBl QMa QMb QMab QMe - F - - --
Blocagem
• O modelo assume que a interação entre blocos e tratamentos é desprezível. Isso foi suposto anteriormente na análise dos planejamentos em bloco aleatorizados. Se essas interações de fato existem, elas não podem ser separadas da componente de erro.
• De fato, o termo de erro nesse modelo consiste das interações entre cada fator principal e bloco e entre os três (fatores A, B e bloco).
Exercício 19
O resultado de um processo químico está sendo estudado. Os dois fatores de interesse são temperatura e pressão. Três níveis de cada fator foram selecionados.
Porém, somente nove realizações podem ser feitas num dia. O experimentador rodou as replicações em dias diferentes. Analise os dados, supondo que os dias são blocos.
y=read.table("e:\\dox\\procquim.txt",header=T) pr=as.factor(y$Pressure) temp=as.factor(y$Temperature) bloco=as.factor(y$dia) modeloB=y$Yield~pr+temp+pr:temp+bloco fitB=aov(modeloB) fv gl SQ QM F value Pr(>F) pr 2 5.508 2.754 5.1838 0.035988 * temp 2 99.854 49.927 93.9807 2.778e-06 *** bloco 1 13.005 13.005 24.4800 0.001124 ** pr:temp 4 4.452 1.113 2.0952 0.173314 Residuals 8 4.250 0.531 -- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Perfis das médias por combinação de níveis apresentam paralelismo, confirmando a não-rejeição da hipótese de ausência de efeito de Interação entre temperatura e pressão. 94 92 90 88 86 84 82 80 Low Medium High 250 260 270
Experimentos fatoriais 2
k
– Cap. 6
• Os experimentos fatoriais são muito usados em experimentos envolvendo vários fatores para os quais é necessário estudar o efeito conjunto dos fatores sobre a resposta.
• No capítulo 5 apresentamos métodos gerais de análise do experimento fatorial.
• Existem casos especiais do experimento fatorial geral.
• O caso mais importante, entre os especiais, é o de
k
um com apenas 2 níveis.
fatores cada • Os níveis podem ser quantitativos ou qualitativos.
• O plano 2
k
é particularmente útil nos estágios iniciais do trabalho experimental, quando muitos fatores deverão ser investigados.
• Ele fornece o menor número de realizações com o qual fatorial completo.
k
fatores podem ser simultaneamente investigados em um planejamento
Experimentos fatoriais 2
k • Consequentemente, esses planos são muito usados em experimentos chamados
factor screening experiments
(filtragem, peneiramento de fatores).
• Como há somente dois níveis para cada fator, supõe-se que a resposta é aproximadamente linear sobre a variedade de níveis dos fatores escolhidos.
• Em muitos experimentos de filtragem de fatores, quando estamos apenas começando a estudar o processo, essa suposição costuma ser razoável.
• Na seção 6.8 são apresentados um método simples para verificar se essa suposição é violada e ações em caso afirmativo.
The Simplest Case: The 2
2
• • • “-” and “+” denote the low and high levels of a factor, respectively Low and high are arbitrary terms Geometrically, the four runs form the corners of a square Factors can be quantitative or qualitative, although their treatment in the final model will be different
Chemical Process Example
A =
reactant concentration,
B
= catalyst amount,
y
= recovery
Analysis Procedure for a Factorial Design
• • Estimate factor
effects Formulate
model – With replication, use full model – With an unreplicated design, use normal probability plots • • Statistical
testing
•
Refine
the model • Analyze
residuals Interpret
results (ANOVA) (graphical)
Estimation of Factor Effects
A
y A
ab
a y A
2
n b
(1) 2
n B
1 2
n
[
ab y B
y B
a
(1) (1)] 2
n
2
n AB
1 2
n ab
[ (1) 1 2
n
2
n
[
ab
(1) 2
n
(1)] ] See textbook, pg. 209-210 For
manual
calculations The effect estimates are:
A
= 8.33,
B
= -5.00,
AB
= 1.67
Practical interpretation?
Design-Expert analysis
Estimation of Factor Effects Form Tentative Model
Model Model Model Model Error Error Term Intercept A Effect SumSqr 8.33333 208.333
% Contribution B AB -5 Lack Of Fit 0 75 1.66667 8.33333 2.57998
0 23.2198
P Error 31.3333 9.70072
64.4995
Lenth's ME Lenth's SME 6.15809
7.95671
Obs.: As somas de quadrados aqui são bastante simples.
Statistical Testing - ANOVA
The F test for the “model” source is testing the significance of the overall model; that is, is either A, B, or AB or some combination of these effects important?
Residuals and Diagnostic Checking
The 2
3
Factorial Design
Effects in The 2
3
Factorial Design
A B C
y y A
B y C
y B
y A
y C
etc, etc, ...
Analysis done via computer
An Example of a 2
3
Factorial Design
A
= gap,
B
= Flow,
C
= Power,
y
= Etch Rate
Table of – and + Signs for the 2 3 Factorial Design (pg. 218)
Properties of the Table
• Except for column
I
, every column has an equal number of + and – signs • The sum of the product of signs in any two columns is zero • Multiplying any column by
I
leaves that column unchanged (identity element) • The product of any two columns yields a column in the table:
AB
2
AB C
AC
•
Orthogonal
design • Orthogonality is an important property shared by all factorial designs
Estimation of Factor Effects
ANOVA Summary – Full Model
Model Interpretation
Cube plots
are often useful visual displays of experimental results
Cube Plot of Ranges
What do the large ranges when gap and power are at the high level tell you?