Penalized Nonlinear Least Squares Estimation of Time-Varying Parameters in Ordinary Differential Equations

Jiguo Cao, Jianhua Z. Huang, Hulin Wu

Research output: Contribution to journalArticlepeer-review

21 Scopus citations

Abstract

Ordinary differential equations (ODEs) are widely used in biomedical research and other scientific areas to model complex dynamic systems. It is an important statistical problem to estimate parameters in ODEs from noisy observations. In this article we propose a method for estimating the time-varying coefficients in an ODE. Our method is a variation of the nonlinear least squares where penalized splines are used to model the functional parameters and the ODE solutions are approximated also using splines. We resort to the implicit function theorem to deal with the nonlinear least squares objective function that is only defined implicitly. The proposed penalized nonlinear least squares method is applied to estimate a HIV dynamic model from a real dataset. Monte Carlo simulations show that the new method can provide much more accurate estimates of functional parameters than the existing two-step local polynomial method which relies on estimation of the derivatives of the state function. Supplemental materials for the article are available online.
Original languageEnglish (US)
Pages (from-to)42-56
Number of pages15
JournalJournal of Computational and Graphical Statistics
Volume21
Issue number1
DOIs
StatePublished - Jan 2012
Externally publishedYes

Fingerprint

Dive into the research topics of 'Penalized Nonlinear Least Squares Estimation of Time-Varying Parameters in Ordinary Differential Equations'. Together they form a unique fingerprint.

Cite this