State-Space Multivariate Autoregressive Models for Estimation of Cortical Connectivity from EEG

We propose using a state-space model to estimate cortical connectivity from scalp-based EEG recordings. A state equation describes the dynamics of the cortical signals and an observation equation describes the manner in which the cortical signals contribute to the scalp measurements. The state equation is based on a multivariate autoregressive (MVAR) process model for the cortical signals. The observation equation describes the physics relating the cortical signals to the scalp EEG measurements and spatially correlated observation noise. An expectation-maximization (EM) algorithm is employed to obtain maximum-likelihood estimates of the MVAR model parameters. The strength of influence between cortical regions is then derived from the MVAR model parameters. Simulation results show that this integrated approach performs significantly better than the two-step approach in which the cortical signals are first estimated from the EEG measurements by attempting to solve the EEG inverse problem and second, an MVAR model is fit to the estimated signals. The method is also applied to data from a subject watching a movie, and confirms that feedforward connections between visual and parietal cortex are generally stronger than feedback connections.

[1]  Katarzyna J. Blinowska,et al.  Determination of EEG activity propagation: pair-wise versus multichannel estimate , 2004, IEEE Transactions on Biomedical Engineering.

[2]  S. Bressler,et al.  Beta oscillations in a large-scale sensorimotor cortical network: directional influences revealed by Granger causality. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Luiz A. Baccalá,et al.  Partial directed coherence: a new concept in neural structure determination , 2001, Biological Cybernetics.

[4]  S. Kosslyn,et al.  Brain areas underlying visual mental imagery and visual perception: an fMRI study. , 2004, Brain research. Cognitive brain research.

[5]  Barry D. Van Veen,et al.  Estimation of cortical multivariate autoregressive models for EEG/MEG using an expectation-maximization algorithm , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[6]  Karin Schwab,et al.  Comparison of linear signal processing techniques to infer directed interactions in multivariate neural systems , 2005, Signal Process..

[7]  C. Striebel,et al.  On the maximum likelihood estimates for linear dynamic systems , 1965 .

[8]  Gregory A. Worrell,et al.  Ictal source analysis: Localization and imaging of causal interactions in humans , 2007, NeuroImage.

[9]  Hualou Liang,et al.  Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment , 2000, Biological Cybernetics.

[10]  Mingzhou Ding,et al.  Evaluating causal relations in neural systems: Granger causality, directed transfer function and statistical assessment of significance , 2001, Biological Cybernetics.

[11]  Barry D. Van Veen,et al.  Cortical patch basis model for spatially extended neural activity , 2006, IEEE Transactions on Biomedical Engineering.

[12]  Helmut Lütkepohl,et al.  Introduction to multiple time series analysis , 1991 .

[13]  Karl J. Friston,et al.  Where bottom-up meets top-down: neuronal interactions during perception and imagery. , 2004, Cerebral cortex.

[14]  Ian T. Jolliffe,et al.  Introduction to Multiple Time Series Analysis , 1993 .

[15]  G. A. Miller,et al.  Comparison of different cortical connectivity estimators for high‐resolution EEG recordings , 2007, Human brain mapping.

[16]  Richard M. Leahy,et al.  Linearly constrained MEG beamformers for MVAR modeling of cortical interactions , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[17]  R. Shumway,et al.  AN APPROACH TO TIME SERIES SMOOTHING AND FORECASTING USING THE EM ALGORITHM , 1982 .