1 STA 536 – Fractional Factorial Experiments at Two Levels Fractional Factorial Experiments at Two Levels Source : sections 5.1-5.5, part of 5.6 Effect.
Post on 29-Mar-2015
236 Views
Preview:
Transcript
CommandButton1CommandButton1CommandButton1CommandButton1 1STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Fractional Factorial Experiments at Two LevelsSource : sections 5.1-5.5, part of 5.6
Effect aliasing, resolution, minimum aberration criteria.
Analysis.
Techniques for resolving ambiguities in aliased effects.
Choice of designs, use of design tables.
Blocking in 2k−p designs.
CommandButton1CommandButton1CommandButton1CommandButton1 2STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Fractional Factorial Experiments at Two Levels
A 2k full factorial design requires 2k runs to be performed.
Full factorial designs are rarely used in practice for large k (say, k > 7).
A fractional factorial design is a subset or fraction of a full factorial design.
For economic reasons, fractional factorial designs are commonly used
CommandButton1CommandButton1CommandButton1CommandButton1 3STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.1 A Leaf Spring Experiment
An experiment to improve a heat treatment process on truck leaf springs.
The heat treatment which forms the camber (or curvature) in leaf springs consists of heating in a high temperature furnace, processing by a forming machine, and quenching in an oil bath.
Response: the height of an unloaded spring, known as free height.
Goal is to make the variation about the target 8.0 as small as possible.
A half fraction of a 25 full factorial design is used to study 5 factors. A 25−1 design
The quench oil temperature (Q) could not be controlled precisely.
CommandButton1CommandButton1CommandButton1CommandButton1 4STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Factors and Levels, Leaf Spring Experiment
FactorLevel
+ -B. high heat temperature (F) 1840 1880
C. heating time (seconds) 23 25
D. transfer time (seconds) 10 12
E. hold down time (seconds) 2 3
Q. quench oil temperature (F) 130-150 150-170
CommandButton1CommandButton1CommandButton1CommandButton1 5STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Leaf Spring Experiment: Design Matrix and Data
CommandButton1CommandButton1CommandButton1CommandButton1 6STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.2 Fractional Factorial Designs: Effect Aliasing and the Criteria of Resolution and Minimum Aberration
If a 25 design is used for the experiment, its 31 degrees of freedom would be allocated as follows:
Using effect hierarchy principle, one would argue that 4fi’s, 5fi’s and even 3fi’s are not likely to be important. There are 10+5+1 = 16 such effects, half of the total runs! Using a 25 design can be wasteful (unless 32 runs cost about the same as 16 runs.)
A half fraction of a 25 full factorial design has 16 runs and can estimate 15 factorial effects.
Ideally, estimate 5 ME’s and 10 2fi’s.
Why Using Fractional Factorial Designs?
CommandButton1CommandButton1CommandButton1CommandButton1 7STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Leaf Spring Experiment: 25-1
Use of a FF design instead of full factorial design is usually done for economic reasons. Since there is no free lunch , what price to pay?
What factorial effects can be estimated in the leaf spring experiment?
Notice that column E equals the product of the columns B, C and D.
E=BCD Consequence: the column for E is used for estimating the
main effect of E and also for the interaction effect between B, C and D, i.e.,
CommandButton1CommandButton1CommandButton1CommandButton1 8STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Leaf Spring Experiment: 25-1
The data from such a design is not capable of distinguishing the estimate of E from the estimate of BCD.
The main effect E is aliased with the B × C × D interaction.
Assuming BCD interaction is negligible, we can estimate E.
A 25−1 has 16 runs and can estimate 15 pairs of aliased factorial effects.
The replicates can be used for estimating error or dispersion effects.
CommandButton1CommandButton1CommandButton1CommandButton1 9STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Leaf Spring Experiment: 25-1
E = BCD or I = BCDE I = column of +’s is the identity element in the group of
multiplications. I = BCDE is the defining relation for the 25−1 design. It
implies all the 15 effect aliasing relations :B =CDE, C = BDE, D = BCE, E = BCD,
BC = DE, BD =CE, BE =CD,Q = BCDEQ, BQ =CDEQ, CQ = BDEQ, DQ = BCEQ,EQ = BCDQ, BCQ = DEQ, BDQ =CEQ, BEQ =CDQ.
CommandButton1CommandButton1CommandButton1CommandButton1 10STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
B =CDE, C = BDE, D = BCE, E = BCD,BC = DE, BD =CE, BE =CD,Q = BCDEQ, BQ =CDEQ, CQ = BDEQ, DQ = BCEQ,EQ = BCDQ, BCQ = DEQ, BDQ =CEQ, BEQ =CDQ.
Each of the four main effects, B, C, D and E, is estimable if the respective three-factor interaction alias is negligible, which is usually a reasonable assumption.
The B×C interaction is estimable only if prior knowledge would suggest that the D × E interaction is negligible.
The Q main effect is estimable under the weak assumption that the B × C × D × E × Q interaction is negligible.
The B ×Q interaction is estimable under the weak assumption that the C × D × E × Q interaction is negligible.
CommandButton1CommandButton1CommandButton1CommandButton1 11STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Definition: Clear Effects
A main effect or two-factor interaction is clear if none of its aliases are main effects or two-factor interactions
and strongly clear if none of its aliases are main effects, two-factor or three-factor interactions.
Therefore a clear effect is estimable under the assumption of negligible 3-factor and higher interactions and a strongly clear effect is estimable under the weaker assumption of negligible 4-factor and higher interactions.
For the leaf spring experiment, Clear effects: B, C, D, E, Q, B × Q, C × Q, D × Q and E ×
Q. Strongly clear effects: Q, B × Q, C × Q, D × Q and E × Q.
CommandButton1CommandButton1CommandButton1CommandButton1 12STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
An Alternate Design with defining relation: I = BCDEQ.
The 15 aliasing relations are:B=CDEQ, C=BDEQ, D=BCEQ, E=BCDQ,
Q=BCDEBC=DEQ, BD=CEQ, BE=CDQ, BQ=CDE,CD=BEQ, CE=BDQ, CQ=BDE, DE=BCQ,
DQ=BCE, EQ=BCD Clear effects: all 5 ME’s and 10 2fi’s. Strongly clear effects: all 5 ME’s.
CommandButton1CommandButton1CommandButton1CommandButton1 13STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Which of the two fractional plans is better?I=BCDE vs I=BCDEQ
The 2nd plan is generally preferred. (Why? Because each of the 15 degrees of freedom is either clear or strongly clear)
The 1st plan is preferred if there is special interest on factor Q and 2fi’s involving Q.
The reason for adopting the 1st plan will be given in Chapter 11 in the context of robust parameter design.
CommandButton1CommandButton1CommandButton1CommandButton1 14STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
A 2k−p fractional factorial design
has k factors, each at 2 levels, consisting of 2k−p runs.
is a 2−p-th fraction of the 2k full factorial designis determined by p defining wordsLetters are the names of the factors: 1, 2, . . ., k
or A,B, . . ..The number of letters in a word is its
wordlength.The group formed by the p defining words is
called the defining contrast subgroup.The group consists of 2p-1 words plus the
identity element I.
CommandButton1CommandButton1CommandButton1CommandButton1 15STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Example: A 26−2 design d
defining relations: 5 = 12 and 6 = 134.
The two independent defining words: I = 125 and I = 1346.
Then I = 125×1346 = 23456.
defining contrast subgroup{I,125,1346,23456} or
I = 125 = 1346 = 23456.
CommandButton1CommandButton1CommandButton1CommandButton1 16STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Example: A 26−2 design d5 = 12 and 6 = 134
I = 125 = 1346 = 23456,1 = 25 = 346 = 123456,2 = 15 = 12346 = 3456,3 = 1235 = 146 = 2456,4 = 1245 = 136 = 2356,5 = 12 = 13456 = 2346,6 = 1256 = 134 = 2345,13 = 235 = 46 = 12456,14 = 245 = 36 = 12356,16 = 256 = 34 = 12345,23 = 135 = 1246 = 456,24 = 145 = 1236 = 356,26 = 156 = 1234 = 345,35 = 123 = 1456 = 246,45 = 124 = 1356 = 236,56 = 126 = 1345 = 234.
clear effects: 3, 4, 6, and 23, 24, 26, 35, 45, 56.
wordlength pattern: W = (1, 1, 1, 0).
resolution: III
CommandButton1CommandButton1CommandButton1CommandButton1 17STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
WordLength Pattern and Resolution
Ai = number of words of length i in its defining contrast subgroup.
wordlength pattern W = (A3, . . . ,Ak)resolution=the length of the shortest word in
the defining contrast subgroup, i.e., smallest r such that Ar≥1.
A design of resolution I or II is useless. (Why? No main effect is allowed to be aliased with another main effect.)
CommandButton1CommandButton1CommandButton1CommandButton1 18STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Resolution criterion
Maximum Resolution criterion, proposed by Box and Hunter (1961),
chooses the 2k−p design with maximum resolution. justified by the hierarchical ordering principle (Section 3.5).Rules for Resolution: Resolution R implies that no effect involving i factors is
aliased with effects involving less than R − i factors. In a resolution III design, some ME’s are aliased with 2fi’s, but
not with other ME’s. In a resolution IV design, some ME’s are aliased with 3fi’s and
some 2fi’s are aliased with other 2fi’s. (all ME’s are clear.) In a resolution V design, some ME’s are aliased with 4fi’s and
some 2fi’s are aliased with 3fi’s. (all ME’s are strongly clear and all 2fi’s are clear.)
CommandButton1CommandButton1CommandButton1CommandButton1 19STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Rules for Resolution IV and V Designs:
In any resolution IV design, the main effects are clear.
In any resolution V design, the main effects are strongly clear and the two-factor interactions are clear.
Among the resolution IV designs with given k and p, those with the largest number of clear two-factor interactions are the best.
CommandButton1CommandButton1CommandButton1CommandButton1 20STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
A Projective Rationale for Resolution
A projection on to any n ( ≤ R − 1) factors is a full 2n factorial design (replicated if n ≤ k − p).
Fig. 4.1: projections of 23−1 designs.
In analysis, if there are at most R−1 important factors, then all factorial effects among the important factors can be estimated.
Projective property of a design
CommandButton1CommandButton1CommandButton1CommandButton1 21STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Minimum Aberration Criterion
Example: two 27−2 designs:
d1 : I = 4567 = 12346 = 12357,
d2 : I = 1236 = 1457 = 234567.
• Both designs have resolution IV.
• wordlength patterns:
W(d1) = (0, 1, 2, 0, 0) and W(d2) = (0, 2, 0, 1, 0).
CommandButton1CommandButton1CommandButton1CommandButton1 22STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
d1 : I = 4567 = 12346 = 12357,d2 : I = 1236 = 1457 = 234567.
d1 has less aberration than d2.d1 has minimum aberration (by complete
search.)aliased 2fi’s:
d1 : 45 = 67, 46 = 57, 47 = 56, d2 : 12 = 36, 13 = 26, 16 = 23, 14 = 57, 15
= 47, 17 = 45d1 is preferred to d2.
CommandButton1CommandButton1CommandButton1CommandButton1 23STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Create FFD in SAS/*FFD 2^(8-3) design using 6=135 7=235 8=1234*/
proc factex;
factors x1-x5;
output out = SavedDesign;
data SavedDesign1; set SavedDesign;
x6=x1*x3*x5; x7=x2*x3*x5; x8=x1*x2*x3*x4;
run;
proc print; run;
/*FFD 2^(8-3) design using 6=12 7=13 8=34*/
data SavedDesign2; set SavedDesign;
x6=x1*x2; x7=x1*x3; x8=x3*x4;
run;
proc print; run;
CommandButton1CommandButton1CommandButton1CommandButton1 24STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
FFD by SAS Macro%macro createFFD;
proc factex;
factors x1-x&k /NLEV=&s;
output out = SavedDesign;
proc IML;
use SavedDesign;
read all into D;
def=&def;
if &s=2 then D=(D+1)/2;
else if &s=3 then D=D+1;
DD=D || mod(D*def`, &s);
create ffddsn from DD;
append from DD;
quit;
%mend;
/*FFD 2^(8-3) with 6=135, 7=235, 8=1234*/
%let k=5; %let q=3; %let s=2;
%let def={1 0 1 0 1, 0 1 1 0 1, 1 1 1 1 0};
%createFFD;
proc print; run;/*FFD 2^(8-3) with 6=12, 7=13, 8=45*/
%let k=5; %let q=3; %let s=2;
%let def={1 1 0 0 0, 1 0 1 0 0, 0 0 0 1 1};
%createFFD;
proc print; run;
CommandButton1CommandButton1CommandButton1CommandButton1 25STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
d1 : I = 4567 = 12346 = 12357,d2 : I = 1236 = 1457 = 234567.
/*d1 : I = 4567 = 12346 = 12357*/
%let k=5; %let s=2;
%let def={1 1 1 1 0, 1 1 1 0 1};
%createFFD;
/*d2 : I = 1236 = 1457 = 234567*/
%let k=5; %let s=2;
%let def={1 1 1 0 0, 1 0 0 1 1};
%createFFD;
CommandButton1CommandButton1CommandButton1CommandButton1 26STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.3 Analysis - the leaf spring experiment
defining relation: I = BCDE aliasing relations:
B = CDE, C = BDE, D = BCE, E = BCD,BC = DE, BD = CE, BE = CD,
Q = BCDEQ, BQ = CDEQ, CQ = BDEQ, DQ = BCEQ,EQ = BCDQ, BCQ = DEQ, BDQ = CEQ, BEQ = CDQ.
9 clear effects: B, C, D, E, Q, BQ, CQ, DQ and EQ 6 pairs of aliased effects: BC = DE, BD = CE, BE = CD,
BCQ =DEQ, BDQ = CEQ, BEQ = CDQ an effect is estimable if its alias is negligible. Same strategy as in full factorial experiments except for
the interpretation and handling of aliased effects.
CommandButton1CommandButton1CommandButton1CommandButton1 27STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
data leaf; input B C D E Q y1-y3 m s2 lns; BQ=B*Q; CQ=C*Q; DQ=D*Q; EQ=E*Q; BC=B*C; BD=B*D; CD=C*D; BCQ=B*C*Q; BDQ=B*D*Q; BEQ=B*E*Q; cards;-1 1 1 -1 -1 7.78 7.78 7.81 7.7900 0.0003 -8.1117 1 1 1 1 -1 8.15 8.18 7.88 8.0700 0.0273 -3.6009-1 -1 1 1 -1 7.50 7.56 7.50 7.5200 0.0012 -6.7254 1 -1 1 -1 -1 7.59 7.56 7.75 7.6333 0.0104 -4.5627-1 1 -1 1 -1 7.94 8.00 7.88 7.9400 0.0036 -5.6268 1 1 -1 -1 -1 7.69 8.09 8.06 7.9467 0.0496 -3.0031-1 -1 -1 -1 -1 7.56 7.62 7.44 7.5400 0.0084 -4.7795 1 -1 -1 1 -1 7.56 7.81 7.69 7.6867 0.0156 -4.1583-1 1 1 -1 1 7.50 7.25 7.12 7.2900 0.0373 -3.2888 1 1 1 1 1 7.88 7.88 7.44 7.7333 0.0645 -2.7406-1 -1 1 1 1 7.50 7.56 7.50 7.5200 0.0012 -6.7254 1 -1 1 -1 1 7.63 7.75 7.56 7.6467 0.0092 -4.6849-1 1 -1 1 1 7.32 7.44 7.44 7.4000 0.0048 -5.3391 1 1 -1 -1 1 7.56 7.69 7.62 7.6233 0.0042 -5.4648-1 -1 -1 -1 1 7.18 7.18 7.25 7.2033 0.0016 -6.4171 1 -1 -1 1 1 7.81 7.50 7.59 7.6333 0.0254 -3.6717 ;
CommandButton1CommandButton1CommandButton1CommandButton1 28STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
SAS code
PROC REG DATA=LEAF outest=effects;
MODEL M=B C D E Q BQ CQ DQ EQ BC BD CD BCQ BDQ BEQ;
RUN;
%halfnormalplot(effects, m);
PROC REG DATA=LEAF outest=effects;
MODEL LNS=B C D E Q BQ CQ DQ EQ BC BD CD BCQ BDQ BEQ;
RUN;
%halfnormalplot(effects, LNS);
CommandButton1CommandButton1CommandButton1CommandButton1 29STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
SAS Macro halfnormalplot%macro halfnormalplot(effectdata, response);/*effectdata is from the output of regression*/data effect2; set &effectdata; drop &response intercept _RMSE_;proc transpose data=effect2 out=effect3;data effect4; set effect3; effect=ABS(col1*2); proc sort data=effect4; by effect;proc print data=effect4;proc rank data=effect4 out=effect5; var effect; ranks nefff; run;proc sql; select * from effect5; quit;data effect5; set effect5;
neff=quantile('NORMAL', 0.5+0.5*(nefff-0.5)/&sqlobs);proc gplot data=effect5; plot effect*neff=_NAME_;run; quit;%mend;
CommandButton1CommandButton1CommandButton1CommandButton1 30STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Factorial Effects, Leaf Spring Experiment
CommandButton1CommandButton1CommandButton1CommandButton1 31STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
For the location effects, visually one may judge that Q,B,C,CQ and possibly E,BQ are significant.
One can apply the studentized maximum modulus test (see chapter 4, not covered in class) to confirm that Q,B,C are significant at 0.05 level.
The B×Q and C×Q plots show that they are synergystic.
CommandButton1CommandButton1CommandButton1CommandButton1 32STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
For illustration, we use the model
For the dispersion effects (based on zi = lns2i
values), the half-normal plot is given in next slide. Visually only effect B stands out. This is confirmed by applying the studentized maximum modulus test. For illustration, we will include B,DQ,BCQ in the following model,
(1)
CommandButton1CommandButton1CommandButton1CommandButton1 33STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
For the dispersion effects
CommandButton1CommandButton1CommandButton1CommandButton1 34STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
CommandButton1CommandButton1CommandButton1CommandButton1 35STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Two-Step Procedure for Optimization
Minimize dispersion. Choose xB = −1, xDxQ = −1 and xBxCxQ = +1. Use interaction plots to determine level settings of C,
D, Q. Should choose (B,C,D,Q) = (−,+,+,−).
Bring location to target (t = 8.0). E is an adjustment factor. After chosen (B,C,D,Q) = (−,+,+,−).
ˆy = 7.6360 + 0.1106(−1) + 0.0519xE + 0.0881(+1) − 0.1298(−1)+0.0423(−1)(−1) − 0.0827(+1)(−1)= 7.8683 + 0.0519xE.
By solving ˆy = 8.0, xE should be set at (8.0 − 7.8683)/0.0519 = 2.54.
Note:The uncoded value for factor E is 2.5 + 2.54(0.5) = 3.77 seconds.Caution: there are no experimental data around 3.77 seconds.
CommandButton1CommandButton1CommandButton1CommandButton1 36STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.4 Techniques for Resolving the Ambiguities in Aliased Effects
Among the three factorial effects that feature in model (1), B is clear and DQ is strongly clear.
However, the term xBxCxQ is aliased with xDxExQ. The following three techniques can be used to resolve the
ambiguities. Subject matter knowledge may suggest some effects in the
alias set are not likely to be significant (or does not have a good physical interpretation).
Or use effect hierarchy principle to assume away some higher order effects.
Or use a follow-up experiment to de-alias these effects. Three methods are given in section 5.4 of WH. Two are considered here.
CommandButton1CommandButton1CommandButton1CommandButton1 37STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Fold-over Technique
Suppose the original experiment is based on a design with generators
d1 : 4 = 12, 5 = 13, 6 = 23, 7 = 123.None of its main effects are clear.
To de-alias them, we can choose another 8 runs (no. 9-16 in Table 4.7) with reversed signs for each of the 7 factors. This follow-up design d2 has the generators
d2 : 4 = −12,5 = −13,6 = −23,7 = 123 With the extra degrees of freedom, we can introduce a
new factor 8 for run number 1-8, and -8 for run number 9-16. See Table 4.7.
The combined design d1+d2 is a design and thus all main effects are clear. (Its defining contrast subgroup is on p.227 of WH).
CommandButton1CommandButton1CommandButton1CommandButton1 38STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Table 5.6: Augmented Design Matrix Using Fold-Over Technique
7=1238=1248=1358=236=>5=138=13*124=2346=238=23*124=1347=1238=124Without 8,
I=2345=1346=1237 =1256=1457=3567 =3567
CommandButton1CommandButton1CommandButton1CommandButton1 39STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Suppose one factor, say 5, is very important. We want to de-alias 5 and all 2fi’s involving 5.
Choose, instead, the following design
Then the combined design d1+d3 is a design with the generators
Since 5 does not appear in (5), 5 is strongly clear and all 2fi’s involving 5 are clear. However, other main effects are not clear (see Table 5.7 of WH for d1+d3).
Choice between d2 and d3 depends on the priority given to the effects
{I=124=236=1237=1346=347=167=2467} Block Factor 8=135. Thus {1358=I=124=236=1237=
1346= 347=167=2467 =23458=12568=2578=4568= 14578=35678=12345678}
Fold-over Technique: Version Two
CommandButton1CommandButton1CommandButton1CommandButton1 40STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Critique of Fold-over Technique
Fold-over technique is not an efficient technique. It requires doubling of the run size and can only de-alias a specific set of effects.
In practice, after analyzing the first experiment, a set of effects will emerge and need to be de-aliased. It will usually require much fewer runs to de-alias a few effects.
A more efficient technique that does not have these deficiencies is the optimum design approach given in Section 5.4.2.
CommandButton1CommandButton1CommandButton1CommandButton1 41STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Optimal Design Approach for Follow-Up Experiments This approach add runs according to a particular optimal
design criterion. The D and Ds criteria shall be discussed. Optimal design criteria depend on the assumed model.
In general, the model should contain: All effects and their aliases (except those judged
unimportant a priori or by the effect hierarchy principle) identified as significant in the initial experiment.
A block variable that accounts for differences in the average value of the response over different time periods.
An intercept. In the leaf spring experiment, we specify the model:
CommandButton1CommandButton1CommandButton1CommandButton1 42STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
D-Criterion In Table 5.8, the columns B, C, D, E, and Q comprise the
design matrix while the columns B, block, BCQ, DEQ, DQ comprise the model matrix. Two runs are to be added to the original 16-run experiment. There are 45 = 1024 possible choices of factor settings for the follow up runs (runs 17 and 18) since each factor can take on either the + or - level in each run.
For each of the 1024 choices of settings for B, C, D, E, Q for runs 17 and 18, denote the corresponding model matrix by Xd, d = 1, ...,1024. We may choose the factor settings d∗ that maximizes the D-criterion, i.e.
Maximizing the D criterion minimizes the volume of the confidence ellipsoid for all model parameters b.
64 choices of d attain the maximum D value of 4,194,304. Two are: d1 : (B,C,D,E,Q) = (+++−+) and (++−+−) d2 : (B,C,D,E,Q) = (+++−+) and (−++++)
CommandButton1CommandButton1CommandButton1CommandButton1 43STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Table 5.8: Augmented Design Matrix and Model Matrix, Leaf Spring Experiment
CommandButton1CommandButton1CommandButton1CommandButton1 44STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
/*Optimal Design Approach for Follow-Up Experiments*/
proc IML;
use leaf;
read all;
x={1 1 -1 -1, 1 -1 1 -1};
res={0 0 0 0 0 0 0}; kk=0;
do i1=1 to 4;
B1=B // x[,i1];
do i2=1 to 4;
C1=C // x[,i2];
do i3=1 to 4;
D1=D // x[,i3];
do i4=1 to 4;
E1=E // x[,i4];
do i5=1 to 4;
Q1=Q // x[,i5];
kk=kk+1;
DD=j(18,1,1)|| (j(16,1,-1)//{1,1}) || B1 || D1#Q1 || B1#C1#Q1 || D1#E1#Q1;
tmp=det(DD`*DD);
res=res // (kk || i1 || i2 || i3 || i4 || i5 || tmp);
end;
end;
end;
end;
end;
call sort(res, {7 2 3 4 5 6}, {7});
print x, res;
quit;
CommandButton1CommandButton1CommandButton1CommandButton1 45STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
d1 : (B,C,D,E,Q) = (+++−+) and (++−+−)
d2 : (B,C,D,E,Q) = (+++−+) and (−++++)
CommandButton1CommandButton1CommandButton1CommandButton1 46STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
CommandButton1CommandButton1CommandButton1CommandButton1 47STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.5 Selection of 2k−p Designs Using Minimum Aberration and Related Criteria Questions:
How to choose a good 2k−p design? What criteria should be used to judge “goodness”?
Example: Two 29−4IV designs from Table 5A.3 (P254)
Design d1 with generators 6=123, 7=124, 8=125, 9=1345
Defining contrast subgroup:I = 1236 = 1247 = 1258 = 13459= 3467 = 3568 = 24569 = 4578 = 23579 = 23489= 12345678 = 15679 = 14689 = 13789 = 26789
Resolution: IV Wordlength pattern: W = (0, 6, 8, 0, 0, 1, 0) Clear effects: all 9 MEs, any 2fis involving 9 Strongly clear effects: 9
CommandButton1CommandButton1CommandButton1CommandButton1 48STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Comparing two 29-4 design
Design d2 with generators 6=123, 7=124, 8=134, 9=2345 Defining contrast subgroup:
I = 1236 = 1247 = 1348 = 23459= 3467 = 2468 = 14569 = 2378 = 13579 = 12589= 1678 = 25679 = 35689 = 45789 = 123456789
Resolution: IV Wordlength pattern: W=(0,7,7,0,0,0,1) Clear effects: all 9 MEs, any 2fis involving 5 or 9 Strongly clear effects: 5 and 9
Which design is better? Answer depends.
CommandButton1CommandButton1CommandButton1CommandButton1 49STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Comparing two 29-4 design
d1 has less aberration than d2 (A4(d1)=6<A4(d2)=7). d1 is preferred to d2 in general. d2 has 15 clear 2fi’s, but d1 has 8 clear 2fi’s. d2 may be superior to d1 if prior information is available
on the importance of 2fi’s.Guideline: Minimum aberration designs are preferred in general (by
hierarchical ordering principal). Among resolution IV designs, those with the largest
number of clear 2fi’s are the best.
CommandButton1CommandButton1CommandButton1CommandButton1 50STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Tables 5A.1-5A.7 2k−p fractional factorial designs.
CommandButton1CommandButton1CommandButton1CommandButton1 51STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Table 5A.2: 16-Run 2k−p FFD (k− p = 4)
CommandButton1CommandButton1CommandButton1CommandButton1 52STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Table 5A.3: 32 Run 2k−p FFD (k−p=5, 6≤k≤11)
CommandButton1CommandButton1CommandButton1CommandButton1 53STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Research problems:
How to construct minimum aberration designs with large run sizes (>128)?
What is the relationship between minimum aberration and maximum #of clear effects?
CommandButton1CommandButton1CommandButton1CommandButton1 54STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Example: Two 26−2 designs in Table 5A.2
design d1 with generators 5 = 123, 6 = 124I = 1235 = 1246 = 3456
resolution IV 6 clear effects: 6 ME’s.
design d2 with generators 5 = 12, 6 = 134I = 125 = 1346 = 23456
resolution III 9 clear effects: 3, 4, 6, 23, 24, 26, 35, 45, 56.
d1 is generally preferred. d2 is useful when some 2fi’s are important.
CommandButton1CommandButton1CommandButton1CommandButton1 55STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Choice of Fractions and Avoidance of Specific Combinations
A 27−4III design with 4=12, 5=13, 6=23 and 7=123.
There are 16 fractions (27−4) by using + or − signs in the 4 generators.
4 = ±12, 5 = ±13, 6 = ±23, 7 = ±123. Most tables present the standard fraction with all +
signs. All 16 fractions are equivalent by changing signs of some
columns. Other fractions are useful to avoid some treatment
combinations (If specific combinations (e.g., (+++) for high pressure, high temperature, high concentration) are deemed undesirable or even disastrous, they can be avoided by choosing a fraction that does not contain them. Example on p.237 of WH)
CommandButton1CommandButton1CommandButton1CommandButton1 56STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Factor assignment:
24−1: 4=12325−2: 4=12, 5=1326−3: 4=12, 5=13,
6=2327−4: 4=12, 5=13,
6=23, 7=123
Table 5A.1 8-run designs
CommandButton1CommandButton1CommandButton1CommandButton1 57STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
5.6 Blocking in Fractional Factorial Designs
Example: To arrange a 26−2 design in 4 blocks in Table 4B.2.
design generators 5 = 123, 6 = 124 treatment defining contrast subgroup:
I = 1235 = 1246 = 3456 block generators B1 = 134,B2 = 234 block defining contrast subgroup:
B1 = 134,B2 = 234,B1B2 = 12
CommandButton1CommandButton1CommandButton1CommandButton1 58STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Aliasing relations:I = 1235 = 1246 = 3456B1 = 134 = 245 = 236 = 156B2 = 234 = 145 = 136 = 256B1B2 = 12 = 35 = 46 = 1234561 = 235 = 246 = 134562 = 135 = 146 = 234563 = 125 = 12346 = 4564 = 12345 = 126 = 3565 = 123 = 12456 = 3466 = 12356 = 124 = 34513 = 25 = 2346 = 145614 = 26 = 2345 = 135615 = 23 = 2456 = 134616 = 24 = 2356 = 134534 = 56 = 1245 = 123636 = 45 = 1256 = 1234
Clear effects: all 6 MEs
Can estimate 3 block effects, 6 clear ME’s and 6 pairs of aliased 2fi’s.
CommandButton1CommandButton1CommandButton1CommandButton1 59STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Definition: A main effect or two-factor interaction is said to be clear
if it is not aliased with any other main effects or two-factor interactions as well as not confounded with any block effects.
The definition of strongly clear effect has the additional requirement that it is not aliased with any three-factor interactions.
Minimum aberration Several different definitions by combining treatment and
block defining contrast subgroups.
CommandButton1CommandButton1CommandButton1CommandButton1 60STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
Table 5B.1-5B.5 (p. 260-264): 2k−p designs in 2q blocks.
Run sizes 8, 16, 32, 64, 128 with k − p = 3, 4, 5, 6, 7
k≤9 and 1≤q≤k−p−1.designs are chosen by # of clear effects.Two 26−2 designs with k = 6, p = 2, q = 1, 2.
Research problem:How to construct optimal blocked 2k−p designs
with large run sizes (≥128)?
CommandButton1CommandButton1CommandButton1CommandButton1 61STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
CommandButton1CommandButton1CommandButton1CommandButton1 62STA 536 – Fractional Factorial Experiments at Two Levels
CommandButton1CommandButton1CommandButton1CommandButton1
top related