Spectral Differential Network Analysis for High-Dimensional Time Series

2024-12-10

AI Summary

Conventional methods for estimating inter-condition network differences in high-dimensional multivariate time series—such as pre- and post-seizure EEG—estimate individual spectral networks separately and subtract them, thereby amplifying estimation error. Performance degrades significantly when the underlying networks are non-sparse but their difference is sparse. Method: We propose a novel frequency-domain approach that directly estimates the difference of inverse spectral density matrices. Our method integrates spectral density modeling, inverse covariance estimation, and frequency-domain statistical inference. Contribution/Results: We introduce the first theoretical framework formalizing “difference sparsity,” imposing L1 regularization solely on the differential structure—bypassing restrictive sparsity assumptions on the original high-dimensional networks. Evaluated on synthetic data and real EEG/micro-ECoG recordings, our method achieves superior consistency, stability, and detection accuracy in differential network estimation, providing an interpretable and robust tool for dynamic brain connectivity analysis.

Technology Category

• Cognitive Modeling & Cognitive Systems: Neural Spike Coding
• Machine Learning: Matrix & Tensor Methods
• Humans and AI: Brain-Sensing and Analysis

Application Category

• Graph Algorithms and Modeling for the Web: Representation, reconstruction, and subgraph or motif discovery in Web-related graphs
• Social Networks and Social Media: Influence propagation, information diffusion, and the prediction on networks
• Web Mining and Content Analysis: Content-based information diffusion

Abstract

Spectral networks derived from multivariate time series data arise in many domains, from brain science to Earth science. Often, it is of interest to study how these networks change under different conditions. For instance, to better understand epilepsy, it would be interesting to capture the changes in the brain connectivity network as a patient experiences a seizure, using electroencephalography data. A common approach relies on estimating the networks in each condition and calculating their difference. Such estimates may behave poorly in high dimensions as the networks themselves may not be sparse in structure while their difference may be. We build upon this observation to develop an estimator of the difference in inverse spectral densities across two conditions. Using an L1 penalty on the difference, consistency is established by only requiring the difference to be sparse. We illustrate the method on synthetic data experiments, on experiments with electroencephalography data, and on experiments with optogentic stimulation and micro-electrocorticography data.

Problem

• Estimating spectral network differences in high-dimensional time series
• Capturing sparse changes in brain connectivity during seizures
• Developing L1-penalized inverse spectral density difference estimator

Innovation

• Estimates inverse spectral density differences
• Uses L1 penalty for sparse differences
• Validated with synthetic and EEG data