A continuous-time adaptive particle filter for estimations under measurement time uncertainties with an application to a plasma-leucine mixed effects model

BackgroundWhen mathematical modelling is applied to many different application areas, a common task is the estimation of states and parameters based on measurements. With this kind of inference making, uncertainties in the time when the measurements have been taken are often neglected, but especially in applications taken from the life sciences, this kind of errors can considerably influence the estimation results. As an example in the context of personalized medicine, the model-based assessment of the effectiveness of drugs is becoming to play an important role. Systems biology may help here by providing good pharmacokinetic and pharmacodynamic (PK/PD) models. Inference on these systems based on data gained from clinical studies with several patient groups becomes a major challenge. Particle filters are a promising approach to tackle these difficulties but are by itself not ready to handle uncertainties in measurement times.ResultsIn this article, we describe a variant of the standard particle filter (PF) algorithm which allows state and parameter estimation with the inclusion of measurement time uncertainties (MTU). The modified particle filter, which we call MTU-PF, also allows the application of an adaptive stepsize choice in the time-continuous case to avoid degeneracy problems. The modification is based on the model assumption of uncertain measurement times. While the assumption of randomness in the measurements themselves is common, the corresponding measurement times are generally taken as deterministic and exactly known. Especially in cases where the data are gained from measurements on blood or tissue samples, a relatively high uncertainty in the true measurement time seems to be a natural assumption. Our method is appropriate in cases where relatively few data are used from a relatively large number of groups or individuals, which introduce mixed effects in the model. This is a typical setting of clinical studies. We demonstrate the method on a small artificial example and apply it to a mixed effects model of plasma-leucine kinetics with data from a clinical study which included 34 patients.ConclusionsComparisons of our MTU-PF with the standard PF and with an alternative Maximum Likelihood estimation method on the small artificial example clearly show that the MTU-PF obtains better estimations. Considering the application to the data from the clinical study, the MTU-PF shows a similar performance with respect to the quality of estimated parameters compared with the standard particle filter, but besides that, the MTU algorithm shows to be less prone to degeneration than the standard particle filter.

[1]  A. Doucet,et al.  Sequential MCMC for Bayesian model selection , 1999, Proceedings of the IEEE Signal Processing Workshop on Higher-Order Statistics. SPW-HOS '99.

[2]  Torsten P. Bohlin,et al.  Practical Grey-box Process Identification: Theory and Applications , 2006 .

[3]  P. Fearnhead.  Markov chain Monte Carlo, Sufficient Statistics, and Particle Filters , 2002 .

[4]  Jon Wakefield,et al.  Statistical methods for population pharmacokinetic modelling , 1998, Statistical methods in medical research.

[5]  A. Doucet,et al.  Particle Markov chain Monte Carlo methods , 2010 .

[6]  Arnaud Doucet,et al.  A survey of convergence results on particle filtering methods for practitioners , 2002, IEEE Trans. Signal Process..

[7]  Fredrik Gustafsson,et al.  On Resampling Algorithms for Particle Filters , 2006, 2006 IEEE Nonlinear Statistical Signal Processing Workshop.

[8]  MartinAdiels,et al.  Overproduction of Very Low–Density Lipoproteins Is the Hallmark of the Dyslipidemia in the Metabolic Syndrome , 2008 .

[9]  Martin Adiels,et al.  Overproduction of Very Low–Density Lipoproteins Is the Hallmark of the Dyslipidemia in the Metabolic Syndrome , 2008, Arteriosclerosis, thrombosis, and vascular biology.

[10]  Nando de Freitas,et al.  Sequential Monte Carlo Methods in Practice , 2001, Statistics for Engineering and Information Science.

[11]  Sophie Donnet,et al.  EM algorithm coupled with particle filter for maximum likelihood parameter estimation of stochastic differential mixed-effects models , 2010 .

[12]  Mikael Sunnåker,et al.  Investigations of a compartmental model for leucine kinetics using non-linear mixed effects models with ordinary and stochastic differential equations. , 2011, Mathematical medicine and biology : a journal of the IMA.

[13]  Timothy J. Robinson,et al.  Sequential Monte Carlo Methods in Practice , 2003 .

[14]  Martin Adiels.  A compartmental model for kinetics of Apolipoprotein B-100 and Triglycerides in VLDL_1 and VLDL_2 in normolipidemic subjects , 2002 .

[15]  Henrik Madsen,et al.  Stochastic Differential Equations in NONMEM®: Implementation, Application, and Comparison with Ordinary Differential Equations , 2005, Pharmaceutical Research.

[16]  D. Talay.  Numerical solution of stochastic differential equations , 1994 .

[17]  Jon Wakefield,et al.  Population modelling in drug development , 1999, Statistical methods in medical research.

[18]  M. Pitt,et al.  Filtering via Simulation: Auxiliary Particle Filters , 1999 .

[19]  Maksat Ashyraliyev,et al.  Systems biology: parameter estimation for biochemical models , 2009, The FEBS journal.

[20]  Nando de Freitas,et al.  Rao-Blackwellised Particle Filtering for Dynamic Bayesian Networks , 2000, UAI.

[21]  Bernt Wennberg,et al.  A new combined multicompartmental model for apolipoprotein B-100 and triglyceride metabolism in VLDL subfractions Published, JLR Papers in Press, October 16, 2004. DOI 10.1194/jlr.M400108-JLR200 , 2005, Journal of Lipid Research.

[22]  Bernt Wennberg,et al.  Overproduction of VLDL1 Driven by Hyperglycemia Is a Dominant Feature of Diabetic Dyslipidemia , 2005, Arteriosclerosis, thrombosis, and vascular biology.

[23]  Simon J. Godsill,et al.  An Overview of Existing Methods and Recent Advances in Sequential Monte Carlo , 2007, Proceedings of the IEEE.

[24]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[25]  Sophie Donnet,et al.  Parametric inference for mixed models defined by stochastic differential equations , 2008 .

[26]  N. Gordon,et al.  Novel approach to nonlinear/non-Gaussian Bayesian state estimation , 1993 .

[27]  Jan Nygaard Nielsen,et al.  Parameter estimation in stochastic differential equations: An overview , 2000 .

[28]  Eric Moulines,et al.  Comparison of resampling schemes for particle filtering , 2005, ISPA 2005. Proceedings of the 4th International Symposium on Image and Signal Processing and Analysis, 2005..

[29]  J Shepherd,et al.  Sensitive methods to study human apolipoprotein B metabolism using stable isotope-labeled amino acids. , 1996, The American journal of physiology.

[30]  C Cobelli,et al.  Compartmental model of leucine kinetics in humans. , 1991, The American journal of physiology.

[31]  Geir Storvik,et al.  Particle filters for state-space models with the presence of unknown static parameters , 2002, IEEE Trans. Signal Process..

[32]  H. Madsen,et al.  Non-Linear Mixed-Effects Models with Stochastic Differential Equations: Implementation of an Estimation Algorithm , 2005, Journal of Pharmacokinetics and Pharmacodynamics.

[33]  Malene Højbjerre,et al.  A population-based Bayesian approach to the minimal model of glucose and insulin homeostasis. , 2005, Statistics in medicine.

[34]  Xiao-Li Hu,et al.  A Basic Convergence Result for Particle Filtering , 2008, IEEE Transactions on Signal Processing.

[35]  I. Chou,et al.  Recent developments in parameter estimation and structure identification of biochemical and genomic systems. , 2009, Mathematical biosciences.

[36]  P. Moral,et al.  Sequential Monte Carlo samplers , 2002, cond-mat/0212648.

[37]  P. Fearnhead.  MCMC, sufficient statistics and particle filters. , 2002 .