Robust Estimation for Ordinary Differential Equation Models
Summary Applied scientists often like to use ordinary differential equations (ODEs) to model complex dynamic processes that arise in biology, engineering, medicine, and many other areas. It is interesting but challenging to estimate ODE parameters from noisy
data, especially when the data have some outliers. We propose a robust method to address this problem. The dynamic process is represented with a nonparametric function, which is a linear combination of basis functions. The nonparametric function is estimated by a robust penalized smoothing
method. The penalty term is defined with the parametric ODE model, which controls the roughness of the nonparametric function and maintains the fidelity of the nonparametric function to the ODE model. The basis coefficients and ODE parameters are estimated in two nested levels of optimization.
The coefficient estimates are treated as an implicit function of ODE parameters, which enables one to derive the analytic gradients for optimization using the implicit function theorem. Simulation studies show that the robust method gives satisfactory estimates for the ODE parameters from
noisy data with outliers. The robust method is demonstrated by estimating a predator–prey ODE model from real ecological data.
Document Type: Research Article
Department of Statistics & Actuarial Science, Simon Fraser University, Burnaby, British Columbia V5A1S6, Canada
Department of Statistics, University of British Columbia, Vancouver, British Columbia V6T1Z2, Canada
Publication date: December 1, 2011