Effectiveness of sparse Bayesian algorithm for MVAR coefficient estimation in MEG/EEG source-space causality analysis

This paper examines the effectiveness of a sparse Bayesian algorithm to estimate multivariate autoregressive coefficients when a large amount of background interference exists. This paper employs computer experiments to compare two methods in the source-space causality analysis: the conventional least-squares method and a sparse Bayesian method. Results of our computer experiments show that the interference affects the least-squares method in a very severe manner. It produces large false-positive results, unless the signal-to-interference ratio is very high. On the other hand, the sparse Bayesian method is relatively insensitive to the existence of interference. However, this robustness of the sparse Bayesian method is attained on the scarifies of the detectability of true causal relationship. Our experiments also show that the surrogate data bootstrapping method tends to give a statistical threshold that are too low for the sparse method. The permutation-test-based method gives a higher (more conservative) threshold and it should be used with the sparse Bayesian method whenever the control period is available.

[1]  Julia P. Owen,et al.  Robust Bayesian estimation of the location, orientation, and time course of multiple correlated neural sources using MEG , 2010, NeuroImage.

[2]  Giulio Tononi,et al.  State-Space Multivariate Autoregressive Models for Estimation of Cortical Connectivity from EEG , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

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

[4]  V. Solo,et al.  Dynamic Granger–Geweke causality modeling with application to interictal spike propagation , 2009, Human brain mapping.

[5]  S. Roberts,et al.  Bayesian multivariate autoregressive models with structured priors , 2002 .

[6]  J. Geweke,et al.  Measurement of Linear Dependence and Feedback between Multiple Time Series , 1982 .

[7]  J. Geweke,et al.  Measures of Conditional Linear Dependence and Feedback between Time Series , 1984 .

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

[9]  S. Bressler,et al.  Granger Causality: Basic Theory and Application to Neuroscience , 2006, q-bio/0608035.

[10]  Markus Neuhäuser,et al.  Permutation Tests , 2011, International Encyclopedia of Statistical Science.

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

[12]  James Theiler,et al.  Testing for nonlinearity in time series: the method of surrogate data , 1992 .

[13]  Katarzyna J. Blinowska,et al.  A new method of the description of the information flow in the brain structures , 1991, Biological Cybernetics.