Statistics and Its Interface

Volume 6 (2013)

Number 4

Coupling the SAEM algorithm and the extended Kalman filter for maximum likelihood estimation in mixed-effects diffusion models

Pages: 519 – 532



Maud Delattre (Inria Saclay & University Paris Sud Orsay & AgroParisTech, UMR 518 MIA, Paris)

Marc Lavielle (Inria Saclay & University Paris Sud Orsay, France)


We consider some general mixed-effects diffusion models, in which the observations are made at discrete time points and include measurement errors. In these models, the observed likelihood is generally not explicit, making maximum likelihood estimation of the parameters particularly complex. We propose a specific inference methodology for these models. In particular, we combine the SAEM algorithm with the extended Kalman filter to estimate the population parameters. We also provide some tools for estimating the individual parameters, for recovering the individual underlying diffusion trajectories and for evaluating the model. The methods are evaluated on simulations and applied to a pharmacokinetics example.


stochastic differential equations, mixed-effects models, SAEM, extended Kalman filter

Full Text (PDF format)