3,737 research outputs found
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
Stochastic MPC Design for a Two-Component Granulation Process
We address the issue of control of a stochastic two-component granulation
process in pharmaceutical applications through using Stochastic Model
Predictive Control (SMPC) and model reduction to obtain the desired particle
distribution. We first use the method of moments to reduce the governing
integro-differential equation down to a nonlinear ordinary differential
equation (ODE). This reduced-order model is employed in the SMPC formulation.
The probabilistic constraints in this formulation keep the variance of
particles' drug concentration in an admissible range. To solve the resulting
stochastic optimization problem, we first employ polynomial chaos expansion to
obtain the Probability Distribution Function (PDF) of the future state
variables using the uncertain variables' distributions. As a result, the
original stochastic optimization problem for a particulate system is converted
to a deterministic dynamic optimization. This approximation lessens the
computation burden of the controller and makes its real time application
possible.Comment: American control Conference, May, 201
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
A minimal HIV-AIDS infection model with general incidence rate and application to Morocco data
We study the global dynamics of a SICA infection model with general incidence
rate. The proposed model is calibrated with cumulative cases of infection by
HIV-AIDS in Morocco from 1986 to 2015. We first prove that our model is
biologically and mathematically well-posed. Stability analysis of different
steady states is performed and threshold parameters are identified where the
model exhibits clearance of infection or maintenance of a chronic infection.
Furthermore, we examine the robustness of the model to some parameter values by
examining the sensitivity of the basic reproduction number. Finally, using
numerical simulations with real data from Morocco, we show that the model
predicts well such reality.Comment: This is a preprint of a paper whose final and definite form is with
'Statistics Opt. Inform. Comput.', Vol. 7, No 2 (2019). See
[http://www.IAPress.org]. Submitted 16/Sept/2018; Revised 10 & 15/Dec/2018;
Accepted 15/Dec/201
Fitting stochastic epidemic models to gene genealogies using linear noise approximation
Phylodynamics is a set of population genetics tools that aim at
reconstructing demographic history of a population based on molecular sequences
of individuals sampled from the population of interest. One important task in
phylodynamics is to estimate changes in (effective) population size. When
applied to infectious disease sequences such estimation of population size
trajectories can provide information about changes in the number of infections.
To model changes in the number of infected individuals, current phylodynamic
methods use non-parametric approaches, parametric approaches, and stochastic
modeling in conjunction with likelihood-free Bayesian methods. The first class
of methods yields results that are hard-to-interpret epidemiologically. The
second class of methods provides estimates of important epidemiological
parameters, such as infection and removal/recovery rates, but ignores variation
in the dynamics of infectious disease spread. The third class of methods is the
most advantageous statistically, but relies on computationally intensive
particle filtering techniques that limits its applications. We propose a
Bayesian model that combines phylodynamic inference and stochastic epidemic
models, and achieves computational tractability by using a linear noise
approximation (LNA) --- a technique that allows us to approximate probability
densities of stochastic epidemic model trajectories. LNA opens the door for
using modern Markov chain Monte Carlo tools to approximate the joint posterior
distribution of the disease transmission parameters and of high dimensional
vectors describing unobserved changes in the stochastic epidemic model
compartment sizes (e.g., numbers of infectious and susceptible individuals). We
apply our estimation technique to Ebola genealogies estimated using viral
genetic data from the 2014 epidemic in Sierra Leone and Liberia.Comment: 43 pages, 6 figures in the main tex
- …