Next Article in Journal
Strength Analysis of Cement Mortar with Carbon Nanotube Dispersion Based on Fractal Dimension of Pore Structure
Previous Article in Journal
Discussion on the Approximate Controllability of Hilfer Fractional Neutral Integro-Differential Inclusions via Almost Sectorial Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of Optimal Split-Plot Designs for Various Design Scenarios

School of Mathematics and Statistics, Shandong Normal University, Jinan 250358, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(10), 608; https://doi.org/10.3390/fractalfract6100608
Submission received: 7 September 2022 / Revised: 16 October 2022 / Accepted: 17 October 2022 / Published: 19 October 2022

Abstract

:
When performing fractional factorial experiments in a completely random order is impractical, fractional factorial split-plot designs are suitable options as an alternative. It is well recognized that the more there are lower order effects of interest at lower order confounding, the better the designs. From this viewpoint, this paper considers the construction of optimal regular two-level fractional factorial split-plot designs. The optimality criteria for two different design scenarios are proposed. Under the newly proposed optimality criteria, the theoretical construction methods of optimal regular two-level fractional factorial split-plot designs are then proposed. In addition, we also explore the theoretical construction methods of some optimal regular two-level fractional factorial split-plot designs under the widely adopted general minimum lower order confounding criterion.

1. Introduction

Regular two-level fractional factorial (FF) designs are commonly used for factorial experiments. When performing an FF design, it is required to perform the experimental runs in a completely random order. However, in some experiments, due to the reasons of being time-consuming or of economic cost, it is impractical or even impossible to perform the FF experimental runs in a completely random manner. For example, consider a modified experiment from [1] in which the purpose is to study the corrosion resistance of steel bars treated with two coatings, say C 1 and C 2 , each at two furnace temperatures, 360 C and 380 C. It takes a long time to reset the furnace and reach a new equilibrium temperature. The factor furnace temperature is called a hard-to-change factor and the factor coating is called an easy-to-change factor. To save experimental time, it is desirable to reduce the times of resetting equilibrium temperature (the hard-to-change factor). To do so, regular two-level fractional factorial split-plot (FFSP) designs are practical design options. For more examples of the experiments which involve hard-to-change factors, one may refer to [1].
For choosing FFSP designs, Ref. [2] proposed the minimum aberration-FFSP (MA-FFSP) criterion by extending the MA criterion proposed in [3] for FF designs. Since then, a large amount of study on MA-FFSP designs has been carried out, including Ref. [4], which discussed the difference between the FF designs and FFSP designs and developed some theories on MA-FFSP designs; Ref. [5], which developed an algorithm for searching optimal MA-FFSP designs; Ref. [6], which studied MA-FFSP designs by developing a finite projective geometric formulation; Ref. [7], which considered the construction of FFSP designs in terms of consulting designs; Ref. [8], which extended the MA criterion to multi-level FFSP designs; Ref. [9], which proposed theoretical construction methods for MA orthogonal split-plot designs; Ref. [10], which considered the design scenario where the whole plot (WP) factors are more important than the sub-plot (SP) factors under the MA criterion; and Ref. [11], which constructed the MA FFSP designs for the design scenario considered in [10] via complementary designs.
According to the effect hierarchy principle and effect sparsity principle (see [12]), main effects and two-factor interactions (2FIs) are always of interest, assuming that the third- and higher-order interactions are negligible. A main effect or 2FI is said to be clear if it is not aliased with any other main effects or 2FIs. Based on the effect hierarchy principle, effect sparsity principle, and the concept of clear effects, some work on choosing optimal FFSP designs were carried out, including Ref. [13], which gave the conditions of an FFSP design to contain clear main effects and 2FIs; Ref. [14], which gave the bounds on the maximum number of clear effects of FFSP designs; Refs. [15,16], which studied the mixed-level FFSP designs with a four-level factor in WP or SP section respectively; Ref. [17] which investigated the conditions for the FFSP designs which involving some two-level factors and an eight-level factor to contain clear effects; Ref. [18], which studied the conditions of FFSP designs with some two-level factors and a 2 t -level factor containing various clear effects; and Ref. [19], which provided the conditions of FFSP designs with some s-level factors and an s t -level factor containing various clear effects.
Apart from the MA and clear effect criterion for the FFSP designs, Ref. [20] extended the general minimum lower order confounding (GMC) criterion for the regular two-level FF designs in [21] to the regular two-level FFSP designs and proposed the GMC-FFSP criterion for assessing the regular two-level FFSP designs. However, the theoretical construction methods of the optimal regular two-level FFSP designs under the GMC-FFSP criterion have not been studied yet.
For a regular two-level FFSP design, the effect involving only WP factors is called a WP effect, and the effect involving at least one SP factor is called an SP effect. The studies on MA orthogonal FFSP designs in [9] were motivated by five different design scenarios; among them, two are presented as follows:
  • Scenario 1: the WP effects and SP effects are equally important.
  • Scenario 2: the SP effects are more important than the WP effects.
In this paper, we investigate the regular two-level FFSP designs for Scenario 1 and Scenario 2 based on a commonly adopted principle that the more there are lower order effects of interest at the lower order confounding, the better the regular two-level FFSP designs. This viewpoint is different from that considered in [9]. In addition, this paper also considers constructing optimal regular two-level FFSP designs under the GMC-FFSP criterion. The contributions of this paper are threefold:
(1)
We develop suitable optimality criteria for choosing regular two-level FFSP designs for Scenarios 1 and Scenario 2 based on the assumption that the effects involving more than two factors are negligible.
(2)
The construction methods of the optimal regular two-level FFSP designs under the newly proposed optimality criteria are provided.
(3)
The construction methods of some optimal regular two-level FFSP designs under the GMC-FFSP criterion are derived.
The rest of the paper is organized as follows. Section 2 includes some useful notation, definitions, and the development of the optimality criteria for designs for Scenario 1 and Scenario 2, respectively. The construction of some optimal regular two-level FFSP designs are provided in Section 3. Conclusions are given in Section 4.

2. Optimality Criteria, Notation and Definitions

Let k 1 = n 1 m 1 , k 2 = n 2 m 2 , k = k 1 + k 2 , and N = 2 k . Throughout the paper, we use the notation 2 ( n 1 + n 2 ) ( m 1 + m 2 ) to denote a regular two-level FFSP design with n 1 WP factors/columns, n 2 SP factors/columns, and N runs. Since the factors are assigned to columns of designs, we do not differentiate between factors and columns. Denote a 1 , a 2 , , a k 1 , b 1 , b 2 , , b k 2 as k independent 2 k × 1 columns at + 1 and 1 levels. The saturated design H = H ( a 1 , a 2 , , a k 1 , b 1 , b 2 , , b k 2 ) with 2 k runs and 2 k 1 columns can be obtained by taking all possible component-wise products among the k independent columns. Let H a = H ( a 1 , a 2 , , a k 1 ) , without special statement; the columns in H and H a are placed one after another in Yates order, i.e.,
H = { a 1 , a 2 , a 1 a 2 , a 3 , a 1 a 3 , , a 1 a 2 a k 1 , b 1 , a 1 b 1 , a 1 a 2 b 1 , , a 1 a k 1 b 1 b k 2 } , H a = { a 1 , a 2 , a 1 a 2 , a 3 , a 1 a 3 , , a 1 a 2 a k 1 } .
Let S H and γ H ; then we denote B i ( S , γ ) = # { d 1 , , d i : d 1 , , d i S , d 1 d i = γ } and g ¯ ( S ) = # { γ : γ H \ S , B 2 ( S , γ ) > 0 } , where # denotes the cardinality of a set, d 1 , , d i are mutually different columns in S and d 1 d i is the column genarated by taking component-wise products of columns d 1 , , d i . Let T = ( T W , T S ) denote a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with T W = { a 1 , a 2 , , a k 1 , a k 1 + 1 , , a n 1 } and T S = { b 1 , b 2 , , b k 2 , b k 2 + 1 , , b n 2 } , where T W and T S denote the WP section and SP section in the 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design, respectively. It is worth noting that we have set T W to contain k 1 independent columns and T S to contain k 2 independent columns here. Given any k independent columns a 1 , a 2 , , a k 1 , b 1 , b 2 , , b k 2 , choosing a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design is equal to choosing m ( = m 1 + m 2 ) more columns a k 1 + 1 ,   , a n 1 , b k 2 + 1 , , b n 2 from H. Certainly, the m columns a k 1 + 1 , , a n 1 , b k 2 + 1 , , b n 2 can be generated by some of the previously stated k independent columns.
Let 1 # C 2 ( k ) ( T ) denote the number of main effects which are aliased with k 2FIs, where k = 0 , 1 , , K with K = n 2 . Let 2 # C 2 ( k ) ( T ) denote the number of 2FIs which are aliased with k 2FIs, where k = 0 , 1 , , K 1 . Let 1 ( s ) # C ( w ) ( 0 ) and 1 ( s ) # C ( w ) ( 1 ) denote the number of SP main effects which are not aliased with any WP effects, and the number of SP main effects which are aliased with at least one WP effect, respectively. Let 2 ( s ) # C ( w ) ( 0 ) and 2 ( s ) # C ( w ) ( 1 ) denote the number of SP 2FIs which are not aliased with any WP effects, and the number of SP 2FIs which are aliased with at least one WP effect, respectively. With these notation, we provide the optimality criteria for choosing 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 1 and Scenario 2, respectively, as follows. The 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs which can sequentially maximize
1 # C ( T ) = ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) , 2 # C 2 ( T ) ) ,
are optimal for Scenario 1, where 1 # C 2 ( T ) = ( 1 # C 2 ( 0 ) ( T ) , 1 # C 2 ( 1 ) ( T ) , , 1 # C 2 ( K ) ( T ) ) and 2 # C 2 ( T ) = ( 2 # C 2 ( 0 ) ( T ) , 2 # C 2 ( 1 ) ( T ) , , 2 # C 2 ( K 1 ) ( T ) ) . The 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs which can sequentially maximize
2 # C ( T ) = ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) , 2 ( s ) # C ( w ) ( 0 ) ( T ) )
are optimal for Scenario 2. By combining (1) and (2), the 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs which can sequentially maximize
3 # C ( T ) = ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) , 2 # C 2 ( T ) , 2 ( s ) # C ( w ) ( 0 ) ( T ) )
are optimal under the GMC-FFSP criterion. Let 2 n m denote a regular two-level FF design with n columns and N = 2 n m runs. For a 2 n m design D, the notation 1 # C 2 ( k ) ( D ) and 2 # C 2 ( k ) ( D ) have the same meanings as 1 # C 2 ( k ) ( T ) and 2 # C 2 ( k ) ( T ) , respectively. A 2 n m design D which can sequentially maximize
( 1 # C 2 ( D ) , 2 # C 2 ( D ) )
is optimal under the GMC criterion. To avoid confusion, hereafter, we use the expression GMC-FF instead of GMC to present the contents relative to the 2 n m designs.
Before introducing the theoretical results of this work, we introduce some more notation. Let F a = F ( a 1 , a 2 , , a k 1 ) be the set of columns which are the component-wise products of all possible odd number of columns among the k 1 independent columns a 1 , a 2 , , a k 1 , i.e., F a = { a 1 , a 2 , a 3 , a 1 a 2 a 3 , a 4 , a 1 a 2 a 4 , a 1 a 3 a 4 , a 2 a 3 a 4 , } . The set F b = F ( b 1 , b 2 , , b k 2 ) and F a b = F ( a 1 , a 2 , , a k 1 , b 1 , b 2 , , b k 2 ) are similarly defined. Denote G a b = F a b \ F a . The columns in F a , F a b and G a b are placed in Yates order, respectively. For any two sets A and B of columns from H, the notation A B denotes the set which consists of all the mutually different columns generated by taking component-wise products between two columns in which one is from A and the other is from B. In [13], it is stated that T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design if and only if
T W H a , T S H \ H a , # ( T W ) = n 1 , # ( T S ) = n 2 ,
where # ( · ) denotes the number of columns in a design.

3. Construction of Optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) Designs

A 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design is said to have resolution R if no c-factor interaction is aliased with any other interaction involving fewer than R c factors. The resolution III 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs have at least one main effect which is aliased with at least one 2FI. In the resolution R=IV 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs, all the main effects are clear but there is at least one 2FI which is aliased with at least one 2FI. In Section 3.1Section 3.3, we provide the construction methods of some optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 1, Scenario 2, and under the GMC-FFSP criterion.

3.1. Construction Methods of Optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) Designs for Scenario 1

We first provide a lemma which generalizes the construction of GMC-FF 2 n m designs for given n and m with 5 N 16 + 1 n N 2 . Theorems 1 and 2 provide the construction methods of some optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 1.
Lemma 1.
For k 1 2 , suppose D is a 2 n m design with respect to
2 k 1 2 + 1 n 1 2 k 1 1 , n 2 = t = s k 2 2 t f o r k 1 1 s k 2 a n d 5 N 16 + 1 n N 2 .
If D consists of the first n 1 columns of F a and the last n 2 columns of G a b , then D is optimal under the GMC-FF criterion.
Proof. 
According to [22,23], a 2 n m design D with D F a b must has resolution at least IV. Therefore, 1 # C 2 ( 0 ) ( D ) = n , and 1 # C 2 ( D ) is sequentially maximized. Next, we prove that 2 # C 2 ( D ) is sequentially maximized among all the 2 n m designs with respect to (6).
Suppose E is a 2 n m design which consists of the first n columns of F a b . According to [24], E is a GMC-FF design which sequentially maximizes 2 # C 2 ( E ) among all the 2 n m designs with respect to (6). Let r = n / 2 k 1 1 . Write D = ( D ¯ 1 , D 1 ) , where D 1 contains the last r × 2 k 1 1 columns of D, and D ¯ 1 = D \ D 1 . Write E = ( E 1 , E ¯ 1 ) , where E 1 contains the first r × 2 k 1 1 columns of E, and E ¯ 1 = E \ E 1 . We can always find γ 1 , γ 2 H \ F a b such that D 1 = γ 1 E 1 and H \ D 1 = γ 2 ( H \ E 1 ) , implying that D ¯ 1 = γ 2 E ¯ 1 . Rewrite D 1 as D 1 = { d 1 F a , d 2 F a , , d r F a } , where d 1 is the grand mean and d 2 , , d r are from H \ F a b . Rewrite E 1 as E 1 = { e 1 F a , e 2 F a , , e r F a } , where e 1 , e 2 , , e r are from H \ F a b . Actually, there exists the facts that
(1)
D 1 D 1 = E 1 E 1 = H \ F a b ,
(2)
D 1 D ¯ 1 = E 1 E ¯ 1 D 1 D 1 ( = E 1 E 1 ) ,
(3)
D ¯ 1 D ¯ 1 = E ¯ 1 E ¯ 1 D 1 D 1 ( = E 1 E 1 ) , and
(4)
( D ¯ 1 D ¯ 1 ) ( D 1 D ¯ 1 ) =
due to the following reasons.
For (1). According to Lemma A.3 in [25], since 2 k 2 + 1 # ( D 1 ) 2 k 1 and D 1 has k independent columns, then D 1 D 1 = H \ F a b . Similarly, we can also obtain E 1 E 1 = H \ F a b .
For (2). Let l 1 denote the first column of D ¯ 1 , then
l 1 ( F a b \ l 1 ) = ( l 1 ( ( F a b \ D 1 ) \ l 1 ) ) ( l 1 D 1 ) = ( l 1 ( ( F a b \ D 1 ) \ l 1 ) ) ( D ¯ 1 D 1 ) = H \ F a b ,
where the second equality is because l 1 D 1 = D ¯ 1 D 1 due to D ¯ 1 F a and the structure of D 1 . Therefore, D ¯ 1 D 1 D 1 D 1 . Similarly, we obtain that
q 1 ( F a b \ q 1 ) = ( q 1 ( ( F a b \ E 1 ) \ q 1 ) ) ( q 1 E 1 ) = ( q 1 ( ( F a b \ E 1 ) \ q 1 ) ) ( E ¯ 1 E 1 ) = H \ F a b
and E ¯ 1 E 1 D 1 D 1 ( = E 1 E 1 ) , where q 1 is the first column in E ¯ 1 . Note that ( D ¯ 1 D 1 ) ( D ¯ 1 ( F a b \ D ) ) = H \ F a b , and ( D ¯ 1 D 1 ) ( D ¯ 1 ( F a b \ D ) ) = . Similarly, there exists ( E ¯ 1 E 1 ) ( E ¯ 1 ( F a b \ E ) ) = H \ F a b and ( E ¯ 1 E 1 ) ( E ¯ 1 ( F a b \ E ) ) = . Since E ¯ 1 ( F a b \ E ) = D ¯ 1 ( F a b \ D ) as F a b \ E = γ 2 ( F a b \ D ) , we have D ¯ 1 D 1 = E ¯ 1 E 1 . This obtains the fact (2).
For (3). Since D ¯ 1 F a and D ¯ 1 = γ 2 E ¯ 1 , it is easy to obtain that D ¯ 1 D ¯ 1 = E ¯ 1 E ¯ 1 H a \ F a . This completes the proof for (3).
For (4). Note that D ¯ 1 D ¯ 1 H a \ F a and any two-column interaction with one column from D ¯ 1 and the other from D 1 is not in H a \ F a . Therefore, ( D ¯ 1 D ¯ 1 ) ( D 1 D ¯ 1 ) = .
Based on the analysis above, the 2FIs of D and E can be classified into three disjoint groups, respectively, as
G 1 :
D 1 D ¯ 1 = E 1 E ¯ 1 ,
G 2 :
D ¯ 1 D ¯ 1 = E ¯ 1 E ¯ 1 and
G 3 :
( D 1 D 1 ) \ ( ( D 1 D ¯ 1 ) ( D ¯ 1 D ¯ 1 ) ) = ( E 1 E 1 ) \ ( ( E 1 E ¯ 1 ) ( E ¯ 1 E ¯ 1 ) ) .
From (1) and (2), for any γ G 1 , there are # ( D ¯ 1 ) two-column pairs ( α 1 , β 1 ) with α 1 D 1 and β 1 D ¯ 1 such that γ = α 1 β 1 , and there are # ( E ¯ 1 ) two-column pairs ( α 2 , β 2 ) with α 2 E 1 and β 2 E ¯ 1 such that γ = α 2 β 2 , where # ( D ¯ 1 ) = # ( E ¯ 1 ) ; if there are t 1 two-column pairs ( α 1 , β 1 ) with α 1 D 1 and β 1 D 1 such that γ = α 1 β 1 , there must be t 1 two-column pairs ( α 2 , β 2 ) with α 2 E 1 and β 2 E 1 such that γ = α 2 β 2 due to D 1 = γ 1 E 1 .
From (1) and (3), for any γ G 2 , if there are t 3 two-column pairs ( α 1 , β 1 ) with α 1 D ¯ 1 and β 1 D ¯ 1 such that γ = α 1 β 1 , there must be t 3 two-column pairs ( α 2 , β 2 ) with α 2 E ¯ 1 and β 2 E ¯ 1 such that γ = α 2 β 2 , due to that D ¯ 1 = γ 2 E ¯ 1 ; if there are t 4 two-column pairs ( α 1 , β 1 ) with α 1 D 1 and β 1 D 1 such that γ = α 1 β 1 , there must be t 4 two-column pairs ( α 2 , β 2 ) with α 2 E 1 and β 2 E 1 such that γ = α 2 β 2 due to that D 1 = γ 1 E 1 .
For any γ G 3 , if there are t 5 two-column pairs ( α 1 , β 1 ) with α 1 D 1 and β 1 D 1 such that γ = α 1 β 1 , there must be t 5 two-column pairs ( α 2 , β 2 ) with α 2 E 1 and β 2 E 1 such that γ = α 2 β 2 due to that D 1 = γ 1 E 1 .
Therefore, we have 2 # C 2 ( D ) = 2 # C 2 ( E ) which is sequentially maximized among all the 2 n m designs with respect to (6) as E is a GMC-FF design according to [24]. This completes the proof. □
Remark 1.
In [24], it is stated that a 2 n m design with 5 N 16 + 1 n N 2 is a GMC-FF design if this design consists of the first (or last) n columns of F a b . Lemma 1 generalizes their construction methods for GMC-FF 2 n m designs with 5 N 16 + 1 n N 2 .
Based on Lemma 1, the following Theorems 1 and 2 provide construction methods of some optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 1.
Theorem 1.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with respect to
2 k 1 2 + 1 n 1 2 k 1 1 , n 2 = t = s k 2 2 t f o r k 1 1 s k 2 a n d 5 N 16 + 1 n N 2 .
If T W consists of the first n 1 columns of F a and T S consists of the last n 2 columns of G a b , then T = ( T W , T S ) is optimal for Scenario 1.
Proof. 
Clearly, T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design as it satisfies (5); thus, 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 . According to Lemma 1, we obtain that T can sequentially maximize ( 1 # C 2 ( T ) , 2 # C 2 ( T ) ) . This completes the proof. □
Example 1 shows the application of Theorem 1.
Example 1.
Consider constructing an optimal 2 ( 6 + 8 ) ( 2 + 7 ) design for Scenario 1. Without loss of generality, let a 1 = 5 , a 2 = 15 , a 3 = 25 , a 4 = 35 and b 1 = 45 , then F a = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 } and G a b = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . Let T W = { 5 , 15 , 25 , 125 , 35 , 135 } and T S = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . According to Theorem 1, T = ( T W , T S ) is optimal for Scenario 1.
Theorem 2.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with n 1 = 2 k 1 1 , n 2 2 k 1 2 k 1 1 and 5 N 16 + 1 n N 2 . Let T W = F a and T S consists of the first n 2 columns of G a b , then T = ( T W , T S ) is optimal for Scenario 1.
Proof. 
Clearly, the design T in this theorem is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design; thus, 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 . Note that T consists of the first n columns of F a b ; thus, T sequentially maximizes ( 1 # C 2 ( T ) , 2 # C 2 ( T ) ) as it is also a GMC-FF design according to [24]. This completes the proof. □
Example 2.
Consider constructing an optimal 2 ( 4 + 8 ) ( 1 + 6 ) design for Scenario 1. Without loss of generality, let a 1 = 5 , a 2 = 15 , a 3 = 25 , b 1 = 35 and b 2 = 45 , then F a = { 5 , 15 , 25 , 125 } and G a b = { 35 , 135 , 235 , 1235 , 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . Let T W = { 5 , 15 , 25 , 125 } and T S = { 35 , 135 , 235 , 1235 , 45 , 145 , 245 , 1245 } . According to Theorem 2, T = ( T W , T S ) is optimal for Scenario 1.
In Theorem 3, we build the connection between GMC-FF 2 n m designs and the optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 1. Before introducing Theorem 3, we first give a useful lemma.
Lemma 2.
Suppose D and B are two 2 n m designs from F a b . If D can be divided into two disjoint parts D 1 and D 2 such that
(i)
B 1 = γ 1 D 1 , B 2 = γ 2 D 2 and B = B 1 B 2 with B 1 B 2 = ,
(ii)
( D 1 D 2 ) ( ( D 1 D 1 ) ( D 2 D 2 ) ) = , and
(iii)
( B 1 B 2 ) ( ( B 1 B 1 ) ( B 2 B 2 ) ) = ,
then ( g ¯ ( D ) , 1 # C 2 ( D ) , 2 # C 2 ( D ) ) = ( g ¯ ( B ) , 1 # C 2 ( B ) , 2 # C 2 ( B ) ) , where each of γ 1 and γ 2 can be the grand mean or any column from H \ F a b , and ∅ denotes the empty set.
Proof. 
Since B 1 = γ 1 D 1 and B 2 = γ 2 D 2 , we have that D 1 D 1 = B 1 B 1 , D 2 D 2 = B 2 B 2 , and γ 1 γ 2 ( D 1 D 2 ) = B 1 B 2 . More specifically, if there are t 1 two-column pairs ( α 1 , α 2 ) with α 1 D 1 and α 2 D 1 such that ν = α 1 α 2 , then there are must be t 1 two-column pairs ( β 1 , β 2 ) with β 1 B 1 and β 2 B 1 such that ν = β 1 β 2 ; for any ν D 2 D 2 , if there are t 2 two-column pairs ( α 1 , α 2 ) with α 1 D 2 and α 2 D 2 such that ν = α 1 α 2 , then there must be t 2 two-column pairs ( β 1 , β 2 ) with β 1 B 2 and β 2 B 2 such that ν = β 1 β 2 ; for any ν D 1 D 2 , if there are t 3 two-column pairs ( α 1 , α 2 ) with α 1 D 1 and α 2 D 2 such that ν = α 1 α 2 , then there must be t 3 two-column pairs ( β 1 , β 2 ) with β 1 B 1 and β 2 B 2 such that γ 1 γ 2 ν = β 1 β 2 .
With the analysis above, we first prove that g ¯ ( D ) = g ¯ ( B ) . Recalling the definition of g ¯ ( D ) , we have
g ¯ ( D ) = # { ν : ν H \ D , B 2 ( D 1 D 2 , ν ) > 0 } = # { ν : ν H \ D , ν ( D 1 D 1 ) ( D 2 D 2 ) ( D 1 D 2 ) } = # { ν : ν H \ D , ν ( D 1 D 1 ) ( D 2 D 2 ) } + # { ν : ν H \ D , ν D 1 D 2 } = # { τ : τ H \ B , τ ( B 1 B 1 ) ( B 2 B 2 ) } + # { τ : τ H \ B , τ B 1 B 2 } , = g ¯ ( B )
where in the fourth equality # { τ : τ H \ B , τ B 1 B 2 } = # { ν : ν H \ D , ν D 1 D 2 } is due to the fact that for any ν 0 H \ D with ν 0 D 1 D 2 we have τ 0 = γ 1 γ 2 ν 0 H \ B . This obtains that g ¯ ( D ) = g ¯ ( B ) .
Since any 2 n m design from F a b has resolution IV, then 1 # C 2 ( D ) = 1 # C 2 ( B ) .
Next, we give the proof that 2 # C 2 ( D ) = 2 # C 2 ( B ) . According to the analysis in the first paragraph, for any ν 0 = α 1 α 2 ( D 1 D 1 ) ( D 2 D 2 ) , we have τ 0 = ν 0 = ( γ i α 1 ) ( γ i α 2 ) ( B 1 B 1 ) ( B 2 B 2 ) . Therefore, we have
# { ν : ν ( D 1 D 1 ) ( D 2 D 2 ) , B 2 ( D , ν ) = k } = # { τ : τ ( B 1 B 1 ) ( B 2 B 2 ) , B 2 ( B , τ ) = k } ,
where k = 0 , 1 , , K . Similarly, for any ν 0 = α 1 α 2 D 1 D 2 , we have τ 0 = γ 1 γ 2 ν 0 = ( γ 1 α 1 ) ( γ 2 α 2 ) B 1 B 2 . Therefore, we have
# { ν : ν D 1 D 2 , B 2 ( D , ν ) = k } = # { τ : τ B 1 B 2 , B 2 ( B , τ ) = k } ,
where k = 0 , 1 , , K . This obtains that 2 # C 2 ( D ) = 2 # C 2 ( B ) and the proof is completed. □
With Lemma 2, we immediately obtain Theorem 3, which connects optimal FFSP designs for Scenario 1 with GMC-FF 2 n m designs.
Theorem 3.
Suppose T = ( T W , T S ) F a b and B F a b are 2 ( n 1 + n 2 ) ( m 1 + m 2 ) and GMC-FF 2 n m designs with 5 N 16 + 1 n N 2 , respectively. For T ¯ = F a b \ T and B ¯ F a b \ B , if T ¯ can be divided into two disjoint parts T 1 and T 2 such that
(i)
B 1 = γ 1 T 1 , B 2 = γ 2 T 2 and B ¯ = B 1 B 2 with B 1 B 2 = ;
(ii)
( T 1 T 2 ) ( ( T 1 T 1 ) ( T 2 T 2 ) ) = , and
(iii)
( B 1 B 2 ) ( ( B 1 B 1 ) ( B 2 B 2 ) ) = ,
then T is optimal for Scenario 1, where each of γ 1 and γ 2 can be the grand mean or any column from H \ F a b .
Proof. 
On one hand, according to Lemma 1 of [24], sequentially maximizing 2 # C 2 ( T ) is equal to sequentially maximizing ( g ¯ ( T ¯ ) , 2 # C 2 ( T ¯ ) ) . On the other hand, according to Lemma 2, we obtain that ( g ¯ ( T ¯ ) , 2 # C 2 ( T ¯ ) ) = ( g ¯ ( B ¯ ) , 2 # C 2 ( B ¯ ) ) indicating that ( g ¯ ( T ¯ ) , 2 # C 2 ( T ¯ ) ) is sequentially maximized. This is because ( g ¯ ( B ¯ ) , 2 # C 2 ( B ¯ ) ) is sequentially maximized among all the 2 n m designs with 5 N 16 + 1 n N 2 . Therefore, we obtain that T can sequentially maximize (1) among all the 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs and thus it is optimal for Scenario 1. □
Theorem 3 provides an approach to conforming that a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design is optimal for Scenario 1. The following example illustrates the application of Theorem 3.
Example 3.
For a given 2 ( 6 + 6 ) ( 2 + 5 ) design T = ( T W , T S ) with T W = { 5 , 15 , 25 , 125 , 35 , 135 } and T S = { 45 , 145 , 245 , 1245 , 345 , 1345 } , we have T ¯ = F a b \ T = { 235 , 1235 , 2345 , 12345 } . Divide T ¯ into two disjoint subsets as T ¯ = T 1 T 2 with T 1 = { 235 , 1235 } and T 2 = { 2345 , 12345 } , then T 1 and T 2 satisfy ( ( T 1 T 1 ) ( T 2 T 2 ) ) ( T 1 T 2 ) = . Let B 1 = 24 T 1 = { 345 , 1345 } , B 2 = T 2 and B ¯ = B 1 B 2 then B = F a b \ B ¯ = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 , 45 , 145 , 245 , 1245 } which is composed of the first 12 columns of F a b . According to Theorem 3, we obtain that T sequentially maximizes (1) among all 2 ( 6 + 6 ) ( 2 + 5 ) FFSP designs. Therefore, design T is optimal for Scenario 1.

3.2. Construction Methods of Optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) Designs for Scenario 2

Lemmas 3 and 4 below derive some properties for 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs which is useful for deriving the construction methods of optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for Scenario 2.
Lemma 3.
For any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design T = ( T W , T S ) , there must be n 1 n 2 2 ( s ) # C ( w ) ( 0 ) ( T ) n 2 2 + n 1 n 2 .
Proof. 
For any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design, the number of 2FIs which have two SP factors and the number of 2FIs which have only one SP factor are n 2 2 and n 1 n 2 , respectively. Therefore, 2 ( s ) # C ( w ) ( 0 ) n 2 2 + n 1 n 2 . As aforementioned, the generator which contains only one SP factor is not allowed, implying that all the 2FIs which have only one SP factor are not aliased with any WP effects. Therefore, we have n 1 n 2 1 ( s ) # C ( w ) ( 0 ) ( T ) . This completes the proof. □
Lemma 4.
For any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design T = ( T W , T S ) with k 2 = 1 , there must be 2 ( s ) # C ( w ) ( 0 ) ( T ) = n 1 n 2 .
Proof. 
The formula k 2 = 1 indicates that there is only one independent SP factor denoted as b 1 . Therefore, the SP dependent factors b 2 , b 3 , , b n 2 can be expressed as b i = b 1 a i 1 a i 2 a i j , where i = 2 , 3 , , n 2 and i 1 , i 2 , , i j = 1 , 2 , , k 1 . Therefore, all of the n 2 2 2FIs which contain two SP factors are aliased with WP effects. As aforementioned, for any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design, the 2FIs which contain only one SP factor are not aliased with any WP effects. Therefore, we have 2 ( s ) # C ( w ) ( 0 ) ( T ) = n 1 n 2 . This completes the proof. □
With Lemma 3, Theorems 4 blow provides construction methods of some FFSP designs which are optimal for Scenario 2.
Theorem 4.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with n 1 2 k 1 1 and n 2 = k 2 , i.e., m 2 = 0 , if T W F a and T S G a b , then T = ( T W , T S ) is optimal for Scenario 2.
Proof. 
Note that T F a b , then T has resolution at least IV. Therefore, T sequentially maximizes ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) ) . The formula m 2 = 0 implies that no SP 2FI is aliased with WP effects meaning that 2 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 2 + n 1 n 2 which is the upper bound of 2 ( s ) # C ( w ) ( 0 ) ( · ) . Therefore, T sequentially maximizes ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) , 2 ( s ) # C ( w ) ( 0 ) ( T ) ) meaning that T is optimal for Scenario 2. □
Example 4.
Consider constructing a 2 ( 4 + 2 ) ( 1 + 0 ) design which is optimal for Scenario 2. Without loss of generality, we set a 1 = 5 , a 2 = 15 , a 3 = 25 , b 1 = 35 and b 2 = 45 . Let T W = { 5 , 15 , 25 , 125 } and T S = { 35 , 45 } . According to Theorem 4, the design T = ( T W , T S ) is an optimal 2 ( 4 + 2 ) ( 1 + 0 ) design for Scenario 2.
With Lemma 3, we obtain Theorem 5 below.
Theorem 5.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with n 1 2 k 1 1 , n 2 2 k 2 1 and m 2 1 . Let T W F a and T S F ( a 1 , b 1 , , b k 2 ) \ a 1 , then T = ( T W , T S ) is optimal for Scenario 2.
Proof. 
Clearly, T is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design as T W H a and T S F a b \ H a . Since any two-column interaction of F ( a 1 , b 1 , , b k 2 ) \ a 1 is not in H a , and any two-column interaction with one column from F a and the other from F ( a 1 , b 1 , , b k 2 ) \ a 1 is not in H a , then T has no SP 2FI which is aliased with any WP effects. Therefore, we have 2 ( s ) # C ( w ) ( 0 ) ( T ) = n 1 n 2 + n 2 2 which is the upper bound for every 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design according to Lemma 3. This completes the proof, noting that T sequentially maximizes 1 # C 2 ( T ) due to its resolution IV. □
Example 5 below illustrates the application of Theorem 5.
Example 5.
Consider constructing a 2 ( 2 + 7 ) ( 0 + 4 ) design which is optimal for Scenario 2. Without loss of generality, we set a 1 = 5 , a 2 = 15 , b 1 = 25 , b 2 = 35 , and b 3 = 45 . Then F a = { 5 , 15 } and F ( a 1 , b 1 , , b k 2 ) \ a 1 = { 25 , 35 , 45 , 235 , 245 , 345 , 2345 } . According to Theorem 5, any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design T = ( T W , T S ) with T W F a and T S F ( a 1 , b 1 , , b k 2 ) \ a 1 is an optimal 2 ( 2 + 7 ) ( 0 + 4 ) design for Scenario 2.
With Theorems 4 and 5, the following corollary is obtained.
Corollary 1.
The 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs constructed by Theorems 4 and 5 have 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( t ) ( T ) = 0 for t = 1 , 2 , , K , and 2 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 2 + n 1 n 2 .
With Lemma 4, we can immediately obtain the results in Theorem 6.
Theorem 6.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with k 1 n 1 2 k 2 , n 2 2 k 2 and m 2 = n 2 1 . Let T W F a and T S contains any n 2 columns of G a b , then T is optimal for Scenario 2.
Example 6.
Consider constructing a 2 ( 5 + 2 ) ( 1 + 1 ) design which is optimal for Scenario 2. Without loss of generality, we set a 1 = 5 , a 2 = 15 , a 3 = 25 , a 4 = 35 and b 1 = 45 . Then F a = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 } and G a b = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . According to Theorem 6, any 2 ( n 1 + n 2 ) ( m 1 + ( n 2 1 ) ) design T = ( T W , T S ) with T W F a and T S G a b is an optimal 2 ( 5 + 2 ) ( 1 + 1 ) design for Scenario 2. Without loss of generality, let T W = { 5 , 15 , 25 , 125 , 35 } and T S = { 145 , 245 } , then T = ( T W , T S ) is optimal for Scenario 2.

3.3. Construction Methods of GMC-FFSP 2 ( n 1 + n 2 ) ( m 1 + m 2 ) Designs

With Theorem 1 and Lemma 4, we immediately obtain Theorem 7 below, which constructs some GMC-FFSP 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs.
Theorem 7.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with 2 k 3 + 1 n 1 2 k 2 , n 2 = 2 k 2 , 5 N 16 + 1 n N 2 and m 2 = n 2 1 . If T W consists of the first n 1 columns of F a and T S = G a b , then T is a GMC-FFSP design.
Example 7.
Consider constructing a 2 ( 5 + 8 ) ( 1 + 7 ) GMC-FFSP design by Thereom 7. Without loss of generality, we set a 1 = 5 , a 2 = 15 , a 3 = 25 , a 4 = 35 and b 1 = 45 . Then F a = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 } and G a b = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . Let T W = { 5 , 15 , 25 , 125 , 35 } and T S = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } , then T = ( T W , T S ) is a 2 ( 5 + 8 ) ( 1 + 7 ) GMC-FFSP design.
With Theorem 2 and Lemma 4, Theorem 8 below provides construction methods of some GMC-FFSP 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs.
Theorem 8.
Suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with n 1 = 2 k 2 , n 2 2 k 2 , 5 N 16 + 1 n N 2 and m 2 = n 2 1 . If T W = F a and T S consists of the first n 2 columns of G a b , then T is a GMC-FFSP design.
Proof. 
The formula n 1 = 2 k 1 indicates that T W consists of k 1 independent columns, i.e., k 1 = k 1 . Therefore, we have k 2 = 1 . In Theorem 2, it is proved that T can sequentially maximize ( 1 ( s ) # C ( w ) ( 0 ) ( T ) = n 2 , 1 # C 2 ( T ) , 2 # C 2 ( T ) ) . According to Lemma 4, for any 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with k 2 = 1 , we have 2 ( s ) # C ( w ) ( 0 ) ( · ) = n 1 n 2 . This completes the proof. □
Example 8.
Consider constructing a 2 ( 8 + 3 ) ( 4 + 2 ) GMC-FFSP design by Theoreom 8. Without loss of generality, we set a 1 = 5 , a 2 = 15 , a 3 = 25 , a 4 = 35 and b 1 = 45 . Then F a = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 } and G a b = { 45 , 145 , 245 , 1245 , 345 , 1345 , 2345 , 12345 } . Let T W = { 5 , 15 , 25 , 125 , 35 , 135 , 235 , 1235 } and T S = { 45 , 145 , 245 } , then T = ( T W , T S ) is a 2 ( 8 + 3 ) ( 4 + 2 ) GMC-FFSP design.
Similar to Theorem 3, the theorem below provides an approach to conforming that some 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs are GMC-FFSP designs.
Theorem 9.
For 5 N 16 + 1 n N 2 , suppose T = ( T W , T S ) is a 2 ( n 1 + n 2 ) ( m 1 + m 2 ) design with T F a b and m 2 = n 2 1 . If there exists a GMC-FF design D F a b such that
(i)
T 1 = γ 1 D 1 , T 2 = γ 2 D 2 and T ¯ = T 1 T 2 with T 1 T 2 = ;
(ii)
( D 1 D 2 ) { ( D 1 D 1 ) ( D 2 D 2 ) } = , and
(iii)
( T 1 T 2 ) { ( T 1 T 1 ) ( T 2 T 2 ) } = ,
then T is a GMC-FFSP design, where T ¯ = F a b \ T , D ¯ = F a b \ D , D ¯ = D 1 D 2 with D 1 D 2 = , each of γ 1 and γ 2 can be the grand mean or any column from H \ F a b , and ∅ denotes the empty set.
Example 9.
For a given 2 ( 12 + 12 ) ( 7 + 11 ) design T = ( T W , T S ) with T W = { 6 , 16 , 26 , 126 , 36 , 136 , 236 , 1236 , 46 , 146 , 246 , 1246 } and T S = { 56 , 156 , 256 , 1256 , 356 , 1356 , 2356 , 12356 , 456 , 1456 , 2456 , 12456 } , we have T ¯ = F a b \ T = { 346 , 1346 , 2346 , 12346 , 3456 , 13456 , 23456 , 123456 } . Divide T ¯ into two disjoint subsets as T ¯ = T 1 T 2 with T 1 = { 346 , 1346 , 2346 , 12346 } and T 2 = { 3456 , 13456 , 23456 , 123456 } , then T 1 and T 2 satisfy ( ( T 1 T 1 ) ( T 2 T 2 ) ) ( T 1 T 2 ) = . Let D 1 = 35 T 1 = { 456 , 1456 , 2456 , 12456 } , D 2 = T 2 and D ¯ = D 1 D 2 then D = F a b \ D ¯ = { 6 , 16 , 26 , 126 , 36 , 136 , 236 , 1236 , 46 , 146 , 246 , 1246 , 346 , 1346 , 2346 , 12346 , 56 , 156 , 256 , 1256 , 356 , 1356 , 2356 , 12356 } which is composed of the first 24 columns of F a b . According to Theorem 9, we obtain that T is a 2 ( 12 + 12 ) ( 7 + 11 ) GMC-FFSP design.

3.4. Some More Illustrative Examples and Further Discussions

In this section, we provide some more examples to illustrate how to recognize the superiority of an FFSP design over another under criteria (1), (2), and (3), respectively.
Consider the following two 2 ( 2 + 7 ) ( 0 + 4 ) designs represented by their independent defining words
D 1 : I = a 1 a 2 b 2 b 3 b 4 = a 1 a 2 b 2 b 5 = a 1 b 1 b 2 b 6 = a 2 b 1 b 2 b 7 a n d D 2 : I = a 1 b 1 b 2 b 4 = a 1 b 2 b 3 b 5 = a 1 b 2 b 3 b 6 = b 1 b 2 b 3 b 7 ,
respectively. With some calculations we obtain that
1 ( s ) # C ( w ) ( 0 ) ( D 1 ) = 7 , 1 # C 2 ( D 1 ) = 9 , 2 # C 2 ( D 1 ) = ( 15 , 0 , 21 ) , 2 ( s ) # C ( w ) ( 0 ) ( D 1 ) = 33 a n d 1 ( s ) # C ( w ) ( 0 ) ( D 2 ) = 7 , 1 # C 2 ( D 2 ) = 9 , 2 # C 2 ( D 2 ) = ( 8 , 0 2 , 28 ) , 2 ( s ) # C ( w ) ( 0 ) ( D 2 ) = 35 .
Under criterion (1), D 1 is better than D 2 due to the following reasons. Note that 2 # C 2 ( 0 ) ( · ) is the first component, in (1), such that 2 # C 2 ( 0 ) ( D 1 ) 2 # C 2 ( 0 ) ( D 2 ) and 2 # C 2 ( 0 ) ( D 1 ) = 15 > 2 # C 2 ( 0 ) ( D 2 ) = 8 . Therefore, D 1 is better than D 2 under criterion (1).
In contrast, the FFSP design D 2 is better than D 1 under criterion (2). Note that criterion (2) prefers FFSP designs with resolution of at least IV, which have more SP 2FIs that are not aliased with any WP effect regardless of 2 # C 2 ( · ) . With this point in mind, since 1 ( s ) # C ( w ) ( 0 ) ( D 1 ) = 1 ( s ) # C ( w ) ( 0 ) ( D 2 ) , 1 # C 2 ( D 1 ) = 1 # C 2 ( D 2 ) and 2 ( s ) # C ( w ) ( 0 ) ( D 2 ) = 35 > 2 ( s ) # C ( w ) ( 0 ) ( D 1 ) = 33 , then design D 2 is better than D 1 under criterion (2).
As for criterion (3), it is clear that, if an FFSP design is better than another under criterion (1), then it is always the case when they are compared under criterion (3), noting that criterion (3) concerns one more component 2 ( s ) # C ( w ) ( 0 ) ( · ) apart from the three common components 1 ( s ) # C ( w ) ( 0 ) ( · ) = n 2 , 1 # C 2 ( · ) and 2 # C 2 ( · ) shared by (1) and (3). Therefore, design D 1 is better than D 2 under criterion (3). To show how to identify a better design under criterion (3), we consider two more examples represented by their independent defining words:
D 3 : I = a 1 b 1 b 2 b 3 b 4 = a 2 b 1 b 2 b 3 b 5 = a 3 b 1 b 2 b 3 b 6 = a 1 a 2 a 3 b 1 b 2 b 3 b 7 = a 1 a 2 b 1 b 8 = a 1 a 3 b 1 b 9 = a 2 a 3 b 1 b 10 = a 1 a 2 b 2 b 11 = a 1 a 3 b 2 b 12 = a 2 a 3 b 2 b 12 = a 1 a 2 b 3 b 14 = a 1 a 3 b 3 b 15 = a 2 a 3 b 3 b 16 = a 1 a 2 a 3 a 4 , D 4 : I = a 1 b 1 b 2 b 3 b 4 = a 1 a 2 b 1 b 2 b 5 = a 1 a 3 b 1 b 2 b 6 = a 1 a 2 a 3 b 1 b 2 b 3 b 7 = a 1 a 2 b 3 b 8 = a 1 a 3 b 3 b 9 = a 1 a 2 a 3 a 4 = a 2 b 1 b 3 b 10 = a 3 b 1 b 3 b 11 = a 2 a 3 b 1 b 12 = a 2 b 2 b 3 b 13 = a 3 b 2 b 3 b 14 = a 2 a 3 b 2 b 15 = a 2 a 3 b 3 b 16 ,
where D 3 and D 4 are two 2 ( 4 + 16 ) ( 1 + 13 ) FFSP designs, respectively. With some calculations, we obtain that
1 ( s ) # C ( w ) ( 0 ) ( D 3 ) = 16 , 1 # C 2 ( D 3 ) = 20 , 2 # C 2 ( D 3 ) = ( 0 3 , 160 , 0 5 , 30 ) , 2 ( s ) # C ( w ) ( 0 ) ( D 3 ) = 160 a n d 1 ( s ) # C ( w ) ( 0 ) ( D 4 ) = 16 , 1 # C 2 ( D 4 ) = 20 , 2 # C 2 ( D 4 ) = ( 0 3 , 160 , 0 5 , 30 ) , 2 ( s ) # C ( w ) ( 0 ) ( D 4 ) = 171 .
Although D 3 and D 4 have equal performance under criterion (1) due to that 1 ( s ) # C ( w ) ( 0 ) ( D 3 ) = 1 ( s ) # C ( w ) ( 0 ) ( D 4 ) , 1 # C 2 ( D 3 ) = 1 # C 2 ( D 4 ) and 2 # C 2 ( D 3 ) = 2 # C 2 ( D 4 ) , design D 4 is better than D 3 under criterion (3) as 2 ( s ) # C ( w ) ( 0 ) ( D 4 ) = 171 > 2 ( s ) # C ( w ) ( 0 ) ( D 3 ) = 160 .
The study of this paper is substantially different from the Refs. [7,8,9,10,11,18,19,26]. More specifically, Ref. [7] considered the regular symmetrical or mixed-level FFSP designs under the minimum secondary aberration criterion, which concerns only the number of SP-factor interactions in the WP alias sets; Ref. [8] studied the matrix presentation for FFSP designs at s levels as well as the maximum resolution and minimum aberration properties for such FFSP designs, where s is a prime number; Ref. [9] proposed generalized minimum aberration criteria for two-level orthogonal FFSP designs in five different design scenarios and tabulated a catalog of optimal 12-, 16-, 20-, and 24-run FFSP designs under their generalized minimum aberration criteria by computer algorithm; Refs. [10,11] both considered construction of FFSP designs under the WP-minimum aberration criterion, which assumes that the whole plot factor are more important. The criteria considered in our paper is different from those in [7,8,9,10,11]. These differences lead to that, for two-level regular FFSP designs, the optimal ones under the criteria considered in [7,8,9,10,11] may not be optimal under criteria (1), (2), and (3), and vice versa. Refs. [18,19] proposed some sufficient and necessary conditions for the asymmetrical split-plot designs to contain various types of clear effects, while our work considers developing theoretical construction methods of regular two-level FFSP designs under the optimality criteria (1), (2), and (3). Ref. [26] mainly focused on the regular two-level FFSP designs with replicated settings of the level combinations for WP factors, while the level combinations for the regular two-level FFSP design in our work are not replicated.
Due to the complex structure of FFSP designs, although we provide a series of theoretical construction methods for optimal FFSP designs under criteria (1), (2), and (3), there are still many optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) FFSP designs which cannot be constructed by our methods. For example, the theoretical construction methods for optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) FFSP designs, under criteria (1), (2), and (3), which satisfy N 4 + 1 n 5 N 16 are not covered in this paper. This is a future research direction worthy of study.

4. Conclusions

The 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs enjoy a wide application when performing a 2 n m design in a completely random order is impractical. A large body of work on choosing 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs under the MA criterion and clear effect criterion was proposed. The GMC-FFSP criterion is a widely used criterion for assessing 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs. This criterion advocates the FFSP designs with more effects at lower order confounding. The FFSP designs chosen under the GMC-FFSP criterion are preferable when we have prior information on the importance ordering of some effects. However, the theoretical construction methods of optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs under the GMC-FFSP criterion have not been studied yet.
This paper investigates theoretical construction methods of GMC-FFSP 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs. In addition, from the angle that the more there are lower order effects of interest at lower order confounding, the better the 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs, we propose optimality criteria for two kinds of design scenarios stated in the Introduction section. Some optimal 2 ( n 1 + n 2 ) ( m 1 + m 2 ) designs for these two kinds of design scenarios are also theoretically constructed under the newly proposed optimality criteria. In the supplementary material, the R code for the proposed designs is provided.

Author Contributions

Conceptualization, B.H. and Y.Z.; methodology, B.H. and Y.Z.; validation, B.H. and Y.Z.; writing—original draft preparation, B.H.; writing—review and editing, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 12171277 and 11801331).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
FFfractional factorial
MAminimum aberration
GMCgeneral minimum lower order confounding
FFSPfractional factorial split-plot
WPwhole plot
SPsubplot
2FItwo-factor interaction

References

  1. Box, G.; Hunter, W.; Hunter, S. Statistics for Experimenters: Design, Innovation, and Discovery, 2nd ed.; Wiley-Interscience: New York, NY, USA, 2005. [Google Scholar]
  2. Huang, P.; Chen, D.; Voelkel, J.O. Minimum-aberration two-level split-plot designs. Technometrics 1998, 40, 314–326. [Google Scholar] [CrossRef]
  3. Box, G.E.P.; Hunter, J.S. The 2k-p fractional factorial designs part I, part II. Technometrics 1961, 3, 311–351 & 449–458. [Google Scholar]
  4. Bingham, D.; Sitter, R.R. Some theoretical results for fractional factorial split-plot designs. Ann. Statist. 1999, 27, 1240–1255. [Google Scholar] [CrossRef]
  5. Bingham, D.; Sitter, R.R. Minimum aberration two-level fractional factorial split-plot designs. Technometrics 1999, 41, 62–70. [Google Scholar] [CrossRef]
  6. Mukerjee, R.; Fang, K.T. Fractional factorial split-plot designs with minimum aberration and maximum estimation capacity. Statist. Sinica 2002, 12, 885–903. [Google Scholar]
  7. Ai, M.Y.; Zhang, R.C. Minimum secondary aberration fractional factorial split-plot designs in terms of consulting designs. Sci. China Ser. A 2006, 49, 494–512. [Google Scholar] [CrossRef]
  8. Yang, J.F.; Liu, M.Q.; Zhang, R.C. Some results on fractional factorial split-plot designs with multi-level factors. Comm. Statist. 2009, 38, 3623–3633. [Google Scholar] [CrossRef]
  9. Tichon, J.G.; Li, W.; Mcleod, R.G. Generalized minimum aberration two-level split-plot designs. J. Statist. Plann. Inference 2012, 142, 1407–1414. [Google Scholar] [CrossRef]
  10. Wang, C.C.; Zhao, Q.Q.; Zhao, S.L. Optimal fractional factorial split-plot designs when the whole plot factors are important. J. Statist. Plann. Inference 2019, 199, 1–13. [Google Scholar] [CrossRef]
  11. Zhao, Q.Q.; Zhao, S.L. Constructing minimum aberration split-plot designs via complementary sets when the whole plot factors are important. J. Statist. Plann. Inference 2020, 209, 123–143. [Google Scholar] [CrossRef]
  12. Wu, C.F.J.; Hamada, M.S. Experiments: Planning, Analysis and Optimization, 2nd ed.; John and Wiley & Sons, Inc.: Hoboken, NJ, USA, 2009. [Google Scholar]
  13. Yang, J.F.; Li, P.F.; Liu, M.Q.; Zhang, R.C. 2(n1+n2)-(k1+k2) fractional factorial split-plot designs containing clear effects. J. Statist. Plann. Inference 2006, 136, 4450–4458. [Google Scholar] [CrossRef]
  14. Zi, X.M.; Zhang, R.C.; Liu, M.Q. Bounds on the maximum numbers of clear two-factor interactions for2(n1+n2)-(k1+k2) fractional factorial split-plot designs. Sci. China Ser. A 2006, 49, 1816–1829. [Google Scholar] [CrossRef]
  15. Zhao, S.L.; Chen, X.F. Mixed two- and four-level fractional factorial split-plot designs with clear effects. J. Statist. Plann. Inference 2012, 142, 1789–1793. [Google Scholar] [CrossRef]
  16. Zhao, S.L.; Chen, X.F. Mixed-level fractional factorial split-plot designs containing clear effects. Metrika 2012, 75, 953–962. [Google Scholar] [CrossRef]
  17. Yuan, Y.; Zhao, S.L. Mixed two- and eight-level fractional factorial split-plot designs containing clear effects. Acta Math. Appl. Sin. 2016, 32, 995–1004. [Google Scholar] [CrossRef]
  18. Han, X.X.; Liu, M.Q.; Yang, J.F.; Zhao, S.L. Mixed 2- and 2r-level fractional factorial split-plot designs with clear effects. J. Statist. Plann. Inference 2020, 204, 206–216. [Google Scholar] [CrossRef]
  19. Han, X.X.; Chen, J.B.; Liu, M.Q.; Zhao, S.L. Asymmetrical split-plot designs with clear effects. Metrika 2020, 83, 779–798. [Google Scholar] [CrossRef]
  20. Wei, J.L.; Yang, J.F.; Li, P.; Zhang, R.C. Split-plot designs with general minimum lower-order confounding. Sci. China Math. 2010, 53, 939–952. [Google Scholar] [CrossRef]
  21. Zhang, R.C.; Li, P.; Zhao, S.L.; Ai, M.Y. A general minimum lower-order confounding criterion for two-level regular designs. Statist. Sinica 2008, 18, 1689–1705. [Google Scholar]
  22. Bruen, A.A.; Haddad, L.; Wehlau, D.L. Binary codes and caps. J. Comb. Des. 1998, 6, 275–284. [Google Scholar] [CrossRef]
  23. Butler, N.A. Results for two-level fractional factorial designs of resolution IV or more. J. Statist. Plann. Inference 2007, 137, 317–323. [Google Scholar] [CrossRef]
  24. Li, P.F.; Zhao, S.L.; Zhang, R.C. A theory on constructing 2n-m designs with general minimum lower order confounding. Statist. Sinica 2011, 21, 1571–1589. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, Y.N.; Zhao, S.L.; Liu, M.Q. On construction of optimal two-level designs with multi block variables. J. Syst. Sci. Complex. 2018, 31, 773–786. [Google Scholar] [CrossRef]
  26. Zhao, Y.N. Split-plot designs with few whole plot factors containing clear effects. Fractal Fract. 2022, 6, 1–14. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Han, B.; Zhao, Y. Construction of Optimal Split-Plot Designs for Various Design Scenarios. Fractal Fract. 2022, 6, 608. https://doi.org/10.3390/fractalfract6100608

AMA Style

Han B, Zhao Y. Construction of Optimal Split-Plot Designs for Various Design Scenarios. Fractal and Fractional. 2022; 6(10):608. https://doi.org/10.3390/fractalfract6100608

Chicago/Turabian Style

Han, Beichen, and Yuna Zhao. 2022. "Construction of Optimal Split-Plot Designs for Various Design Scenarios" Fractal and Fractional 6, no. 10: 608. https://doi.org/10.3390/fractalfract6100608

APA Style

Han, B., & Zhao, Y. (2022). Construction of Optimal Split-Plot Designs for Various Design Scenarios. Fractal and Fractional, 6(10), 608. https://doi.org/10.3390/fractalfract6100608

Article Metrics

Back to TopTop