Multitaper Analysis of Evolutionary Spectral Density Matrix From Multivariate Spiking Observations

Extracting the spectral representations of the neural processes that underlie spiking activity is key to understanding how the brain rhythms mediate cognitive functions. While spectral estimation of continuous time-series is well studied, inferring the spectral representation of latent processes from spiking observations is a challenging problem. In this paper, we address this issue by proposing a spectral estimation methodology that can be directly applied to multivariate spiking observations in order to extract the evolutionary spectral representation of the latent non-stationary processes. We compare the performance of our proposed technique with several existing methods using simulated data, which reveals significant gains in terms of the bias-variance trade-off.

[1]  W. Singer,et al.  Neural Synchrony in Brain Disorders: Relevance for Cognitive Dysfunctions and Pathophysiology , 2006, Neuron.

[2]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[3]  J. R. Rosenberg,et al.  Time and Frequency Domain Analysis of Spike Train and Time Series Data , 1999 .

[4]  E. Brown.  Course 14 - Theory of Point Processes for Neural Systems , 2005 .

[5]  D. Slepian.  Prolate spheroidal wave functions, fourier analysis, and uncertainty — V: the discrete case , 1978, The Bell System Technical Journal.

[6]  A. Schnitzler,et al.  Normal and pathological oscillatory communication in the brain , 2005, Nature Reviews Neuroscience.

[7]  Piet de Jong,et al.  Covariances for smoothed estimates in state space models , 1988 .

[8]  Sergey L. Gratiy,et al.  Fully integrated silicon probes for high-density recording of neural activity , 2017, Nature.

[9]  C. Striebel,et al.  On the maximum likelihood estimates for linear dynamic systems , 1965 .

[10]  Jiangang Du,et al.  Multiplexed, High Density Electrophysiology with Nanofabricated Neural Probes , 2011, PloS one.

[11]  N. Huang,et al.  The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[12]  M. Priestley.  Evolutionary Spectra and Non‐Stationary Processes , 1965 .

[13]  Emery N. Brown,et al.  Robust Estimation of Sparse Narrowband Spectra from Neuronal Spiking Data , 2017, IEEE Transactions on Biomedical Engineering.

[14]  Donald B. Percival,et al.  Spectral Analysis for Physical Applications , 1993 .

[15]  Laura D. Lewis,et al.  Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness , 2012, Proceedings of the National Academy of Sciences.

[16]  Sonja Grün,et al.  Unitary Events in Multiple Single-Neuron Spiking Activity: II. Nonstationary Data , 2002, Neural Computation.

[17]  L. M. Ward,et al.  Synchronous neural oscillations and cognitive processes , 2003, Trends in Cognitive Sciences.

[18]  Behtash Babadi,et al.  Dynamic Bayesian Multitaper Spectral Analysis , 2017, IEEE Transactions on Signal Processing.

[19]  R. Kass,et al.  Multiple neural spike train data analysis: state-of-the-art and future challenges , 2004, Nature Neuroscience.

[20]  L. Paninski.  Maximum likelihood estimation of cascade point-process neural encoding models , 2004, Network.

[21]  Emery N. Brown,et al.  A Review of Multitaper Spectral Analysis , 2014, IEEE Transactions on Biomedical Engineering.

[22]  A. Walden.  A unified view of multitaper multivariate spectral estimation , 2000 .

[23]  Emery N. Brown,et al.  Estimating a State-space Model from Point Process Observations Emery N. Brown , 2022 .

[24]  Patrick Flandrin,et al.  Wigner-Ville spectral analysis of nonstationary processes , 1985, IEEE Trans. Acoust. Speech Signal Process..

[25]  D. Thomson,et al.  Spectrum estimation and harmonic analysis , 1982, Proceedings of the IEEE.