The generalized shrinkage estimator for the analysis of functional connectivity of brain signals

Abstract

We develop a new statistical method for estimating functional connectivity between neurophysiological signals represented by a multivariate time series. We use partial coherence as the measure of functional connectivity. Partial coherence identifies the frequency bands that drive the direct linear association between any pair of channels. To estimate partial coherence, one would first need an estimate of the spectral density matrix of the multivariate time series. Parametric estimators of the spectral density matrix provide good frequency resolution but could be sensitive when the parametric model is misspecified. Smoothing-based nonparametric estimators are robust to model misspecification and are consistent but may have poor frequency resolution. In this work, we develop the generalized shrinkage estimator, which is a weighted average of a parametric estimator and a nonparametric estimator. The optimal weights are frequency-specific and derived under the quadratic risk criterion so that the estimator, either the parametric estimator or the nonparametric estimator, that performs better at a particular frequency receives heavier weight. We validate the proposed estimator in a simulation study and apply it on electroencephalogram recordings from a visual-motor experiment.

Published in at http://dx.doi.org/10.1214/10-AOAS396 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (1)

• Structural shrinkage of nonparametric spectral estimators for multivariate time series

Cited by in corpus (15)

• Analyzing complex functional brain networks: fusing statistics and network science to understand the brain
• Graph Autoencoders for Embedding Learning in Brain Networks and Major Depressive Disorder Identification
• An Efficient Approach to Graphical Modelling of Time Series
• Multi-Scale Factor Analysis of High-Dimensional Brain Signals
• Estimating whole brain dynamics using spectral clustering
• Random Matrix Derived Shrinkage of Spectral Precision Matrices
• Spectral Analysis of High-dimensional Time Series
• Modeling Brain Connectivity with Graphical Models on Frequency Domain
• A Multi-Resolution Spatio-Temporal Model for Brain Activation and Connectivity in fMRI Data
• Statistical Challenges in Modeling Big Brain Signals
• Sparse Symmetric Tensor Regression for Functional Connectivity Analysis
• Exploratory Analysis of High Dimensional Time Series with Applications to Multichannel Electroencephalograms
• Modeling high dimensional multichannel brain signals
• Smooth Online Parameter Estimation for time varying VAR models with application to rat's LFP data
• Learning Multi-Frequency Partial Correlation Graphs