A-optimal designs for heteroscedastic multifactor regression models

Abstract

This paper searches for A-optimal designs for Kronecker product and additive regression models when the errors are heteroscedastic. Sufficient conditions are given so that A-optimal designs for the multifactor models can be built from A-optimal designs for their sub-models with a single factor. The results of an efficiency study carried out to check the adequacy of the products of optimal designs for uni-factor marginal models when these are used to estimate different multi-factor models are also reported.     

E-optimal designs for regression models with quantitative factors— a reasonable choice?

For regression models with quantitative factors it is illustrated that the E-optimal design can be extremely inefficient in the sense that it degenerates to a design which takes all observations at only one point. This phenomenon is caused by the different size of the elements in the covariance matrix of the least-squares estimator for the unknown parameters. For these reasons we propose to replace the E-criterion by a corresponding standardized version. The advantage of this approach is demonstrated for the polynomial regression on a nonnegative interval, where the classical and standardized E-optimal designs can be found explicitly. The described phenomena are not restricted to the E-criterion but appear for nearly all optimality criteria proposed in the literature. Therefore standardization is recommended for optimal experimental design in regression models with quantitative factors. The optimal designs with respect to the new standardized criteria satisfy a similar invariance property as the famous D-optimal designs, which allows an easy calculation of standardized optimal designs on many linearly transformed design spaces.     

General orthogonal designs for parameter estimation of ANOVA models under weighted sum-to-zero constraints

ABSTRACT

Traditional studies on optimal designs for ANOVA parameter estimation are based on the framework of equal probabilities of appearance for each factor's levels. However, this premise does not hold in a variety of experimental problems, and it is of theoretical and practical interest to investigate optimal designs for parameters with unequal appearing odds. In this paper, we propose a general orthogonal design via matrix image, in which all columns’ matrix images are orthogonal with each other. Our main results show that such designs have A- and E-optimalities on the estimation of ANOVA parameters which have unequal appearing odds. In addition, we develop two simple methods to construct the proposed designs. The optimality of the design is also validated by a simulation study.     

T-optimal designs for multi-factor polynomial regression models via a semidefinite relaxation method

Statistics and Computing - We consider T-optimal experiment design problems for discriminating multi-factor polynomial regression models where the design space is defined by polynomial inequalities...     

A semidefinite programming study of the Elfving theorem

The theorem of Elfving is one of the most important and earliest results which have led to the theory of optimal design of experiments. This paper presents a fresh study of it from the viewpoint of modern semidefinite programming. There is one-to-one correspondence between solutions of the derived semidefinite programming problem (SDP) and c-optimal designs. We also derive a uniqueness theorem which ensures a unique optimal design without assuming the linear independence property over the largest set of supporting points. The SDP can also be cast as an ?₁‐convex program which has recently been extensively studied and often yields sparse solutions. Our numerical experiments on the trigonometric regression model confirm that the SDP does produce a sparse optimal design.     

Comparison of design for quadratic regression on cubes

Designs for quadratic regression are considered when the possible choices of the controllable variable are points x=(x₁,x₂,…,x_q) in the q-dimensional cube of side 2. The designs that are optimum with respect to such criteria as those of D-, A-, and E-optimality are compared in their performance relative to these and other criteria. Some of the results are developed algebraically; others, numerically. The possible supports of E-optimum designs are much more numerous than the D-optimum supports characterized earlier. The A-optimum design appears to be fairly robust in its efficiency, under variation of criterion.     

Extrapolation designs and Φp-optimum designs for cubic regression on the q-ball

This paper continues earlier work of the authors in carrying out the program discussed in Kiefer (1975), of comparing the performance of designs under various optimality criteria. Designs for extrapolation problems are also obtained. The setting is that in which the controllable variable takes on values in the q-dimensional unit ball, and the regression is cubic. Thus, the ideas of comparison are tested for a model more complex than the quadratic models discussed previously. The E-optimum design performs well in terms of other criteria, as well as for extrapolation to larger balls. A method of simplifying the calculations to obtain approximately optimum designs, is illustrated.     

Optimal orthogonal block designs for four-mixture components in two blocks based on F-squares for Becker's models and K-model

Orthogonal block designs in mixture experiments have been extensively studied by various authors. Aggarwal et al. [M.L. Aggarwal, P. Singh, V. Sarin, and B. Husain, Mixture designs in orthogonal blocks using F-squares, METRON – Int. J. Statist. LXVII(2) (2009), pp. 105–128] considered the case of components assuming the same volume fractions and obtained mixture designs in orthogonal blocks using F-squares. In this paper, we have used the class of designs presented by Aggarwal et al. and have obtained D-, A- and E-optimal orthogonal block designs for four components in two blocks for Becker's mixture models and K-model, respectively. Orthogonality conditions for the considered models are also given.     

Bayesian D-optimal designs for exponential regression models

We consider the Bayesian D-optimal design problem for exponential growth models with one, two or three parameters. For the one-parameter model conditions on the shape of the density of the prior distribution and on the range of its support are given guaranteeing that a one-point design is also Bayesian D-optimal within the class of all designs. In the case of two parameters the best two-point designs are determined and for special prior distributions it is proved that these designs are Bayesian D-optimal. Finally, the exponential growth model with three parameters is investigated. The best three-point designs are characterized by a nonlinear equation. The global optimality of these designs cannot be proved analytically and it is demonstrated that these designs are also Bayesian D-optimal within the class of all designs if gamma-distributions are used as prior distributions.     

Restricted minimax robust designs for misspecified regression models

The authors propose and explore new regression designs. Within a particular parametric class, these designs are minimax robust against bias caused by model misspecification while attaining reasonable levels of efficiency as well. The introduction of this restricted class of designs is motivated by a desire to avoid the mathematical and numerical intractability found in the unrestricted minimax theory. Robustness is provided against a family of model departures sufficiently broad that the minimax design measures are necessarily absolutely continuous. Examples of implementation involve approximate polynomial and second order multiple regression.     

D-optimal designs for beta regression models with single predictor

The problem of construction of D-optimal designs for beta regression models involving one predictor is considered for the mean-precision parameterization suggested by Ferrari and Cribari-Neto [Beta regression for modelling rates and proportions. J Appl Stat. 2004;31:799–815]. Here we use the logit link function for the mean sub-model. These designs are presented and discussed for unrestricted as well as restricted design regions by considering the precision parameter as (1) a known constant and (2) an unknown constant. Efficiency comparison of obtained designs with commonly used equi-weighted, equi-spaced designs is made to recommend designs for practical use. Real-life applications are given to show the usefulness of these designs.     

Bayesian and maximin optimal designs for heteroscedastic regression models

The authors consider the problem of constructing standardized maximin D‐optimal designs for weighted polynomial regression models. In particular they show that by following the approach to the construction of maximin designs introduced recently by Dette, Haines & Imhof (2003), such designs can be obtained as weak limits of the corresponding Bayesian _q‐optimal designs. They further demonstrate that the results are more broadly applicable to certain families of nonlinear models. The authors examine two specific weighted polynomial models in some detail and illustrate their results by means of a weighted quadratic regression model and the Bleasdale–Nelder model. They also present a capstone example involving a generalized exponential growth model.     

D-optimal designs for combined polynomial and trigonometric regression on a partial circle

Consider the D-optimal designs for a combined polynomial and trigonometric regression on a partial circle. It is shown that the optimal design is equally supported and the structure of the optimal design depends only on the length of the design interval and the support points are analytic functions of this parameter. Moreover, the Taylor expansion of the optimal support points can be determined efficiently by a recursive procedure. Examples are presented to illustrate the procedures for computing the optimal designs.     

Optimal designs based on exact confidence regions for parameter estimation of a nonlinear regression model

This paper is concerned with the proposal of optimality criteria, referred to as X  - and X^′$X^{\prime }$-optimality criteria, and the construction of X  - and X^′$X^{\prime }$-optimal designs, for nonlinear regression models. These optimal designs aim at improving the estimation of parameters of this class of models. The principle of these criteria is the minimization, with respect to the design, of the expected volume of a particular exact parametric confidence region. In this paper we give detailed definitions, properties, and computation methods of X  - and X^′$X^{\prime }$-optimal designs. We also compare these designs with the classic local D-optimal designs, with regard to robustness and efficiency, for two very well-known academic models (Box–Lucas and Michaelis–Menten models).     

Minimax designs for 2k factorial experiments for generalized linear models

ABSTRACT

Formulas for A- and C-optimal allocations for binary factorial experiments in the context of generalized linear models are derived. Since the optimal allocations depend on GLM weights, which often are unknown, a minimax strategy is considered. This is shown to be simple to apply to factorial experiments. Efficiency is used to evaluate the resulting design. In some cases, the minimax design equals the optimal design. For other cases no general conclusion can be drawn. An example of a two-factor logit model suggests that the minimax design performs well, and often better than a uniform allocation.     

Computing A-optimal Designs for Weighted Polynomial Regression by Taylor Expansion

This article is concerned with the problem of constructing A-optimal design for polynomial regression with analytic weight function on the interval [m ? a, m + a], m, a > 0. It is shown that the structure of the optimal design depends on a and weight function only, as a close to 0. Moreover, if the weight function is an analytic function a, then a scaled version of optimal support points, and weights are analytic functions of a at a = 0. We make use of a Taylor expansion to provide a recursive procedure for calculating the A-optimal designs. Examples are presented to illustrate the procedures for computing the optimal designs.     

The effect of small sample on optimal designs for logistic regression models

Asymptotic methods are commonly used in statistical inference for unknown parameters in binary data models. These methods are based on large sample theory, a condition which may be in conflict with small sample size and hence leads to poor results in the optimal designs theory. In this paper, we apply the second order expansions of the maximum likelihood estimator and derive a matrix formula for the mean square error (MSE) to obtain more precise optimal designs based on the MSE. Numerical results indicate the new optimal designs are more efficient than the optimal designs based on the information matrix.     

Third-order inference for autocorrelation in nonlinear regression models

We propose third-order likelihood-based methods to derive highly accurate p-value approximations for testing autocorrelated disturbances in nonlinear regression models. The proposed methods are particularly accurate for small- and medium-sized samples whereas commonly used first-order methods like the signed log-likelihood ratio test, the Kobayashi (1991) test, and the standardized test can be seriously misleading in these cases. Two Monte Carlo simulations are provided to show how the proposed methods outperform the above first-order methods. An empirical example applied to US population census data is also provided to illustrate the implementation of the proposed method and its usefulness in practice.     

Minimax A-, c-, and I-optimal regression designs for models with heteroscedastic errors

It is well known that it is difficult to construct minimax optimal designs. Furthermore, since in practice we never know the true error variance, it is important to allow small deviations and construct robust optimal designs. We investigate a class of minimax optimal regression designs for models with heteroscedastic errors that are robust against possible misspecification of the error variance. Commonly used A-, c-, and I-optimality criteria are included in this class of minimax optimal designs. Several theoretical results are obtained, including a necessary condition and a reflection symmetry for these minimax optimal designs. In this article, we focus mainly on linear models and assume that an approximate error variance function is available. However, we also briefly discuss how the methodology works for nonlinear models. We then propose an effective algorithm to solve challenging nonconvex optimization problems to find minimax designs on discrete design spaces. Examples are given to illustrate minimax optimal designs and their properties.     

Some robust designs for polynomial regression models

The author identifies static optimal designs for polynomial regression models with or without intercept. His optimality criterion is an average between the D‐optimality criterion for the estimation of low‐degree terms and the D₈‐optimality criterion for testing the significance of higher degree terms. His work relies on classical results concerning canonical moments and the theory of continued fractions.