Families of Hadamard group divisible designs

By a family of designs we mean a set of designs whose parameters can be represented as functions of an auxiliary variable t where the design will exist for infinitely many values of t. The best known family is probably the family of finite projective planes with υ = b = t² + t + 1, r = k = t + 1, and λ = 1. In some instances, notably coding theory, the existence of families is essential to provide the degree of precision required which can well vary from one coding problem to another. A natural vehicle for developing binary codes is the class of Hadamard matrices. Bush (1977) introduced the idea of constructing semi-regular designs using Hadamard matrices whereas the present study is concerned mostly with construction of regular designs using Hadamard matrices. While codes constructed from these designs are not optimal in the usual sense, it is possible that they may still have substantial value since, with different values of λ₁ and λ₂, there are different error correcting capabilities.     

On the orthogonal designs of order 40

Using a new method we construct all 17 remaining (unresolved for over 20 years) full orthogonal designs of order 40 in three variables. This implies that all full orthogonal designs OD(2^t5;x,y, 2^t5−x−y) exist for all t⩾3. The last two remaining orthogonal designs of order 40 in 2 variables are obtained as a special case of two of these designs.     

Efficient crossover designs in the presence of interactions between direct and carry-over treatment effects

Crossover designs, or repeated measurements designs, are used for experiments in which t treatments are applied to each of n experimental units successively over p time periods. Such experiments are widely used in areas such as clinical trials, experimental psychology and agricultural field trials. In addition to the direct effect on the response of the treatment in the period of application, there is also the possible presence of a residual, or carry-over, effect of a treatment from one or more previous periods. We use a model in which the residual effect from a treatment depends upon the treatment applied in the succeeding period; that is, a model which includes interactions between the treatment direct and residual effects. We assume that residual effects do not persist further than one succeeding period.A particular class of strongly balanced repeated measurements designs with n=t² units and which are uniform on the periods is examined. A lower bound for the A-efficiency of the designs for estimating the direct effects is derived and it is shown that such designs are highly efficient for any number of periods p=2,…,2t.     

Optimal designs for approximating the path of a stochastic process

We consider a centered stochastic process {X(t):t ∈ T} with known and continuous covariance function. On the basis of observations X(t₁), …, X(t_n) we approximate the whole path by orthogonal projection and measure the performance of the chosen design d = (t₁, …, t_n)′ by the corresponding mean squared L²-distance. For covariance functions on T² = [0, 1]², which satisfy a generalized Sacks-Ylvisaker regularity condition of order zero, we construct asymptotically optimal sequences of designs. Moreover, we characterize the achievement of a lower error bound, given by Micchelli and Wahba (1981), and study the question of whether this bound can be attained.     

Construction of orthogonal Latin hypercube designs with flexible run sizes

Latin hypercube designs (LHDs) have recently found wide applications in computer experiments. A number of methods have been proposed to construct LHDs with orthogonality among main-effects. When second-order effects are present, it is desirable that an orthogonal LHD satisfies the property that the sum of elementwise products of any three columns (whether distinct or not) is 0. The orthogonal LHDs constructed by Ye (1998), Cioppa and Lucas (2007), Sun et al. (2009) and Georgiou (2009) all have this property. However, the run size n of any design in the former three references must be a power of two (n=2^c) or a power of two plus one (n=2^c+1), which is a rather severe restriction. In this paper, we construct orthogonal LHDs with more flexible run sizes which also have the property that the sum of elementwise product of any three columns is 0. Further, we compare the proposed designs with some existing orthogonal LHDs, and prove that any orthogonal LHD with this property, including the proposed orthogonal LHD, is optimal in the sense of having the minimum values of ave(|t|), t_max, ave(|q|) and q_max.     

On computer generation of balanced arrays

This paper provides an algebraic (and hence computing) procedure for generation of balanced arrays having two symbols, m rows, specified minimum and maximum column weights, arbitrary strength t≦m, and index set parameters μ^t₁, μ^t₂,…, μ^t_t. μ^t₀ is unspecified, and calculated as part of the algorithm, although the procedure for specifying it is straightforward and can be used if desired. Array generation is herein reduced to finding integral solutions to a linear programming problem. It is shown that the integral solutions of the system of equations comprise all balanced arrays with the given set of parameters.A computing algorithm is provided which constructs the system of equations to be solved; it has been interfaced with a standard linear programming package to provide some preliminary results.Additional algorithms whose development should result in substantial decreases in computing costs are discussed.     

Optimal Designs for Approximating a Stochastic Process with Respect to a Minimax Criterion

We study the problem of approximating a stochastic process Y = {Y(t: t ∈ T} with known and continuous covariance function R on the basis of finitely many observations Y(t ₁,), …, Y(t _n ). Dependent on the knowledge about the mean function, we use different approximations ? and measure their performance by the corresponding maximum mean squared error sub _t∈T E(Y(t) ? ?(t))². For a compact T ? ? ^p we prove sufficient conditions for the existence of optimal designs. For the class of covariance functions on T ² = [0, 1]² which satisfy generalized Sacks/Ylvisaker regularity conditions of order zero or are of product type, we construct sequences of designs for which the proposed approximations perform asymptotically optimal.     

Nested Hadamard difference sets

A Hadamard difference set (HDS) has the parameters (4N², 2N² − N, N² − N). In the abelian case it is equivalent to a perfect binary array, which is a multidimensional matrix with elements ±1 such that all out-of-phase periodic autocorrelation coefficients are zero. We show that if a group of the form H × Z²_p^r contains a (hp^2r, √hp^r(2√hp^r − 1), √hp^r(√hp^r − 1)) HDS (HDS), p a prime not dividing |H| = h and p^j ≡ −1 (mod exp(H)) for some j, then H × Z²_p^t has a (hp^2t, √hp^t(2√hp^t − 1), √hp^t(√hp^t − 1)) HDS for every 0⩽t⩽r. Thus, if these families do not exist, we simply need to show that H × Z²_p does not support a HDS. We give two examples of families that are ruled out by this procedure.     

One-factor-at-a-time designs of resolution V

Two classes of designs of resolution V are constructed using one-factor-at-a-time techniques. They facilitate sequential learning and are more economical in run size than regular 2_V^k−p designs. Comparisons of D efficiency with other designs are given to assess the suitability of the proposed designs. Their D_s efficiencies can be dramatically improved with the addition of a few runs.     

Minimum aberration majorization in non-isomorphic saturated designs

In this paper we propose a new criterion, minimum aberration majorization, for comparing non-isomorphic saturated designs. This criterion is based on the generalized word-length pattern proposed by Ma and Fang (Metrika 53 (2001) 85) and Xu and Wu (Ann. Statist. 29 (2001) 1066) and majorization theory. The criterion has been successfully applied to check non-isomorphism and rank order saturated designs. Examples are given through five non-isomorphic L₁₆(2¹⁵) designs and two L₂₇(3¹³) designs.     

Some determinant optimality aspects of main effect fold-over designs

Saha and Mohanty (1970) presented a main effect fold-over design consisting of 14 treatment combinations of the 24×33 factorial, which had the nice property of being even balanced. Calling this design DSM, this paper establishes the following specific results: (i) DSM is not d-optimal in the subclass Δe of all 14 point even balanced main effect fold-over designs of the 24×33 factorial; (ii) DSM is not d-optimal in the subclass $\text{Δ}{}^1\text{?Δ}{}^{\mathrm{e}}$ of all 14 point even and odd balanced main effect fold-over designs of the 24×33 factorial; (iii) DSM is even optimal in $\text{Δ}{}^1$ and Δe. In addition to these results two 14 point designs in $\text{Δ}{}^1$ are presented which are d-optimal and via a counter example it is shown that these designs are not odd optimal. Finally, several general matrix algebra results are given which should be useful in resolving d-optimality problems of fold-over designs of the kn11×kn22 factorial.     

Note on Draper and Lin's article “capacity considerations for two-level fractional factorial designs”

In Table 4 of Draper and Lin (1990) the authors present eleven 2^{k − p}_R designs for which doubt is expressed as to whether they are saturated. This note removes doubt on nine of these designs indicating that they are indeed saturated. Some observations are made in terms of the resolution of a design for the remaining two designs among the 11. The conclusions are derived by properly interpreting an updated table of binary linear codes of Verhoeff (1987).     

Balanced fractional factorial designs of resolution 2l+1 for interesting effects orthogonal to some nuisance parameters: 2m1+m2 series

For 2^m₁+m₂ factorial designs, this paper investigates balanced fractional 2^m₁ factorial designs of resolution 2l+1 with some nuisance parameters concerning the second factors. They are derivable from partially balanced arrays and further permit estimation of the effects up to the l-factor interactions concerning the first factors orthogonally to the nuisance parameters.     

Generalized Bhaskar Rao designs

Generalized Bhaskar Rao designs with non-zero elements from an abelian group G are constructed. In particular this paper shows that the necessary conditions are sufficient for the existence of generalized Bhaskar Rao designs with k=3 for the following groups: ?G? is odd, G=Z^r₂, and G=Z^r₂×H where 3? ?H? and r?1. It also constructs generalized Bhaskar Rao designs with υ=k, which is equivalent to υ rows of a generalized Hadamard matrix of order n where υ?n.     

Theory of factorial designs of the parallel flats type. I: the coefficient matrix

Let GF(s) be the finite field with s elements.(Thus, when s=3, the elements of GF(s) are 0, 1 and 2.)Let A(r×n), of rank r, and c_i(i=1,…,f), (r×1), be matrices over GF(s). (Thus, for n=4, r=2, f=2, we could have A=[¹¹¹⁰₀₁₂₁], c₁=[¹₀], c₂=[⁰₂].) Let T_i (i=1,…,f) be the flat in EG(n, s) consisting of the set of all the s^n?r solutions of the equations At=c_i, wheret′=(t₁,…,t_n) is a vector of variables.(Thus, EG(4, 3) consists of the 3⁴=81 points of the form (t₁,t₂,t₃,t₄), where t's take the values 0,1,2 (in GF(3)). The number of solutions of the equations At=c_i is s^n?r, where r=Rank(A), and the set of such solutions is said to form an (n?r)-flat, i.e. a flat of (n?r) dimensions. In our example, both T₁ and T₂ are 2-flats consisting of 3^4?2=9 points each. The flats T₁,T₂,…,T_f are said to be parallel since, clearly, no two of them can have a common point. In the example, the points of T₁ are (1000), (0011), (2022), (0102), (2110), (1121), (2201), (1212) and (0220). Also, T₂ consists of (0002), (2010), (1021), (2101), (1112), (0120), (1200), (0211) and (2222).) Let T be the fractional design for a sⁿ symmetric factorial experiment obtained by taking T₁,T₂,…,T_f together. (Thus, in the example, 3⁴=81 treatments of the 3⁴ factorial experiment correspond one-one with the points of EG(4,3), and T will be the design (i.e. a subset of the 81 treatments) consisting of the 18 points of T₁ and T₂ enumerated above.)In this paper, we lay the foundation of the general theory of such ‘parallel’ types of designs. We define certain functions of A called the alias component matrices, and use these to partition the coefficient matrix X (n×v), occuring in the corresponding linear model, into components X._j(j=0,1,…,g), such that the information matrix X is the direct sum of the X′._jX._j. Here, v is the total number of parameters, which consist of (possibly μ), and a (general) set of (geometric) factorial effects (each carrying (s?1) degrees of freedom as usual). For j≠0, we show that the spectrum of X′._jX._j does not change if we change (in a certain important way) the usual definition of the effects. Assuming that such change has been adopted, we consider the partition of the X._j into the X_ij (i=1,…,f). Furthermore, the X_ij are in turn partitioned into smaller matrices (which we shall here call the) X_ijh. We show that each X_ijh can be factored into a product of 3 matrices J, ζ (not depending on i,j, and h) and Q(j,h,i)where both the Kronecker and ordinary product are used. We introduce a ring R using the additive groups of the rational field and GF(s), and show that the Q(j,h,i) belong to a ring isomorphic to R. When s is a prime number, we show that R is the cyclotomic field. Finally, we show that the study of the X._j and X′._jX._j can be done in a much simpler manner, in terms of certain relatively small sized matrices over R.     

Comparison of rotatable designs for regression on balls,I (quadratic)

Designs for quadratic and cubic regression are considered when the possible choices of the controlable variable are points x=( x₁,x₂,…,x_q) in the q-dimensional. Full of radius R, B_q(R) ={x:Σ⁴_ix²_i?R²}. The designs that are optimum among rotatable designs with respect to the D-, A-, and E-optimality criteria are compared in their performance relative to these and other criteria, including extrapolation. Additionally, the performance of a design optimum for one value of R, when it is implemented for a different value of R, is investigated. Some of the results are developed algebraically; others, numerically. For example, in quadratic regression the A-optimum design appears to be fairly robust in its efficiency, under variation of criterion.     

A-optimal partially balanced fractional 2m1+m2 factorial designs of resolution V,with 4≤m1>+m2≤6

By use of the algebraic structure, we obtain an explicit expression for the characteristic polynomial of the information matrix of a partially balanced fractional 2^m₁+m₂ factorial design of resolution V derived from a partially balanced array. For 4≤m₁+m₂≤6, A-optimal designs considered here are also presented for reasonable number of assemblies.     

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.     

Construction of optimal supersaturated design with large number of levels

Inspired by the ideas of column and row juxtaposition in Liu and Lin (2009) and level transformation in Yamada and Lin (1999), this paper presents a new method for constructing optimal supersaturated designs (SSDs). This method provides a convenient way to construct mixed-level designs with relatively large numbers of levels, avoiding the blind search and numerous calculations by computers. The goodness of the resulting SSDs is judged by the χ² (Yamada and Lin, 1999 and Yamada and Matsui, 2002) and J₂ (Xu, 2002) criteria. Some nice properties of the new designs are also provided.     

Weighted A-optimality for fractional 2m factorial designs of resolution V

A weighted A-optimality (WA-optimality) criterion is discussed for selecting a fractional 2^m factorial design of resolution V. A WA-optimality criterion having one weight may be considered for designs. It is shown that designs derived from orthogonal arrays are WA-optimal for any weight. From a WA-optimal design, a procedure for finding WA-optimal designs for various weights is given. WA-optimal balanced designs are presented for 4 ⩽ m ⩽ 7 and for the values of n assemblies in certain ranges. It is pointed out that designs for m = 7 and for n = 41, 42 given in Chopra and Srivastava (1973a) or in the corrected paper by Chopra et al. (1986), are not A-optimal.