An identifiable model to assess frequency-domain granger causality in the presence of significant instantaneous interactions

We present a new approach for the investigation of Granger causality in the frequency domain by means of the partial directed coherence (PDC). The approach is based on the utilization of an extended multivariate autoregressive (MVAR) model, including instantaneous effects in addition to the lagged effects traditionally studied, to fit the observed multiple time series prior to PDC computation. Model identification is performed combining standard MVAR coefficient estimation with a recent technique for instantaneous causal modeling based on independent component analysis. The approach is first validated on simulated MVAR processes showing that, in the presence of instantaneous effects, only the extended model is able to interpret the imposed Granger causality patterns, while the traditional MVAR approach may yield strongly biased PDC estimates. The subsequent application to multichannel EEG time series confirms the potentiality of the approach in real data applications, as the importance of instantaneous effects led to significant differences in the PDC estimated after traditional and extended MVAR identification.

[1]  C. Granger Investigating causal relations by econometric models and cross-spectral methods , 1969 .

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

[3]  Luca Faes,et al.  Multivariate Autoregressive Model with Instantaneous Effects to Improve Brain Connectivity Estimation , 2009 .

[4]  M. Kaminski,et al.  Topographic analysis of coherence and propagation of EEG activity during sleep and wakefulness. , 1997, Electroencephalography and clinical neurophysiology.

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

[6]  Aapo Hyvärinen,et al.  Causal modelling combining instantaneous and lagged effects: an identifiable model based on non-Gaussianity , 2008, ICML '08.

[7]  L Faes,et al.  Assessing frequency domain causality in cardiovascular time series with instantaneous interactions. , 2010, Methods of information in medicine.

[8]  Aapo Hyvärinen,et al.  A Linear Non-Gaussian Acyclic Model for Causal Discovery , 2006, J. Mach. Learn. Res..

[9]  L.A. Baccald,et al.  Generalized Partial Directed Coherence , 2007, 2007 15th International Conference on Digital Signal Processing.

[10]  Aapo Hyvärinen,et al.  Fast and robust fixed-point algorithms for independent component analysis , 1999, IEEE Trans. Neural Networks.

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

[12]  H. Akaike A new look at the statistical model identification , 1974 .

[13]  Luca Faes,et al.  Testing Frequency-Domain Causality in Multivariate Time Series , 2010, IEEE Transactions on Biomedical Engineering.

[14]  G. Baselli,et al.  Spectral decomposition in multichannel recordings based on multivariate parametric identification , 1997, IEEE Transactions on Biomedical Engineering.