Open Access
Translator Disclaimer
October 2015 Computing exact $D$-optimal designs by mixed integer second-order cone programming
Guillaume Sagnol, Radoslav Harman
Ann. Statist. 43(5): 2198-2224 (October 2015). DOI: 10.1214/15-AOS1339


Let the design of an experiment be represented by an $s$-dimensional vector $\mathbf{w}$ of weights with nonnegative components. Let the quality of $\mathbf{w}$ for the estimation of the parameters of the statistical model be measured by the criterion of $D$-optimality, defined as the $m$th root of the determinant of the information matrix $M(\mathbf{w} )=\sum_{i=1}^{s}w_{i}A_{i}A_{i}^{T}$, where $A_{i},i=1,\ldots,s$ are known matrices with $m$ rows.

In this paper, we show that the criterion of $D$-optimality is second-order cone representable. As a result, the method of second-order cone programming can be used to compute an approximate $D$-optimal design with any system of linear constraints on the vector of weights. More importantly, the proposed characterization allows us to compute an exact $D$-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 $D_{K}$-optimality, which measures the quality of $\mathbf{w}$ for the estimation of a linear parameter subsystem defined by a full-rank coefficient matrix $K$.

We prove that some other widely used criteria are also second-order cone representable, for instance, the criteria of $A$-, $A_{K}$-, $G$- and $I$-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.


Download Citation

Guillaume Sagnol. Radoslav Harman. "Computing exact $D$-optimal designs by mixed integer second-order cone programming." Ann. Statist. 43 (5) 2198 - 2224, October 2015.


Received: 1 December 2013; Revised: 1 January 2015; Published: October 2015
First available in Project Euclid: 16 September 2015

zbMATH: 1331.62384
MathSciNet: MR3396983
Digital Object Identifier: 10.1214/15-AOS1339

Primary: 62K05
Secondary: 65K05

Keywords: $D$-criterion , exact optimal designs , mixed integer programming , optimal experimental design , second-order cone programming

Rights: Copyright © 2015 Institute of Mathematical Statistics


Vol.43 • No. 5 • October 2015
Back to Top