1. Introduction
1.1. Research Background
Experimental design refers to the scientific and rational arrangement of experiments after clarifying the factors to be examined and the research objectives, in order to achieve the best experimental results. This process is also called experimental design. The purpose of experimental design is to explore the relationships between variables through systematic methods, so as to optimize processes and improve product quality.
The development of experimental design can be traced back to the early 20th century. British statistician R.A. Fisher published the first example of experimental design in collaboration with W.A. Mackenzie in 1923, and proposed the basic ideas of experimental design in 1926. In 1935, Fisher [1] published his famous book “The Design of Experiments”, in which he proposed three principles that experimental design should follow: randomization, local control, and replication. These principles aim to reduce the impact of accidental factors, so that experimental data has an appropriate mathematical model, which can be used for data analysis with the method of variance analysis.
In China, research on experimental design began in the 1950s, and there were new insights in the viewpoints, theories, and methods of orthogonal experimental design. The famous mathematician Professor Hua Luogeng actively advocated and popularized the “optimization method” in China, thus popularizing the concept of experimental design. In 1978, mathematicians Wang Yuan [2] and Fang Kaitai [3] proposed the uniform design, which considers how to scatter the design points evenly within the experimental range, so as to obtain the most information with fewer experimental points.
In the field of experimental design, with the continuous in-depth of scientific research and industrial applications, traditional orthogonal parameterization design methods have gradually shown limitations in some specific scenarios. Especially in experimental situations where it is necessary to screen key factors from many factors and hope to minimize changes to the existing process, baseline design has emerged, bringing new ideas and methods to the field of experimental design.
Baseline design is an experimental design in which each factor has a specified baseline level (default or preferred level), measuring the impact of changes in one or more factor levels on the response while other factors remain at the baseline level. This design parameterization method is different from orthogonal parameterization, which measures the impact of changes in one or more factor levels on the response, averaged over all possible level combinations of all other factors. The goal of baseline design is to identify designs with the minimum
-aberration.
1.2. Research Content
Based on the results of Miller and Tang [4], this paper further develops the relationship between K-aberration and word length pattern. Use
to represent a regular two-level design with
factors and
runs, and the two levels in the design are represented by 0, 1 respectively. For any regular
design
, the set of all
column submatrices of
is denoted as
, where
. Let
represent the total number of rows in
where all elements are 1, where
is a subdesign composed of some columns of
. Under the baseline parameterization model, for a regular
design
,
represents the total deviation caused by all
th-order factor interactions on the estimation of main effects in design
. Mukerjee and Hunter [5] proved that,
where,
,
.
Given the number of rows and columns of the design, the two-level orthogonal design that minimizes the sequence
in lexicographic order is called the
-aberration optimal design.
Based on the
expression as the theoretical foundation, this paper has thoroughly explored the research process of the
expression under various circumstances, analyzed the relationship between
and the word length pattern for regular two-level designs with resolution
, listed all possible defining word scenarios, and provided the expression results for
, offering corresponding theoretical results for the design of such experiments. That is, by sequentially minimizing the
values up to
, the design method can be obtained more quickly.
2. The Relationship between K-Criterion and
Word Length Pattern
In the study of experimental design, we find that
and
are closely related to the subarrays
or
of the regular
design. According to the definition of
in the formula, when the subarray
composed of any
columns in design
does not contain defining words, when calculating
, the contribution of
to
is
. When the subarray
composed of any
columns in design
contains defining words, that is, several columns in
form defining words, the situation becomes complicated. The existence of defining words means that some column combinations in
have special structures, which will affect the contribution of
to
. In addition, similar situations need to be analyzed when calculating
. Therefore, investigating whether the subarray
contains defining words is crucial for simplifying the calculation of
. For the convenience of the subsequent analysis of the
formula, the following lemmas are given to provide theoretical explanations for all possible cases of defining words.
2.1. Three Basic Lemmas
Lemma 1. Suppose
is a regular
design with resolution
, let
, then
(i) When
,
contains at most three independent defining words;
(ii) When
,
contains at most two independent defining words;
(iii) When
,
contains at most one defining word,
where
.
Proof:
(i) When
,
,
, then the number of independent defining words contained in
is as follows:
(a1)
does not contain independent defining words;
(a2)
contains one independent defining word, the length of which may be 3 or 4 or 5 or 6;
(a3)
contains two independent defining words, their lengths may be 3 and 3 or 3 and 4;
(a4)
contains three independent defining words, all of whose lengths are 3.
Situations (a1) and (a2) are obvious. In situation (a3), when the lengths of the two independent defining words are 3 and 3, there are two cases for their composition, that is, they contain a common column (for example,
) and do not contain a common column (for example,
); when the lengths of the two independent defining words are 3 and 4, their composition is that they contain a common column (for example,
). In situation (a4), if the length of one of the independent defining words exceeds 3, then it is impossible to have two other independent defining words with lengths ≤3. Therefore, when
contains three independent defining words, the length of the independent defining words can only be 3. Since 6 factors form 3 independent defining words, there must be a common column between every two independent defining words. Let
be the common column factors between the three independent defining words,
represents an unknown factor, then the relationship between the three independent defining words can be obtained as:
, such arrangement can supplement the
place with the other three factors. Since these three independent defining words will generate a non-independent defining word with a length of 3, so when
,
cannot contain four independent defining words.
(ii) When
,
,
, then the number of independent defining words contained in
is as follows:
(a1)
does not contain independent defining words;
(a2)
contains one independent defining word, the length of which may be 4 or 5 or 6 or 7;
(a3)
contains two independent defining words, their lengths may be 4 and 4 or 4 and 5.
Situations (a1) and (a2) are obvious. In situation (a3), the independent defining words with lengths 4 and 4 may have one common column or two common columns. The independent defining words with lengths 4 and 5 have two common columns. Consider whether
may contain three independent defining words. Suppose
contains three independent defining words, all of which have a length of 4, and two of the independent defining words have two common columns, denoted as
, and there is the following defining relationship, that is,
, which can generate a non-independent defining word
, so that the other independent defining word has at most two common columns with these three defining words, thus two additional factors need to be added, exceeding the range of seven factors. Therefore, when
,
contains at most two independent defining words. The cases of
and
can be proved in the same way.
(iii) When
,
cannot have two independent defining words. Selecting two independent defining words with a length of 7 from 10 factors will inevitably obtain a non-independent defining word with a length
, which does not conform to the condition of resolution 7, so when the resolution is 7,
contains at most one independent defining word. The cases of
can be proved in the same way.
Lemma 2. Consider a regular matrix
, which is an orthogonal matrix of strength 2. If this design is used to create a baseline design, then the values of
and
in the expression satisfy the following conditions:
(i) When
,
,
(ii) When
,
,
(iii) When
and
is odd
(iv) When
and
is even
Proof:
(i) For any regular two-level design, let
be the submatrix composed of any
columns of the design matrix. Let
be the resolution of the design, then when
,
. Substituting it into the expressions of
and
, we can get
,
, proved. (ii) When
,
still satisfies the above condition, that is,
, so
.
But for
,
that is,
, then it is necessary to consider whether this
columns contain defining words. Let
, where
represents the value of the sum of the rows of the
matrix modulo 2 and then added together. When
is a non-defining word matrix,
,
. When
is a defining word matrix,
or
,
. Therefore, when this
columns do not contain defining words,
,
,
, it can be concluded that
when this
columns contain defining words,
. Among them, this submatrix contains
designs and
half of the
designs, that is, 0
designs and
half of the
designs. Whether the
half of the
designs contain all 1 rows needs to consider the parity of
and the
value of the defining words. The following four situations are divided:
(a1) When
is odd,
, that is, there is no all 1 row,
,
,
;
(a2) When
is odd,
, that is, there is an all 1 row,
,
,
;
(a3) When
is even,
, that is, there is an all 1 row,
,
,
;
(a4) When
is even,
, that is, there is no all 1 row,
,
,
.
The
obtained in the four situations are all
, so we can get:
(iii) The calculation of
needs to consider three situations:
Table 1. Calculation situations of
.
Defining word situation |
|
|
This
columns do not contain any defining words |
|
|
This
columns form defining words and
|
0 |
|
This
columns form defining words and
|
|
|
: The number of W in each situation.
According to Table 1, the calculation result of
can be obtained.
The calculation of
needs to consider four situations:
(a1)
is a defining word with a length of
and
;
(a2)
is a defining word with a length of
and
;
(a3)
contains a defining word with a length of
and
;
(a4)
contains a defining word with a length of
and
.
Table 2. Calculation situations of
.
Situation |
|
:
|
|
(a1) |
|
|
|
(a2) |
0 |
|
|
(a3) |
0 |
and
|
|
(a4) |
|
and
|
|
: The number of each
value corresponding to
.
: The number of
in (a1)-(a4).
According to Table 2, the calculation result of
can be obtained.
(iii) is proved, and the result of (iv) can be obtained in the same way.
Lemma 3. The
value of the defining word with a length of 5 obtained by the combination of independent defining words with lengths 3 and 4 containing a common column is determined. Let the
values of the above defining words with lengths 3 and 4 be
and
respectively, and the
value of the defining word with a length of 5 obtained by the combination be
, the following results can be obtained:
(i) If
and
are both
, then
is
;
(ii) If
and
are both
, then
is
;
(iii) If one of
and
is
and the other is
, then
is
.
Proof:
Let
, the defining relationship is
, let
,
, where the three factors in
form a defining word with a length of 3,
, the four factors in
form a defining word with a length of 4,
, their common column is
.
indicates that there are an even number of 1 s in each row of the three columns
,
indicates that there are an even number of 1 s in each row of the three columns
. When the column
is 1, the corresponding two columns
have an odd number of 1 s, and
have an even number of 1 s. When the column
is 0, the corresponding two columns
have an even number of 1 s, and
have an even number of 1 s, that is,
have an even number of 1 s, so when
and
are both
, then
is
. Thus, (i) is proved, and (ii) and (iii) can be proved in the same way.
2.2. The Relationship between
and WLP When the Resolution
This section mainly studies the relationship between
and WLP in regular
designs with resolution
. By discussing the defining words contained in
and
, the analysis process of
and
is given, and finally the relationship between
and WLP in regular
designs with resolution
is obtained.
Theorem 1. Suppose
is a regular
design with resolution
, then
where
Proof: First, calculate the expression of
, let
, then there are ten possible situations for the columns in
:
(a1)
contains two independent defining words with a length of 3 and
;
(a2)
contains two independent defining words with a length of 3 and
;
(a3)
contains two independent defining words with a length of 3, one of which has
and the other has
;
(a4)
only contains one defining word with a length of 3 and
;
(a5)
only contains one defining word with a length of 3 and
;
(a6)
only contains one defining word with a length of 4 and
;
(a7)
only contains one defining word with a length of 4 and
;
(a8)
only contains one defining word with a length of 5 and
;
(a9)
only contains one defining word with a length of 5 and
;
(a10)
contains 5 independent columns.
The reason why (a1), (a2), (a3) appear independently from other situations is that two defining words can combine to form other situations, and to avoid repeated calculations, they are listed separately. The defining words in (a1) and (a2) can combine to form an independent defining word with a length of 4 and
. The defining words in (a3) can combine to form an independent defining word with a length of 4 and
.
Below are the proofs of the values of
and the number of each
in situations (a1)-(a10).
For (a1), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0, the number is
.
For (a2), two independent defining words with a length of 3 and
can combine to form an independent defining word with a length of 4 and
. The cross column of the two independent defining words is affected by the other two columns of the defining words. When the other two columns of the defining words are all 1, the cross column must be 1, so there is an all-one row in this
,
is
, the number is
.
For (a3) and (a4), the situation is the same as (a1),
, the numbers are
and
respectively. In situation (a4), each
has and only has one defining word with a length of 3 and
. The number of such defining words is
,
, that is, selecting 5 factors from
factors that meet the conditions of (a4), which is
. To avoid the appearance of repeated situations, it is necessary to remove the situations where
in (a1) and (a2). Since (a1) contains two
defining words in
, the final number in situation (a4) is
.
For (a5), each
contains a defining word with a length of 3 and
, so the number of all-one rows in the defining word, combined with the other two independent columns, is
. In this situation, each
contains 6
, two of which are combinations of the defining word and one independent column, their
is
, and the other four
are independent columns, their
is
. The number is based on
and needs to remove the non-independent defining words with a length of 5 and
generated in (a9), (a10), (a15), (a16), (a18). Thus, the number is
.
For (a6), each
contains a defining word with a length of 4 and
, so the number of all-one rows in the defining word is
, and the number of all-one rows combined with another independent column is
, that is,
is
, the number is
.
For (a7), since there exists a defining word with a length of 4 and
, there cannot be an all-one row in
,
is 0, the number is
.
For (a8), since there exists a defining word with a length of 5 and
, there cannot be an all-one row in
,
is 0, the number is
.
For (a9),
contains an all-one row, five factors form a defining word, so
is
, the number is
.
For (a10), since
contains 5 independent columns,
is
. Selecting 5 factors from
factors has
possibilities, removing all previous situations, the number is
, denoted as
.
Summarizing the above situations, we obtain Table 3:
Table 3. Calculation of
in Theorem 1.
Situation |
|
|
(a1) |
0 |
|
(a2) |
|
|
(a3) |
0 |
|
(a4) |
0 |
|
(a5) |
|
|
(a6) |
|
|
(a7) |
0 |
|
(a8) |
0 |
|
(a9) |
|
|
(a10) |
|
|
: The number of
in each situation.
Thus, the expression of
is obtained.
Next, calculate the expression of
. When
,
has the following twenty situations:
(a1)
contains only one defining word with a length of 6 and
;
(a2)
contains only one defining word with a length of 6 and
;
(a3)
contains only one defining word with a length of 5 and
;
(a4)
contains only one defining word with a length of 5 and
;
(a5)
contains only one defining word with a length of 4 and
;
(a6)
contains only one defining word with a length of 4 and
;
(a7)
contains only one defining word with a length of 3 and
;
(a8)
contains only one defining word with a length of 3 and
;
(a9)
contains two independent defining words with a length of 3 and
, sharing one common column;
(a10)
contains two independent defining words with a length of 3 and
, sharing one common column;
(a11)
contains two independent defining words with a length of 3, one with
and the other with
, sharing one common column;
(a12)
contains two independent defining words with a length of 3 and
, without sharing any common column;
(a13)
contains two independent defining words with a length of 3 and
, without sharing any common column;
(a14)
contains two independent defining words with a length of 3, one with
and the other with
, without sharing any common column;
(a15)
contains three independent defining words with a length of 3, each pair sharing one common column, divided into the following four cases:
(1)
does not contain any defining word with
;
(2)
contains only one defining word with
;
(3)
contains only two defining words with
;
(4)
contains only three defining words with
.
(a16)
contains one independent defining word with a length of 3 and
and one independent defining word with a length of 4 and
, sharing one common column;
(a17)
contains one independent defining word with a length of 3 and
and one independent defining word with a length of 4 and
, sharing one common column;
(a18)
contains one independent defining word with a length of 3 and
and one independent defining word with a length of 4 and
, sharing one common column;
(a19)
contains one independent defining word with a length of 3 and
and one independent defining word with a length of 4 and
, sharing one common column;
(a20)
contains 6 independent columns.
Next, using Lemma 2’s calculation method, the values of
and
and the number of each
in situations (a1)-(a20) are proved.
For (a1), each
contains a defining word with a length of 6 and
, so the number of all-one rows in the defining word is
. In this situation, each
contains 6
which are all independent columns, that is,
is
. The number is based on
and needs to remove the non-independent defining words with a length of 6 and
generated by two independent defining words with a length of 3 in (a12) and (a14). Thus, the number is
;
For (a2), since there exists a defining word with a length of 6 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
which are all independent columns, that is,
is
. The number is based on
and needs to remove the non-independent defining words with a length of 6 and
generated by two independent defining words with a length of 3 in (a13). Thus, the number is
;
For (a3), since there exists a defining word with a length of 5 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, one of which is a defining word, its
is 0, and the other five
are independent columns, their
is
. The number is based on
and needs to remove the non-independent defining words with a length of 5 and
generated by two independent defining words with lengths of 3 and 4 in (a16) and (a19). Thus, the number is
;
For (a4), each
contains a defining word with a length of 5 and
, so the number of all-one rows in the defining word, combined with another independent column, is
. In this situation, each
contains 6
, one of which is a defining word, its
is
, and the other five
are independent columns, their
is
. The number is based on
and needs to remove the non-independent defining words with a length of 5 and
generated by two independent defining words with lengths of 3 and 4 in (a16) and (a19). Thus, the number is
;
For (a5), each
contains a defining word with a length of 4 and
, so the number of all-one rows in the defining word, combined with the other two independent columns, is
. In this situation, each
contains 6
, two of which are combinations of defining words with a length of 4 and one independent column, their
is
, and the other four
are independent columns, their
is
. The number is based on
and needs to remove the non-independent defining words with a length of 4 and
generated in (a9), (a10), (a15), (a16), (a18). Thus, the number is
, denoted as
;
For (a6), since there exists a defining word with a length of 4 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, two of which are combinations of defining words with a length of 4 and one independent column, their
is 0, and the other four
are independent columns, their
is
. The number is similar to the analysis in (a5), which is
, denoted as
;
For (a7), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, three of which are combinations of defining words with a length of 3 and two independent columns, their
is 0, and the other three
are independent columns, their
is
. The number needs to remove (a9), (a11), (a12), (a14), (a15), (a16), (a17), which is
, denoted as
;
For (a8), each
contains a defining word with a length of 3 and
, so the number of all-one rows in the defining word, combined with the other three independent columns, is
. In this situation, each
contains 6
, three of which are combinations of defining words with a length of 3 and two independent columns, their
is
, and the other three
are independent columns, their
is
. The number is similar to the situation of (a8), which is
, denoted as
;
For (a9), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, their
are all 0, and there is no all-one column in the case without common columns. The number is
;
For (a10),
contains an all-one row,
is
. In this situation, each
contains 6
, the
in the case without common columns is
, and the
in the other five cases is
. The number is
;
For (a11), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, two of which have
as
, and four have
as 0. The number is
;
For (a12), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
which are all 0. The number is
;
For (a13),
contains an all-one row,
is
. In this situation, each
contains 6
which are all
. The number is
;
For (a14), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, three of which are 0, and the other three are
. The number is
;
For (a15), three independent defining words with a length of 3 will generate a non-independent defining word with a length of 3, the
value of which is determined by the
values of the three independent defining words. As long as there exists a defining word with a length of 3 and
, there cannot be an all-one row in
, so the
values of situations (1), (2), and (3) are all 0, and the
value of situation (4) is
. When there are two or more defining words with
, each
contains 6
which are all 0, so the
of situations (1) and (2) are all 0. In situation (3), each
contains 6
, one of which is
, and the other five are 0. In situation (4), each
contains 6
whose
are all
. Their numbers are denoted as
,
,
,
respectively;
For (a16), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, three of which are 0, and the other three are
. The number is
;
For (a17), since there exists a defining word with a length of 3 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
whose
are all
. The number is
;
For (a18), according to Lemma 3, a defining word with a length of 3 and
and a defining word with a length of 4 and
will generate a defining word with a length of 5 and
,
is
. In this situation, each
contains 6
, five of which are 0, and the other one is
. The number is
;
For (a19), since there exists a defining word with a length of 4 and
, there cannot be an all-one row in
,
is 0. In this situation, each
contains 6
, three of which are 0, and the other three are
. The number is
;
Table 4. Calculation of
in Theorem 1.
Situation |
|
:
|
|
(a1) |
|
|
|
(a2) |
0 |
|
|
(a3) |
0 |
and
|
|
(a4) |
|
and
|
|
(a5) |
|
and
|
|
(a6) |
0 |
and
|
|
(a7) |
0 |
and
|
|
(a8) |
|
and
|
|
(a9) |
0 |
|
|
(a10) |
|
and
|
|
(a11) |
0 |
and
|
|
(a12) |
0 |
|
|
(a13) |
|
|
|
(a14) |
0 |
and
|
|
(a15-1) |
0 |
|
|
(a15-2) |
0 |
|
|
(a15-3) |
0 |
and
|
|
(a15-4) |
|
|
|
(a16) |
0 |
and
|
|
(a17) |
0 |
and
|
|
(a18) |
|
|
|
(a19) |
0 |
and
|
|
(a20) |
|
|
|
: The number of
in each situation.
For (a20), since
contains 6 independent columns,
is
. In this situation, each
contains 6
whose
are all
. Selecting 6 factors from
factors has
possibilities, removing all previous situations, the number is denoted as
.
Table 4 summarizes the above discussions.
According to Table 4, the calculation result of
can be obtained.
The above completes the proof of Theorem 1.
3. Applications
To verify the effectiveness of the
expression in orthogonal two-level fractional factorial designs, a resolution 3 experimental scheme can be designed to assess its computational advantages and practical application value through simulation studies.
Suppose we are studying a chemical reaction process and need to consider 9 factors (
) that affect the reaction efficiency. Due to the limited number of experimental runs, we choose a fractional factorial design, specifically a
fractional factorial design, which means 9 factors and
experimental runs. According to the design catalog proposed by Xu [6], we select one type of experiment with a word length pattern of (4, 14, 8, 0, 4, 1, 0), and its defining relations are as follows:
Based on the defining relations, the
expression is calculated, and the
value is minimized in sequence by using level permutation to seek the optimal design. According to Lemma 2, we can obtain
, where
is a constant.
, by which the
values of the defining words with length 3 can be set to
to minimize
. Assuming that
has reached the minimum value, we next consider
. Using C to represent the constant, the expression for
is obtained as
, by which the
values of the defining words with length 4 can be set to
to minimize
. After verification, it can be concluded that this design is the optimal design under the
-aberration criterion.