Description and evaluation of a nested cube experimental design

An incomplete factorial design based on an extension of the Fawiliar 2^kfactorial called a nested cube is proposed for use in response surface investigationso The simplicity and general efficiency of the nested cube suggest its suitability to many areas of research, especially that repeated at many locations orconductea over a long period. Comparisons to potentially competing designs are provided for bias in response estination due to fitting an ioappropriate model and for profiles of variance. merits of the nested cube are (1) a level of relative bias and variance judged to be favorable though not optimal, (2) an ability to utilize a minimum blag estimator not available to competing designs, and (3) a simplicity associated with use of equal spacing

and nearly equal replication on the margin for each factor level.     

A survey of nested designs

A nest with parameters (r,k,λ)→(r′,k′,λ′) is a BIBD on (b,v,r,k,λ) where each block has a distinguished sublock of cardinality k, the sublocks forming a (b,v,r,k,λ)-design.These designs are ‘nested’ in the sense of W.T. Federer (1972), who recommended the use of these designs for the sequential addition of periods in marketing experiments in order to retain Youden design properties as rows are added. Note that for a Youden design, the b columns and v treatments are in an SBIBD arrangement with parameters v=b, k=r, and λ.     

Generalized staggered nested designs for variance components estimation

Staggered nested experimental designs are the most popular class of unbalanced nested designs in practical fields. The most important features of the staggered nested design are that it has a very simple open-ended structure and each sum of squares in the analysis of variance has almost the same degrees of freedom. Based on the features, a class of unbalanced nested designs that is a generalization of the staggered nested design is proposed in this paper. Formulae for the estimation of variance components and their sums are provided. Comparing the variances of the estimators to the staggered nested designs, it is found that some of the generalized staggered nested designs are more efficient than the traditional staggered nested design in estimating some of the variance components and their sums. An example is provided for illustration.     

Comparison of designs for the three-fold nested random model

The quality of estimation of variance components depends on the design used as well as on the unknown values of the variance components. In this article, three designs are compared, namely, the balanced, staggered, and inverted nested designs for the three-fold nested random model. The comparison is based on the so-called quantile dispersion graphs using analysis of variance (ANOVA) and maximum likelihood (ML) estimates of the variance components. It is demonstrated that the staggered nested design gives more stable estimates of the variance component for the highest nesting factor than the balanced design. The reverse, however, is true in case of lower nested factors. A comparison between ANOVA and ML estimation of the variance components is also made using each of the aforementioned designs.     

Optimal blocking and foldover plans for nonregular two-level designs

This paper discusses the issue of choosing optimal designs when both blocking and foldover techniques are simultaneously employed to nonregular two-level fractional factorial designs. By using the indicator function, the treatment and block generalized wordlength patterns of the combined blocked design under a general foldover plan are defined. Some general properties of combined block designs are also obtained. Our results extend the findings of Ai et al. (2010) from regular designs to nonregular designs. Based on these theoretical results, a catalog of optimal blocking and foldover plans in terms of the generalized aberration criterion for nonregular initial design with 12, 16 and 20 runs is tabulated, respectively.     

Projection designs for mixture experiments in orthogonal blocks

Mixture experiments are often carried out in the presence of process variables, such as days of the week or different machines in a manufacturing process, or different ovens in bread and cake making. In such experiments it is particularly useful to be able to arrange the design in orthogonal blocks, so that the model in tue mixture vanauies may ue iitteu inucpenuentiy or tne UIOCK enects mtrouuceu to take account of the changes in the process variables. It is possible in some situations that some of the ingredients in the mixture, such as additives or flavourings, are present in soian quantities, pernaps as iuw a.s 5% ur even !%, resulting in the design space being restricted to only part of the mixture simplex. Hau and Box (1990) discussed the construction of experimental designs for situations where constraints are placed on the design variables. They considered projecting standard response surface designs, including factorial designs and central composite designs, into the restricted design space, and showed that the desirable property of block orthogonality is preserved by the projections considered. Here we present a number of examples of projection designs and illustrate their use when some of the ingredients are restricted to small values, such that the design space is restricted to a sub-region within the usual simplex in the mixture variables.     

Statistical tests for random effects in staggered nested designs

In this paper, we present certain statistical tests under a staggered nested design set-up, for the hypotheses that certain variance components are zero. To do so, the particular variance-covariance, structure induced by the staggering is exploited and certain results of multivariate analysis are used. In most problems, the test statistics can be easily computed. An example is provided for illustration and some power computations for comparison of test statistics are shown.     

Parametric bootstrap tests for unbalanced nested designs under heteroscedasticity

In this article, we consider the two-factor unbalanced nested design model without the assumption of equal error variance. For the problem of testing ‘main effects’ of both factors, we propose a parametric bootstrap (PB) approach and compare it with the existing generalized F (GF) test. The Type I error rates of the tests are evaluated using Monte Carlo simulation. Our studies show that the PB test performs better than the GF test. The PB test performs very satisfactorily even for small samples while the GF test exhibit poor Type I error properties when the number of factorial combinations or treatments goes up. It is also noted that the same tests can be used to test the significance of the random effect variance component in a two-factor mixed effects nested model under unequal error variances.     

Minimum aberration blocking schemes for 128-run designs

Several criteria have been proposed for ranking blocked fractional factorial designs. For large fractional factorial designs, the most appropriate minimum aberration criterion was one proposed by Cheng and Wu (2002). We justify this assertion and propose a novel construction method to overcome the computational challenge encountered in large fractional factorial designs. Tables of minimum aberration blocked designs are presented for N=128 runs and n=8–64 factors.     

Bayesian optimal blocking of factorial designs

The presence of block effects makes the optimal selection of fractional factorial designs a difficult task. The existing frequentist methods try to combine treatment and block wordlength patterns and apply minimum aberration criterion to find the optimal design. However, ambiguities exist in combining the two wordlength patterns and therefore, the optimality of such designs can be challenged. Here we propose a Bayesian approach to overcome this problem. The main technique is to postulate a model and a prior distribution to satisfy the common assumptions in blocking and then, to develop an optimal design criterion for the efficient estimation of treatment effects. We apply our method to develop regular, nonregular, and mixed-level blocked designs. Several examples are presented to illustrate the advantages of the proposed method.     

Minimax designs for estimating slopes in a trigonometric regression model

Abstract

Designs for the first order trigonometric regression model over an interval on the real line are considered for the situation where estimation of the slope of the response surface at various points in the factor space is of primary interest. Minimization of the variance of the estimated slope at a point maximized over all points in the region of interest is taken as the design criterion. Optimal designs under the minimax criterion are derived for the situation where the design region and the region of interest are identical and a symmetric “partial cycle”. Some comparisons of the minimax designs with the traditional D- and A-optimal designs are provided. Efficiencies of some exact designs under the minimax criterion are also investigated.     

Existence of nested designs with block size five

The block designs considered here are nested in the unified sense of Preece (Biometrika 54 (1967) 479–486) and Federer (in: T.A. Baneroft (Ed.), Statistical Papers in Honor of George W. Snedecor, 1972, pp. 91–114), that is, each block of the larger balanced incomplete block design contains several distinguished families of mutually disjoint sub-blocks, the sub-blocks of the same family belonging to the same system, such that each system forms the collection of blocks of some balanced incomplete block design. In this paper, it is shown that the necessary conditions for the existence of such designs are also sufficient for block size 5 and sub-block sizes 2 and 3. This, together with known results, implies the entire existence of such designs with block size 5 in general.     

Estimation of variance components in staggered nested designs

SUMMARY Variance components are estimated by two different methods for a general p stage random-effects staggered nested design. In addition to estimation from an analysis of variance, a new approach is introduced. The main features of this new technique are its simplicity and its ability to yield non-negative estimates of the variance components. The performances of the two procedures are compared using simulation and the meansquared-error criterion.     

Construction of nested space-filling designs using difference matrices

Experiments that study complex real world systems in business, engineering and sciences can be conducted at different levels of accuracy or sophistication. Nested space-filling designs are suitable for such multi-fidelity experiments. In this paper, we propose a systematic method to construct nested space-filling designs for experiments with two levels of accuracy. The method that makes use of nested difference matrices can be easily performed, many nested space-filling designs for experiments with two levels of accuracy can thus be constructed, and the resulting designs achieve stratification in low dimensions. In addition, the proposed method can also be used to obtain sliced space-filling designs for conducting computer experiments with both qualitative and quantitative factors.     

Optimum designs for estimation of optimum point under cost constraint

In this paper, we consider the estimation of the optimum factor combination in a response surface model. Assuming that the response function is quadratic concave and there is a linear cost constraint on the factor combination, we attempt to find the optimum design using the trace optimality criterion. As the criterion function involves the unknown parameters, we adopt a pseudo-Bayesian approach to resolve the problem.     

Optimal designs for non-competitive enzyme inhibition kinetic models

In this paper, we present a new method for determining optimal designs for enzyme inhibition kinetic models, which are used to model the influence of the concentration of a substrate and an inhibition on the velocity of a reaction. The approach uses a nonlinear transformation of the vector of predictors such that the model in the new coordinates is given by an incomplete response surface model. Although there exist no explicit solutions of the optimal design problem for incomplete response surface models so far, the corresponding design problem in the new coordinates is substantially more transparent, such that explicit or numerical solutions can be determined more easily. The designs for the original problem can finally be found by an inverse transformation of the optimal designs determined for the response surface model. We illustrate the method determining explicit solutions for the D-optimal design and for the optimal design problem for estimating the individual coefficients in a non-competitive enzyme inhibition kinetic model.     

Orthogonal blocking of response surface split-plot designs

When all experimental runs cannot be performed under homogeneous conditions, blocking can be used to increase the power for testing the treatment effects. Orthogonal blocking provides the same estimator of the polynomial effects as the one that would be obtained by ignoring the blocks. In many real-life design scenarios, there is at least one factor that is hard to change, leading to a split-plot structure. This paper shows that for a balanced ordinary least square–generalized least square equivalent split-plot design, orthogonal blocking can be achieved. Orthogonally blocked split-plot central composite designs are constructed and a catalog is provided.     

Minimum aberration blocking schemes for two- and three-level fractional factorial designs

The concept of minimum aberration has been extended to choose blocked fractional factorial designs (FFDs). The minimum aberration criterion ranks blocked FFDs according to their treatment and block wordlength patterns, which are often obtained by counting words in the treatment defining contrast subgroups and alias sets. When the number of factors is large, there are a huge number of words to be counted, causing some difficulties in computation. Based on coding theory, the concept of minimum moment aberration, proposed by Xu [Statist. Sinica, 13 (2003) 691–708] for unblocked FFDs, is extended to blocked FFDs. A method is then proposed for constructing minimum aberration blocked FFDs without using defining contrast subgroups and alias sets. Minimum aberration blocked FFDs for all 32 runs, 64 runs up to 32 factors, and all 81 runs are given with respect to three combined wordlength patterns.