Sieve estimation of constant and time-varying coefficients in nonlinear ordinary differential equation models by considering both numerical error and measurement error

Abstract

This article considers estimation of constant and time-varying coefficients in nonlinear ordinary differential equation (ODE) models where analytic closed-form solutions are not available. The numerical solution-based nonlinear least squares (NLS) estimator is investigated in this study. A numerical algorithm such as the Runge--Kutta method is used to approximate the ODE solution. The asymptotic properties are established for the proposed estimators considering both numerical error and measurement error. The B-spline is used to approximate the time-varying coefficients, and the corresponding asymptotic theories in this case are investigated under the framework of the sieve approach. Our results show that if the maximum step size of the $p$-order numerical algorithm goes to zero at a rate faster than $n^{-1/(p\wedge 4)}$, the numerical error is negligible compared to the measurement error. This result provides a theoretical guidance in selection of the step size for numerical evaluations of ODEs. Moreover, we have shown that the numerical solution-based NLS estimator and the sieve NLS estimator are strongly consistent. The sieve estimator of constant parameters is asymptotically normal with the same asymptotic co-variance as that of the case where the true ODE solution is exactly known, while the estimator of the time-varying parameter has the optimal convergence rate under some regularity conditions. The theoretical results are also developed for the case when the step size of the ODE numerical solver does not go to zero fast enough or the numerical error is comparable to the measurement error. We illustrate our approach with both simulation studies and clinical data on HIV viral dynamics.

Published in at http://dx.doi.org/10.1214/09-AOS784 the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (2)

• Two likelihood-based semiparametric estimation methods for panel count data with covariates
• Parameter estimation of ODE's via nonparametric estimators

Cited by in corpus (13)

• $\sqrt{n}$-consistent parameter estimation for systems of ordinary differential equations: bypassing numerical integration via smoothing
• Parametric Estimation of Ordinary Differential Equations with Orthogonality Conditions
• Bayesian Analysis of ODE's: solver optimal accuracy and Bayes factors
• Laplace based approximate posterior inference for differential equation models
• Estimation of ordinary differential equation models with discretization error quantification
• Estimation of nonlinear differential equation model for glucose-insulin dynamics in type I diabetic patients using generalized smoothing
• A regularization method for the parameter estimation problem in ordinary differential equations via discrete optimal control theory
• Bayesian estimation in differential equation models
• Constrained Nonlinear and Mixed Effects Differential Equation Models for Dynamic Cell Polarity Signaling
• Nonparametric Estimation of Isotropic Covariance Function
• Optimal Rate of Direct Estimators in Systems of Ordinary Differential Equations Linear in Functions of the Parameters
• Ordinary Differential Equation Models and their Computation Methods
• The Influence of Numerical Error on an Inverse Problem Methodology in PDE Models