Direct Estimation of Parameters in ODE Models Using WENDy: Weak-form Estimation of Nonlinear Dynamics

Abstract

We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of C-infinity bump functions of varying support sizes. We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at (https://github.com/MathBioCU/WENDy).

arXiv comment28 pages, 16 figures

References in corpus (12)

• Discovering governing equations from data: Sparse identification of nonlinear dynamical systems
• Ensemble-SINDy: Robust sparse model discovery in the low-data, high-noise limit, with active learning and control
• Weak SINDy For Partial Differential Equations
• PySINDy: A comprehensive Python package for robust sparse system identification
• Using Noisy or Incomplete Data to Discover Models of Spatiotemporal Dynamics
• Parameter estimation of ODE's via nonparametric estimators
• Variational system identification of the partial differential equations governing pattern-forming physics: Inference under varying fidelity and noise
• Robust and optimal sparse regression for nonlinear PDE models
• Learning Mean-Field Equations from Particle Data Using WSINDy
• Parametric Estimation of Ordinary Differential Equations with Orthogonality Conditions
• Data-driven discovery and extrapolation of parameterized pattern-forming dynamics
• Learning Anisotropic Interaction Rules from Individual Trajectories in a Heterogeneous Cellular Population

Cited by in corpus (4)

• Physics-informed active learning with simultaneous weak-form latent space dynamics identification
• Influence of initial conditions on data-driven model identification and information entropy for ideal mhd problems
• Statistical Mechanics of Dynamical System Identification
• A Constrained Optimisation Framework for Parameter Identification of the SIRD Model