Slope rotatability over all directions designs for k = 4

Slope rotatability over all directions (SROAD) is a useful concept when the slope of a second-order response is to be studied. SROAD designs ensure that knowledge jof the slope is acquired symmetrically, whatever direction later becomes of more interest as the data are analyzed. In a prior paper, we explored designs for k=2 and 3 dimensions, which do not have the full symmetries of second-order designs but which still possess the SROAD property.Here, we discuss designs in higher dimensions.The introductory sections 1 and 2 are essentially identical to those of the prior paper.     

Modified Robust Second-Order Slope-Rotatable Designs

Generally it is very difficult to construct robust slope-rotatable designs along axial directions. Present paper focuses on modified second-order slope-rotatable designs (SOSRDs) with correlated errors. Modified robust second-order slope-rotatability conditions are derived for a general variance–covariance structure of errors. These conditions get simplified for intraclass correlation structure. A few robust second-order slope-rotatable designs (over all directions, or with equal maximum directional variance slope, or D-optimal slope) are examined with respect to modified robust slope-rotatability. It is observed that robust second-order slope-rotatable designs over all directions, or with equal maximum directional variance slope, or D-optimal slope are not generally modified robust second-order slope-rotatable designs.     

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.     

Modified second-order slope rotatable designs with equi-spaced levels

In this paper, a new method of construction of modified second-order slope rotatable designs with equi-spaced levels using central composite designs and balanced incomplete block designs is suggested.     

Minimax designs for estimating the slope of a second-order response surface in a cubic region

Minimization of the variance of the estimated slope of a response surface maximized over all points in the factor space is taken as the design criterion. Optimal designs under the criterion are derived for second-order polynomial regression over hypercubic regions.     

Augmented Box-Behnken Designs for Fitting Third-Order Response Surfaces

Box-Behnken designs are popular with experimenters who wish to estimate a second-order model, due to their having three levels, their simplicity and their high efficiency for the second-order model. However, there are situations in which the model is inadequate due to lack of fit caused by higher-order terms. These designs have little ability to estimate third-order terms. Using combinations of factorial points, axial points, and complementary design points, we augment these designs and develop catalogues of third-order designs for 3–12 factors. These augmented designs can be used to estimate the parameters of a third-order response surface model. Since the aim is to make the most of a situation in which the experiment was designed for an inadequate model, the designs are clearly suboptimal and not rotatable for the third-order model, but can still provide useful information.     

Space-filling orthogonal arrays of strength two

This paper considers the use of orthogonal arrays of strength two as experimental designs for fitting a surrogate model. Contrary to standard space-filling designs or Latin hypercube designs, the points of an orthogonal array of strength two are well distributed when they are projected on the two-dimensional faces of the unit cube. The aim is to determine if this property allows one to fit an accurate surrogate model when the computer response is governed by second-order interactions of some input variables. The first part of the paper is devoted to the construction of orthogonal arrays with space-filling properties. In the second part, orthogonal arrays are compared with standard designs for fitting a Gaussian process model.     

On D-optimal robust second order slope-rotatable designs

Das and Park (2006) introduced slope-rotatable designs overall directions for correlated observations which is known as A-optimal robust slope-rotatable designs. This article focuses D-optimal slope-rotatable designs for second-order response surface model with correlated observations. It has been established that robust second-order rotatable designs are also D-optimal robust slope-rotatable designs. A class of D-optimal robust second-order slope-rotatable designs has been derived for special correlation structures of errors.     

Central composite designs for estimating the optimum conditions for a second-order model

Central composite designs which maximize both the precision and the accuracy of estimates of the extremal point of a second-order response surface for fixed values of the model parameters are constructed. Two optimality criteria are developed, the one relating to precision and based on the sum of the first-order approximations to the asymptotic variances and the other to accuracy and based on the sum of squares of the second-order approximations to the asymptotic biases of the estimates of the coordinates of the extremal point. Exact and continuous central composite designs are introduced and in particular designs which place no restriction on the pattern of the weights, termed benchmark designs, and designs which comprise equally weighted factorial and equally weighted axial points, termed axial-factorial designs, are explored. Algebraic results proved somewhat elusive and the requisite designs are obtained by a mix of algebra and numeric calculation or simply numerically. An illustrative example is presented and some interesting features which emerge from that example are discussed.     

D-optimal designs for Poisson regression models

Most of the current research on optimal experimental designs for generalized linear models focuses on logistic regression models. In this paper, D-optimal designs for Poisson regression models are discussed. For the one-variable first-order Poisson regression model, it has been found that the D-optimal design, in terms of effective dose levels, is independent of the model parameters. However, it is not the case for more complicated models. We investigate how the D-optimal designs depend on the model parameters for the one-variable second-order model and two-variable interaction model. The performance of some “standard” designs that appeal to practitioners is also studied.     

Optimal designs for estimating the slope of a regression

In the common linear model with quantitative predictors we consider the problem of designing experiments for estimating the slope of the expected response in a regression. We discuss locally optimal designs, where the experimenter is only interested in the slope at a particular point, and standardized minimax optimal designs, which could be used if precise estimation of the slope over a given region is required. General results on the number of support points of locally optimal designs are derived if the regression functions form a Chebyshev system. For polynomial regression and Fourier regression models of arbitrary degree the optimal designs for estimating the slope of the regression are determined explicitly for many cases of practical interest.     

Efficient three-level screening designs using weighing matrices

Screening is the first stage of many industrial experiments and is used to determine efficiently and effectively a small number of potential factors among a large number of factors which may affect a particular response. In a recent paper, Jones and Nachtsheim [A class of three-level designs for definitive screening in the presence of second-order effects. J. Qual. Technol. 2011;43:1–15] have given a class of three-level designs for screening in the presence of second-order effects using a variant of the coordinate exchange algorithm as it was given by Meyer and Nachtsheim [The coordinate-exchange algorithm for constructing exact optimal experimental designs. Technometrics 1995;37:60–69]. Xiao et al. [Constructing definitive screening designs using conference matrices. J. Qual. Technol. 2012;44:2–8] have used conference matrices to construct definitive screening designs with good properties. In this paper, we propose a method for the construction of efficient three-level screening designs based on weighing matrices and their complete foldover. This method can be considered as a generalization of the method proposed by Xiao et al. [Constructing definitive screening designs using conference matrices. J. Qual. Technol. 2012;44:2–8]. Many new orthogonal three-level screening designs are constructed and their properties are explored. These designs are highly D-efficient and provide uncorrelated estimates of main effects that are unbiased by any second-order effect. Our approach is relatively straightforward and no computer search is needed since our designs are constructed using known weighing matrices.     

Experimental designs when there are one or more factor constraints

In response surface methodology, designs of orders one or two are often needed such that some or all the factor levels satisfy one or more linear constraints. A method is discussed for obtaining such designs by projection of a standard design onto the constraint hyperplane. It is shown that a projected design obtained from a rotatable design is also rotatable, and for a rotatable design that is also orthogonal (in particular any orthogonal first-order design) a least squares analysis carried out on the generating design supplies a least squares solution for the constrained design subject to the constraints. Some useful properties of the generating design, such as orthogonal blocking and fractionation are retained in the projected design. Some second-order mixture designs generated by two-level factorials are discussed.     

On d- and e- minimax optimal designs for estimating the axial slopes of a second-order response surface over hypercubic regions

Design of experiments for estimating the slopes of a response surface is considered. Design criteria analogous to the traditional ones but based upon the variance-covariance matrix of the estimated slopes along factor axes are proposed. Optimal designs under the proposed criteria are derived for second-order polynomial regression over hypercubic regions. Best de¬signs within some commonly used classes of designs are also obtained and their efficiencies are investigated.     

Slope-rotatable designs with equal maximum directional variance for second order response surface models

The concept of sloperotaiability with equal maximum directional vari ance for second order response surface models is introduced as a new design property. This requires that the maximum variance of the estimated slope over all possible directions be only a function of p, which is the distance from the design originif is shown that a rotatable design satisfies this property Also, minimization of tiie maximum variance of the estimated slope over all possible directions is proposed as a new design optirnality criterion, and op¬timal designs are called slope-directional minirnax designs. For the class of cquiradial designs, the slope-directional minirnax designs are compared with D— optimal designs.     

Second-Order Response Surface Model with Neighbor Effects

This article considers the second-order response surface model in which the experimental units, i.e., plots experience the neighbor effects from immediate left and right neighboring plots assuming the plots to be placed adjacent linearly with no gaps. Conditions have been derived for the estimation of coefficients of second-order response surface model. A method of constructing designs for fitting second-order response surface in the presence of neighbor effects has been developed. The designs so obtained are found to be rotatable.     

A measure of robust slope-rotatability for second-order response surface experimental designs

In response surface methodology, rotatability and slope-rotatability are natural and highly desirable properties for second-order regression models. In this paper a measure of robust slope-rotatability for second-order response surface designs with a general correlated error structure is developed and illustrated with different examples for autocorrelated error structure.     

New 3-level response surface designs constructed from incomplete block designs

Box and Behnken [1958. Some new three level second-order designs for surface fitting. Statistical Technical Research Group Technical Report No. 26. Princeton University, Princeton, NJ; 1960. Some new three level designs for the study of quantitative variables. Technometrics 2, 455–475.] introduced a class of 3-level second-order designs for fitting the second-order response surface model. These 17 Box–Behnken designs (BB designs) are available for 3–12 and 16 factors. Although BB designs were developed nearly 50 years ago, they and the central-composite designs of Box and Wilson [1951. On the experimental attainment of optimum conditions. J. Royal Statist. Soc., Ser. B 13, 1–45.] are still the most often recommended response surface designs. Of the 17 aforementioned BB designs, 10 were constructed from balanced incomplete block designs (BIBDs) and seven were constructed from partially BIBDs (PBIBDs). In this paper we show that these seven BB designs constructed from PBIBDs can be improved in terms of rotatability as well as average prediction variance, D$D$- and G$G$-efficiency. In addition, we also report new orthogonally blocked solutions for 5, 8, 9, 11 and 13 factors. Note that an 11-factor BB design is available but cannot be orthogonally blocked. All new designs can be found at http://www.math.montana.edu/∼$\sim$jobo/bbd/.     

Robustness of some designs against missing data

This article studies the robustness of several types of designs against missing data. The robustness of orthogonal resolution III fractional factorial designs and second-order rotatable designs is studied when a single observation is missing. We also study the robustness of balanced incomplete block designs when a block is missing and of Youden square designs when a column is missing.     

A Comparison of Three-Level Orthogonal Arrays in the Presence of Different Correlation Structures in Observations

When the experimenter suspects that there might be a quadratic relation between the response variable and the explanatory parameters, a design with at least three points must be employed to establish and explore this relation (second-order design). Orthogonal arrays (OAs) with three levels are often used as second-order response surface designs. Generally, we assume that the data are independent observations; however, there are many situations where this assumption may not be sustainable. In this paper, we want to compare three-level OAs with 18, 27, and 36 runs under the presence of three specific forms of correlation in observations. The aim is to derive the best designs that can be efficiently used for response surface modeling.