PuSH - Publication Server of Helmholtz Zentrum München

Optimization and profile calculation of ODE models using second order adjoint sensitivity analysis.

Bioinformatics 34, 151-159 (2018)
Publishers Version Research data DOI PMC
Open Access Gold (Paid Option)
Creative Commons Lizenzvertrag
as soon as is submitted to ZB.
Motivation: Parameter estimation methods for ordinary differential equation (ODE) models of biological processes can exploit gradients and Hessians of objective functions to achieve convergence and computational efficiency. However, the computational complexity of established methods to evaluate the Hessian scales linearly with the number of state variables and quadratically with the number of parameters. This limits their application to low-dimensional problems. Results: We introduce second order adjoint sensitivity analysis for the computation of Hessians and a hybrid optimization-integration-based approach for profile likelihood computation. Second order adjoint sensitivity analysis scales linearly with the number of parameters and state variables. The Hessians are effectively exploited by the proposed profile likelihood computation approach. We evaluate our approaches on published biological models with real measurement data. Our study reveals an improved computational efficiency and robustness of optimization compared to established approaches, when using Hessians computed with adjoint sensitivity analysis. The hybrid computation method was more than 2-fold faster than the best competitor. Thus, the proposed methods and implemented algorithms allow for the improvement of parameter estimation for medium and large scale ODE models.
Altmetric
Additional Metrics?
Edit extra informations Login
Publication type Article: Journal article
Document type Scientific Article
Keywords Trust-region Method; Confidence-intervals; Dynamical-systems; Likelihood; Minimization; Algorithm; Approximation; Information; Gradient; Suite
Reviewing status