Results for two-level fractional factorial designs of resolution IV or more

The main theorem of this paper shows that foldover designs are the only (regular or nonregular) two-level factorial designs of resolution IV (strength 3) or more for n   runs and n/3?m?n/2$n/3?m?n/2$ factors. This theorem is a generalization of a coding theory result of Davydov and Tombak [1990. Quasiperfect linear binary codes with distance 4 and complete caps in projective geometry. Problems Inform. Transmission 25, 265–275] which, under translation, effectively states that foldover (or even) designs are the only regular two-level factorial designs of resolution IV or more for n   runs and 5n/16?m?n/2$5n/16?m?n/2$ factors. This paper also contains other theorems including an alternative proof of Davydov and Tombak's result.     

Classification of orthogonal arrays by integer programming

The problem of classifying all isomorphism classes of OA(N,k,s,t)$\mathrm{OA}(N,k,s,t)$'s is shown to be equivalent to finding all isomorphism classes of non-negative integer solutions to a system of linear equations under the symmetry group of the system of equations. A branch-and-cut algorithm developed by Margot [2002. Pruning by isomorphism in branch-and-cut. Math. Programming Ser. A 94, 71–90; 2003a. Exploiting orbits in symmetric ILP. Math. Programming Ser. B 98, 3–21; 2003b. Small covering designs by branch-and-cut. Math. Programming Ser. B 94, 207–220; 2007. Symmetric ILP: coloring and small integers. Discrete Optim., 4, 40–62] for solving integer programming problems with large symmetry groups is used to find all non-isomorphic OA(24,7,2,2)$\mathrm{OA}(24,7,2,2)$'s, OA(24,k,2,3)$\mathrm{OA}(24,k,2,3)$'s for 6?k?11$6?k?11$, OA(32,k,2,3)$\mathrm{OA}(32,k,2,3)$'s for 6?k?11$6?k?11$, OA(40,k,2,3)$\mathrm{OA}(40,k,2,3)$'s for 6?k?10$6?k?10$, OA(48,k,2,3)$\mathrm{OA}(48,k,2,3)$'s for 6?k?8$6?k?8$, OA(56,k,2,3)$\mathrm{OA}(56,k,2,3)$'s, OA(80,k,2,4)$\mathrm{OA}(80,k,2,4)$'s, OA(112,k,2,4)$\mathrm{OA}(112,k,2,4)$'s, for k=6,7$k=6,7$, OA(64,k,2,4)$\mathrm{OA}(64,k,2,4)$'s, OA(96,k,2,4)$\mathrm{OA}(96,k,2,4)$'s for k=7,8$k=7,8$, and OA(144,k,2,4)$\mathrm{OA}(144,k,2,4)$'s for k=8,9$k=8,9$. Further applications to classifying covering arrays with the minimum number of runs and packing arrays with the maximum number of runs are presented.     

Super-simple cyclic designs with small values

Super-simple cyclic designs are useful on constructing codes and designs such as superimposed codes, perfect hash families and optical orthogonal codes with index two. In this paper, we show that there exists a super-simple cyclic (v,4,λ)$(v,4,\lambda )$ for 7?v?41$7?v?41$ and all admissible λ$\lambda$ with two definite exceptions of (v,λ)=(9,3),(13,5)$(v,\lambda )=(9,3),(13,5)$ and one possible exception of (v,λ)=(39,18)$(v,\lambda )=(39,18)$. Some useful algorithms are explained for computer search and new designs are displayed.     

Tools for constructing optimal two-level factorial designs for a linear model containing main effects and one two-factor interaction

In Hedayat and Pesotan [1992, Two-level factorial designs for main effects and selected two-factor interactions. Statist. Sinica 2, 453–464.] the concepts of a g(n,e)$g(n,e)$-design and a g(n,e)$g(n,e)$-matrix are introduced to study designs of n$n$ factor two-level experiments which can unbiasedly estimate the mean, the n$n$ main effects and e$e$ specified two-factor interactions appearing in an orthogonal polynomial model and it is observed that the construction of a g-design is equivalent to the construction of a g  -matrix. This paper deals with the construction of D-optimal g(n,1)$g(n,1)$-matrices. A standard form for a g(n,1)$g(n,1)$-matrix is introduced and some lower and upper bounds on the absolute determinant value of a D-optimal g(n,1)$g(n,1)$-matrix in the class of all g(n,1)$g(n,1)$-matrices are obtained and an approach to construct D-optimal g(n,1)$g(n,1)$-matrices is given for 2?n?8$2?n?8$. For two specific subclasses, namely a certain class of g(n,1)$g(n,1)$-matrices within the class of g(n,1)$g(n,1)$-matrices of index one and the class C(H)$C(H)$ of g(8t+2,1)$g(8t+2,1)$-matrices constructed from a normalized Hadamard matrix H   of order 8t+4(t?1)$8t+4(t?1)$ two techniques for the construction of the restricted D-optimal matrices are given.     

An isomorphism check for two-level fractional factorial designs

Two fractional factorial designs are isomorphic if one can be obtained from the other by reordering the treatment combinations, relabelling the factor levels and relabelling the factors. By defining a word-pattern matrix, we are able to create a new isomorphism check which is much faster than existing checks for certain situations. We combine this with a new, extremely fast, sufficient condition for non-isomorphism to avoid checking certain cases. We then create a faster search algorithm by combining the Bingham and Sitter [1999. Minimum aberration fractional factorial split-plot designs. Technometrics 41, 62–70] search algorithm, the isomorphism check algorithm of Clark and Dean [2001. Equivalence of fractional factorial designs. Statist. Sinica 11, 537–547] with our proposed isomorphism check. The algorithm is used to extend the known set of existing non-isomorphic 128-run two-level regular designs with resolution ?4$?4$ to situations with 12, 13, 14, 15 and 16 factors, 256- and 512-run designs with resolution ?5$?5$ and ?17$?17$ factors and 1024-run even designs with resolution ?6$?6$ and ?18$?18$ factors.     

Optimal block designs for three treatments when observations are correlated

E$E$-optimal designs for comparing three treatments in blocks of size three are identified, where intrablock observations are correlated according to a first order autoregressive error process with parameter ρ∈(0,1)$\rho \in (0,1)$. For number of blocks b   of the form b=3n+1$b=3n+1$, there are two distinct optimal designs depending on the value of ρ$\rho$, with the best design being unequally replicated for large ρ$\rho$. For other values of b$b$, binary, equireplicate designs with specified within-block assignment patterns are best. In many cases, the stronger majorization optimality is established.     

Improved WLP and GWP lower bounds based on exact integer programming

By using exact integer programming (IP) (integer programming in infinite precision) bounds on the word-length patterns (WLPs) and generalized word-length patterns (GWPs) for fractional factorial designs are improved. In the literature, bounds on WLPs are formulated as linear programming (LP) problems. Although the solutions to such problems must be integral, the optimization is performed without the integrality constraints. Two examples of this approach are bounds on the number of words of length four for resolution IV regular designs, and a lower bound for the GWP of two-level orthogonal arrays. We reformulate these optimization problems as IP problems with additional valid constraints in the literature and improve the bounds in many cases. We compare the improved bound to the enumeration results in the literature to find many cases for which our bounds are achieved. By using the constraints in our integer programs we prove that f(16λ,2,4)?9$f(16\lambda ,2,4)?9$ if λ$\lambda$ is odd where f(2^tλ,2,t)$f(2^t\lambda ,2,t)$ is the maximum n   for which an OA(N,n,2,t)$\mathrm{OA}(N,n,2,t)$ exists. We also present a theorem for constructing GMA OA(N,N/2-u,2,3)$\mathrm{OA}(N,N/2-u,2,3)$ for u=1,…,5$u=1,\dots ,5$.     

Comparing several treatments with a control

Consider the partially balanced one-way layout for comparing k   treatments _μi,${\mu }_i,$1?i?k,$1?i?k,$ with a control _μ0${\mu }_0$. We propose a new test which is similar to the test statistics of Marcus [1976. The powers of some tests of the equality of normal means against an ordered alternative. Biometrika 63, 177–183]. By simulation we find that the proposed test has a good power performance when compared with other tests. Moreover, it can produce confidence intervals for _μi-_μ0,1?i?k.${\mu }_i-{\mu }_0,1?i?k.$     

A counterexample to Beder's conjectures about Hadamard matrices

In this note we provide a counterexample which resolves conjectures about Hadamard matrices made in this journal. Beder [1998. Conjectures about Hadamard matrices. Journal of Statistical Planning and Inference 72, 7–14] conjectured that if H$\mathit{H}$ is a maximal m×n$m\times n$ row-Hadamard matrix then m is a multiple of 4; and that if n   is a power of 2 then every row-Hadamard matrix can be extended to a Hadamard matrix. Using binary integer programming we obtain a maximal 13×32$13\times 32$ row-Hadamard matrix, which disproves both conjectures. Additionally for n being a multiple of 4 up to 64, we tabulate values of m   for which we have found a maximal row-Hadamard matrix. Based on the tabulated results we conjecture that a m×n$m\times n$ row-Hadamard matrix with m?n-7$m?n-7$ can be extended to a Hadamard matrix.     

Limit laws for the cumulative number of ties for the maximum in a random sequence

Let {_Xn,n?1}$\{X_n,n?1\}$ be a sequence of independent identically distributed random variables, taking nonnegative integer values. An observation _Xn$X_n$ is a tie for the maximum if _Xn=max{_X1,…,_Xn-1}$X_n=\max \{X_1,\dots ,X_{n-1}\}$. In this paper, we obtain weak and strong laws of large numbers and central limit theorems for the cumulative number of ties for the maximum among the first n$n$ observations.     

A sequential design for a clinical trial with a linear prognostic factor

We study a randomized adaptive design to assign one of the L$L$ treatments to patients who arrive sequentially by means of an urn model. At each stage n$n$, a reward is distributed between treatments. The treatment applied is rewarded according to its response, 0?_Yn?1$0?Y_n?1$, and 1-_Yn$1-Y_n$ is distributed among the other treatments according to their performance until stage n-1$n-1$. Patients can be classified in K+1$K+1$ levels and we assume that the effect of this level in the response to the treatments is linear. We study the asymptotic behavior of the design when the ordinary least square estimators are used as a measure of performance until stage n-1$n-1$.     

General minimum lower order confounding designs: An overview and a construction theory

For fractional factorial (FF) designs, Zhang et al. (2008) introduced a new pattern for assessing regular designs, called aliased effect-number pattern (AENP), and based on the AENP, proposed a general minimum lower order confounding (denoted by GMC for short) criterion for selecting design. In this paper, we first have an overview of the existing optimality criteria of FF designs, and then propose a construction theory for 2^n−m$2^{n-m}$ GMC designs with 33N/128≤n≤5N/16$33N/128\le n\le 5N/16$, where N=2^n−m$N=2^{n-m}$ is the run size and n is the number of factors, for all N's and n  's, via the doubling theory and SOS resolution IV designs. The doubling theory is extended with a new approach. By introducing a notion of rechanged (RC) Yates order for the regular saturated design, the construction result turns out to be quite transparent: every GMC 2^n−m$2^{n-m}$ design simply consists of the last n columns of the saturated design with a specific RC Yates order. This can be very conveniently applied in practice.     

On the E-optimality of complete block designs under a mixed interference model

In experiments in which the response to a treatment can be affected by other treatments, the interference model with neighbor effects is usually used. It is known that circular neighbor balanced designs (CNBDs) are universally optimal under such a model if the neighbor effects are fixed (Druilhet, 1999) or random (4 and 7). However, such designs cannot exist for every combination of design parameters. In the class of block designs with the same number of treatments as experimental units per block, a CNBD cannot exist if the number of blocks, b  , is equal to p(t−1)±1$p(t-1)\pm 1$, where p is a positive integer and t is the number of treatments. Filipiak et al. (2008) gave the structure of the left-neighboring matrix of E-optimal complete block designs with p  =1 under the model with fixed neighbor effects. The purpose of this paper is to generalize E-optimality results for designs with p∈N$p\in \mathbb{N}$ assuming random neighbor effects.     

Model-robust supersaturated and partially supersaturated designs

Supersaturated designs are an increasingly popular tool for screening factors in the presence of effect sparsity. The advantage of this class of designs over resolution III factorial designs or Plackett–Burman designs is that n, the number of runs, can be substantially smaller than the number of factors, m. A limitation associated with most supersaturated designs produced thus far is that the capability of these designs for estimating g active effects has not been discussed. In addition to exploring this capability, we develop a new class of model-robust supersaturated designs that, for a given n and m, maximizes the number g   of active effects that can be estimated simultaneously. The capabilities of model-robust supersaturated designs for model discrimination are assessed using a model-discrimination criterion, the subspace angle. Finally, we introduce the class of partially supersaturated designs, intended for use when we require a specific subset of _m1$m_1$ core factors to be estimable, and the sparsity of effects principle applies to the remaining (m-_m1$m-m_1$) factors.     

Generalizing Clatworthy group divisible designs

In this paper a neat construction is provided for three new families of group divisible designs that generalize some designs from Clatworthy's table of the only 11 designs with two associate classes that have block size four, three groups, and replication numbers at most 10. In each case (namely, _λ1=4${\lambda }_1=4$ and _λ2=5${\lambda }_2=5$, _λ1=4${\lambda }_1=4$ and _λ2=2${\lambda }_2=2$, and _λ1=8${\lambda }_1=8$ and _λ2=4${\lambda }_2=4$), we have proved that the necessary conditions found are also sufficient for the existence of such GDD's with block size four and three groups, with one possible exception.