A survey and evaluation of methods for determination of combinatorial equivalence of factorial designs

Equivalent factorial designs have identical statistical properties for estimation of factorial contrasts and for model fitting. Non-equivalent designs, however, may have the same statistical properties under one particular model but different properties under a different model. In this paper, we describe known methods for the determination of equivalence or non-equivalence of two-level factorial designs, whether they be regular factorial designs, non-regular orthogonal arrays, or have no particular structure. In addition, we evaluate a number of potential fast screening methods for detecting non-equivalence of designs. Although the paper concentrates mainly on symmetric designs with factors at two levels, we also evaluate methods of determining combinatorial equivalence and non-equivalence of three-level designs and indicate extensions to larger numbers of levels and to asymmetric designs.     

On equivalence of fractional factorial designs based on singular value decomposition

The singular value decomposition of a real matrix always exists and is essentially unique. Based on the singular value decomposition of the design matrices of two general 2-level fractional factorial designs, new necessary and sufficient conditions for the determination of combinatorial equivalence or non-equivalence of the corresponding designs are derived. Equivalent fractional factorial designs have identical statistical properties for estimation of factorial contrasts and for model fitting. Non-equivalent designs, however, may have the same statistical properties under one particular model but different properties under a different model. Results extend to designs with factors at larger number of levels.     

Approximate theory of multiway block designs

Optimality properties of multiway block designs are deduced from the general results of J. Kiefer's approximate-design theory. In the model with additive effects these optimality properties solely depend on the two-dimensional marginals of the designs. Uniform designs, and designs whose two-dimensional marginals are products of the one-dimensional marginals, are shown to be optimal. Approximate Youden designs are introduced for the case when the support sets of the two-dimensional marginals are prescribed in advance. They are optimal in a relatively small class of competing designs only.     

The analysis of mixed and random effects models for equireplicate variance-balanced block designs

Equirephcate variance-balanced block designs possess many desirable properties and include several important subclasses of designs, such as balanced incomplete block designs. However, while attention has been paid to the derivation of properties and the construction of these designs, methods of analyzing experiments using these designs have not been developed completely. Thus, we consider models appropriate to these designs, construct test statistics for basic hypotheses and derive the distributions for these statistics. Four additive models that include all of the possible combinations of fixed and random effects are examined. The analyses are generalizations of the analyses for balanced incomplete block designs, and include the intra-block, intra-treatment and inter-block analyses, Estimators of fixed treatment effects and their distributions also are given.     

On algorithmic construction of maximin distance designs

Maximin distance designs are useful for conducting expensive computer experiments. In this article, we compare some global optimization algorithms for constructing such designs. We also introduce several related space-filling designs, including nested maximin distance designs, sliced maximin distance designs, and general maximin distance designs with better projection properties. These designs possess more flexible structures than their analogs in the literature. Examples of these designs constructed by the algorithms are presented.     

Block-circulant matrices for constructing optimal Latin hypercube designs

Computer simulations are usually needed to study a complex physical process. In this paper, we propose new procedures for constructing orthogonal or low-correlation block-circulant Latin hypercube designs. The basic concept of these methods is to use vectors with a constant periodic autocorrelation function to obtain suitable block-circulant Latin hypercube designs. A general procedure for constructing orthogonal Latin hypercube designs with favorable properties and allowing run sizes being different from a power of 2 (or a power of 2 plus 1), is presented here for the first time. In addition, an expansion of the method is given for constructing Latin hypercube designs with low correlation. This expansion is useful when orthogonal Latin hypercube designs do not exist. The properties of the generated designs are further investigated. Some examples of the new designs, as generated by the proposed procedures, are tabulated. In addition, a brief comparison with the designs that appear in the literature is given.     

αn–Designs

This paper defines a broad class of resolvable incomplete block designs called α_n–designs, of which the original α–designs are a special case with n = 1. The statistical and mathematical properties of α–designs extend naturally to these n –dimensional designs. They are a flexible class of resolvable designs appropriate for use in factorial experiments, in constructing efficient t –latinized resolvable block designs, and for enhancing the existing class of α–designs for a single treatment factor.     

Response surface design evaluation and comparison

Designing an experiment to fit a response surface model typically involves selecting among several candidate designs. There are often many competing criteria that could be considered in selecting the design, and practitioners are typically forced to make trade-offs between these objectives when choosing the final design. Traditional alphabetic optimality criteria are often used in evaluating and comparing competing designs. These optimality criteria are single-number summaries for quality properties of the design such as the precision with which the model parameters are estimated or the uncertainty associated with prediction. Other important considerations include the robustness of the design to model misspecification and potential problems arising from spurious or missing data. Several qualitative and quantitative properties of good response surface designs are discussed, and some of their important trade-offs are considered. Graphical methods for evaluating design performance for several important response surface problems are discussed and we show how these techniques can be used to compare competing designs. These graphical methods are generally superior to the simplistic summaries of alphabetic optimality criteria. Several special cases are considered, including robust parameter designs, split-plot designs, mixture experiment designs, and designs for generalized linear models.     

Orthogonal Latin hypercube designs from generalized orthogonal designs

Latin hypercube designs is a class of experimental designs that is important when computer simulations are needed to study a physical process. In this paper, we proposed some general criteria for evaluating Latin hypercube designs through their alias matrices. Moreover, a general method is proposed for constructing orthogonal Latin hypercube designs. In particular, links between orthogonal designs (ODs), generalized orthogonal designs (GODs) and orthogonal Latin hypercube designs are established. The generated Latin hypercube designs have some favorable properties such as uniformity, orthogonality of the first and some second order terms, and optimality under the defined criteria.     

Orthogonal designs for computer experiments

Computer experiments are becoming increasingly popular surrogates for physical experiments in recent years. In this paper, we propose a general procedure for constructing orthogonal designs with many levels and distinct runs. These designs are suitable for computer experiments. The basic idea of the method is to fold-over suitable orthogonal matrices. The properties of the derived designs are studied and a brief comparison with the designs in the literature is given.     

Cost-cautious designs for confirmatory bioassay

Confirmatory bioassay experiments take place in late stages of the drug discovery process when a small number of compounds have to be compared with respect to their properties. As the cost of the observations may differ considerably, the design problem is well specified by the cost of compound used rather than by the number of observations. We show that cost-efficient designs can be constructed using useful properties of the minimum support designs. These designs are particularly suited for studies where the parameters of the model to be estimated are known with high accuracy prior to the experiment, although they prove to be robust against typical inaccuracies of these values. When the parameters of the model can only be specified with ranges of values or by a probability distribution, we use a Bayesian criterion of optimality to construct the required designs. Typically, the number of their support points depends on the prior knowledge for the model parameters. In all cases we recommend identifying a set of designs with good statistical properties but different potential costs to choose from.     

On Exact K-optimal Designs Minimizing the Condition Number

A new design criterion based on the condition number of an information matrix is proposed to construct optimal designs for linear models, and the resulting designs are called K-optimal designs. The relationship between exact and asymptotic K-optimal designs is derived. Since it is usually hard to find exact optimal designs analytically, we apply a simulated annealing algorithm to compute K-optimal design points on continuous design spaces. Specific issues are addressed to make the algorithm effective. Through exact designs, we can examine some properties of the K-optimal designs such as symmetry and the number of support points. Examples and results are given for polynomial regression models and linear models for fractional factorial experiments. In addition, K-optimal designs are compared with A-optimal and D-optimal designs for polynomial regression models, showing that K-optimal designs are quite similar to A-optimal designs.     

Regular graph designs

In this survey paper, measures of optimality of experimental designs are discussed from a combinatorial viewpoint. The main vehicle for this discussion is the regular graph design. This family of designs was introduced in the hope that they would be ‘good’ designs in cases where no balanced incomplete block design can exist. Various properties of these designs are discussed.     

Some methods for constructing efficiency-balanced block designs

Three methods are presented for constructing connected efficiency-balanced block designs from other block designs with the same properties. The resulting designs differ from the original ones in the number of blocks and/or in the number of experimental units and their arrangement, while the number of treatments remains unaltered. Some remarks on the proposed methods of construction refer also to variance-balanced block designs.     

Classification of efficient two-level fractional factorial designs of resolution IV or more

Two-level fractional factorial designs that are efficient in terms of aberration or other aliasing properties are classified into four types of designs of resolution IV or more: the half-fraction designs, the even designs, the five-column designs and the join designs. The designs are shown to have concise grid representations which provide simple interpretations of their aliasing structure. New efficient 128-run designs are presented and blocking of the designs is considered.     

Construction of effect-wise orthogonal factorial designs

Generalizing the concept of Kronecker products of designs, two distinct methods have been suggested for the construction of effect-wise orthogonal factorial designs. The methods described ensure desirable properties with respect to main effects, cover almost all cases of factorial designs and require, in most cases, a smaller number of replications than any of the existing methods.     

ADJUSTED ORTHOGONAL NESTED ROW-COLUMN DESIGNS

Adjusted orthogonality in nested row-column designs is defined and a sufficient condition established for its existence. It is shown that the properties of an adjusted orthogonal nested row-column design are directly related to those of its separate row and column component designs. A method for constructing efficient adjusted orthogonal designs involving a single replicate of every treatment in each of two blocks is given.     

A generating algorithm for linear trend-free and nearly linear trend-free block designs

The theory and properties of trend-free (TF) and nearly trend-free (NTF) block designs are wel1 developed. Applications have been hampered because a methodology for design construction has not been available.

This article begins with a short review of concepts and properties of TF and NTF block designs. The major contribution is provision of an algorithm for the construction of linear and nearly linear TF block designs. The algorithm is incorporated in a computer program in FORTRAN 77 provided in an appendix for the IBM PC or compatible microcomputer, a program adaptable also to other computers. Three sets of block designs generated by the program are given as examples.

A numerical example of analysis of a linear trend-free balanced incomplete block design is provided.     

Minimax Designs for the Stability of Slope Estimation on Second-order Response Surfaces

In this paper, designs for the stability of the slope estimation on a second-order response surface are considered. Minimization of the point dispersion measure, which is maximized over all points in the region of interest is taken as the optimality criterion, and the minimax properties in some class of designs are derived in spherical and cubic regions of interest. We study the efficiencies of the minimax designs relative to other optimal designs with various criteria.     

Optimal orthogonal-array-based latin hypercubes

The use of optimal orthogonal array latin hypercube designs is proposed. Orthogonal arrays were proposed for constructing latin hypercube designs by Tang (1993). Such designs generally have better space filling properties than random latin hypercube designs. Even so, these designs do not necessarily fill the space particularly well. As a result, we consider orthogonal-array-based latin hypercube designs that try to achieve optimality in some sense. Optimization is performed by adapting strategies found in Morris & Mitchell (1995) and Ye et al. (2000). The strategies here search only orthogonal-array-based latin hypercube designs and, as a result, optimal designs are found in a more efficient fashion. The designs found are in general agreement with existing optimal designs reported elsewhere.