Sensitivity approximation by the Peano-Baker series

Abstract In this paper we develop a new method for numerically approximating sensitivities in parameter-dependent ordinary differential equations (ODEs). Our approach, intended for situations where the standard forward and adjoint sensitivity analyses become too computationally costly for practical purposes, is based on the Peano-Baker series from control theory. Using this series, we construct a representation of the sensitivity matrix $$\textbf{S}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>S</mml:mi> </mml:math> and, from this representation, a numerical method for approximating $$\textbf{S}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>S</mml:mi> </mml:math> . We prove that, under standard regularity assumptions, the error of our method scales as $$\mathcal {O}(\Delta t^2_{\textrm{max}})$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>O</mml:mi> <mml:mo>(</mml:mo> <mml:mi>Δ</mml:mi> <mml:msubsup> <mml:mi>t</mml:mi> <mml:mtext>max</mml:mtext> <mml:mn>2</mml:mn> </mml:msubsup> <mml:mo>)</mml:mo> </mml:mrow> </mml:math> , where $$\Delta t_{\textrm{max}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>Δ</mml:mi> <mml:msub> <mml:mi>t</mml:mi> <mml:mtext>max</mml:mtext> </mml:msub> </mml:mrow> </mml:math> is the largest time step used when numerically solving the ODE. We illustrate the performance of the method in several numerical experiments, taken from both the systems biology setting and more classical dynamical systems. The experiments show the sought-after improvement in running time of our method compared to the forward sensitivity approach. In experiments involving a random linear system, the forward approach requires roughly $$\sqrt{n}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msqrt> <mml:mi>n</mml:mi> </mml:msqrt> </mml:math> longer computational time, where n is the dimension of the parameter space, than our proposed method.