Bayesian smoothing of dipoles in magneto-/electroencephalography

We describe a novel method for dynamic estimation of multi-dipole states from magneto-/electroencephalography (M/EEG) time series. The new approach builds on the recent development of particle filters for M/EEG; these algorithms approximate, with samples and weights, the posterior distribution of the neural sources at time t given the data up to time t. However, for off-line inference purposes it is preferable to work with the smoothing distribution, i.e. the distribution for the neural sources at time t conditioned on the whole time series. In this study, we use a Monte Carlo algorithm to approximate the smoothing distribution for a time-varying set of current dipoles. We show, using numerical simulations, that the estimates provided by the smoothing distribution are more accurate than those provided by the filtering distribution, particularly at the appearance of the source. We validate the proposed algorithm using an experimental data set recorded from an epileptic patient. Improved localization of the source onset can be particularly relevant in source modeling of epileptic patients, where the source onset brings information on the epileptogenic zone.

[1]  E. Somersalo,et al.  Non-stationary magnetoencephalography by Bayesian filtering of dipole models , 2003 .

[2]  Alberto Sorrentino,et al.  Bayesian multi-dipole modelling of a single topography in MEG by adaptive sequential Monte Carlo samplers , 2013, 1305.4511.

[3]  L. Dousset.  Understanding Human Relations (Kinship Systems) , 2011 .

[4]  Emery N. Brown,et al.  A spatiotemporal dynamic distributed solution to the MEG inverse problem , 2011, NeuroImage.

[5]  Polina Golland,et al.  A distributed spatio-temporal EEG/MEG inverse solver , 2009, NeuroImage.

[6]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[7]  Thomas E. Nichols,et al.  Dynamic filtering of static dipoles in magnetoencephalography , 2013 .

[8]  Chaitali Chakrabarti,et al.  Multi-source Neural Activity Estimation and Sensor Scheduling: Algorithms and Hardware Implementation , 2013, J. Signal Process. Syst..

[9]  Simona Temereanca,et al.  STATE-SPACE SOLUTIONS TO THE DYNAMIC MAGNETOENCEPHALOGRAPHY INVERSE PROBLEM USING HIGH PERFORMANCE COMPUTING. , 2011, The annals of applied statistics.

[10]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[11]  Eric Moulines,et al.  Inference in hidden Markov models , 2010, Springer series in statistics.

[12]  P. Fearnhead,et al.  A sequential smoothing algorithm with linear computational cost. , 2010 .

[13]  Xi Chen,et al.  A Bayesian particle filtering method for brain source localisation , 2014, Digit. Signal Process..

[14]  A. Doucet,et al.  Smoothing algorithms for state–space models , 2010 .

[15]  Simon J. Godsill,et al.  On sequential Monte Carlo sampling methods for Bayesian filtering , 2000, Stat. Comput..

[16]  A. Doucet,et al.  A Tutorial on Particle Filtering and Smoothing: Fifteen years later , 2008 .

[17]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[18]  R. Salmelin,et al.  Global optimization in the localization of neuromagnetic sources , 1998, IEEE Transactions on Biomedical Engineering.

[19]  A. Gramfort,et al.  Mixed-norm estimates for the M/EEG inverse problem using accelerated gradient methods , 2012, Physics in medicine and biology.

[20]  Lauri Parkkonen,et al.  Particle filters and RAP-MUSIC in MEG source modelling: A comparison , 2007 .

[21]  S. Sanei,et al.  Sequential monte carlo techniques for EEG dipole placing and tracking , 2008, 2008 5th IEEE Sensor Array and Multichannel Signal Processing Workshop.

[22]  Ba-Ngu Vo,et al.  A Consistent Metric for Performance Evaluation of Multi-Object Filters , 2008, IEEE Transactions on Signal Processing.

[23]  Lauri Parkkonen,et al.  Dynamical MEG source modeling with multi‐target Bayesian filtering , 2009, Human brain mapping.

[24]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[25]  J. Sarvas.  Basic mathematical and electromagnetic concepts of the biomagnetic inverse problem. , 1987, Physics in medicine and biology.

[26]  Haipeng Shen,et al.  Two‐way regularization for MEG source reconstruction via multilevel coordinate descent , 2013, Stat. Anal. Data Min..

[27]  Simon J. Godsill,et al.  Rao-Blackwellized Particle Smoothers for Conditionally Linear Gaussian Models , 2015, IEEE Journal of Selected Topics in Signal Processing.

[28]  D. M. Schmidt,et al.  Bayesian inference applied to the electromagnetic inverse problem , 1998, Human brain mapping.

[29]  Alberto Sorrentino,et al.  Sequential Monte Carlo samplers for semi-linear inverse problems and application to magnetoencephalography , 2014, 1409.8109.

[30]  Lauri Parkkonen,et al.  Particle filters: a new method for reconstructing multiple current dipoles from meg data , 2007 .

[31]  Nando de Freitas,et al.  Fast particle smoothing: if I had a million particles , 2006, ICML.