Statistical inference in mechanistic models: time warping for improved gradient matching

Comput Stat. 2018;33(2):1091-1123. doi: 10.1007/s00180-017-0753-z. Epub 2017 Aug 9.

Abstract

Inference in mechanistic models of non-linear differential equations is a challenging problem in current computational statistics. Due to the high computational costs of numerically solving the differential equations in every step of an iterative parameter adaptation scheme, approximate methods based on gradient matching have become popular. However, these methods critically depend on the smoothing scheme for function interpolation. The present article adapts an idea from manifold learning and demonstrates that a time warping approach aiming to homogenize intrinsic length scales can lead to a significant improvement in parameter estimation accuracy. We demonstrate the effectiveness of this scheme on noisy data from two dynamical systems with periodic limit cycle, a biopathway, and an application from soft-tissue mechanics. Our study also provides a comparative evaluation on a wide range of signal-to-noise ratios.

Keywords: Differential equations; Dynamical systems; Objective function; Reproducing kernel Hilbert space.