Designs balanced for circular residual effects

Magda (1980) introduced a model for repeated measurements designs with a circular structure of the residual effects. He proved the universal optimality of circular balanced uniform designs over a subclass of the possible designs. We strengthen his result to optimality over the set of all designs with the same number of experimental units, periods and treatments.     

Minimal balanced repeated measurements designs

Experimental designs in which treatments are applied to the experimental units, one at a time, in sequences over a number of periods, have been used in several scientific investigations and are known as repeated measurements designs. Besides direct effects, these designs allow estimation of residual effects of treatments along with adjustment for them. Assuming the existence of first-order residual effects of treatments, Hedayat & Afsarinejad (1975) gave a method of constructing minimal balanced repeated measurements [RM(v,n,p)] design for v treatments using n=2v experimental units for p [=(v+1)/2] periods when v is a prime or prime power. Here, a general method of construction of these designs for all odd v has been given along with an outline for their analysis. In terms of variances of estimated elementary contrasts between treatment effects (direct and residual), these designs are seen to be partially variance balanced based on the circular association scheme.     

ALL-BIAS DESIGNS FOR POLYNOMIAL SPLINE REGRESSION MODELS

Polynomial spline regression models of low degree have proved useful in modeling responses from designed experiments in science and engineering when simple polynomial models are inadequate. Where there is uncertainty in the number and location of the knots, or breakpoints, of the spline, then designs that minimize the systematic errors resulting from model misspecification may be appropriate. This paper gives a method for constructing such all‐bias designs for a single variable spline when the distinct knots in the assumed and true models come from some specified set. A class of designs is defined in terms of the inter‐knot intervals and sufficient conditions are obtained for a design within this class to be all‐bias under linear, quadratic and cubic spline models. An example of the construction of all‐bias designs is given.     

Robust designs for linear mixed effects models

Summary.  In health sciences, medicine and social sciences linear mixed effects models are often used to analyse time-structured data. The search for optimal designs for these models is often hampered by two problems. The first problem is that these designs are only locally optimal. The second problem is that an optimal design for one model may not be optimal for other models. In this paper the maximin principle is adopted to handle both problems, simultaneously. The maximin criterion is formulated by means of a relative efficiency measure, which gives an indication of how much efficiency is lost when the uncertainty about the models over a prior domain of parameters is taken into account. The procedure is illustrated by means of three growth studies. Results are presented for a vocabulary growth study from education, a bone gain study from medical research and an epidemiological decline in height study. It is shown that, for the mixed effects polynomial models that are applied to these studies, the maximin designs remain highly efficient for different sets of models and combinations of parameter values.     

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.     

Minimax designs for approximately linear models with AR(1) errors

We obtain designs for linear regression models under two main departures from the classical assumptions: (1) the response is taken to be only approximately linear, and (2) the errors are not assumed to be independent, but to instead follow a first-order autoregressive process. These designs have the property that they minimize (a modification of) the maximum integrated mean squared error of the estimated response, with the maximum taken over a class of departures from strict linearity and over all autoregression parameters ρ,|ρ,| < 1, of fixed sign. Specific methods of implementation are discussed. We find that an asymptotically optimal procedure for AR(1) models consists of choosing points from that design measure which is optimal for uncorrelated errors, and then implementing them in an appropriate order.     

Robust integer-valued designs for linear random intercept models

This article considers the robust design problem for linear random intercept models with both departures from fixed effects and correlated errors on a finite design space. Two strategies are proposed. One is a worst-case method minimizing the maximum value of the MSE of estimates for the fixed effects over the departure. The other is an average-case method minimizing the average value of the MSE with respect to some priors for the class of departure functions and correlation structures of random errors. Two examples are given to show robust designs for two polynomial models.     

Optimal cross-over designs for nonlinear mixed models using a first-order expansion

We consider the construction of optimal cross-over designs for nonlinear mixed effect models based on the first-order expansion. We show that for AB/BA designs a balanced subject allocation is optimal when the parameters depend on treatments only. For multiple period, multiple sequence designs, uniform designs are optimal among dual balanced designs under the same conditions. As a by-product, the same results hold for multivariate linear mixed models with variances depending on treatments.     

Robustness of designed experiments against missing data

This paper investigates the robustness of designed experiments for estimating linear functions of a subset of parameters in a general linear model against the loss of any t( ≥1) observations. Necessary and sufficient conditions for robustness of a design under a homoscedastic model are derived. It is shown that a design robust under a homoscedastic model is also robust under a general heteroscedastic model with correlated observations. As a particular case, necessary and sufficient conditions are obtained for the robustness of block designs against the loss of data. Simple sufficient conditions are also provided for the binary block designs to be robust against the loss of data. Some classes of designs, robust up to three missing observations, are identified. A-efficiency of the residual design is evaluated for certain block designs for several patterns of two missing observations. The efficiency of the residual design has also been worked out when all the observations in any two blocks, not necessarily disjoint, are lost. The lower bound to A-efficiency has also been obtained for the loss of t observations. Finally, a general expression is obtained for the efficiency of the residual design when all the observations of m ( ≥1) disjoint blocks are lost.     

Minimal circular balanced repeated measurement designs in unequal period sizes

Abstract

Repeated measurement designs (RMDs) are widely used in medicine, pharmacology, animal sciences and psychology. In these fields, there are several situations where these designs should be used in periods of different sizes. With the use of RMD, residual effects or carry over effects may arise and balanced RMDs are solution to this problem. In this article, therefore, some infinite series are developed through method of cyclic shifts to obtain circular balanced repeated measurements designs in periods of two different sizes.     

Construction of circular partially balanced-Repeated measurement designs using cyclic shifts

Repeated measurement designs are widely used in medicine, pharmacology, animal sciences and psychology. These designs balance out the residual effects. The situations where balanced repeated measurements designs require a large number of the subjects, partially-balanced repeated measurements designs should be used. In this paper some infinite series are developed which provide circular partially-balanced repeated measurement designs for p (periods) even. Catalogues of circular partially-balanced repeated measurement designs are also presented for v (treatments) ≤ 100 with p = 5, 7 & 9.     

Considerations on group-wise identical designs for linear mixed models

We consider a general class of mixed models, where the individual parameter vector is composed of a linear function of the population parameter vector plus an individual random effects vector. The linear function can vary for the different individuals. We show that the search for optimal designs for the estimation of the population parameter vector can be restricted to the class of group-wise identical designs, i.e., for each of the groups defined by the different linear functions only one individual elementary design has to be optimized. A way to apply the result to non-linear mixed models is described.     

The optimal Lp norm estimator in linear regression models

The least squares estimator is usually applied when estimating the parameters in linear regression models. As this estimator is sensitive to departures from normality in the residual distribution, several alternatives have been proposed. The L_p norm estimators is one class of such alternatives. It has been proposed that the kurtosis of the residual distribution be taken into account when a choice of estimator in the L_p norm class is made (i.e. the choice of p). In this paper, the asymtotic variance of the estimators is used as the criterion in the choice of p. It is shown that when this criterion is applied, other characteristics of the residual distribution than the kurtosis (namely moments of order p-2 and 2p-2) are important.     

A clustering-based coordinate exchange algorithm for generating G-optimal experimental designs

In the optimal experimental design literature, the G-optimality is defined as minimizing the maximum prediction variance over the entire experimental design space. Although the G-optimality is a highly desirable property in many applications, there are few computer algorithms developed for constructing G-optimal designs. Some existing methods employ an exhaustive search over all candidate designs, which is time-consuming and inefficient. In this paper, a new algorithm for constructing G-optimal experimental designs is developed for both linear and generalized linear models. The new algorithm is made based on the clustering of candidate or evaluation points over the design space and it is a combination of point exchange algorithm and coordinate exchange algorithm. In addition, a robust design algorithm is proposed for generalized linear models with modification of an existing method. The proposed algorithm are compared with the methods proposed by Rodriguez et al. [Generating and assessing exact G-optimal designs. J. Qual. Technol. 2010;42(1):3–20] and Borkowski [Using a genetic algorithm to generate small exact response surface designs. J. Prob. Stat. Sci. 2003;1(1):65–88] for linear models and with the simulated annealing method and the genetic algorithm for generalized linear models through several examples in terms of the G-efficiency and computation time. The result shows that the proposed algorithm can obtain a design with higher G-efficiency in a much shorter time. Moreover, the computation time of the proposed algorithm only increases polynomially when the size of model increases.     

On selecting the best treatment in a generalized linear model

The problem of selecting the best treatment is studied under generalized linear models. For certain balanced designs, it is shown that simple rules are Bayes with respect to any non-informative prior on the treatment effects under any monotone invariant loss. When the nuisance parameters such as block effects are assumed to follow a uniform (improper) prior or a normal prior, Bayes rules are obtained for the normal linear model under more suitable balanced designs, keeping the generality of the loss and the generality of the non-informativeness on the prior of the treatment effects. These results are extended to certain types of informative priors on the treatment effects. When the designs are unbalanced, algorithms based on the Gibbs sampler and the Laplace method are provided to compute the Bayes rules.     

Construction and analysis of edge designs from skew-symmetric supplementary difference sets

The purpose of screening experiments is to identify the dominant variables from a set of many potentially active variables which may affect some characteristic y. Edge designs were recently introduced in the literature and are constructed by using conferences matrices and were proved to be robust. We introduce a new class of edge designs which are constructed from skew-symmetric supplementary difference sets. These designs are particularly useful since they can be applied for experiments with an even number of factors and they may exist for orders where conference matrices do not exist. Using this methodology, examples of new edge designs for 6, 14, 22, 26, 38, 42, 46, 58, and 62 factors are constructed. Of special interest are the new edge designs for studying 22 and 58 factors because edge designs with these parameters have not been constructed in the literature since conference matrices of the corresponding order do not exist. The suggested new edge designs achieve the same model-robustness as the traditional edge designs. We also suggest the use of a mirror edge method as a test for the linearity of the true underlying model. We give the details of the methodology and provide some illustrating examples for this new approach. We also show that the new designs have good D-efficiencies when applied to first order models.     

Analyzing Binomial Data in a Split-Plot Design: Classical Approach or Modern Techniques?

Modeling data that are non-normally distributed with random effects is the major challenge in analyzing binomial data in split-plot designs. Seven methods for analyzing such data using mixed, generalized linear, or generalized linear mixed models are compared for the size and power of the tests. This study shows that analyzing random effects properly is more important than adjusting the analysis for non-normality. Methods based on mixed and generalized linear mixed models hold Type I error rates better than generalized linear models. Mixed model methods tend to have higher power than generalized linear mixed models when the sample size is small.     

Optimum Mixture Design for the Detection of Synergistic Effects

In mixture experiments, optimal designs for the estimation of parameters, both linear and non-linear, have been discussed by several authors. Optimal designs for the estimation of a subset of parameters have also been investigated. However, designs for testing the effects of certain factors and interactions have been studied only in the context of response surface models. In this article, we attempt to find the optimum design for testing the presence of synergistic effects in a mixture model. The classical F-test has been considered and the optimum design has been obtained so as to maximize the power of the test. It is observed that the barycenters are necessarily the support points of the trace-optimal design.     

LEVERAGE ADJUSTMENTS FOR DISPERSION MODELLING IN GENERALIZED NONLINEAR MODELS

For normal linear models, it is generally accepted that residual maximum likelihood estimation is appropriate when covariance components require estimation. This paper considers generalized linear models in which both the mean and the dispersion are allowed to depend on unknown parameters and on covariates. For these models there is no closed form equivalent to residual maximum likelihood except in very special cases. Using a modified profile likelihood for the dispersion parameters, an adjusted score vector and adjusted information matrix are found under an asymptotic development that holds as the leverages in the mean model become small. Subsequently, the expectation of the fitted deviances is obtained directly to show that the adjusted score vector is unbiased at least to O(1/n) . Exact results are obtained in the single‐sample case. The results reduce to residual maximum likelihood estimation in the normal linear case.     

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.