EM algorithm coupled with particle filter for maximum likelihood parameter estimation of stochastic differential mixed-effects models
نویسندگان
چکیده
Biological processes measured repeatedly among a series of individuals are standardly analyzed by mixed models. These biological processes can be adequately modeled by parametric Stochastic Differential Equations (SDEs). We focus on the parametric maximum likelihood estimation of this mixed-effects model defined by SDE. As the likelihood is not explicit, we propose a stochastic version of the Expectation-Maximization algorithm combined with the Particle Markov Chain Monte Carlo method. When the transition density of the SDE is explicit, we prove the convergence of the SAEM-PMCMC algorithm towards the maximum likelihood estimator. Two simulated examples are considered: an Ornstein-Uhlenbeck process with two random parameters and a time-inhomogeneous SDE (Gompertz SDE) with a stochastic volatility error model and three random parameters. When the transition density is unknown, we prove the convergence of a different version of the algorithm based on the Euler approximation of the SDE towards the maximum likelihood estimator.
منابع مشابه
Using PMCMC in EM algorithm for stochastic mixed models: theoretical and practical issues
Biological processes measured repeatedly among a series of individuals are standardly analyzed by mixed models. Recently, stochastic processes have been introduced to model the variability along time for each subject. Although the likelihood of these stochastic mixed models is intractable, various estimation methods have been proposed when the latent stochastic process is a discrete time finite...
متن کاملThe Development of Maximum Likelihood Estimation Approaches for Adaptive Estimation of Free Speed and Critical Density in Vehicle Freeways
The performance of many traffic control strategies depends on how much the traffic flow models have been accurately calibrated. One of the most applicable traffic flow model in traffic control and management is LWR or METANET model. Practically, key parameters in LWR model, including free flow speed and critical density, are parameterized using flow and speed measurements gathered by inductive ...
متن کاملA Regularized Particle Filter EM Algorithm Based on Gaussian Randomization with an Application to Plant Growth Modeling
Parameter estimation in complex models arising in real data applications is a topic which still attracts a lot of interest. In this article, we study a specific data and parameter augmentation method which gives us the opportunity to estimate more easily the parameters of the initial model. For this reason, the notion of Gaussian randomization of a model with respect to some of its parameters i...
متن کاملParameter Estimation in Spatial Generalized Linear Mixed Models with Skew Gaussian Random Effects using Laplace Approximation
Spatial generalized linear mixed models are used commonly for modelling non-Gaussian discrete spatial responses. We present an algorithm for parameter estimation of the models using Laplace approximation of likelihood function. In these models, the spatial correlation structure of data is carried out by random effects or latent variables. In most spatial analysis, it is assumed that rando...
متن کاملParameter Estimation in Stochastic Differential Mixed-Effects Models
Stochastic differential equation (SDE) models have shown useful to describe continuous time processes, e.g. a physiological process evolving in an individual. Biomedical experiments often imply repeated measurements on a series of individuals or experimental units and individual differences can be represented by incorporating random effects into the model. When both system noise and individual ...
متن کامل