Deterministic global optimization algorithm based on outer approximation for the parameter estimation of nonlinear dynamic biological systems
OPEN ACCESS
Author / Producer
Date
2012-05-12
Publication Type
Journal Article
ETH Bibliography
no
Citations
Altmetric
OPEN ACCESS
Data
Rights / License
Abstract
Background
The estimation of parameter values for mathematical models of biological systems is an optimization problem that is particularly challenging due to the nonlinearities involved. One major difficulty is the existence of multiple minima in which standard optimization methods may fall during the search. Deterministic global optimization methods overcome this limitation, ensuring convergence to the global optimum within a desired tolerance. Global optimization techniques are usually classified into stochastic and deterministic. The former typically lead to lower CPU times but offer no guarantee of convergence to the global minimum in a finite number of iterations. In contrast, deterministic methods provide solutions of a given quality (i.e., optimality gap), but tend to lead to large computational burdens.
Results
This work presents a deterministic outer approximation-based algorithm for the global optimization of dynamic problems arising in the parameter estimation of models of biological systems. Our approach, which offers a theoretical guarantee of convergence to global minimum, is based on reformulating the set of ordinary differential equations into an equivalent set of algebraic equations through the use of orthogonal collocation methods, giving rise to a nonconvex nonlinear programming (NLP) problem. This nonconvex NLP is decomposed into two hierarchical levels: a master mixed-integer linear programming problem (MILP) that provides a rigorous lower bound on the optimal solution, and a reduced-space slave NLP that yields an upper bound. The algorithm iterates between these two levels until a termination criterion is satisfied.
Conclusion
The capabilities of our approach were tested in two benchmark problems, in which the performance of our algorithm was compared with that of the commercial global optimization package BARON. The proposed strategy produced near optimal solutions (i.e., within a desired tolerance) in a fraction of the CPU time required by BARON.
Permanent link
Publication status
published
External links
Editor
Book title
Journal / series
Volume
13 (1)
Pages / Article No.
90
Publisher
BioMed Central
Event
Edition / version
Methods
Software
Geographic location
Date collected
Date created
Subject
Organisational unit
09655 - Guillén Gosálbez, Gonzalo / Guillén Gosálbez, Gonzalo