Two-Level Fractional Factorial Designs

2y ago
30 Views
2 Downloads
211.45 KB
14 Pages
Last View : 1m ago
Last Download : 10m ago
Upload by : Kelvin Chao
Transcription

5Two-Level Fractional Factorial Designs Because the number of runs in a 2k factorial design increases rapidly as the number of factorsincreases, it is often impossible to run the full factorial design given available resources. If the experimenter can reasonably assume that certain high-order interactions (often 3-wayand higher order) are negligible, then it is often possible to estimate and test for the significanceof main effects and low-order interactions from only a fraction of the full factorial design. A design which contains a subset of factor level combinations from a full factorial design iscalled a fractional factorial design. A fractional factorial design is often used as a screening experiment involving many factorswith the goal of identifying only those factors having large effects. Once specific factors areidentified as important, they are investigated in greater detail in subsequent experiments. The successful use of two-level fractional factorial designs is based on three ideas:1. The sparsity of effects principle: When there are many variables under consideration,it is typical for the system or process to be dominated by main effects and low-orderinteractions.2. The projective property: A fractional factorial design can be projected into strongerdesigns in a subset of the significant factors.3. Sequential experimentation: If is possible to combine the runs from two or morefractional factorial designs to sequentially form a larger design that allows estimation ofall interaction effects of interest.5.12k 1 Fractional Factorial Designs Situation: There are k factors of interest each having 2 levels, but there are only enoughresources to run 1/2 of the full factorial 2k design. Thus, we say we want to run a 1/2fraction of a 2k design. This design is called a 2k 1 fractional factorial design.– Suppose there are 5 factors of interest (A, B, C, D, and E), and there are only enoughresources for 16 experimental runs. Thus, we want to run a 1/2 fraction of a 25 design.This design is called a 25 1 fractional factorial design. A 1/2 fraction can be generated from any interaction, but using the highest-order interactionis the standard. The interaction used to generate the 1/2 fraction is called the generator ofthe fractional factorial design.– When there are 3 factors, use ABC as the generator of the 23 1 design, when there are4 factors, use ABCD as the generator of the 24 1 design, . . ., when there are k factors,use ABC . . . K as the generator of the 2k 1 design. To generate a 2k 1 fractional factorial design from the highest-order interaction:(i) Select only those treatment combinations that have a plus ( ) sign in the ABC · · · Kcolumn of the 2k design, OR,(ii) Select only those treatment combinations that have a minus ( ) sign in the ABC · · · Kcolumn of the 2k design.Because the I column also contains all plus signs, we refer to I ABC · · · K orI ABC · · · Kas a defining relation for the 2k 1 design.71

Example When there are 3 factors, we have the following table of pluses and minuses.Treatmenttotalabcabcabacbc(1)I A Factorial EffectC AB AC B BC ABC AliasedeffectsI A B C I A B C Thus, for defining relation I ABC , the treatment combinations with ABC 1 formthe 23 1 design, or, the 1/2 fraction of the 23 design. The remaining 4 treatment combinationswith ABC 1 are discarded (i.e, discard the ab, ac, bc, and (1) rows). Because only 1/2 of the full factorial design is run, each of the 2k effects (including theintercept) is aliased with one other effect. That is, estimation of aliased effects are calculatedidentically and, therefore, cannot be separated from each other. To find the estimate of any model effect, use the difference in means between the and rows as we did with the 2k design. Example In the 23 1 design with defining relation I ABC, let lA , lB , and lC representestimates of the main effects. ThenlA 1[(a abc) (b c)]21[(b abc) (a c)]21 [(c abc) (a b)]2lB lCNow, let lAB , lAC , and lBC represent estimates of the two-factor interaction effects. ThenlBC 1[2]1[21 [2lAC ]lAB] Note that lA lBC , lB lAC , and lC lAB . Hence, it is impossible to separate the effectof A from BC, the effect of B from AC, and the effect of C from AB. Thus, the A and BCeffects are aliased, the B and AC effects are aliased, and the C and AB effects are aliased. In addition, note that the intercept I is aliased with ABC. If two effects E1 and E2 are aliased, then we cannot estimate each separately. That is, eachof the two columns produces an estimate of the sum of the two effects. For effect i, this isdenoted li E1 E2 for i 1, 2.72

Example: In the 23 1 design, we havelA A BClBC A BClB B ACandlAC B AClC C ABlAB C AB The alias structure for any 2k 1 design can be determined by taking the defining relationI ABC · · · K and multiplying it by any effect. The resulting effect is the aliased effect.– In the 23 1 design, to find the aliases of the main effects A, B, and C, we multiplyI ABC by A, B, and C, and reduce:A A · I A · ABC A2 BC BC A BCB B·I B ACC C ·I C AB Now consider the other 1/2 fraction having defining relation I ABC. This design containsthe treatment combinations that have minus signs for ABC · · · K. This is known as thealternate or complementary 1/2 fraction.– In the 23 1 design, the alternate 1/2 fraction are the 4 treatment combinations corresponding to rows (1), ab, ac, and bc.Treatmenttotalabacbc(1)I A B Factorial EffectC AB AC BC ABC AliasedeffectsI ABCA BCB ACC AB– The defining relation is I ABC and the alias structure using the multiplication ruleisA A · I A · ( ABC) A2 BC BC A BCB B·I B ACC C ·I C BCThus, the A and BC effects are aliased, the B and AC effects are aliased, and the C00and AB effects are aliased. The alternate fraction yields estimates lA, lBand lC0 where0lA 0lB lC0 Creating a 2k 1 fractional-factorial design is equivalent to forming 2 blocks from a 2k designand then selecting one block to run. In practice, it does not matter which fraction we use(although it is common to use the ( , ) case to avoid using the negative of an estimate. Because of aliased pairs of effects, only one effect from each pair can be included in the model.73

Example: A 24 1 design (a 1/2 fraction of a 24 design) Using generator I ABCD, the following table contains the 8 combinations ABCD 1.This is one possible 24 1 design.I A B C D ABCD Notation Using generator I ABCD, the following table contains the 8 combinations ABCD 1.This is a second possible 24 1 design.I A B C D ABCD Notation In the 24 1 design, to find the aliases of the A, B, C, D, AB, AC, . . . we multiply I ABCDby A, B, C, D, AB, AC, . . . and then reduce:A A · I A · ABCD A2 BCD BCDlA lBCD A BCDB B·I lB lACD B ACDC C ·I lC lABD C ABDD D·I lD lABC D ABCAB AB · I AB · ABCD A2 B 2 CD CDlAB lCD AB CDAC AC · I AC · ABCD A2 BC 2 D BDlAC lBD AC BDAD AD · I lAD lBC 74

5.2Design ResolutionA design is resolution R if no p-factor effect is aliased with another effect containing less thanR p factors. Resolution III Designs: No main effect is aliased with any other main effect, but at least onemain effect is aliased with a two-factor interaction. Resolution IV Designs: No main effect is aliased with any other main effect or two-factor interaction, but at least one two-factor interaction is aliased with another two-factor interaction. Resolution V Designs: No main effect or two-factor interaction is aliased with any other maineffect or two-factor interaction, but at least one two-factor interaction is aliased with a threefactor interaction. For example:– The 23 1 design with defining relation I ABC is resolution III and is denoted 23 1III .– The 24 1 design with defining relation I ABCD is resolution IV and is denoted 24 1IV .– The 25 1 design with defining relation I ABCDE is resolution V and is denoted 25 1V . When choosing a design, we want the highest resolution possible to minimize the number ofinteractions that must be considered negligible in order to obtain a unique interpretation ofthe data.5.32k 2 Fractional Factorial Designs There are k factors of interest each having 2 levels. There are only enough resources to run1/4 of the full factorial 2k design. Thus, we say we want to run a 1/4 fraction of a 2kdesign. This design is called a 2k 2 fractional factorial design.– Example: There are 6 factors of interest (A, B, C, D, E, F ). There are only enoughresources for 16 experimental runs which is 1/4 of a 26 design. Thus, we want to run a1/4 fraction of a 26 design. This design is called a 26 2 fractional factorial design. When selecting a 1/4 fraction, we want to be sure that we select design points that will enableus to estimate effects of interest. Generation of such a design (if it exists) requires carefully choosing two interactions to generatethe design of maximum resolution and then decide on the sign of each generator. These two signed interactions are called the generators of the 2k 2 fractional factorial design,and with their generalized interactions form the complete defining relation for the design.– Suppose there are 6 factors and we choose ABCE and BCDF to be the generators ofthe 26 2 design. Then there are 4 ways to assign signs to these:( , ) I ABCE and I BCDF( , ) I ABCE and I BCDF( , ) I ABCE and I BCDF( , ) I ABCE and I BCDFEach of these four assignments will generate a unique 1/4 fraction of a 26 design.– Suppose we choose the ( , ) case: I ABCE and I BCDF . The generalizedinteraction ABCE · BCDF . Thus, the complete definingrelation is given by I ABCE BCDF .75

The resolution of any design equals the length of the shortest ‘word’ (excluding I) in thecomplete defining relation.– In the 26 2 design with complete defining relation I ABCE BCDF ADEF , thelength of the shortest ‘word’ is four. Therefore, this is a 26 2IV fractional factorial design. A 2k 2 fractional factorial design is formed by selecting one of the four sign pairs ( ( , ),( , ), ( , ) or ( , )) and the selecting only those treatment combinations that matchthe pair of signs in the columns corresponding to the two generating effects. This can beaccomplished in two ways:(i) List all 2k combinations and selecting the rows with plus signs in the two columns corresponding to the two generators.(ii) Or, more simply, list all 2k 2 combinations for a factorial design having k 2 factors.Then create columns for the remaining two factors based on the effects in the completedefining relation. Example: Generate the 26 2 design with defining relation I ABCE BCDF ADEFusing method (ii).– Generate a 24 design in the first four factors (A, B, C, and D). Then, choose aliases forE and F that can be expressed only in terms of A, B, C, and D. Use these to determinethe columns for E and F .E ABC BCDEF ADFF – Now, multiply the appropriate columns to form columns for E and F .A 24 DesignB C D E F 01653752 Because only 1/4 of the full factorial design is run, each of the 2k effects (including the intercept) is aliased with three other effects. That is, estimation of aliased effects are calculatedidentically and, therefore, cannot be separated from each other. The alias structure for any 2k 2 design can be determined by taking the defining relation andmultiplying it by any effect. The resulting four effects are all aliased.76

Example revisited: In the 26 2 design with defining relation I ABCE BCDF ADEF ,the effects are aliased as follows:A BCE ABCDF DEFB ACE CDF ABDEFC ABE BDF ACDEFD E F ABCEF BCD ADEAB CE ACDF BDEFAC BE ABDF CDEFAD BCDE ABCF EFAE AF BD ACDE CF ABEFBF ACEF CD ABDEABD ACD BDE ABF CEFI ABCE BCDF ADEF Thus, when determining a model for analysis, only one effect in each alias class is allowedto be in the model. The complete estimation structure is given below with the main effects and two-factor interactions highlighted in bold-face.l1l2l3l4l5l6l7l8 A ABCE ABCDF DEF B ACE CDF ABDEF C ABE BDF ACDEF D ABCDE BCF AEF E ABC BCDEF ADF F ABCEF BCD ADE AB CE ACDF BDEF AC BE ABDF CDEFl9 AD BCDE ABCF EFl10 AE BC ABCDEF DFl11 AF BCEF ABCD DEl12 BD ACDE CF ABEFl13 BF ACEF CD ABDEl14 ABD CDE ACF BEFl15 ACD BDE ABF CEFl16 I ABCE BCDF ADEF If we restrict consideration to only main effects and two-factor interactions (i.e., assume thatall 3, 4, 5, and 6 factor interactions are negligible) then when estimating effects, we get thesimplified estimation structure:l1l2l3l4l5 A B C D El6 Fl7 AB CEl8 AC BEl9 AD EFl10 AE BC DFl11l12l13l16 AF DE BD CF BF CD I Thus, in l7 to l13 , we can only include one of the aliased effects in the model.Example of a Data Analysis We will now analyze the data from this 26 2 design. The model contains only main effects and a subset of the two-factor interactions. This leavesonly 2 df for the M SE . You can only put one effect from each alias class in the model. For example,– AB and CE are aliased we can put AB or CE in the model (but not both).– AE, DF , and BC are aliased we can put only one of these 3 effects in the model. Pool together the aliased two factor interactions that have very large p-values of .8729, .9744,and .9744 for F A DE, BD CF , and F B CD, respectively. Reanalyze the data. You now have 5 df for the new M SE .77

ANOVA FOR THE 2**(6-2) DESIGNGeneral Linear Models ProcedureDependent Variable: YSourceDFSum ofSquaresMeanSquareModelErrorCorrected R-Square0.985716C.V.25.25055Root MSE6.8966F ValuePr F10.620.0893 -- Only 2 df for ErrorY Mean27.313ANOVA FOR THE 2**(6-2) DESIGNSourceABCDEFA*BA*CA*DA*EF*AB*DF*BDF( CE)( BE)( BF)( DF BC)( DE)( CF)( CD)Type III SS1111111111111Mean SquareF ValuePr 0.562510.5625115.5625115.562514.062514.06251.5625 -1.56250.0625 -0.06250.0625 .07500.68390.25940.64110.87290.97440.9744 -- A -- B -- AB CE - Pool these - first After pooling the three terms into the M SE , we get the following ANOVA. The A and B main effects are highly significant. The other significant effect is the AB CE effect. Because A and B are significant, weconclude that the AB effect is the dominant effect in AB CE.ANOVA FOR THE 2**(6-2) DESIGN WITH POOLED TERMSThe GLM ProcedureDependent Variable: YSourceDFSum ofSquaresModelErrorCorrected .985462Coeff Var16.11088Mean SquareF ValuePr F656.26250019.36250033.890.0006Root MSE4.40028478Y Mean27.31250

SourceABCDEFA*BA*CA*DA*E( CE)( BE)( BF)( DF BC)DFType III SSMean SquareF ValuePr 2.190.160.390.030.0329.130.555.970.730.0015 330 -- A -- B -- AB -- AD or BF DM ’LOG;CLEAR;OUT;CLEAR;’;OPTIONS NODATE NONUMBER PS 60 LS ***;*** 2**(6-2) DESIGN WITH I ABCE BCDF ADEF ****;DATADODODODOIN;D -1 TO 1C -1 TO 1B -1 TO 1A -1 TO 1E A*B*C;F B*C*D;INPUT Y @@;END; END; END;CARDS;6 10 32 60 4 15;BYBYBYBY2;2;2;2;OUTPUT;END;26 60 8 12 34 60 16 5 37 52PROC GLM DATA IN;CLASS A B C D E;MODEL Y A B C D E FA*B A*C A*D A*E A*F B*D B*F / SS3;TITLE ’ANOVA FOR THE 2**(6-2) DESIGN’;PROC GLM DATA IN;CLASS A B C D E;MODEL Y A B C D E F A*B A*C A*D A*E / SS3;TITLE ’ANOVA FOR THE 2**(6-2) DESIGN WITH POOLED TERMS’;RUN;79

80

5.4Using SAS to Generate a Fractional Factorial Design Use SAS to generate a 1/4 fraction of the 26 design (a 26 2 design). SAS used defining relation I BCDE ACDF ABEF .********************************;*** GENERATE A 2**6-2 DESIGN ***;********************************;PROC FACTEX;FACTORS X1 X2 X3 X4 X5 X6 / NLEV 2;SIZE FRACTION 4;EXAMINE ALIASING(3) CONFOUNDING DESIGN;MODEL RESOLUTION 4;OUTPUT OUT FF62X1 NVALS (100 200)X2 NVALS (5 10)X3 CVALS (’ON’ ’OFF’)X4 CVALS (’DAY’ ’NIGHT’)X5 NVALS (0 1)X6 NVALS (50 55);TITLE ’GENERATION OF A 2**(6-2) RESOLUTION IV DESIGN’;PROC PRINT DATA FF62;RUN;GENERATION OF A 2**(6-2) RESOLUTION IV DESIGNDesign -1-1-11611111181

Factor Confounding RulesX5 X2*X3*X4X6 X1*X3*X4So - --I BCDE ACDF ABEFGENERATION OF A 2**(6-2) RESOLUTION IV DESIGNAliasing StructureX1 X2*X5*X6 X3*X4*X6X2 X1*X5*X6 X3*X4*X5X3 X1*X4*X6 X2*X4*X5X4 X1*X3*X6 X2*X3*X5X5 X1*X2*X6 X2*X3*X4X6 X1*X2*X5 X1*X3*X4X1*X2 X5*X6X1*X3 X4*X6X1*X4 X3*X6X1*X5 X2*X6X1*X6 X2*X5 X3*X4X2*X3 X4*X5X2*X4 X3*X5X1*X2*X3 X1*X4*X5 X2*X4*X6 X3*X5*X6X1*X2*X4 X1*X3*X5 X2*X3*X6 5550555050555550505582E BCDF ACD

5.5Using SAS to Generate a Fractional Factorial Design with *************;*** GENERATE A 2**6-1 DESIGN WITH 2 BLOCKS OF SIZE 16 ***********;PROC FACTEX;FACTORS X1 X2 X3 X4 X5 X6 / NLEV 2;BLOCKS NBLOCKS 2;SIZE FRACTION 2;EXAMINE ALIASING(6) CONFOUNDING DESIGN;MODEL RESOLUTION 5;OUTPUT OUT FF61B2X1 NVALS (100 200)X2 NVALS (5 10)X3 CVALS (’ON’ ’OFF’)X4 CVALS (’DAY’ ’NIGHT’)X5 NVALS (0 1)X6 NVALS (50 55)BLOCKNAME SHIFT CVALS (’A.M.’ ’P.M.’);TITLE ’GENERATION OF A 2**(6-1) DESIGN IN 2 BLOCKS’;PROC PRINT DATA FF61B2; RUN; GENERATION OF A 2**(6-1) DESIGN IN 2 BLOCKSDesign 111-1-1132111111283

Factor Confounding RulesX6 X1*X2*X3*X4*X5Block Pseudo-factor Confounding Rules[B1] X3*X4*X5Aliasing Structure0 X1*X2*X3*X4*X5*X6X1 X2*X3*X4*X5*X6X2 X1*X3*X4*X5*X6X3 X1*X2*X4*X5*X6X4 X1*X2*X3*X5*X6X5 X1*X2*X3*X4*X6X6 X1*X2*X3*X4*X5X1*X2 X3*X4*X5*X6X1*X3 X2*X4*X5*X6X1*X4 X2*X3*X5*X6X1*X5 X2*X3*X4*X6X1*X6 X2*X3*X4*X5X2*X3 X1*X4*X5*X6X2*X4 X1*X3*X5*X6X2*X5 X1*X3*X4*X6X2*X6 X1*X3*X4*X5X3*X4 X1*X2*X5*X6X3*X5 X1*X2*X4*X6X3*X6 X1*X2*X4*X5X4*X5 X1*X2*X3*X6X4*X6 X1*X2*X3*X5X5*X6 X1*X2*X3*X4X1*X2*X3 X4*X5*X6X1*X2*X4 X3*X5*X6X1*X2*X5 X3*X4*X6[B] X1*X2*X6 X3*X4*X5X1*X3*X4 X2*X5*X6X1*X3*X5 X2*X4*X6X1*X3*X6 X2*X4*X5X1*X4*X5 X2*X3*X6X1*X4*X6 X2*X3*X5X1*X5*X6 X2*X3*X4GENERATION OF A 2**(6-1) DESIGN IN 2 555555584

5 Two-Level Fractional Factorial Designs Because the number of runs in a 2k factorial design increases rapidly as the number of factors increases, it is often impossible to run the full factorial design given available resources. If the experimenter can reasonably assume that certain high-order interactions (often 3-way

Related Documents:

Types of Experimental Designs! Reducing Cost of Full Factorial Design: " Reduce the no. of levels of each factor. If all factors have 2 levels, we have a 2 k factorial design. " Reduce the number of factors. " Use fractional factorial designs. 4 Types of Experimental Designs! Fractional Factorial Design: " Use a fraction of the full factorial .

2. Nonparametric factorial designs and hypotheses We describe the idea of the nonparametric marginal model and its connection to di erent types of commonly arising factorial designs for longitudinal data. To classify common factorial designs, we introduce a notational s

Abstract: The aim of this paper, to provide the detail survey of historical and recent developments of fractional calculus. Recently fractional calculus has been attracted much attention since it plays an important role in many . Fractional calculus has been 300 years old history, the development of fractional calculus is mainly focused on .

Three-Factor Factorial Designs: Fixed Factors A, B, C 175 Three Factor Factorial Example In a paper production process, the e ects of percentage of hardwood concentration in raw wood pulp, the vat pressure, and the cooking time on the paper strength were studied. There were a 3 levels of hardwood concentration (CONC 2%, 4%, 8%).

Ito formula for the two-parameter fractional Brownian motion 1 Itô formula for the two-parameter fractional Brownian motion . to study the problem of stochastic calculus for two-parameter Gaussian processes. The canonical example of such processes is the the fractional . time, including a Tanaka formula. Section 4 describes the extension .

1 ur Approach4 O 5 CRP-10 Designs 7 Shaker Designs 9 Solid Panel Designs 11 Engineered Panel Designs 13 Mitered Designs 15 Applied Moulding Designs 17 MDF Designs 19 Mullion Designs 21 Slab Designs 22 Decorative Laminate Veneers 23 Thermo Structured Surfaces 25 High Gloss Surfaces 27 Interior Access 28 Drawer Construction 29 Range Hoods 30 Mouldings 31 Incomparable Integration 32 Architectural .

This paper provides a basic introduction to fractional calculus, a branch of mathematical analysis that studies the possibility of taking any real power of the di erentiation operator. We introduce two di erent def-initions of the fractional derivative, namely the Riemann-Liouville and

Trading on SIX Swiss Exchange Introduction 7 159 Sensitivity: Public 1. Introduction SIX Swiss Exchange AG’s (SIX Swiss Exchange) trader training and testing programmes set high