University, Wuhan, China c College of Physical Science and Technology, Central China Normal. University, Wuhan, China Published online: 25 Apr 2014.

Similar documents
Park, Pennsylvania, USA. Full terms and conditions of use:

Gilles Bourgeois a, Richard A. Cunjak a, Daniel Caissie a & Nassir El-Jabi b a Science Brunch, Department of Fisheries and Oceans, Box

Use and Abuse of Regression

Testing Goodness-of-Fit for Exponential Distribution Based on Cumulative Residual Entropy

Published online: 17 May 2012.

Nacional de La Pampa, Santa Rosa, La Pampa, Argentina b Instituto de Matemática Aplicada San Luis, Consejo Nacional de Investigaciones Científicas

The American Statistician Publication details, including instructions for authors and subscription information:

Guangzhou, P.R. China

Dissipation Function in Hyperbolic Thermoelasticity

Published online: 05 Oct 2006.

Precise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions

University, Tempe, Arizona, USA b Department of Mathematics and Statistics, University of New. Mexico, Albuquerque, New Mexico, USA

Communications in Algebra Publication details, including instructions for authors and subscription information:

Full terms and conditions of use:

Open problems. Christian Berg a a Department of Mathematical Sciences, University of. Copenhagen, Copenhagen, Denmark Published online: 07 Nov 2014.

Characterizations of Student's t-distribution via regressions of order statistics George P. Yanev a ; M. Ahsanullah b a

A note on optimal foldover design

Ankara, Turkey Published online: 20 Sep 2013.

Online publication date: 22 March 2010

The Homogeneous Markov System (HMS) as an Elastic Medium. The Three-Dimensional Case

George L. Fischer a, Thomas R. Moore b c & Robert W. Boyd b a Department of Physics and The Institute of Optics,

MSE Performance and Minimax Regret Significance Points for a HPT Estimator when each Individual Regression Coefficient is Estimated

Version of record first published: 01 Sep 2006.

Tong University, Shanghai , China Published online: 27 May 2014.

Geometric View of Measurement Errors

Online publication date: 30 March 2011

To cite this article: Edward E. Roskam & Jules Ellis (1992) Reaction to Other Commentaries, Multivariate Behavioral Research, 27:2,

Online publication date: 01 March 2010 PLEASE SCROLL DOWN FOR ARTICLE

CCSM: Cross correlogram spectral matching F. Van Der Meer & W. Bakker Published online: 25 Nov 2010.

Full terms and conditions of use:

Dresden, Dresden, Germany Published online: 09 Jan 2009.

PLEASE SCROLL DOWN FOR ARTICLE

Derivation of SPDEs for Correlated Random Walk Transport Models in One and Two Dimensions

OF SCIENCE AND TECHNOLOGY, TAEJON, KOREA

A Simple Approximate Procedure for Constructing Binomial and Poisson Tolerance Intervals

Diatom Research Publication details, including instructions for authors and subscription information:

The Fourier transform of the unit step function B. L. Burrows a ; D. J. Colwell a a

University of Thessaloniki, Thessaloniki, Greece

Erciyes University, Kayseri, Turkey

PLEASE SCROLL DOWN FOR ARTICLE

Discussion on Change-Points: From Sequential Detection to Biology and Back by David Siegmund

Tore Henriksen a & Geir Ulfstein b a Faculty of Law, University of Tromsø, Tromsø, Norway. Available online: 18 Feb 2011

Optimal Foldover Plans for Two-Level Fractional Factorial Designs

G. S. Denisov a, G. V. Gusakova b & A. L. Smolyansky b a Institute of Physics, Leningrad State University, Leningrad, B-

Published online: 10 Apr 2012.

E(s 2 )-OPTIMALITY AND MINIMUM DISCREPANCY IN 2-LEVEL SUPERSATURATED DESIGNS

France. Published online: 23 Apr 2007.

FB 4, University of Osnabrück, Osnabrück

Geometrical optics and blackbody radiation Pablo BenÍTez ab ; Roland Winston a ;Juan C. Miñano b a

PLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:

Communications in Algebra Publication details, including instructions for authors and subscription information:

Acyclic, Cyclic and Polycyclic P n

PLEASE SCROLL DOWN FOR ARTICLE

Online publication date: 12 January 2010

Quan Wang a, Yuanxin Xi a, Zhiliang Wan a, Yanqing Lu a, Yongyuan Zhu a, Yanfeng Chen a & Naiben Ming a a National Laboratory of Solid State

PLEASE SCROLL DOWN FOR ARTICLE

Online publication date: 29 April 2011

Global Existence of Large BV Solutions in a Model of Granular Flow

A. H. Abuzaid a, A. G. Hussin b & I. B. Mohamed a a Institute of Mathematical Sciences, University of Malaya, 50603,

Acta Mathematica Sinica, English Series Springer-Verlag Berlin Heidelberg & The Editorial Office of AMS 2015

Projective spaces in flag varieties

Published online: 04 Jun 2015.

T. Runka a, M. Kozielski a, M. Drozdowski a & L. Szczepańska b a Institute of Physics, Poznan University of

Classification of three-word indicator functions of two-level factorial designs

A NEW CLASS OF NESTED (NEARLY) ORTHOGONAL LATIN HYPERCUBE DESIGNS

Interaction balance in symmetrical factorial designs with generalized minimum aberration

PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE. Full terms and conditions of use:

35-959, Rzeszów, Poland b Institute of Computer Science, Jagiellonian University,

A Strongly Convergent Method for Nonsmooth Convex Minimization in Hilbert Spaces

Fractal Dimension of Turbulent Premixed Flames Alan R. Kerstein a a

Analysis of Leakage Current Mechanisms in BiFeO 3. Thin Films P. Pipinys a ; A. Rimeika a ; V. Lapeika a a

Some Nonregular Designs From the Nordstrom and Robinson Code and Their Statistical Properties

The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features

A note on adaptation in garch models Gloria González-Rivera a a

MINIMUM MOMENT ABERRATION FOR NONREGULAR DESIGNS AND SUPERSATURATED DESIGNS

PLEASE SCROLL DOWN FOR ARTICLE

C.K. Li a a Department of Mathematics and

Yan-Qing Lu a, Quan Wang a, Yuan-Xin Xi a, Zhi- Liang Wan a, Xue-Jing Zhang a & Nai-Ben Ming a a National Laboratory of Solid State

Moment Aberration Projection for Nonregular Fractional Factorial Designs

Móstoles, Madrid. Spain. Universidad de Málaga. Málaga. Spain

Melbourne, Victoria, 3010, Australia. To link to this article:

Avenue G. Pompidou, BP 56, La Valette du Var cédex, 83162, France

P. C. Wason a & P. N. Johnson-Laird a a Psycholinguistics Research Unit and Department of

András István Fazekas a b & Éva V. Nagy c a Hungarian Power Companies Ltd., Budapest, Hungary. Available online: 29 Jun 2011

Connections between the resolutions of general two-level factorial designs

Marko A. A. Boon a, John H. J. Einmahl b & Ian W. McKeague c a Department of Mathematics and Computer Science, Eindhoven

Miami, Florida, USA. Engineering, University of California, Irvine, California, USA. Science, Rehovot, Israel

PLEASE SCROLL DOWN FOR ARTICLE

Sports Technology Publication details, including instructions for authors and subscription information:

INTELLIGENT SEARCH FOR AND MINIMUM ABERRATION DESIGNS

Published online: 04 Oct 2006.

Accepted author version posted online: 28 Nov 2012.

An Approach to Constructing Good Two-level Orthogonal Factorial Designs with Large Run Sizes

Minimax design criterion for fractional factorial designs

Anisotropic Hall effect in Al 13 TM 4

PLEASE SCROLL DOWN FOR ARTICLE

Maximal Rank - Minimum Aberration Regular Two-Level Split-Plot Fractional Factorial Designs

Fokker-Planck Solution for a Neuronal Spiking Model Derek J. Daniel a a

Transcription:

This article was downloaded by: [0.9.78.106] On: 0 April 01, At: 16:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 10795 Registered office: Mortimer House, 7-1 Mortimer Street, London W1T JH, UK Communications in Statistics - Theory and Methods Publication details, including instructions for authors and subscription information: http://www.tandfonline.com/loi/lsta0 A Lower Bound for the Wrap-around L - discrepancy on Combined Designs of Mixed Two- and Three-level Factorials Zujun Ou ac, Hong Qin b & Xu Cai c a College of Mathematics and Statistics, Jishou University, Jishou, China b Faculty of Mathematics and Statistics, Central China Normal University, Wuhan, China c College of Physical Science and Technology, Central China Normal University, Wuhan, China Published online: 5 Apr 01. To cite this article: Zujun Ou, Hong Qin & Xu Cai (01) A Lower Bound for the Wrap-around L - discrepancy on Combined Designs of Mixed Two- and Three-level Factorials, Communications in Statistics - Theory and Methods, :10-1, 7-85 To link to this article: http://dx.doi.org/10.1080/061096.01.77608 PLEASE SCROLL DOWN FOR ARTICLE Taylor & Francis makes every effort to ensure the accuracy of all the information (the Content ) contained in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden. Terms &

Downloaded by [0.9.78.106] at 16:7 0 April 01 Conditions of access and use can be found at http://www.tandfonline.com/page/termsand-conditions

Communications in Statistics Theory and Methods, : 7 85, 01 Copyright Taylor & Francis Group, LLC ISSN: 061-096 print / 15-15X online DOI: 10.1080/061096.01.77608 Regression and Experimental Designs A Lower Bound for the Wrap-around L -discrepancy on Combined Designs of Mixed Two- and Three-level Factorials ZUJUN OU, 1, HONG QIN, AND XU CAI Downloaded by [0.9.78.106] at 16:7 0 April 01 1 College of Mathematics and Statistics, Jishou University, Jishou, China Faculty of Mathematics and Statistics, Central China Normal University, Wuhan, China College of Physical Science and Technology, Central China Normal University, Wuhan, China The objective of this article is to study the issue of employing the uniformity criterion measured by wrap-around L -discrepancy to assess the optimal foldover plans. For mixed two- and three-level fractional factorials as the original designs, general foldover plan and combined design under a foldover plan are defined, and the equivalence between any foldover plan and its complementary foldover plan is investigated. A lower bound of wrap-around L -discrepancy of combined designs under a general foldover plan is obtained, which can be used as a benchmark for searching optimal foldover plans. Moreover, it also provides a theoretical justification for optimal foldover plans in terms of uniformity criterion. Keywords Mixed level factorials; Wrap-around L -discrepancy; Combined design; Foldover plan; Lower bound. Mathematics Subject Classification 6K15; 6K99. 1. Introduction The foldover is a useful technique in construction of two-level factorial designs. Many works on optimal foldover plans in terms of aberration criterion or clear effects criterion have been published. We can refer to Box et al. (1978), Montgomery and Runger (1996), Li and Mee (00), Li and Lin (00), Li et al. (00), Ye and Li (00), Li and Jacroux (007) and Wang et al. (009). Recently, there has been considerable interest in employing the discrepancy to assess the optimal foldover plans. The foldover plan such that the combined design (defined below) Received August 0, 01; Accepted February 11, 01. Address correspondence to Hong Qin, Faculty of Mathematics and Statistics, Central China Normal University, Wuhan 0079, China. E-mail: qinhong@mail.ccnu.edu.cn 7

Lower Bound of Discrepancy on Combined Designs 75 Downloaded by [0.9.78.106] at 16:7 0 April 01 has the smallest discrepancy value over all foldover plans is called the optimal foldover plan. Fang et al. (00) first used the uniformity criterion measured by centered L -discrepancy (Hickernell, 1998) to search the optimal foldover plan. Lei et al. (010, 01) obtained some lower bounds of centered L -discrepancy of combined designs when all Hamming distances between any distinct pair of runs of the initial design are equal. Subsequently, Ou et al. (011) obtained lower bounds of various discrepancies, including wrap-around L -discrepancy, of combined designs for general case. Qin et al. (01) extended the ones of Ou et al. (011) for symmetric designs to a set of asymmetric designs under a special case, where the last p-level column of foldover design is retained as the same in the original design. There are many measures to assess the uniformity of various designs. Among them, the wrap-around L -discrepancy has nice properties, such as it is invariant under reordering the runs, relabeling coordinates, and coordinate shift, we can refer to Fang et al. (005) and Zhou and Ning (008). In this article, we take the wrap-around L -discrepancy as the uniformity measure, and try to find a lower bound of wrap-around L -discrepancy of combined designs. To find a good lower bound of combined designs for a kind of discrepancy is an important issue, because this lower bound can serve as a benchmark for searching optimal foldover plans. While the work of Ou et al. (011) gave a lower bound of wrap-around L -discrepancy, of combined designs under a general foldover plan, the present paper aims at obtaining further results. We extend the findings in Ou et al. (011) from two-level symmetrical factorials to mixed two- and three-level fractional factorials since they are most popular among mixed factorials, and extend the ones in Qin et al. (01) from a kind of special foldover plan to the full foldover plan space if the last factor has three levels. First, the original design is a mixed two- and three-level factorial, which is different from that in Ou et al. (011). Second, the foldover plan considered in Qin et al. (01) is a special case of this article when the last factor has three levels. This article is organized as follows. Some notations and preliminaries are described in Sec.. Section gives the equivalence between any foldover plan and its complementary foldover plan, and a lower bound of wrap-around L -discrepancy of combined design under a general foldover plan also is presented in Sec.. Section gives some examples to illustrate our theoretical results and lends further support to our theoretical results.. Notations and Preliminaries Let U(n; q 1 q s ) denote an U-type design with n runs and s factors, which corresponds to an n s matrix d with entries 0, 1,...,q j 1inthejth column such that all the entries in this column appear equally often. The set of all such designs is denoted by U(n; q 1 q s ). When some q i s are equal, the notations U(n; qs 1 1 qs l l ) and U(n; q s 1 1 qs l l ) are used to express the U-type design and the corresponding set, where s 1 + +s l = s. For any design d U(n; q 1 q s ), note that the rows and columns of d are identified with the runs and factors, respectively. In this article, we mainly focus on the design in U(n; s 1 s ) since two- and three-level mixed factorials are most popular among mixed factorials. Consider a design in U(n; s 1 s ) with the design matrix d = (x 1,x,...,x s ), where x j is the jth column of d,j = 1,...,s. Define Ɣ ={γ = (γ 1,...,γ s1,γ s1 +1,...,γ s ) γ j1 =

76 Ou et al. 0, 1, 1 j 1 s 1 ; γ j = 0, 1,,s 1 j s}, then for any γ = (γ 1,...,γ s ) Ɣ,itdefines a foldover plan for design d U(n; s 1 s ), for 1 j 1 s 1,thej 1 th factor is mapped to x j 1 (γ j1 ) = (x j 1 γ j1 ) = (x 1j1 + γ j1,...,x nj1 + γ j1 ) (mod ) and for s 1 j s, thej th factor is mapped to x j (γ j ) = (x j γ j ) = (x 1j + γ j,...,x nj + γ j ) (mod ). Downloaded by [0.9.78.106] at 16:7 0 April 01 For any design d U(n; s 1 s ) and a foldover plan γ Ɣ, the foldover design, denoted by d γ, is obtained by mapping the columns in d defined above according to the foldover plan γ. Thus, it is to be noted that each foldover design is generated by a foldover plan. The design obtained by augmenting the runs of the foldover design d γ to those of the design d is called the combined design, denoted by d(γ ), that is, d(γ ) = ( d ). d γ Let 1 s be the 1 s vector with all elements unity, for any γ = (γ 1,...,γ s1,γ s1 +1,...,γ s ) Ɣ, define γ = (γ (), 1 s γ () ), where γ () = (γ 1,...,γ s1 ), γ () = (γ s1 +1,...,γ s ) and 1 s γ () = ( γ s1 +1,..., γ s )(mod), and γ is called as the complementary foldover plan of γ in this article. Example.1. Consider a design d 1 U(1; ) with n = 1, s 1 = and s = 1asthe original design, these runs are shown in the first part of Table 1 (Run No. 1 1). If we take γ 1 = (1, 0, 1, 0, ) Ɣ, then we can obtain the foldover design d 1γ1 (Run No. 1 ) by the above foldover method. The original design d and the foldover design d 1γ1 forms the combined design d 1 (γ 1 ) (Run No. 1 ). Furthermore, the complementary foldover plan of γ 1 is γ 1 = (1, 0, 1, 0, 1). It is to be noted that any treatment combination (x 1,x,...,x s ) can be mapped to (u 1,u,...,u s ), where, u j1 = x j 1, 1 j 1 s 1, and u j = x j,s 1 j s. 6 In this article, we consider the wrap-around L -discrepancy, due to Hickernell (1998), given by [WD(d)] = ( ) s n n s i,j=1 l=1 [ ] u il u jl (1 u il u jl ). (1)

Lower Bound of Discrepancy on Combined Designs 77 Table 1 Combined design d 1 (γ 1 ) corresponding to γ 1 d 1 Run No. x 1 x x x x 5 Downloaded by [0.9.78.106] at 16:7 0 April 01 1 0 0 0 0 0 0 1 0 1 0 1 0 1 0 0 1 1 1 1 0 5 0 0 1 1 1 6 0 1 1 0 1 7 1 0 0 0 1 8 1 1 0 1 1 9 0 0 1 1 10 0 1 0 0 11 1 0 0 1 1 1 1 1 0 d 1γ1 Run No. γ 1 = (1 0 1 0 ) 1 1 0 1 0 1 1 1 1 1 15 0 0 0 0 16 0 1 0 1 17 1 0 0 1 0 18 1 1 0 0 0 19 0 0 1 0 0 0 0 1 1 1 0 1 1 0 0 1 1 1 1 1 0 1 0 0 1 1 1 0 1 0 0 1 Under a foldover plan γ, the wrap-around L -discrepancy value of the combined design d(γ ), denoted by WD(d(γ )), can be calculated by the following formula: [WD(d(γ ))] = ( ) s n n s l=s 1 +1 i,j=1 l=1 n s i,j=1 l=1 [ ] u il u jl (1 u il u jl ) ( n s1 ( ) ) u il + ( 1) γl+1 u jl γ l (1 u il + ( 1) γl+1 u jl γ l ) ( ) u il u (γ l) jl (1 u il u (γ l) jl ), ()

78 Ou et al. where u (γ l) jl = xl j (γ l)+1 and x l 6 j (γ l) = (x jl + γ l )(mod). For γ = (γ 1,γ,...,γ s ) Ɣ, lete 1 ={v γ v = 1,v = 1,...,s 1 } and E ={v γ v 0,v = s 1,...,s} of cardinality t 1 and t, respectively, let E 1 ={1,,...,s 1 } E 1 and E ={s 1,s 1 +,...,s} E.LetƔ t1 t be the set of all elements γ Ɣ such that γ exactly have t 1 and t non zero components among {γ 1,...,γ s1 } and {γ s1 +1,...,γ s }, respectively. Given d, s 1, s, t 1, and t,letf(n, s 1,s,t 1,t ) = min{[wd(d(γ ))] : γ Ɣ t1 t }. The optimal foldover plan in this paper is defined as the foldover plan γ such that [WD(d(γ )] = f (n, s 1,s,t 1,t ). The next section attempts to derive a lower bound to [WD(d(γ))].. Lower Bound to the Discrepancy Measure of the Combined Design Downloaded by [0.9.78.106] at 16:7 0 April 01.1. The Equivalence Between γ and γ in Terms of Uniformity Theorem.1. Let d be an U-type design in U(n; s 1 s ), for any γ Ɣ, γ is the complementary foldover plan of γ, then γ is equivalent to γ in terms of uniformity measured by wrap-around L -discrepancy, that is, [WD(d(γ ))] = [WD(d( γ ))]. Proof. For any γ Ɣ, leth ij be the Hamming distance between the ith and jth runs of d(γ ) in algebraic coding theory, h ij (γ ) be the Hamming distance between the ith run of d and jth run of the foldover design d γ.form = 0, 1,...,s, define H m (γ ) = 1 n {(i, j): h ij (γ ) = m, i, j = 1,...,n}, where denotes the cardinality of a set. It is to be note that H m (γ ) = H m ( γ )for m = 0, 1,...,s. Since the third term of the right side of () is decided by {H m (γ )} s m=0, and the first two terms of the right side of () is free from the foldover plan γ, this yields [WD(d(γ))] = [WD(d( γ ))], and γ is equivalent to γ, which completes the proof. Corollary.1. Let d be an U-type design in U(n; s 1 s ), γ Ɣ is the optimal foldover plan of d in the sense of uniformity measured by wrap-aroud L -discrepancy, then the complementary foldover plan γ of γ also is an optimal foldover plan of d. Remark.1. Theorem.1 is very useful to search the optimal foldover plans in the full foldover plan space Ɣ. Specifically, it only requires s 1 ( s )/ times search for a design d U(n; s 1 s ) based on Theorem.1, comparing to exhaustive search method which requires s 1 s times search... A Lower Bound of [WD(d(γ))] Based on the Row Distance It is to be noted that for l E, u il u (γ l) il (1 u il u (γ l) il ) = 9,

Lower Bound of Discrepancy on Combined Designs 79 then we can rewrite Eq. () as follows: [WD(d(γ ))] = ( ) s n n n 1 i j n l=1 1 i j n ( ) s n s ( l E 1 ( ) s1 t 1 ( 5 ) t1 ( [ ] u il u jl (1 u il u jl ) ) s t ( ) t 18 ( ) u il + u jl 1 (1 u il + u jl 1 ) Downloaded by [0.9.78.106] at 16:7 0 April 01 l E 1 ( l E l E ( ) u il u jl (1 u il u jl ) ( ) u il u (γ l) jl (1 u il u (γ l) jl ) ( ) u il u jl (1 u il u jl ). () From the analytical expression of Eq. (), it is easy to see that [WD(d(γ))] is only a function of products of αij l u il u jl (1 u il u jl ), βij l u il+u jl 1 (1 u il +u jl 1 ) and λ l ij u il u (γ l) jl (1 u il u (γ l) jl ). However, for an U-type design, its α-, β-, and λ-values can only be limited to a specific set. More precisely, for an U-type design d U(n; s 1 s ), when 1 l s 1, both α- and β-values can only take possible values, i.e., 0, 1/, when s 1 l s, both α- and λ-values can only take two possible values, i.e., 0, /9. Table gives the distribution of α-values over the set {αij l :1 i j n} for l Ē i,i = 1,, 1 l s 1 and s 1 l s, respectively. Table gives the distribution of β-values over the set {βij l :1 i j n, l E 1} and the distribution of λ-values over the set {λ l ij :1 i j n, l E }. Note that given (n, s 1,s,t 1,t ), these distributions are the same for each design in d U(n; s 1 s ) and each foldover plan γ Ɣ t1 t. The following theorem provides a lower bound to [WD(d(γ))]. Table Distributions of α-values of a d U(n; s 1 s ) and a γ Ɣ t1 t l E 1 1 l s 1 l E s 1 l s α-values Number α-values Number α-values Number α-values Number 0 1 (s 1 t 1 )n(n ) 0 (s 1 t 1 )n 1 s 1 n(n ) 0 s 1 n 9 (s t )n(n ) 0 (s t )n 9 s n(n ) s n

80 Ou et al. Table Distributions of β-values and λ-values of a d U(n; s 1 s ) and a γ Ɣ t1 t l E 1 l E β-values Number λ-values Number 0 t 1 n t n 1 0 t 1 n(n ) 9 t n(n ) Theorem.. Consider a design d U(n; s 1 s ), for any 0 t i s i,i = 1,, and for any γ = (γ 1,...,γ s ) Ɣ t1 t, we have Downloaded by [0.9.78.106] at 16:7 0 April 01 [WD(d(γ ))] ( ) s n ( ) s1n 5 (n 1) ( 18 ( ) 5 s1n t1 (n 1) ( ) s ( 1 + ( 5 6 ) sn (n 1) + n 1 n ( ) sn t (n 1). 18 ) t1 ( ) t ) + n 1 7 n ( ) s1n s1+t1 + sn s+t (n 1) (n 1) ( ) s1(n ) (n 1) + s(n ) (n 1) Moreover, [WD(d(γ))] achieves the above lower bound if all s l=1 ( αl ij ) and all l E 1 ( βl ij ) l E 1 E ( αl ij ) l E ( λl ij ) are the same for 1 i j n. Proof. By Eq. (), to minimize [WD(d(γ ))] over Ɣ for a given design d U(n; s 1 s ) is equivalent to minimizing with respect to α l ij and = 1 i j n l E 1 1 = s 1 i j n l=1 ( ) βl ij l E 1 E ( ) αl ij ( ) ( ) αl ij λl ij l E with respect to α l ij, βl ij, and λl ij. From Tables and, we know that for given (n, s 1,s,t 1,t ), the distributions of α-, β- and λ-values are the same, so both s 1 i j n l=1 ( ) αl ij and 1 i j n l E 1 ( ) βl ij l E 1 E ( ) ( ) αl ij λl ij l E are constants on Ɣ t1 t. Moreover, αl ij > 0, βl ij > 0 and λl ij > 0. Based on the geometric and arithmetic mean inequality, [WD(d(γ))] arrives at its minimum if all s l=1 ( αl ij ) and all l E 1 ( βl ij ) l E 1 E ( αl ij ) l E ( λl ij ) are the same for 1 i j n. The expression of the lower bound of [WD(d(γ))] is straightforward according to Table and Table for 1 and, which completes the proof.

Lower Bound of Discrepancy on Combined Designs 81 Remark.. If we ignore the foldover structure of d(γ ), then d(γ ) can be regarded as a design in U(n; s 1 s ). Following Zhou and Ning (008), a unify lower bound of [WD(d(γ))] for any t 1 = 0, 1,...,s 1 and t = 0, 1,...,s can be obtained as follows: where LB(n, s 1,s ) = [WD(d(γ ))] LB(n, s 1,s ), ( ) s n + n 1 n ( ( ) 1 n 1 ) s [ s 1 (n 1)+ s (n ) ] ( ) 5 s 1 n n 1 ( 18 ) s n (n 1). () Downloaded by [0.9.78.106] at 16:7 0 April 01 Numerical examples (see Examples.1. in Sec. ) show that the lower bounds in Theorem. are tighter than LB(n, s 1,s ) for any t 1 = 0, 1,...,s 1 and t = 0, 1,...,s, since the lower bounds in Theorem. make full use of the foldover structure of d(γ ). Corollary.. If s 1 = 0 and consider a design d U(n; s ), for any 0 t s and for any γ Ɣ 0t, we have [WD(d(γ ))] ( ) s n + n 1 n ( ) sn s+t (n 1) ( ) s ( 1 + ( ) t ) + n 1 7 n ( ) s n t (n 1). 18 ( ) s(n ) (n 1) ( 18 ) sn (n 1) Remark.. If s = 1 and consider a design d U(n; s 1 ), for any 0 t 1 s 1, t = 0, 1, and for any γ Ɣ t1 t, one can obtain a lower bound of [WD(d(γ ))] by Theorem.. For any design d U(n; s 1 p)(p ), Qin et al. (01) also obtained a lower bound of centered L -discrepancy of combined design of d for any 0 t 1 s 1, t = 0, that is, for any γ Ɣ t1 0 which the levels of the last factor are retained as they are in the design d.. Numerical Examples In this section, we give some examples to show our theoretical results. For convenience, we denote the lower bound of [WD(d(γ))] which given in Theorem. as LB(t 1,t ) for any 0 t i s i, i = 1,, and γ Ɣ t1 t in the following Tables. One of the optimal foldover plan over the full foldover plan set Ɣ of d, the corresponding minimum squared wrap-around L -discrepancy value and its lower bound are marked by bold font. Given design d U(n; s 1 s ), one can find that the lower bounds LB(t 1,t ) in Theorem. are tighter than the lower bound LB(n, s 1,s ) given in Remark. for any t 1 = 0, 1,...,s 1 and t = 0, 1,...,s in the following four examples. Example.1. (Continued) Consider the design d 1 U(1; ) as shown in Example.1. Table shows the numerical results of d 1.From(),wehaveLB(1,, 1) = 0.5516.

8 Ou et al. Table Numerical results for d 1 t = 0 t = 1 t 1 [WD(d 1 (γ ))] LB(t 1,t ) γ [WD(d 1 (γ ))] LB(t 1,t ) γ 0 0.690 0.5806 (00000) 0.6058 0.561 (0 000) 1 0.6006 0.5650 (10000) 0.601 0.5550 (0 0101) 0.61909 0.5591 (11000) 0.61981 0.559 (1 010) 0.6095 0.55 (11010) 0.6089 0.558 (0 1111) 0.6111 0.5597 (11110) 0.61905 0.5691 (1 111) Downloaded by [0.9.78.106] at 16:7 0 April 01 Example.. Consider a design d U(18; 7 ), given below, with n = 18,s 1 = 1,s = 7. 000 000 000 111 111 111 000 111 000 111 01 01 01 01 01 01 01 01 10 01 10 01 d =. 01 10 01 01 01 10 01 10 01 10 01 01 01 01 10 10 01 01 01 01 01 01 10 10 Table 5 shows the numerical results of d.from(),wehavelb(18, 1, 7) = 1.1. Table 5 Numerical results for d t 1 = 0 t 1 = 1 t [WD(d (γ ))] LB(t 1,t ) γ [WD(d (γ ))] LB(t 1,t ) γ 0 1.0 1.870 (0 0 000000) 1.011 1.0 (1 0 000000) 1 1.00 1.86 (0 0 00000) 1.89 1.817 (1 1 000000) 1.895 1.86 (0 0 0000) 1.86 1.550 (1 1 00000) 1.86 1.579 (0 0 01001) 1.86 1.98 (1 1 000011) 1.807 1.1 (0 000) 1.89 1.5 (1 1 000) 5 1.88 1.7 (0 10110) 1.86 1.77 (1 1 00) 6 1.88 1.68 (0 011) 1.86 1.8 (1 1 0) 7 1.807 1.6 (0 1 111) 1.86 1.651 (1 1 11)

Lower Bound of Discrepancy on Combined Designs 8 Table 6 Numerical results for d t [WD(d (γ ))] LB(t 1,t ) γ 0 0.1867 0.1867 (0000) 1 0.1810 0.17016 (1000) 0.1810 0.160 (1100) 0.1810 0.167 (1110) 0.1810 0.16661 (1111) Example.. Consider a design d U(9; ), given below, with n = 9,s 1 = 0,s =. Downloaded by [0.9.78.106] at 16:7 0 April 01 000 111 01 01 01 d =. 01 10 01 01 10 10 Table 6 shows the numerical results of d.from(),wehavelb(9, 0, ) = 0.1611. One can find that the lower bound LB(t 1,t ) given in Theorem. is tight when t = 0. Example. Consider a design d U(1; 11 11 ), given below, with n = 1,s 1 = 11,s = 11. 11011100010 00001001111 01101110001 000101100 10110111000 0010000 01011011100 0000000 00101101110 11010110 00010110111 1101100 d =. 10001011011 110101 11000101101 11010111 11100010110 0111 01110001011 111001 10111000101 10111001 00000000000 1100110 All Hamming distances between any distinct pair of runs of d are equal to 1. Table 7 shows the partial numerical results of d only for t = 0 in order to save space. From (), we have LB(1, 11, 11) = 89.. One can find that the lower bounds LB(t 1,t )given in Theorem. are tight when (t 1,t ) = (0, 0) and (t 1,t ) = (11, 0).

8 Ou et al. Table 7 Partial numerical results for d t = 0 t 1 [WD(d (γ ))] LB(t 1,t ) γ Downloaded by [0.9.78.106] at 16:7 0 April 01 0 699.77 699.77 (0 0 0 0000000000000000000) 1 658.0 65.1119 (0 1 0 0000000000000000000) 65.015 615.51 (1 1 0 0000000000000000000) 598.090 58.650 (0 1 1 1000000000000000000) 576.1515 560.170 (0 1 1 0110000000000000000) 5 558.61 50.8005 (0 1 1 0010100100000000000) 6 5.56 55.7018 (1 1 0 1011100000000000000) 7 51.10 51.1798 (1 0 0 1111100100000000000) 8 50.801 505.6561 (1 1 1 0111100100000000000) 9 510.85 99.690 (1 0 1 1111011100000000000) 10 50.086 95.757 (1 1 1 1111111000000000000) 11 9.67 9.67 (1 1 1 1111111100000000000) Acknowledgments The authors would like to thank the referee and Chief Editor for their valuable comments and suggestions that lead to improve the presentation of this article. Funding This work was partially supported by the National Natural Science Foundation of China (Grant Nos. 117117, 1101177), SRFDP (No. 0090111000), China Postdoctoral Science Foundation (No. 01M51716) and the Research Project Funded by Jishou University (No. jsdxxcfxbskyxm0111). References Box, G. E. P., Hunter, W. G., Hunter, J. S. (1978). Statistics for Experiments. New York: John Wiley and Sons. Chatterjee, K., Fang, K. T., Qin, H. (006). A lower bound for the centered L -discrepancy on asymmetric factorials and its application. Metrika 6: 55. Fang, K. T., Lin, D. K. J., Qin, H. (00). A note on optimal foldover design. Statist. Probab. Lett. 6:5 50. Fang, K. T., Lu, X., Winker, P. (00). Lower bounds for centered and wrap-around L -discrepancies and construction of uniform design by threshold accepting. J. Complex. 19:69 711. Fang, K. T., Tang, Y., Yin, J. X. (005). Lower bounds for wrap-around L -discrepancy and constructions of symmetrical uniform designs. J. Complex. 1:757 771. Hickernell, F. J. (1998). A generalized discrepancy and quadrature error bound. Math. Computat. 67:99. Lei, Y. J., Ou, Z. J., Qin, H., Zou, N. (01). A note on lower bound of centered L -discrepancy on combined designs. Acta Mathematica Sinica. 8():79 800. Lei, Y. J., Qin, H., Zou, N. (010). Some lower bounds of centered L -discrepancy on foldover designs. Acta Math Sci. 0A(6):1555 1561.

Lower Bound of Discrepancy on Combined Designs 85 Downloaded by [0.9.78.106] at 16:7 0 April 01 Li, F., Jacroux, M. (007). Optimal foldover plans for blocked m k fractional factorial designs. J. Statist. Plann. Infer. 17:9 5. Li, H., Mee, R. W. (00). Better foldover fractions for resolution III k p designs. Technometrics :78 8. Li, W., Lin, D. K. J. (00). Optimal foldover plans for two-level fractional factorial designs. Technometrics 5:1 19. Li, W., Lin, D. K. J., Ye, K. Q. (00). Optimal foldover plans for non-regular orthogonal designs. Technometrics 5:7 51. Montgomery, D. C., Runger, G. C. (1996). Foldover of k p resolution IV experimental designs. J. Qual. Technol. 8:6 50. Ou, Z. J., Chatterjee, K., Qin, H. (011). Lower bounds of various discrepancies on combined designs. Metrika 7:109 119. Qin, H., Chatterjee, K. (009). Lower bounds for the uniformity pattern of asymmetric fractional factorials. Commun. Statist. Theor. Meth. 8:18 19. Qin, H., Chatterjee, K., Ou, Z. J. (01). A lower bound for the centered L -discrepancy on combined designs under the asymmetric factorials. Statistics. 7:99 100. Wang, B., Robert, G. M., John, F. B. (010). A note on the selection of optimal foldover plans for 16- and -run fractional factorial designs. J. Statist. Plann. Infer. 10:197 1500. Ye, K. Q., Li, W. (00). Some properties of blocked and unblocked foldover of k p designs. Statistica Sinica 1:0 08. Zhang, A. J., Fang, K. T., Li, R., Sudjianto, A. (005). Majorization framework for balanced lattice designs. Ann. Statist. :87 85. Zhou, Y. D., Ning, J. H. (008). Lower bounds of the wrap-around L -discrepancy and relationships between MLHD and uniform design with a large size. J. Statist. Plann. Infer. 18:0 9.