Set 6
Set 6
Big Picture
• CRD: 1 treatment factor with p levels + no blocking factors
• RBD: 1 treatment factor with p levels + 1 blocking factor, b blocks of size p each
• BIBD: 1 treatment factor with p levels + 1 blocking factor, b blocks of size k (< p) each
• LSD: 1 treatment factor with p levels + 2 blocking factors, p blocks of size p each
• GLSD: 1 treatment factor with p levels + 3 blocking factors, p blocks of size p each
Big Picture
• CRD: 1 treatment factor with p levels + no blocking factors
• RBD: 1 treatment factor with p levels + 1 blocking factor, b blocks of size p each
• BIBD: 1 treatment factor with p levels + 1 blocking factor, b blocks of size k (< p) each
• LSD: 1 treatment factor with p levels + 2 blocking factors, p blocks of size p each
• GLSD: 1 treatment factor with p levels + 3 blocking factors, p blocks of size p each
Factorial design
Simplest 2 – treatment factor design
Factorial design
Simplest 2 – treatment factor design
Factorial design
Simplest 2 – treatment factor design
Example
2! Factorial design
Example:
2! Factorial design
Example:
Main effect of A:
40 + 52 20 + 30
𝑦'%& − 𝑦'%' = − = 21
2 2
2! Factorial design
Example:
Main effect of A:
40 + 52 20 + 30
𝑦'%& − 𝑦'%' = − = 21
2 2
Main effect of B:
30 + 52 20 + 40
𝑦'(& − 𝑦'(' = − = 11
2 2
2! Factorial design
Example:
Main effect of A:
40 + 52 20 + 30
𝑦'%& − 𝑦'%' = − = 21
2 2
Main effect of B:
30 + 52 20 + 40
𝑦'(& − 𝑦'(' = − = 11
2 2
2! Factorial design
Example:
Main effect of A: 21
Main effect of B: 11
Interaction effect of AB: 1
Main effect of A:
Main effect of B:
Interaction effect of AB:
2! Factorial design
Example:
Main effect of A: 21
Main effect of B: 11
Interaction effect of AB: 1
Main effect of A:
Main effect of B:
Interaction effect of AB:
Regression model
2! Factorial design
Example:
Assuming treatment factors are quantitative, define
𝑥" = −1 𝐴 𝑙𝑜𝑤 ; 𝑥" = +1 (𝐴 ℎ𝑖𝑔ℎ).
𝑥# = −1 𝐵 𝑙𝑜𝑤 ; 𝑥# = +1 (𝐵 ℎ𝑖𝑔ℎ).
2! Factorial design
Example:
Assuming treatment factors are quantitative, define
𝑥" = −1 𝐴 𝑙𝑜𝑤 ; 𝑥" = +1 (𝐴 ℎ𝑖𝑔ℎ).
𝑥# = −1 𝐵 𝑙𝑜𝑤 ; 𝑥# = +1 (𝐵 ℎ𝑖𝑔ℎ).
Example
3! Factorial design
Example: an engineer is designing a battery that will be subjected to extreme
variations in temperature and has three choices for the plate material.
Client’s question: Is there a choice of material that would give uniformly long life
regardless of temperature?
3! Factorial design
ANOVA model : CRD with two treatment factors
3! Factorial design
Residual analysis
3! Factorial design
ANOVA model
3! Factorial design
Multiple comparison
One approach to this situation is to fix factor B at a specific level and apply Tukey’s
test to the means of factor A at that level.
Exercise
Exercise
Exercise
Exercise
#
More on 2 factorial designs
Setup
• No blocking factors – all experimental units are homogenous
• Errors are normally distributed
$ '("' (" $
Main effect of A: 𝑦!!" − 𝑦!!# = %!"# & 𝐴, 𝑇𝑟 = %&
− %&
$
Main effect of B: 𝑦!)" − 𝑦!)# = %!"# & ⟨𝐵, 𝑇𝑟⟩
$
Main effect of AB: 𝑦!!)" − 𝑦!!)# = ⟨𝐴𝐵, 𝑇𝑟⟩
%!"# &
Prof. Pritam Ranjan / OM&QT Area 30
Factorial Designs
$ '"'("'*"'(* $ "("*"(*
Main effect of A: 𝑦!!" − 𝑦!!# = %!"# & 𝐴, 𝑇𝑟 = +&
− +&
$
Main effect of B: 𝑦!)" − 𝑦!)# = %!"# & ⟨𝐵, 𝑇𝑟⟩
$
Main effect of AB: 𝑦!!)" − 𝑦!!)# = ⟨𝐴𝐵, 𝑇𝑟⟩
%!"# &
Prof. Pritam Ranjan / OM&QT Area 31
Factorial Designs
Example – 6.1
Example – 6.1
Example – 6.1
Ø lm_fit=lm(etch~A+C+A:C, data=dat)
Ø summary(lm_fit)
# Prediction
yhat = predict(lm_fit)
res = rstandard(lm_fit)
# Diagnostic plots
plot(etch,yhat)
plot(yhat,res); abline(h=0)
qqnorm(res);qqline(res)
Example – 6.1
# Diagnostic plots
plot(etch,yhat)
plot(yhat,res); abline(h=0)
qqnorm(res);qqline(res)
#
Unreplicated 2 factorial
designs
Example – 6.2
• A factorial experiment is carried out in a chemical plant to
study the factors thought to influence the filtration rate of a
chemical product. The factors are temperature (A), pressure
(B), concentration of formaldehyde (C), and stirring rate (D).
Example – 6.2
• Data analysis results
Example – 6.2
• How to test the significance of
factors and factor interactions?
Half-normal plot
• How to test the significance of factors and factor interactions?
1. Half normal plot - a plot of the absolute value of the effect estimates
against their cumulative normal probabilities
2. Length’s test
3. Permutation test
Ø library(DoE.base)
Ø halfnormal(lm_fit, ME.partial = TRUE)
Example-6.2 - diagnostics
Example-6.2
• Effect Hierarchy
• Lower order effects are more important than higher order effects
• Effect Heredity
• AB is included only if both A and B are already in the model (Strong)
• AB is included only if at least one of A and B are in the model (Weak)
Example-6.2
• Find a good final response surface model (R – illustration).
Example-6.2
? ?
• Look at both 𝑅<=> and 𝑅@ABCC
Example – 6.6
• An experiment is conducted to test new ideas to increase
direct mail sales by the credit card division of a financial
services company The factors tested in the experiment are
Practice Problems
Exercises
Exercises
Exercises
Exercises
• (6.22) The effect estimates from a 24 factorial design are as follows: ABCD =
−1.5138, ABC = −1.2661, ABD =−0.9852, ACD = −0.7566, BCD = −0.4842, CD =
−0.0795, BD = −0.0793, AD = 0.5988, BC = 0.9216, AC = 1.1616, AB = 1.3266, D =
4.6744, C = 5.1458, B = 8.2469, and A = 12.7151. Are you comfortable with the
conclusions that all main effects are active?
Exercises
Thank You
Questions ?
Reference – Chapters 5, 6
Next Topic
• Blocking and Confounding in Factorial Designs