93
72
v1v2v3 (latest)

Computing exact DD-optimal designs by mixed integer second-order cone programming

Abstract

Let the design of an experiment be represented by an ss-dimensional vector w\mathbf {w} of weights with nonnegative components. Let the quality of w\mathbf {w} for the estimation of the parameters of the statistical model be measured by the criterion of DD-optimality, defined as the mmth root of the determinant of the information matrix M(w)=i=1swiAiAiTM(\mathbf {w})=\sum_{i=1}^sw_iA_iA_i^T, where Ai,i=1,,sA_i,i=1,\ldots,s are known matrices with mm rows. In this paper, we show that the criterion of DD-optimality is second-order cone representable. As a result, the method of second-order cone programming can be used to compute an approximate DD-optimal design with any system of linear constraints on the vector of weights. More importantly, the proposed characterization allows us to compute an exact DD-optimal design, which is possible thanks to high-quality branch-and-cut solvers specialized to solve mixed integer second-order cone programming problems. Our results extend to the case of the criterion of DKD_K-optimality, which measures the quality of w\mathbf {w} for the estimation of a linear parameter subsystem defined by a full-rank coefficient matrix KK. We prove that some other widely used criteria are also second-order cone representable, for instance, the criteria of AA-, AKA_K-, GG- and II-optimality. We present several numerical examples demonstrating the efficiency and general applicability of the proposed method. We show that in many cases the mixed integer second-order cone programming approach allows us to find a provably optimal exact design, while the standard heuristics systematically miss the optimum.

View on arXiv
Comments on this paper