Multitaper Analysis of Evolutionary Spectra 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 non-stationary processes based on spiking observations is a challenging problem. In this paper, we address this issue by developing a multitaper spectral estimation methodology that can be directly applied to multivariate spiking observations in order to extract the evolutionary spectral density of the latent non-stationary processes that drive spiking activity, based on point process theory. We establish theoretical bounds on the bias-variance trade-off of the proposed estimator. Finally, we compare the performance of our proposed technique with existing methods using simulation studies and application to real data, which reveal significant gains in terms of the bias-variance trade-off.

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

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

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

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

[5]  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.

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

[7]  D. B. Preston.  Spectral Analysis and Time Series , 1983 .

[8]  Demba Ba,et al.  State-space multitaper time-frequency analysis , 2017, Proceedings of the National Academy of Sciences.

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

[10]  Behtash Babadi,et al.  Multitaper Spectral Analysis of Neuronal Spiking Activity Driven by Latent Stationary Processes , 2019, Signal Process..

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

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

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

[14]  Paul R. White,et al.  THE ANALYSIS OF NON-STATIONARY SIGNALS USING TIME-FREQUENCY METHODS , 1996 .

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

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

[17]  E. Brown,et al.  A point-process model of human heartbeat intervals: new definitions of heart rate and heart rate variability. , 2005, American journal of physiology. Heart and circulatory physiology.

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

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

[20]  Thomas P. Bronez,et al.  On the performance advantage of multitaper spectral analysis , 1992, IEEE Trans. Signal Process..

[21]  Behtash Babadi,et al.  Multitaper Analysis of Evolutionary Spectral Density Matrix From Multivariate Spiking Observations , 2019, 2019 IEEE Data Science Workshop (DSW).

[22]  Leon Cohen,et al.  Time Frequency Analysis: Theory and Applications , 1994 .

[23]  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.

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

[25]  Gerald Matz,et al.  Generalized evolutionary spectral analysis and the Weyl spectrum of nonstationary random processes , 1997, IEEE Trans. Signal Process..

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

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

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

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

[30]  Murray Rosenblatt,et al.  Prolate spheroidal spectral estimates , 2008 .

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

[32]  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.

[33]  Gerald Matz,et al.  Nonstationary spectral analysis based on time-frequency operator symbols and underspread approximations , 2006, IEEE Transactions on Information Theory.

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

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

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

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

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

[39]  George Casella,et al.  Explaining the Saddlepoint Approximation , 1999 .