Self-Consistent MUSIC: An approach to the localization of true brain interactions from EEG/MEG data

MUltiple SIgnal Classification (MUSIC) is a standard localization method which is based on the idea of dividing the vector space of the data into two subspaces: signal subspace and noise subspace. The brain, divided into several grid points, is scanned entirely and the grid point with the maximum consistency with the signal subspace is considered as the source location. In one of the MUSIC variants called Recursively Applied and Projected MUSIC (RAP-MUSIC), multiple iterations are proposed in order to decrease the location estimation uncertainties introduced by subspace estimation errors. In this paper, we suggest a new method called Self-Consistent MUSIC (SC-MUSIC) which extends RAP-MUSIC to a self-consistent algorithm. This method, SC-MUSIC, is based on the idea that the presence of several sources has a bias on the localization of each source. This bias can be reduced by projecting out all other sources mutually rather than iteratively. While the new method is applicable in all situations when MUSIC is applicable we will study here the localization of interacting sources using the imaginary part of the cross-spectrum due to the robustness of this measure to the artifacts of volume conduction. For an odd number of sources this matrix is rank deficient similar to covariance matrices of fully correlated sources. In such cases MUSIC and RAP-MUSIC fail completely while the new method accurately localizes all sources. We present results of the method using simulations of odd and even number of interacting sources in the presence of different noise levels. We compare the method with three other source localization methods: RAP-MUSIC, dipole fit and MOCA (combined with minimum norm estimate) through simulations. SC-MUSIC shows substantial improvement in the localization accuracy compared to these methods. We also show results for real MEG data of a single subject in the resting state. Four sources are localized in the sensorimotor area at f=11Hz which is the expected region for the idle rhythm.

[1]  G. Nolte,et al.  Analytic expansion of the EEG lead field for realistic volume conductors , 2005, Physics in medicine and biology.

[2]  Guido Nolte,et al.  Understanding brain connectivity from EEG data by identifying systems composed of interacting sources , 2008, NeuroImage.

[3]  Andreas Ziehe,et al.  Constructing Surrogate Data to Control for Artifacts of Volume Conduction for Functional Connectivity Measures , 2010 .

[4]  Richard M. Leahy,et al.  Source localization using recursively applied and projected (RAP) MUSIC , 1997, Conference Record of the Thirty-First Asilomar Conference on Signals, Systems and Computers (Cat. No.97CB36136).

[5]  R. Ilmoniemi,et al.  Interpreting magnetic fields of the brain: minimum norm estimates , 2006, Medical and Biological Engineering and Computing.

[6]  Richard M. Leahy,et al.  Identifying true cortical interactions in MEG using the nulling beamformer , 2010, NeuroImage.

[7]  D. Tucker,et al.  EEG coherency II: experimental comparisons of multiple measures , 1999, Clinical Neurophysiology.

[8]  Wolf Singer,et al.  Neuronal Synchrony: A Versatile Code for the Definition of Relations? , 1999, Neuron.

[9]  F. Varela,et al.  Measuring phase synchrony in brain signals , 1999, Human brain mapping.

[10]  K. Müller,et al.  Robustly estimating the flow direction of information in complex physical systems. , 2007, Physical review letters.

[11]  M. V. D. Heuvel,et al.  Exploring the brain network: A review on resting-state fMRI functional connectivity , 2010, European Neuropsychopharmacology.

[12]  R. Desimone,et al.  Modulation of Oscillatory Neuronal Synchronization by Selective Visual Attention , 2001, Science.

[13]  J.C. Mosher,et al.  Multiple dipole modeling and localization from spatio-temporal MEG data , 1992, IEEE Transactions on Biomedical Engineering.

[14]  R. O. Schmidt,et al.  Multiple emitter location and signal Parameter estimation , 1986 .

[15]  Petre Stoica,et al.  MUSIC, maximum likelihood, and Cramer-Rao bound , 1989, IEEE Transactions on Acoustics, Speech, and Signal Processing.

[16]  T. Holroyd,et al.  Localizing brain interactions from rhythmic EEG/MEG data , 2004, The 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[17]  Heidi E Kirsch,et al.  Resting functional connectivity in patients with brain tumors in eloquent areas , 2011, Annals of neurology.

[18]  James P. Reilly,et al.  Detection of the number of signals: a predicted eigen-threshold approach , 1991, IEEE Trans. Signal Process..

[19]  M. Hallett,et al.  Identifying true brain interaction from EEG data using the imaginary part of coherency , 2004, Clinical Neurophysiology.

[20]  J. Martinerie,et al.  The brainweb: Phase synchronization and large-scale integration , 2001, Nature Reviews Neuroscience.

[21]  Aapo Hyvärinen,et al.  Independent Component Analysis: A Tutorial , 1999 .

[22]  W. Singer,et al.  Dynamic predictions: Oscillations and synchrony in top–down processing , 2001, Nature Reviews Neuroscience.

[23]  Matthew J. Brookes,et al.  Measuring functional connectivity using MEG: Methodology and comparison with fcMRI , 2011, NeuroImage.

[24]  Ernesto Pereda,et al.  Assessment of electroencephalographic functional connectivity in term and preterm neonates , 2011, Clinical Neurophysiology.

[25]  Guido Nolte,et al.  Wedge MUSIC: A novel approach to examine experimental differences of brain source connectivity patterns from EEG/MEG data , 2014, NeuroImage.

[26]  Andreas Ziehe,et al.  TDSEP { an e(cid:14)cient algorithm for blind separation using time structure , 1998 .

[27]  P. Fries.  Neuronal gamma-band synchronization as a fundamental process in cortical computation. , 2009, Annual review of neuroscience.

[28]  Julia P. Owen,et al.  Removal of Spurious Coherence in MEG Source-Space Coherence Analysis , 2011, IEEE Transactions on Biomedical Engineering.

[29]  Eric Maris.  A resampling method for estimating the signal subspace of spatio-temporal EEG/MEG data , 2003, IEEE Transactions on Biomedical Engineering.

[30]  Steven L. Bressler,et al.  Wiener–Granger Causality: A well established methodology , 2011, NeuroImage.

[31]  Guido Nolte,et al.  Identifying causal networks of neuronal sources from EEG/MEG data with the phase slope index: a simulation study , 2013, Biomedizinische Technik. Biomedical engineering.

[32]  James V. Stone.  Independent Component Analysis: A Tutorial Introduction , 2007 .

[33]  Guido Nolte,et al.  Localizing True Brain Interactions from EEG and MEG Data with Subspace Methods and Modified Beamformers , 2012, Comput. Math. Methods Medicine.

[34]  M. Berger,et al.  Mapping functional connectivity in patients with brain lesions , 2008, Annals of neurology.

[35]  J. Schoffelen,et al.  Source connectivity analysis with MEG and EEG , 2009, Human brain mapping.

[36]  G. Nolte.  The magnetic lead field theorem in the quasi-static approximation and its use for magnetoencephalography forward calculation in realistic volume conductors. , 2003, Physics in medicine and biology.

[37]  P. Fries.  A mechanism for cognitive dynamics: neuronal communication through neuronal coherence , 2005, Trends in Cognitive Sciences.

[38]  Karen O. Egiazarian,et al.  Measuring directional coupling between EEG sources , 2008, NeuroImage.

[39]  K. Miller,et al.  Direct electrophysiological measurement of human default network areas , 2009, Proceedings of the National Academy of Sciences.

[40]  Huawei Chen,et al.  Coherent signal-subspace processing of acoustic vector sensor array for DOA estimation of wideband sources , 2005, Signal Process..

[41]  Guido Nolte,et al.  Estimating true brain connectivity from EEG/MEG data invariant to linear and static transformations in sensor space , 2012, NeuroImage.

[42]  Hong Wang,et al.  Coherent signal-subspace processing for the detection and estimation of angles of arrival of multiple wide-band sources , 1985, IEEE Trans. Acoust. Speech Signal Process..

[43]  C. Almli,et al.  Unbiased nonlinear average age-appropriate brain templates from birth to adulthood , 2009, NeuroImage.

[44]  Erkki Oja,et al.  Independent component analysis: algorithms and applications , 2000, Neural Networks.

[45]  M. Peters,et al.  The volume conductor may act as a temporal filter on the ECG and EEG , 1998, Medical and Biological Engineering and Computing.