Next Article in Journal
TOPSIS Method Based on Hamacher Choquet-Integral Aggregation Operators for Atanassov-Intuitionistic Fuzzy Sets and Their Applications in Decision-Making
Next Article in Special Issue
Discrete Single-Factor Extension of the Exponential Distribution: Features and Modeling
Previous Article in Journal
Relationship Assessment of COVID-19, Air Pollution, and Copper Demand from the Perspective of Copper Price
Previous Article in Special Issue
Representative Points Based on Power Exponential Kernel Discrepancy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projection Uniformity of Asymmetric Fractional Factorials

1
School of Statistics and Mathematics, Zhongnan University of Economics and Law, Wuhan 430073, China
2
School of Mathematics and Statistics, Jishou University, Jishou 416000, China
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(12), 716; https://doi.org/10.3390/axioms11120716
Submission received: 4 November 2022 / Revised: 29 November 2022 / Accepted: 6 December 2022 / Published: 10 December 2022
(This article belongs to the Special Issue Computational Statistics & Data Analysis)

Abstract

:
The objective of this paper is to study the issue of the projection uniformity of asymmetric fractional factorials. On the basis of level permutation and mixture discrepancy, the average projection mixture discrepancy to measure the uniformity for low-dimensional projection designs is defined, the uniformity pattern and minimum projection uniformity criterion are presented for evaluating and comparing any asymmetric factorials. Moreover, lower bounds to uniformity pattern have been obtained, and some illustrative examples are also provided.

1. Introduction

Many criteria were proposed for comparing U-type designs, but none of these criteria can directly distinguish non-isomorphic saturated designs. A special criterion can measure all these subdesigns, and the related values are called its projection pattern. We can use the distribution or the vector of these projection values as a tool to distinguish the underlying designs. Ref. [1] firstly defined the projection discrepancy pattern and proposed the minimum projection uniformity (MPU) criterion, which is equivalent to generalized minimum aberration criterion (GMA [2]). Ref. [3] studied the projection discrepancies of two-level fractional factorials in terms of the centered L 2 -discrepancy (CD [4]). Subsequently, ref. [5] discussed the relationships among criteria of MPU proposed in [1] and minimum generalized aberration [6]. Following this projection discrepancy, [7] studied the projection properties of two-level factorials in view of geometry and proposed the uniformity pattern and MPU criterion to assess and compare two-level factorials. The relations between MPU and minimum aberration, and GMA and orthogonality are clarified; this close relationship raises the hope of improving the connection between uniform design theory and factorial design theory.
Following the uniform pattern and MPU, projection uniformity of asymmetric design based on CD and wrap-around L 2 -discrepancy (WD [8]) has been studied, respectively. As a measure of uniformity, CD does not have fewer cursed dimensions and WD is not sensitive to a shift for one or more dimensions, Mixture discrepancy (MD [9]) retains the good properties of CD and WD and overcomes the shortcomings of both. Aided by the level permutation technique in [10,11], ref. [12] obtained the relationship between the mean of mixture discrepancies and the generalized word–length pattern for multi-level designs. Ref. [13] defined the MPU criterion for two- and three-level factorials under MD. Refs. [14,15] generalize the findings in [13] to q-level and mixed two- and three-level factorials, respectively. Moreover, ref. [16] proposed the uniform projection design that have the smallest average CD values of all two-dimensional projections and are shown to have good-filling properties over all sub-spaces in terms of the distance, uniformity, and orthogonality. Based on the findings of [16], many applications and studies on uniform projection designs have emerged [17,18,19,20,21,22].
While the work of [13,14,15] discussed the projection uniformity for two-level, three-level, q-level, and mixed two- and three-level designs under MD, respectively, the present paper aims at obtaining further results. We extend the findings in [13,14,15] to general asymmetrical factorials. First, the uniformity pattern and MPU criterion are proposed for selecting asymmetrical designs. Second, we build some analytic linkages between uniformity pattern, orthogonality, and generalized word–length pattern. Third, we integrate two lower bound methods in [23], which can be served as a benchmark for searching MPU designs. Finally, the results of [13,14,15] can be used as our special cases, and some numerical examples are provided to illustrate our theoretical results.
This paper is organized as follows: Section 2 describes some notations and basic concepts such as distance distribution and generalized word–length pattern, which are useful throughout in this paper. Section 3 defines the average projection mixture discrepancy and related uniformity pattern, presents a statistical justification of MPU criterion, and establishes a connection between MPU and GMA. Section 4 provides a lower bound of the uniformity pattern. Some illustrative examples to verify our theoretical results are presented in Section 5.

2. Notations and Preliminaries

Consider a class of U-type designs, denoted by U ( n ; q 1 s 1 × q 2 s 2 ) , of mixed q 1 - and q 2 -level factorials in n runs and s ( = s 1 + s 2 ) factors, where each factor of the first s 1 factors takes values from a set of { 0 , 1 , , q 1 1 } equally often and each factor of the last s 2 factors takes values from a set of { 0 , 1 , , q 2 1 } equally often. For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , a typical treatment combination (or run) of design d is defined by w = ( w ( 1 ) , w ( 2 ) ) , where, for i = 1 , 2 , w ( i ) = ( w 1 ( i ) , , w s i ( i ) ) , w j ( 1 ) { 0 , 1 , , q 1 1 } and w j ( 2 ) { 0 , 1 , , q 2 1 } . Denote d = ( d ( 1 ) , d ( 2 ) ) , where w ( i ) d ( i ) , i = 1 , 2 . If all the possible q 1 t 1 × q 2 t 2 level combinations corresponding to any t ( = t 1 + t 2 ) columns of design d appear equally often, 0 t 1 s 1 , 0 t 2 s 2 , design d is called to be an orthogonal array of strength t and denoted by O A ( n ; q 1 s 1 × q 2 s 2 , t ) .
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , its distance distribution is defined by
E j 1 j 2 ( d ) = 1 n { ( i , k ) : H i 1 k 1 ( 1 ) = j 1 , H i 2 k 2 ( 2 ) = j 2 } ,
where | u | is the cardinality of the set | u | , H i k t is the Hamming distance between two runs i and k of design d ( t ) , t = 1 , 2 , 0 j 1 s 1 , 0 j 2 s 2 .
The MacWilliams transforms of the { E j 1 j 2 ( d ) } of any design d U ( n ; q 1 s 1 × q 2 s 2 ) are defined as
E i 1 i 2 ( d ) = 1 n j 1 = 0 s 1 j 2 = 0 s 2 P i 1 ( j 1 ; s 1 , q 1 ) P i 2 ( j 2 ; s 2 , q 2 ) E j 1 j 2 ( d ) , i 1 = 0 , , s 1 , i 2 = 0 , , s 2 ,
where P i ( j ; s , q ) = r = 0 i ( 1 ) r ( q 1 ) i r j r s j i r is the Krawtchouk polynomial, m k = m ( m 1 ) · · · ( m k + 1 ) / k ! and m k = 0 for m < k .
Ref. [2] showed that the generalized word–length pattern is the MacWilliams transform of the distance distribution, that is,
A i ( d ) = i 1 + i 2 = i E i 1 i 2 ( d ) ,
where the vector ( A 1 ( d ) , , A s ( d ) ) is called the generalized word–length pattern. For any two designs d 1 and d 2 in U ( n ; q 1 s 1 × q 2 s 2 ) , d 1 is said to have less aberration than d 2 if there exists a positive integer t s , such that A t ( d 1 ) < A t ( d 2 ) and A i ( d 1 ) = A i ( d 2 ) for i = 1 , , t 1 . The design d 1 has generalized minimum aberration if there is no other design with less aberration than d 1 .
For any positive integer g s , defined C g = { ( g 1 , g 2 ) : g 1 = 0 , , s 1 , g 2 = 0 , , s 2 , g 1 + g 2 = g } , and for any ( g 1 , g 2 ) C g , let S g 1 g 2 be the set of all nonempty subsets of { 1 , , s } with the first g 1 elements from { 1 , 2 , , s 1 } and the next g 2 elements from { s 1 + 1 , , s 1 + s 2 } . For any g, 1 g s , let S g be the set of all nonempty subsets of { 1 , 2 , , s } with cardinality g, it is to be noted that S g = ( g 1 , g 2 ) C g S g 1 g 2 .
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , define the nonempty set u = u 1 u 2 = { u 11 , , u 1 g 1 } { u 21 , , u 2 g 2 } S g 1 g 2 and g = g 1 + g 2 , let d u be the corresponding projection design of d onto factors with indexes from u. A typical treatment combination of d u is represented as w u = ( w u ( 1 ) , w u ( 2 ) ) , where w u ( i ) = ( w u i 1 ( i ) , , w u i g i ( i ) ) , w u i g i ( i ) { 0 , 1 , , q i 1 } , i = 1 , 2 . Let H i k u be the Hamming distance between two runs i u and k u of the projection design d u , denote δ i k u = g H i k u as the coincide number between two runs i u and k u , where i u = ( i 1 u , i 2 u ) and k u = ( k 1 u , k 2 u ) .

3. Projection Uniformity of U ( n ; q 1 s 1 × q 2 s 2 )

For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , g ( = g 1 + g 2 ) s and u S g , let M D u ( d ) be the mixture discrepancy value of the corresponding projection design d u ; following [9], we can derive the below formula for M D u ( d ) ,
[ M D u ( d ) ] 2 = 7 12 g 2 n i = 1 n j u f 1 ( x i j ) + 1 n 2 i = 1 n k = 1 n j u f ( x i j , x k j ) ,
where f 1 ( x i j ) = 2 3 1 4 | x i j 1 2 | 1 4 | x i j 1 2 | 2 , f ( x i j , x k j ) = 7 8 1 4 | x i j 1 2 | 1 4 | x k j 1 2 | 3 4 | x i j x k j | + 1 2 | x i j x k j | 2 , i, k = 1 , , n .
When considering all q 1 ! × q 2 ! possible level permutations for every factor of d U ( n ; q 1 s 1 × q 2 s 2 ) , there are ( q 1 ! ) s 1 × ( q 2 ! ) s 2 combinatorially isomorphic designs of d that can be obtained, and denote the set of these designs as P ( d ) . Similarly, for any positive integer g ( = g 1 + g 2 ) s and u S g , we can obtain ( q 1 ! ) g 1 × ( q 2 ! ) g 2 combinatorially isomorphic designs of d u ; the corresponding set of these combinatorially isomorphic designs d u is denoted by P ( d u ) . The mean of projection mixture discrepancies of all the designs in P ( d u ) is denoted by A M D u ( d ) , that is,
A M D u ( d ) = 1 ( q 1 ! ) g 1 ( q 2 ! ) g 2 d u P ( d u ) [ M D u ( d ) ] 2 .
The following lemma, which can be proved similarly as [14,15], gives the expression for A M D u ( d ) .
Lemma 1. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , u S g and 1 g s ,
( i ) when both q 1 and q 2 are even,
A M D u ( d ) = 7 12 g 2 28 q 1 2 + 1 48 q 1 2 g 1 28 q 2 2 + 1 48 q 2 2 g 2 + 1 n 3 4 g 1 3 4 g 2 i 1 = 0 g 1 i 2 = 0 g 2 7 q 1 2 9 q 1 i 1 7 q 2 2 9 q 2 i 2 E i 1 i 2 ( d u ) ;
( i i ) when both q 1 and q 2 are odd,
A M D u ( d ) = 7 12 g 2 7 q 1 2 + 1 12 q 1 2 g 1 7 q 2 2 + 1 12 q 2 2 g 2 + 1 n 6 q 1 2 + 1 8 q 1 2 g 1 6 q 2 2 + 1 8 q 2 2 g 2 × i 1 = 0 g 1 i 2 = 0 g 2 14 q 1 2 4 q 1 + 3 18 q 1 2 + 3 i 1 14 q 2 2 4 q 2 + 3 18 q 2 2 + 3 i 2 E i 1 i 2 ( d u ) ;
( i i i ) when q 1 is even and q 2 is odd,
A M D u ( d ) = 7 12 g 2 28 q 1 2 + 1 48 q 1 2 g 1 7 q 2 2 + 1 12 q 2 2 g 2 + 1 n 3 4 g 1 6 q 2 2 + 1 8 q 2 2 g 2 × i 1 = 0 g 1 i 2 = 0 g 2 7 q 1 2 9 q 1 i 1 14 q 2 2 4 q 2 + 3 18 q 2 2 + 3 i 2 E i 1 i 2 ( d u ) .
We can obtain the following lemma when the design d is an orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) .
Lemma 2. 
Suppose design d is an orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) , then
A M D u ( d ) = Φ u ,
where | u | = g 1 + g 2 , 1 g 1 + g 2 t , Φ u is a constant only depending on q 1 , q 2 , g 1 and g 2 . In particular,
( i ) when both q 1 and q 2 are even,
Φ u = 7 12 g 2 28 q 1 2 + 1 48 q 1 2 g 1 28 q 2 2 + 1 48 q 2 2 g 2 + 7 q 1 2 + 2 12 q 1 2 g 1 7 q 2 2 + 2 12 q 2 2 g 2 ;
( i i ) when both q 1 and q 2 are odd,
Φ u = 7 12 g 2 7 q 1 2 + 1 12 q 1 2 g 1 7 q 2 2 + 1 12 q 2 2 g 2 + 14 q 1 2 + 7 24 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ;
( i i i ) when q 1 is even and q 2 is odd,
Φ u = 7 12 g 2 28 q 1 2 + 1 48 q 1 2 g 1 7 q 2 2 + 1 12 q 2 2 g 2 + 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 .
It is well known that strength is an important measure of orthogonality. For comparing the difference between design d U ( n ; q 1 s 1 × q 2 s 2 ) and orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) of strength t, the definition of uniformity pattern of design d is given as follows, which provides a measure of the projection uniformity of d onto different dimensions.
Definition 1. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , any positive integer g ( = g 1 + g 2 ) s and u S g , define
M I g ( d ) = | u | = g A M D u ( d ) Φ u ,
where Φ u is shown in Lemma 2. The vector ( M I 1 ( d ) , , M I s ( d ) ) is called the uniformity pattern of design d.
We now state the above discussion as the following theorem, which gives a relationship between the uniformity pattern ( M I 1 ( d ) , , M I s ( d ) ) of design d and the strength t of orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) .
Theorem 1. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , design d is an orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) if and only if M I k ( d ) = 0 for k = 1 , , t and M I t + 1 ( d ) 0 .
Theorem 1 indicates that there is a close relationship between M I t ( d ) and strength t for a design d U ( n ; q 1 s 1 × q 2 s 2 ) , that is, the smaller the value of M I t ( d ) , the design d will be closer to an orthogonal array of strength t. Based on Theorem 1, { M I k ( d ) } may be used as a measure for evaluating designs; it suggests to define some similar criteria, such as MPU.
Definition 2. 
For two designs d 1 , d 2 U ( n ; q 1 s 1 × q 2 s 2 ) , there is an integer t such that M I t ( d 1 ) M I t ( d 2 ) and M I k ( d 1 ) = M I k ( d 2 ) for k = 1 , , t 1 ; then, d 1 is said to have less MPU than d 2 . If there is no other design in U ( n ; q 1 s 1 × q 2 s 2 ) that has less MPU than d 1 , then d 1 is said to have MPU, or d 1 is an MPU design.
Here, we mainly establish the connections between projection uniformity and orthogonality, and some relationships between criteria of MPU and GMA will also be included.
Theorem 2. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , any positive integer g ( = g 1 + g 2 ) s and u S g , we have
M I g ( d ) = | u | = g α g 1 g 2 ( r 1 , r 2 ) R β r 1 r 2 s 1 r 1 s 1 g 1 s 2 r 2 s 2 g 2 A r 1 + r 2 ( d ) ,
where R = { ( r 1 , r 2 ) : r 1 = 0 , , g 1 , r 2 = 0 , , g 2 , ( r 1 , r 2 ) ( 0 , 0 ) } , and
( i ) when both q 1 and q 2 are even,
α g 1 g 2 = 7 q 1 2 + 2 12 q 1 2 g 1 7 q 2 2 + 2 12 q 2 2 g 2 , β r 1 r 2 = 2 q 1 + 2 7 q 1 2 + 2 r 1 2 q 2 + 2 7 q 2 2 + 2 r 2 ;
( i i ) when both q 1 and q 2 are odd,
α g 1 g 2 = 14 q 1 2 + 7 24 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 , β r 1 r 2 = 4 q 1 + 4 14 q 1 2 + 7 r 1 4 q 2 + 4 14 q 2 2 + 7 r 2 ;
( i i i ) when q 1 is even and q 2 is odd,
α g 1 g 2 = 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 , β r 1 r 2 = 2 q 1 + 2 7 q 1 2 + 2 r 1 4 q 2 + 4 14 q 2 2 + 7 r 2 .

4. A Lower Bound of Uniformity Pattern

This section provides a lower bound of uniformity pattern defined in Definition 1. It is very important that the lower bounds of uniformity pattern can be served as a benchmark not only in searching for uniform designs with minimum projection uniformity but also in helping to validate that some good designs are in fact uniform.
Define Δ u = q 1 e v 1 + p 1 e v 3 + q 2 ( e v 2 e v 3 ) when p 1 > q 2 , and Δ u = p 2 e v 1 + q 2 e v 4 + p 1 ( e v 2 e v 4 ) when p 1 q 2 .
Theorem 3. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) and positive integer g ( = g 1 + g 2 ) s , we have
M I g ( d ) L M I g ( d ) ,
( i ) when both q 1 and q 2 are even,
L M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 7 q 1 2 12 q 1 g 1 7 q 2 2 12 q 2 g 2 u S g Δ u ,
where Ψ g 1 g 2 = 1 n ( 3 4 ) g 1 ( 3 4 ) g 2 ( 7 q 1 2 + 2 12 q 1 2 ) g 1 ( 7 q 2 2 + 2 12 q 2 2 ) g 2 ;
( i i ) when both q 1 and q 2 are odd,
L M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 14 q 1 2 4 q 1 + 3 24 q 1 2 g 1 14 q 2 2 4 q 2 + 3 24 q 2 2 g 2 u S g Δ u ,
where Ψ g 1 g 2 = 1 n ( 6 q 1 2 + 1 8 q 1 2 ) g 1 ( 6 q 2 2 + 1 8 q 2 2 ) g 2 ( 14 q 1 2 + 7 24 q 1 2 ) g 1 ( 14 q 2 2 + 7 24 q 2 2 ) g 2 ;
( i i i ) when q 1 is even and q 2 is odd,
L M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 7 q 1 2 12 q 1 g 1 14 q 2 2 4 q 2 + 3 24 q 2 2 g 2 u S g Δ u ,
where Ψ g 1 g 2 = 1 n ( 3 4 ) g 1 ( 6 q 2 2 + 1 8 q 2 2 ) g 2 ( 7 q 1 2 + 2 12 q 1 2 ) g 1 ( 14 q 2 2 + 7 24 q 2 2 ) g 2 .
Theorem 4. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) and positive integer g ( = g 1 + g 2 ) s ,
M I g ( d ) L M I g ( d ) ,
( i ) when both q 1 and q 2 are even,
L M I g ( d ) = | u | = g [ 1 n 2 7 q 1 2 12 q 1 g 1 7 q 2 2 12 q 2 g 2 g 1 i 1 g 2 i 2 × 2 q 1 + 2 7 q 1 2 i 1 2 q 2 + 2 7 q 2 2 i 2 θ i 1 i 2 7 q 1 + 2 12 q 1 g 1 7 q 2 + 2 12 q 2 g 2 ] ;
( i i ) when both q 1 and q 2 are odd,
L M I g ( d ) = | u | = g [ 1 n 2 14 q 1 2 4 q 1 + 3 24 q 1 g 1 14 q 2 2 4 q 2 + 3 24 q 2 g 2 g 1 i 1 g 2 i 2 × 4 q 1 2 + 4 q 1 14 q 1 2 4 q 1 + 3 i 1 4 q 2 2 + 4 q 2 14 q 2 2 4 q 2 + 3 i 2 θ i 1 i 2 14 q 1 2 + 7 24 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ] ;
( i i i ) when q 1 is even and q 2 is odd,
L M I g ( d ) = | u | = g [ 1 n 2 7 q 1 2 12 q 1 g 1 14 q 2 2 4 q 2 + 3 24 q 2 g 2 g 1 i 1 g 2 i 2 × 2 q 1 + 2 7 q 1 2 i 1 4 q 2 2 + 4 q 2 14 q 2 2 4 q 2 + 3 i 2 θ i 1 i 2 7 q 1 + 2 12 q 1 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ] ,
where θ i 1 i 2 = n λ i 1 i 2 + μ i 1 i 2 ( 1 + λ i 1 i 2 ) , μ i 1 i 2 = n q 1 i 1 q 2 i 2 λ i 1 i 2 , λ i 1 i 2 be the largest integer contained in n / ( q 1 i 1 q 2 i 2 ) .
Note that Theorem 3 is based on Hamming distances between any two runs of d, but Theorem 4 comes from the quadratic form y d T D y d in Appendix A Equation (A1). Some numerical examples show that these two lower bounds are not tight simultaneously. Therefore, we give another lower bound of uniformity pattern as the following theorem:
Theorem 5. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) and positive integer g ( = g 1 + g 2 ) s , we have
M I g ( d ) L M I g * ( d ) ,
where L M I g * ( d ) = max { L M I g ( d ) , L M I g ( d ) } .

5. Illustrative Examples

In this section, some numerical examples are provided to illustrate our theoretical results.
Example 1. 
Consider a design d 1 U ( 4 ; 2 3 × 4 3 ) , which are given below:
d 1 = 0 0 0 0 3 2 1 0 1 2 0 1 0 1 1 1 2 0 1 1 0 3 1 3 .
The number of columns in design d 1 is greater than the number of rows, its uniformity pattern in Definition 1, and its lower bound values in Theorems 3–5 are listed in Table 1.
It is clear that d 1 is an orthogonal array of strength 1 and attains the lower bounds in Theorem 3.
Example 2. 
Consider design d 2 U ( 20 ; 2 3 × 5 ) and d 3 U ( 48 ; 2 5 × 3 ) , which are given below,
d 2 = 0 0 1 1 1 1 0 0 0 0 1 1 1 1 0 0 0 0 1 1 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1 0 0 1 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1 0 0 1 0 1 1 0 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 T , d 3 = 1111111100000000 0000000011111111 1111111100000000 1111000011110000 0000111100001111 1111000011110000 1100110011001100 0011001100110011 1100110011001100 1010101010101010 0101010101010101 1010101010101010 1001011001101001 0110100110010110 1001011001101001 0000000000000000 1111111111111111 2222222222222222 T .
The number of rows in designs d 2 and d 3 are greater than the number of columns, and the numerical results of both are shown in Table 2.
As can be seen from Table 2, designs d 2 and d 3 are an orthogonal array with strengths of 2 and 4, respectively, and both reach the lower bound in Theorem 4.
It can be seen from Table 1 and Table 2 that the lower bounds of uniformity pattern of designs d 1 , d 2 , and d 3 are achieved, so d 1 , d 2 , and d 3 are all MPU designs. We can also see that L M I g ( d ) is better than L M I g ( d ) for large n and smaller s. Similar to the findings of Fang et al. (2018) [24], none of the lower bounds in Theorems 3 and 4 are absolutely dominant for all combinations of the number of runs n and of factors s. Therefore, we choose the maximum value of Theorems 3–5.

6. Conclusions

In this paper, the projection uniformity and related properties under mixture discrepancy of asymmetric factorials are explored. The relationship between uniformity pattern and generalized minimum aberration is established. A lower bound of uniformity pattern is also obtained, which can be served as a benchmark for searching minimum projection uniformity designs. These results provide a theoretical basis for searching optimal asymmetric designs with minimum projection uniformity measured by average projection mixture discrepancy. Overall, this paper extends the results of [13,14,15] to the asymmetric case, which makes the corresponding theory more flexible.
The results in this paper can be extended to any asymmetric designs d U ( N ; q 1 s 1 × · · · × q n s n ) . Taking the first t factors as even and the last n t factors as odd, and using some simple calculation of tired multiplication, similar definition and results of uniformity pattern and lower bounds can be obtained.

Author Contributions

Conceptualization, Z.O. and K.W.; methodology, N.Z.; software, K.W. and N.Z.; validation, N.Z. and K.W.; resources, H.Q.; data curation, N.Z. and K.W.; writing—original draft preparation, K.W.; writing—review and editing, H.Q., N.Z., and Z.O.; funding acquisition, H.Q. and Z.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 11871237, 11961027, 12161040, 11701213), the Natural Science Foundation of Hunan Provincial (Grant Nos. 2021JJ30550, 2020JJ4497), the Scientific Research Plan Item of Hunan Provincial Department of Education (Grant No. 22A0355), and the Discipline Coordination Construction Project of Zhongnan University of Economics and Law (Grant No. XKHJ202125).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors have no conflict of interest regarding this paper.

Appendix A

Proof of Lemma 2.
If a design d U ( n ; q 1 s 1 × q 2 s 2 ) is an orthogonal array O A ( n ; q 1 s 1 × q 2 s 2 , t ) of strength t; then, for any nonnegative integer g ( = g 1 + g 2 ) s and u = u 1 u 2 S g 1 g 2 , all possible q 1 g 1 × q 2 g 2 level combinations among any g columns of projection design d u appear equally often. Given row i u 0 = ( i u 1 0 , i u 2 0 ) d u , it is easy to obtain that | { ( i u 0 , k u ) : H i 0 k u 1 = j 1 , H i 0 k u 2 = j 2 , k u d } | = g 1 j 1 g 2 j 2 n ( q 1 1 ) j 1 ( q 2 1 ) j 2 q 1 g 1 q 2 g 2 .
Therefore, the third term in the right side of Formula (4) can be expressed as
1 n 3 4 g 1 6 q 2 2 + 1 8 q 2 2 g 2 i 1 = 0 g 1 i 2 = 0 g 2 7 q 1 2 9 q 1 i 1 14 q 2 2 4 q 2 + 3 18 q 2 2 + 3 i 2 E i 1 i 2 ( d u ) = 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ,
which completes the proof. □
Proof of Theorem 2.
From formulas (1), (3), (4) and Definition 1, we have
M I g ( d ) = | u | = g [ 1 n 3 4 g 1 6 q 2 2 + 1 8 q 2 2 g 2 i 1 = 0 g 1 i 2 = 0 g 2 7 q 1 2 9 q 1 i 1 14 q 2 2 4 q 2 + 3 18 q 2 2 + 3 i 2 E i 1 i 2 ( d u ) 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ] = | u | = g [ 3 4 q 1 g 1 6 q 2 2 + 1 8 q 2 3 g 2 r 1 = 0 g 1 r 2 = 0 g 2 i 1 = 0 g 1 i 2 = 0 g 2 7 q 1 2 9 q 1 i 1 14 q 2 2 4 q 2 + 3 18 q 2 2 + 3 i 2 × P i 1 ( r 1 ; g 1 , q 1 ) P i 2 ( r 2 ; g 2 , q 2 ) E r 1 r 2 ( d u ) 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ] = | u | = g [ 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 r 1 = 0 g 1 r 2 = 0 g 2 2 q 1 + 2 7 q 1 2 + 2 r 1 4 q 2 + 4 14 q 2 2 + 7 r 2 E r 1 r 2 ( d u ) 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ] = | u | = g 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ( r 1 , r 2 ) R 2 q 1 + 2 7 q 1 2 + 2 r 1 4 q 2 + 4 14 q 2 2 + 7 r 2 E r 1 r 2 ( d u ) = | u | = g 7 q 1 2 + 2 12 q 1 2 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ( r 1 , r 2 ) R 2 q 1 + 2 7 q 1 2 + 2 r 1 4 q 2 + 4 14 q 2 2 + 7 r 2 × s 1 r 1 s 1 g 1 s 2 r 2 s 2 g 2 A r ( d ) ,
which completes the proof. □
In order to prove Theorem 3, we need to know Lemmas A1–A3, where Lemma A1 can be obtained from Lemma 1 and Definition 1.
Lemma A1. 
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , positive integer g ( = g 1 + g 2 ) s and u = u 1 u 2 S g ,
( i ) when both q 1 and q 2 are even,
M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 | u | = g 7 q 1 2 12 q 1 g 1 7 q 2 2 12 q 2 g 2 i = 1 n k ( i ) = 1 n e θ i k u ,
where Ψ g 1 g 2 is shown in Theorem 3, θ i k u = ln ( 9 q 1 7 q 1 2 ) · δ i k u 1 + ln ( 9 q 2 7 q 2 2 ) · δ i k u 2 ;
( i i ) when both q 1 and q 2 are odd,
M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 | u | = g 14 q 1 2 4 q 1 + 3 24 q 1 2 g 1 14 q 2 2 4 q 2 + 3 24 q 2 2 g 2 i = 1 n k ( i ) = 1 n e θ i k u ,
where Ψ g 1 g 2 is shown in Theorem 3, θ i k u = ln ( 18 q 1 2 + 3 14 q 1 2 4 q 1 + 3 ) · δ i k u 1 + ln ( 18 q 2 2 + 3 14 q 2 2 4 q 2 + 3 ) · δ i k u 2 ;
( i i i ) when q 1 is even and q 2 is odd,
M I g ( d ) = | u | = g Ψ g 1 g 2 + 1 n 2 | u | = g 7 q 1 2 12 q 1 g 1 14 q 2 2 4 q 2 + 3 24 q 2 2 g 2 i = 1 n k ( i ) = 1 n e θ i k u ,
where Ψ g 1 g 2 is shown in Theorem 3, θ i k u = ln ( 9 q 1 7 q 1 2 ) · δ i k u 1 + ln ( 18 q 2 2 + 3 14 q 2 2 4 q 2 + 3 ) · δ i k u 2 .
The proof of Lemma A1 is similar to [14], so it is omitted.
Lemma A2 
([25]). For any design d U ( n ; q s ) and positive integer t, we have
i = 1 n k ( i ) = 1 n ( δ i k ) t = P w t + Q ( w + 1 ) t .
where w = ( n q ) s q ( n 1 ) , P and Q are integers such that P + Q = n ( n 1 ) , and A means the largest integer contained in A.
Lemma A3 
([26]). For any design d U ( n ; q 1 s 1 × q 2 s 2 ) and positive integer t, we have
i = 1 n k ( i ) = 1 n θ i k = α 1 n ( n q 1 ) g 1 q 1 + α 2 n ( n q 2 ) g 2 q 2 , a n d i = 1 n k ( i ) = 1 n ( θ i k ) t Q 1 v 1 t + Q 2 v 2 t + ( P 1 Q 2 ) v 3 t , w h e n P 1 > Q 2 ; P 2 v 1 t + P 1 v 2 t + ( Q 2 P 1 ) v 4 t , w h e n P 1 Q 2 .
where α 1 > 0 and α 2 > 0 are weights, P 1 and Q 1 are integers such that P 1 + Q 1 = n ( n 1 ) and P 1 w 1 + Q 1 ( w 1 + 1 ) = n ( n q 1 ) s 1 / q 1 , P 2 and Q 2 are integers such that P 2 + Q 2 = n ( n 1 ) and P 2 w 2 + Q 2 ( w 2 + 1 ) = n ( n q 2 ) s 2 / q 2 . Let v 1 = α 1 ( w 1 + 1 ) + α 2 w 2 , v 2 = α 1 w 1 + α 2 ( w 2 + 1 ) , v 3 = α 1 w 1 + α 2 w 2 , v 4 = α 1 ( w 1 + 1 ) + α 2 ( w 2 + 1 ) , w 1 = ( n q 1 ) s 1 q 1 ( n 1 ) , w 2 = ( n q 2 ) s 2 q 2 ( n 1 ) .
Proof of Theorem 4.
According to [23,24], let I q and 1 q respectively be the q × q identity matrix and the q × 1 vector with all elements unity, define
L ( 0 ) = 1 q T , L ( 1 ) = I q , J q = 1 q 1 q T .
Let D g 1 ( 1 ) and D g 2 ( 2 ) be the g 1 -fold and g 2 -fold Kronecker products of D 0 ( 1 ) and D 0 ( 2 ) , respectively. Let Ω be the set of all binary ( q 1 + q 2 ) tuples, Ω i 1 i 2 be the set of Ω consisting of those binary ( g 1 + g 2 ) -tuples with exactly i 1 elements of x 1 unity and i 2 elements of x 2 unity, respectively, where Ω = { x = ( x ( 1 ) , x ( 2 ) ) : x ( 1 ) = ( x 1 ( 1 ) , , x g 1 ( 1 ) ) Ω ( 1 ) , x ( 2 ) = ( x 1 ( 2 ) , , x g 2 ( 2 ) ) Ω ( 2 ) } .
D = D g 1 ( 1 ) D g 2 ( 2 ) , D g 1 ( 1 ) = i 1 = 1 g 1 D 0 ( 1 ) , D g 2 ( 2 ) = i 2 = 1 g 2 D 0 ( 2 ) .
For any design d U ( n ; q 1 s 1 × q 2 s 2 ) , Lemma A1 gives an expression between the uniformity pattern and the number of coincide. Based on this, we can obtain
( i ) when both q 1 and q 2 are even,
D 0 ( 1 ) = q 1 + 1 6 q 1 I q 1 + 7 q 1 2 12 q 1 J q 1 , D 0 ( 2 ) = q 2 + 1 6 q 2 I q 2 + 7 q 2 2 12 q 2 J q 2 ;
( i i ) when both q 1 and q 2 are odd,
D 0 ( 1 ) = q 1 + 1 6 q 1 I q 1 + 14 q 1 2 4 q 1 + 3 24 q 1 2 J q 1 , D 0 ( 2 ) = q 2 + 1 6 q 2 I q 2 + 14 q 2 2 4 q 2 + 3 24 q 2 2 J q 2 ;
( i i i ) when q 1 is even and q 2 is odd,
D 0 ( 1 ) = q 1 + 1 6 q 1 I q 1 + 7 q 1 2 12 q 1 J q 1 , D 0 ( 2 ) = q 2 + 1 6 q 2 I q 2 + 14 q 2 2 4 q 2 + 3 24 q 2 2 J q 2 .
Considering the case ( i i i ) where q 1 is even and q 2 is odd, we have
M I g ( d ) = | u | = g 1 n 2 y d T D y d 7 q 1 + 2 12 q 1 g 1 14 q 2 2 + 7 24 q 2 2 g 2 ,
where
D = γ g 1 g 2 x ( 1 ) Ω ( 1 ) x ( 2 ) Ω ( 2 ) 2 q 1 + 2 7 q 1 2 x i ( 1 ) 4 q 2 2 + 4 q 2 14 q 2 2 4 q 2 + 3 x i ( 2 ) H ( x ) H ( x ) , y d D y d = γ g 1 g 2 i 1 = 0 g 1 i 2 = 0 g 2 2 q 1 + 2 7 q 1 2 i 1 4 q 2 2 + 4 q 2 14 q 2 2 4 q 2 + 3 i 2 x Ω i 1 i 2 y d H ( x ) H ( x ) y d ,
and γ g 1 g 2 = ( q 1 + 1 6 q 1 ) g 1 ( q 2 + 1 6 q 2 ) g 2 .
Let y d ( x ) be the number of times the treatment combination x occurs in d and y d be the n × 1 vector with elements y d ( x ) arranged in the lexicographic order. For any x Ω i 1 i 2 , the elements of the q 1 i 1 q 2 i 2 × 1 vector H ( x ) y d are nonnegative integers with sum n; then, by [24], we have
y d H ( x ) H ( x ) y d λ i 1 i 2 2 ( q 1 i 1 q 2 i 2 μ i 1 i 2 ) + ( λ i 1 i 2 + 1 ) 2 μ i 1 i 2 = n λ i 1 i 2 + μ i 1 i 2 ( λ i 1 i 2 + 1 ) ,
which completes the proof of Case ( i i i ) .
The proof of Case ( i ) and Case ( i i ) are similar to Case ( i i i ) . □

References

  1. Hickernell, F.J.; Liu, M.Q. Uniformity designs limit aliasing. Biometrika 2002, 89, 893–904. [Google Scholar] [CrossRef]
  2. Xu, H.Q.; Wu, C.F.J. Generalized minimum aberration for asymmetrical fractional factorial designs. Ann. Stat. 2001, 29, 1066–1077. [Google Scholar]
  3. Ma, C.X.; Fang, K.T.; Lin, D.K.J. A note on uniformity and orthogonality. J. Stat. Plan. Inference 2003, 113, 323–334. [Google Scholar] [CrossRef]
  4. Hickernell, F.J. A generalized discrepancy and quadrature error bound. Math. Comput. 1998, 67, 299–322. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, M.Q.; Fang, K.T.; Hickernell, F.J. Connections among different criteria for asymmetrical fractional factorial designs. Stat. Sin. 2006, 16, 1285–1297. [Google Scholar]
  6. Ma, C.X.; Fang, K.T. A note on generalized aberration factorial designs. Metrika 2001, 53, 85–93. [Google Scholar] [CrossRef]
  7. Fang, K.T.; Qin, H. Uniformity pattern and related criteria for two-level factorials. Sci. China Ser. Math. 2005, 48, 1–11. [Google Scholar] [CrossRef]
  8. Hickernell, F.J. Lattice rules: How well do they measure up? In Random and Quasi-Random Point Sets; Hellekalek, P., Larcher, G., Eds.; Springer: New York, NY, USA, 1998; pp. 106–166. [Google Scholar]
  9. Zhou, Y.D.; Fang, K.T.; Ning, J.H. Mixture discrepancy for quasi-random point sets. J. Complex. 2013, 29, 283–301. [Google Scholar] [CrossRef]
  10. Tang, Y.; Xu, H.Q.; Lin, D.K.J. Uniform fractional factorial designs. Ann. Stat. 2012, 40, 891–907. [Google Scholar] [CrossRef] [Green Version]
  11. Zhou, Y.D.; Xu, H.Q. Space-filling fractional factorial designs. J. Am. Stat. Assoc. 2014, 109, 1134–1144. [Google Scholar] [CrossRef]
  12. Chen, W.; Qi, Z.F.; Zhou, Y.D. Constructing uniformity designs under mixture discrepancy. Stat. Probab. Lett. 2015, 97, 76–82. [Google Scholar] [CrossRef]
  13. Yi, S.Y.; Zhou, Y.D. Projection uniformity under mixture discrepancy. Stat. Probab. Lett. 2018, 140, 96–105. [Google Scholar] [CrossRef]
  14. Wang, K.; Ou, Z.J.; Liu, J.Q.; Li, H.Y. Uniformity pattern of q-level factorials under mixture discrepancy. Stat. Pap. 2021, 62, 1777–1793. [Google Scholar] [CrossRef]
  15. Wang, K.; Qin, H.; Ou, Z.J. Uniformity pattern of mixed two-and three-level factorials under average projection mixture discrepancy. Statistics 2022, 56, 121–133. [Google Scholar] [CrossRef]
  16. Sun, F.S.; Wang, Y.P.; Xu, H.Q. Uniform projection designs. Ann. Stat. 2019, 47, 641–661. [Google Scholar] [CrossRef]
  17. Chen, H.; Zhang, Y.; Yang, X. Uniform projection nested Latin hypercube designs. Stat. Pap. 2021, 62, 2031–2045. [Google Scholar] [CrossRef]
  18. Liu, S.X.; Wang, Y.P.; Sun, F.S. Two-dimensional projection uniformity for space-filling designs. Can. J. Stat. 2022. [Google Scholar] [CrossRef]
  19. Li, W.L.; Liu, M.Q.; Yang, J.F. Construction of column-orthogonal strong orthogonal arrays. Stat. Pap. 2022, 63, 515–530. [Google Scholar] [CrossRef]
  20. Sun, C.Y.; Tang, B.X. Uniform projection designs and strong orthogonal arrays. J. Am. Stat. Assoc. 2021. [Google Scholar] [CrossRef]
  21. Wang, Y.P.; Sun, F.S.; Xu, H.Q. On Design Orthogonality, Maximin Distance, and Projection Uniformity for Computer Experiments. J. Am. Stat. Assoc. 2022, 117, 375–385. [Google Scholar] [CrossRef]
  22. Zhou, Y.S.; Xiao, Q.; Sun, F.S. Construction of uniform projection designs via level permutation and expansion. J. Stat. Plan. Inference 2023, 222, 209–225. [Google Scholar] [CrossRef]
  23. Chatterjee, K.; Li, Z.H.; Qin, H. Some new lower bounds to centered and wrap-round L2-discrepancies. Stat. Probab. Lett. 2012, 82, 1367–1373. [Google Scholar] [CrossRef]
  24. Fang, K.T.; Liu, M.Q.; Qin, H.; Zhou, Y.D. Theory and Application of Uniform Experimental Designs; Springer: Singapore, 2018. [Google Scholar]
  25. Zhang, Q.H.; Wang, Z.H.; Hu, J.W.; Qin, H. A new lower bound for wrap-around L2-discrepancy on two and three mixed level factorials. Stat. Probab. Lett. 2015, 96, 133–140. [Google Scholar] [CrossRef]
  26. Fang, K.T.; Ma, C.X.; Mukerjee, R. Uniformity in Fractional Factorials. In Monte Carlo and Quasi-Monte Carlo Methods; Springer: Berlin/Heidelberg, Germany, 2002; pp. 232–241. [Google Scholar]
Table 1. Numerical results of designs d 1 .
Table 1. Numerical results of designs d 1 .
g123456
M I g ( d 1 ) 00.08300.21930.21700.09540.0157
L M I g ( d 1 ) 00.08300.21930.21700.09540.0157
L M I g ( d 1 ) 00.01460.03970.04290.03180.0157
L M I g * ( d 1 ) 00.08300.21930.21700.09540.0157
Table 2. Numerical results of designs d 2 and d 3 .
Table 2. Numerical results of designs d 2 and d 3 .
g123456
M I g ( d 2 ) 00 7.8125 × 10 5 1.2148 × 10 4
L M I g ( d 2 ) 0−0.0450−0.0282−0.0040
L M I g ( d 2 ) 00 7.8125 × 10 5 1.2148 × 10 4
L M I g * ( d 2 ) 00 7.8125 × 10 5 1.2148 × 10 4
M I g ( d 3 ) 0000 3.3908 × 10 6 4.0973 × 10 6
L M I g ( d 3 ) 0−0.1837−0.1742−0.1300−0.0365−0.0044
L M I g ( d 3 ) 0000 3.3908 × 10 6 4.0973 × 10 6
L M I g * ( d 3 ) 0000 3.3908 × 10 6 4.0973 × 10 6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, K.; Ou, Z.; Qin, H.; Zou, N. Projection Uniformity of Asymmetric Fractional Factorials. Axioms 2022, 11, 716. https://doi.org/10.3390/axioms11120716

AMA Style

Wang K, Ou Z, Qin H, Zou N. Projection Uniformity of Asymmetric Fractional Factorials. Axioms. 2022; 11(12):716. https://doi.org/10.3390/axioms11120716

Chicago/Turabian Style

Wang, Kang, Zujun Ou, Hong Qin, and Na Zou. 2022. "Projection Uniformity of Asymmetric Fractional Factorials" Axioms 11, no. 12: 716. https://doi.org/10.3390/axioms11120716

APA Style

Wang, K., Ou, Z., Qin, H., & Zou, N. (2022). Projection Uniformity of Asymmetric Fractional Factorials. Axioms, 11(12), 716. https://doi.org/10.3390/axioms11120716

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop