PuSH - Publication Server of Helmholtz Zentrum München

Lakrisenko, P. ; Pathirana, D.* ; Weindl, D. ; Hasenauer, J.

Benchmarking methods for computing local sensitivities in ordinary differential equation models at dynamic and steady states.

PLoS ONE 19:e0312148 (2024)
Publ. Version/Full Text DOI PMC
Open Access Gold
Creative Commons Lizenzvertrag
Estimating parameters of dynamic models from experimental data is a challenging, and often computationally-demanding task. It requires a large number of model simulations and objective function gradient computations, if gradient-based optimization is used. In many cases, steady-state computation is a part of model simulation, either due to steady-state data or an assumption that the system is at steady state at the initial time point. Various methods are available for steady-state and gradient computation. Yet, the most efficient pair of methods (one for steady states, one for gradients) for a particular model is often not clear. In order to facilitate the selection of methods, we explore six method pairs for computing the steady state and sensitivities at steady state using six real-world problems. The method pairs involve numerical integration or Newton's method to compute the steady-state, and-for both forward and adjoint sensitivity analysis-numerical integration or a tailored method to compute the sensitivities at steady-state. Our evaluation shows that all method pairs provide accurate steady-state and gradient values, and that the two method pairs that combine numerical integration for the steady-state with a tailored method for the sensitivities at steady-state were the most robust, and amongst the most computationally-efficient. We also observed that while Newton's method for steady-state computation yields a substantial speedup compared to numerical integration, it may lead to a large number of simulation failures. Overall, our study provides a concise overview across current methods for computing sensitivities at steady state. While our study shows that there is no universally-best method pair, it also provides guidance to modelers in choosing the right methods for a problem at hand.
Altmetric
Additional Metrics?
Edit extra informations Login
Publication type Article: Journal article
Document type Scientific Article
Corresponding Author
Keywords Patterns; Systems; Mass
ISSN (print) / ISBN 1932-6203
Journal PLoS ONE
Quellenangaben Volume: 19, Issue: 10, Pages: , Article Number: e0312148 Supplement: ,
Publisher Public Library of Science (PLoS)
Publishing Place Lawrence, Kan.
Non-patent literature Publications
Reviewing status Peer reviewed
Grants University of Bonn
German Federal Ministry of Education and Research (BMBF)
Deutsche Forschungsgemeinschaft (DFG, German Research foundation)