Causality analysis between brain areas based on multivariate autoregressive models of MEG sensor data
© Michalareas et al; licensee BioMed Central Ltd. 2009
Published: 13 July 2009
In recent years, significant efforts have been aimed towards the development of methods for identifying causality in information flow between brain areas from MEG/EEG data. As such multichannel time-series correspond to inherently multivariate processes, bivariate causality measures such as the Granger Causality Test can only be inconsistently applied. Two recent methods, the Directed Transfer Function (DTF) and the Partial Directed Coherence (PDC) , have been widely used in intra-brain connectivity studies. These methods study the frequency characteristics of multivariate autoregressive models (MVAR) built on time-series of activated brain sources. These time-series are reconstructed from MEG/EEG sensor data through spatial filters derived commonly by beamforming algorithms . The large number of potential activation sources produced by such algorithms corresponds to a large number of activation time series, too many for the derivation of robust multivariate autoregressive models. Typical methods of overcoming this obstacle make a priori assumptions about certain brain locations of potential activity or eliminate areas with low neural activity index (NAI). This is a significant drawback, as neglected areas or areas with low NAI might contain important functional connectivity paths and causality information. In the main part of this work, we investigate the derivation of the multivariate autoregressive model directly on MEG sensor data and subsequently its projection in the source space within the brain through the spatial filters. Here, the modeling process is performed on the sensor space which has moderate dimensionality as compared to the high-dimensional source space. This leads to greater model robustness as well as significantly reduced computation times. Simulated brain activation signals were used and MEG pseudo-data was produced from the model of a 4D-Neuroimaging 248 channel system. Activation signals were modeled as synaptic potentials in 10,20 and 40 Hz with added white process noise.
The investigation of causality analysis through multivariate autoregressive models, built directly in the sensor data space, showed that not only is it a feasible process but also it offers the ability to perform causality analysis between a large number of sources in a computationally efficient manner. These promising results suggest the use of the proposed method for effective derivation of high-dimensional models in the source space that can then in turn be used for exhaustive search for causal effects in functional brain connections. Further improvements are investigated for increasing the efficiency of such a search. One under validation is the reduction of the source space through the use of bivariate causality measures which, although inaccurate, can potentially provide a coarse topology of the activated areas which can then be isolated and investigated more in detail through DTF and PDC.
- Gourévitch B, Bouquin-Jeannes RL, Faucon G: Linear and nonlinear causality between signals: methods, examples and neurophysiological applications. Biological Cybernetics. 2006, 95: 349-369. 10.1007/s00422-006-0098-0.PubMedView ArticleGoogle Scholar
- Hui HB, Leahy RM: Linearly constrained beamformers for MVAR modelling of cortical interactions. Biomedical Imaging: Nano to Macro. 2006, 237-240. doi: 10.1109/ISBI.2006.1624896.Google Scholar
- Baccalá LA, Sameshima K: Partial Directed Coherence: a new concept in neural structure determination. Biological Cybernetics. 2001, 84: 435-439. 10.1007/PL00007990.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd.