Approximate and exact optimal designs for paired comparisons of partial profiles when there are two groups of factors

For paired comparison experiments involving pairs of multifactor options differing in a specified number of factors the problem of finding optimal designs is considered, when only main effects are to be estimated. It is presumed that the set of factors can be partitioned into two groups such that the number of levels is constant within each group. The optimal designs for this frequently encountered case are also optimal for the corresponding choice experiments under the hypothesis that the parameters in the multinomial logit model are equal to zero.     

Simulation-based sequential Bayesian design

We consider simulation-based methods for exploration and maximization of expected utility in sequential decision problems. We consider problems which require backward induction with analytically intractable expected utility integrals at each stage. We propose to use forward simulation to approximate the integral expressions, and a reduction of the allowable action space to avoid problems related to an increasing number of possible trajectories in the backward induction. The artificially reduced action space allows strategies to depend on the full history of earlier observations and decisions only indirectly through a low dimensional summary statistic. The proposed rule provides a finite-dimensional approximation to the unrestricted infinite-dimensional optimal decision rule. We illustrate the proposed approach with an application to an optimal stopping problem in a clinical trial.     

Equivalence theorem for Schur optimality of experimental designs

An experimental design is said to be Schur optimal, if it is optimal with respect to the class of all Schur isotonic criteria, which includes Kiefer's criteria of _Φp${\Phi }_p$-optimality, distance optimality criteria and many others. In the paper we formulate an easily verifiable necessary and sufficient condition for Schur optimality in the set of all approximate designs of a linear regression experiment with uncorrelated errors. We also show that several common models admit a Schur optimal design, for example the trigonometric model, the first-degree model on the Euclidean ball, and the Berman's model.     

An adaptive rule for stopping a discovery process under time censorship

Consider a finite population of large but unknown size of hidden objects. Consider searching for these objects for a period of time, at a certain cost, and receiving a reward depending on the sizes of the objects found. Suppose that the size and discovery time of the objects both have unknown distributions, but the conditional distribution of time given size is exponential with an unknown non-negative and non-decreasing function of the size as failure rate. The goal is to find an optimal way to stop the discovery process. Assuming that the above parameters are known, an optimal stopping time is derived and its asymptotic properties are studied. Then, an adaptive rule based on order restricted estimates of the distributions from truncated data is presented. This adaptive rule is shown to perform nearly as well as the optimal stopping time for large population size.     

Asymptotic properties of kernel estimates of a regression function

The consistency of the Parzen kernel-type as well as recursive kernel estimates of a regression function is shown. The rates of the convergence are studied and compared. Moreover, the problem of selecting asymptotically optimal kernels is discussed.     

Optimal discrimination designs for exponential regression models

We investigate optimal designs for discriminating between exponential regression models of different complexity, which are widely used in the biological sciences; see, e.g., Landaw [1995. Robust sampling designs for compartmental models under large prior eigenvalue uncertainties. Math. Comput. Biomed. Appl. 181–187] or Gibaldi and Perrier [1982. Pharmacokinetics. Marcel Dekker, New York]. We discuss different approaches for the construction of appropriate optimality criteria, and find sharper upper bounds on the number of support points of locally optimal discrimination designs than those given by Caratheodory's Theorem. These results greatly facilitate the numerical construction of optimal designs. Various examples of optimal designs are then presented and compared to different other designs. Moreover, to protect the experiment against misspecifications of the nonlinear model parameters, we adapt the design criteria such that the resulting designs are robust with respect to such misspecifications and, again, provide several examples, which demonstrate the advantages of our approach.     

Robust Linear Calibration

We regard the simple linear calibration problem where only the response y of the regression line y = β₀ + β₁ t is observed with errors. The experimental conditions t are observed without error. For the errors of the observations y we assume that there may be some gross errors providing outlying observations. This situation can be modeled by a conditionally contaminated regression model. In this model the classical calibration estimator based on the least squares estimator has an unbounded asymptotic bias. Therefore we introduce calibration estimators based on robust one-step-M-estimators which have a bounded asymptotic bias. For this class of estimators we discuss two problems: The optimal estimators and their corresponding optimal designs. We derive the locally optimal solutions and show that the maximin efficient designs for non-robust estimation and robust estimation coincide.     

Optimal crossover designs when carryover effects are proportional to direct effects

Crossover designs are used for a variety of different applications. While these designs have a number of attractive features, they also induce a number of special problems and concerns. One of these is the possible presence of carryover effects. Even with the use of washout periods, which are for many applications widely accepted as an indispensable component, the effect of a treatment from a previous period may not be completely eliminated. A model that has recently received renewed attention in the literature is the model in which first-order carryover effects are assumed to be proportional to direct treatment effects. Under this model, assuming that the constant of proportionality is known, we identify optimal and efficient designs for the direct effects for different values of the constant of proportionality. We also consider the implication of these results for the case that the constant of proportionality is not known.     

MV-optimality standardized through the coefficient of variation

Many of the usual criteria for optimal experimental designs do not take into account the different scale of the variance of the parameters. Dette [1997. Designing experiments with respect to “standardized” optimality criteria. J. Roy. Statist. Soc. Ser. B Stat. Methodol. 59(1), 97–110] provided a standardization based on the efficiencies for estimating each of the parameters. This approach provides designs with similar efficiencies for all of the parameters.     

D–optimality of a class of saturated main–effect plans and allied results

This paper applies the theory of unimodular matrices to prove that all saturated main effect plans of an s₁ × s₂ factorial are equivalent from the point of view of D–optimality and are hence all D–optimal. The A– and E–optimal plans in this context have also been derived. An application in sequential experimentation has been considered     

Efficient experimental designs for sigmoidal growth models

For the Weibull- and Richards-regression model robust designs are determined by maximizing a minimum of D  - or _D1$D_1$-efficiencies, taken over a certain range of the non-linear parameters. It is demonstrated that the derived designs yield a satisfactory solution of the optimal design problem for this type of model in the sense that these designs are efficient and robust with respect to misspecification of the unknown parameters. Moreover, the designs can also be used for testing the postulated form of the regression model against a simplified sub-model.     

On representing maximin efficient designs by Taylor series

This paper considers exponential and rational regression models that are nonlinear in some parameters. Recently, locally D-optimal designs for such models were investigated in [Melas, V. B., 2005. On the functional approach to optimal designs for nonlinear models. J. Statist. Plann. Inference 132, 93–116] based upon a functional approach. In this article a similar method is applied to construct maximin efficient D-optimal designs. This approach allows one to represent the support points of the designs by Taylor series, which gives us the opportunity to construct the designs by hand using tables of the coefficients of the series. Such tables are provided here for models with two nonlinear parameters. Furthermore, the recurrent formulas for constructing the tables for arbitrary numbers of parameters are introduced.     

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.     

Addition of runs to a two-level supersaturated design

The purpose of this article is to introduce a new class of extended E(s²)-optimal two level supersaturated designs obtained by adding runs to an existing E(s²)-optimal two level supersaturated design. The extended design is a union of two optimal SSDs belonging to different classes. New lower bound to E(s²) has been obtained for the extended supersaturated designs. Some examples and a small catalogue of E(s²)-optimal SSDs are also included.     

Some sufficient conditions for establishing (M,S)-optimality

Two sufficient conditions are given for an incomplete block design to be (M,S- optimal. For binary designs the conditions are (i) that the elements in each row, excluding the diagonal element, of the association matrix differ by at most one, and (ii) that the off-diagonal elements of the block characteristic matrix differ by at most one. It is also shown how the conditions can be utilized for nonbinary designs and that for blocks of size two the sufficient condition in terms of the association matrix can be attained.     

On robustness of optimal balanced resolution v plans

In this paper, we present Srivastava-Chopra optimal balanced resolution V plans for 2^m factorials (4?m?8) which are robust in the sense that, when any observation is missing, each of these designs will remain as a resolution V plan.     

Robust prediction and extrapolation designs for censored data

In this paper we present the construction of robust designs for a possibly misspecified generalized linear regression model when the data are censored. The minimax designs and unbiased designs are found for maximum likelihood estimation in the context of both prediction and extrapolation problems. This paper extends preceding work of robust designs for complete data by incorporating censoring and maximum likelihood estimation. It also broadens former work of robust designs for censored data from others by considering both nonlinearity and much more arbitrary uncertainty in the fitted regression response and by dropping all restrictions on the structure of the regressors. Solutions are derived by a nonsmooth optimization technique analytically and given in full generality. A typical example in accelerated life testing is also demonstrated. We also investigate implementation schemes which are utilized to approximate a robust design having a density. Some exact designs are obtained using an optimal implementation scheme.     

A Jackknife Method for Estimation of Variance Components

This paper concerns a method of estimation of variance components in a random effect linear model. It is mainly a resampling method and relies on the Jackknife principle. The derived estimators are presented as least squares estimators in an appropriate linear model, and one of them appears as a MINQUE (Minimum Norm Quadratic Unbiased Estimation) estimator. Our resampling method is illustrated by an example given by C. R. Rao [7] and some optimal properties of our estimator are derived for this example. In the last part, this method is used to derive an estimation of variance components in a random effect linear model when one of the components is assumed to be known.     

An extension of the ranked set sampling theory

This paper is concerned with ranked set sampling theory which is useful to estimate the population mean when the order of a sample of small size can be found without measurements or with rough methods. Consider n sets of elements each set having size m. All elements of each set are ranked but only one is selected and quantified. The average of the quantified elements is adopted as the estimator. In this paper we introduce the notion of selective probability which is a generalization of a notion from Yanagawa and Shirahata (1976). Uniformly optimal unbiased procedures are found for some (n,m). Furthermore, procedures which are unbiased for all distributions and are good for symmetric distributions are studied for (n,m) which do not allow uniformly optimal unbiased procedures.     

Robust designs for one-point extrapolation

We consider the construction of designs for the extrapolation of a regression response to one point outside of the design space. The response function is an only approximately known function of a specified linear function. As well, we allow for variance heterogeneity. We find minimax designs and corresponding optimal regression weights in the context of the following problems: (P1) for nonlinear least squares estimation with homoscedasticity, determine a design to minimize the maximum value of the mean squared extrapolation error (MSEE), with the maximum being evaluated over the possible departures from the response function; (P2) for nonlinear least squares estimation with heteroscedasticity, determine a design to minimize the maximum value of MSEE, with the maximum being evaluated over both types of departures; (P3) for nonlinear weighted least squares estimation, determine both weights and a design to minimize the maximum MSEE; (P4) choose weights and design points to minimize the maximum MSEE, subject to a side condition of unbiasedness. Solutions to (P1)–(P4) are given in complete generality. Numerical comparisons indicate that our designs and weights perform well in combining robustness and efficiency. Applications to accelerated life testing are highlighted.