9,277 research outputs found
Parameter estimation of ODE's via nonparametric estimators
Ordinary differential equations (ODE's) are widespread models in physics,
chemistry and biology. In particular, this mathematical formalism is used for
describing the evolution of complex systems and it might consist of
high-dimensional sets of coupled nonlinear differential equations. In this
setting, we propose a general method for estimating the parameters indexing
ODE's from times series. Our method is able to alleviate the computational
difficulties encountered by the classical parametric methods. These
difficulties are due to the implicit definition of the model. We propose the
use of a nonparametric estimator of regression functions as a first-step in the
construction of an M-estimator, and we show the consistency of the derived
estimator under general conditions. In the case of spline estimators, we prove
asymptotic normality, and that the rate of convergence is the usual
-rate for parametric estimators. Some perspectives of refinements of
this new family of parametric estimators are given.Comment: Published in at http://dx.doi.org/10.1214/07-EJS132 the Electronic
Journal of Statistics (http://www.i-journals.org/ejs/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Estimation of constant and time-varying dynamic parameters of HIV infection in a nonlinear differential equation model
Modeling viral dynamics in HIV/AIDS studies has resulted in a deep
understanding of pathogenesis of HIV infection from which novel antiviral
treatment guidance and strategies have been derived. Viral dynamics models
based on nonlinear differential equations have been proposed and well developed
over the past few decades. However, it is quite challenging to use experimental
or clinical data to estimate the unknown parameters (both constant and
time-varying parameters) in complex nonlinear differential equation models.
Therefore, investigators usually fix some parameter values, from the literature
or by experience, to obtain only parameter estimates of interest from clinical
or experimental data. However, when such prior information is not available, it
is desirable to determine all the parameter estimates from data. In this paper
we intend to combine the newly developed approaches, a multi-stage
smoothing-based (MSSB) method and the spline-enhanced nonlinear least squares
(SNLS) approach, to estimate all HIV viral dynamic parameters in a nonlinear
differential equation model. In particular, to the best of our knowledge, this
is the first attempt to propose a comparatively thorough procedure, accounting
for both efficiency and accuracy, to rigorously estimate all key kinetic
parameters in a nonlinear differential equation model of HIV dynamics from
clinical data. These parameters include the proliferation rate and death rate
of uninfected HIV-targeted cells, the average number of virions produced by an
infected cell, and the infection rate which is related to the antiviral
treatment effect and is time-varying. To validate the estimation methods, we
verified the identifiability of the HIV viral dynamic model and performed
simulation studies.Comment: Published in at http://dx.doi.org/10.1214/09-AOAS290 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Statistical modelling of conidial discharge of entomophthoralean fungi using a newly discovered Pandora species
Entomophthoralean fungi are insect pathogenic fungi and are characterized by
their active discharge of infective conidia that infect insects. Our aim was to
study the effects of temperature on the discharge and to characterize the
variation in the associated temporal pattern of a newly discovered Pandora
species with focus on peak location and shape of the discharge. Mycelia were
incubated at various temperatures in darkness, and conidial discharge was
measured over time. We used a novel modification of a statistical model
(pavpop), that simultaneously estimates phase and amplitude effects, into a
setting of generalized linear models. This model is used to test hypotheses of
peak location and discharge of conidia. The statistical analysis showed that
high temperature leads to an early and fast decreasing peak, whereas there were
no significant differences in total number of discharged conidia. Using the
proposed model we also quantified the biological variation in the timing of the
peak location at a fixed temperature.Comment: 23 pages including supplementary materia
Approximation techniques for parameter estimation and feedback control for distributed models of large flexible structures
Approximation ideas are discussed that can be used in parameter estimation and feedback control for Euler-Bernoulli models of elastic systems. Focusing on parameter estimation problems, ways by which one can obtain convergence results for cubic spline based schemes for hybrid models involving an elastic cantilevered beam with tip mass and base acceleration are outlined. Sample numerical findings are also presented
Computational methods for estimation of parameters in hyperbolic systems
Approximation techniques for estimating spatially varying coefficients and unknown boundary parameters in second order hyperbolic systems are discussed. Methods for state approximation (cubic splines, tau-Legendre) and approximation of function space parameters (interpolatory splines) are outlined and numerical findings for use of the resulting schemes in model "one dimensional seismic inversion' problems are summarized
Fast inference in nonlinear dynamical systems using gradient matching
Parameter inference in mechanistic models of
coupled differential equations is a topical problem.
We propose a new method based on kernel
ridge regression and gradient matching, and
an objective function that simultaneously encourages
goodness of fit and penalises inconsistencies
with the differential equations. Fast minimisation
is achieved by exploiting partial convexity
inherent in this function, and setting up an iterative
algorithm in the vein of the EM algorithm.
An evaluation of the proposed method on various
benchmark data suggests that it compares
favourably with state-of-the-art alternatives
Sieve estimation of constant and time-varying coefficients in nonlinear ordinary differential equation models by considering both numerical error and measurement error
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 -order
numerical algorithm goes to zero at a rate faster than , 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.Comment: 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
- …