Parameter estimation and change-point detection from Dynamic Contrast Enhanced MRI data using stochastic differential equations

Dynamic Contrast Enhanced imaging (DCE-imaging) following a contrast agent bolus allows the extraction of information on tissue micro-vascularization. The dynamic signals obtained from DCE-imaging are modeled by pharmacokinetic compartmental models which integrate the Arterial Input Function. These models use ordinary differential equations (ODEs) to describe the exchanges between the arterial and capillary plasma and the extravascular-extracellular space. Their least squares fitting takes into account measurement noises but fails to deal with unpredictable fluctuations due to external/internal sources of variations (patients’ anxiety, time-varying parameters, measurement errors in the input function, etc.). Adding Brownian components to the ODEs leads to stochastic differential equations (SDEs). In DCE-imaging, SDEs are discretely observed with an additional measurement noise. We propose to estimate the parameters of these noisy SDEs by maximum likelihood, using the Kalman filter. In DCE-imaging, the contrast agent injected in vein arrives in plasma with an unknown time delay. The delay parameter induces a change-point in the drift of the SDE and ODE models, which is estimated also. Estimations based on the SDE and ODE pharmacokinetic models are compared to real DCE-MRI data. They show that the use of SDE provides robustness in the estimation results. A simulation study confirms these results.     

A comparison of three different stochastic population models with regard to persistence time

Results are summarized from the literature on three commonly used stochastic population models with regard to persistence time. In addition, several new results are introduced to clearly illustrate similarities between the models. Specifically, the relations between the mean persistence time and higher-order moments for discrete-time Markov chain models, continuous-time Markov chain models, and stochastic differential equation models are compared for populations experiencing demographic variability. Similarities between the models are demonstrated analytically, and computational results are provided to show that estimated persistence times for the three stochastic models are generally in good agreement when the models are consistently formulated. As an example, the three stochastic models are applied to a population satisfying logistic growth. Logistic growth is interesting as different birth and death rates can yield the same logistic differential equation. However, the persistence behavior of the population is strongly dependent on the explicit forms for the birth and death rates. Computational results demonstrate how dramatically the mean persistence time can vary for different populations that experience the same logistic growth.     

Bayesian inference for stochastic kinetic models using a diffusion approximation

This article is concerned with the Bayesian estimation of stochastic rate constants in the context of dynamic models of intracellular processes. The underlying discrete stochastic kinetic model is replaced by a diffusion approximation (or stochastic differential equation approach) where a white noise term models stochastic behavior and the model is identified using equispaced time course data. The estimation framework involves the introduction of m- 1 latent data points between every pair of observations. MCMC methods are then used to sample the posterior distribution of the latent process and the model parameters. The methodology is applied to the estimation of parameters in a prokaryotic autoregulatory gene network.     

Mean-square stability of a stochastic model for bacteriophage infection with time delays

We consider the stability properties of the positive equilibrium of a stochastic model for bacteriophage infection with discrete time delay. Conditions for mean-square stability of the trivial solution of the linearized system around the equilibrium are given by the construction of suitable Lyapunov functionals. The numerical simulations of the strong solutions of the arising stochastic delay differential system suggest that, even for the original non-linear model, the longer the incubation time the more the phage and bacteria populations can coexist on a stable equilibrium in a noisy environment for very long time.     

Testing for dependence between failure time and visit compliance with interval-censored data

Interval-censored failure-time data arise when subjects miss prescheduled visits at which the failure is to be assessed. The resulting intervals in which the failure is known to have occurred are overlapping. Most approaches to the analysis of these data assume that the visit-compliance process is ignorable with respect to likelihood analysis of the failure-time distribution. While this assumption offers considerable simplification, it is not always plausible. Here we test for dependence between the failure- and visit-compliance processes, applicable to studies in which data collection continues after the occurrence of the failure. We do not make any of the assumptions made by previous authors about the joint distribution of the visit-compliance process, a covariate process, and the failure time. Instead, we consider conditional models of the true failure history given the current visit compliance at each visit time, allowing for correlation across visit times. Because failure status is not known at some visit times due to missed visits, only models of the observed failure history given current visit compliance are estimable. We describe how the parameters from these models can be used to test for a negative association and how bounds on unestimable parameters provided by the observed data are needed additionally to infer a positive association. We illustrate the method with data from an AIDS study and we investigate the power of the test through a simulation study.