Multitaper Analysis of Semi-Stationary Spectra From Multivariate Neuronal Spiking Observations

Extracting the spectral representations of neural processes that underlie spiking activity is key to understanding how 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 challenging due to the underlying nonlinearities that limit the spectrotemporal resolution of existing methods. 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 semi-stationary spectral density of the latent non-stationary processes that govern spiking activity. We establish theoretical bounds on the bias-variance trade-off of our proposed estimator. Finally, application of our proposed technique to simulated and real data reveals significant performance gains over existing methods.

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

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

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

[4]  Emery N. Brown,et al.  Dynamic Analysis of Neural Encoding by Point Process Adaptive Filtering , 2004, Neural Computation.

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

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

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

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

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

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

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

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

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

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

[15]  Daniel Levenstein,et al.  Network Homeostasis and State Dynamics of Neocortical Sleep , 2016, Neuron.

[16]  Liam Paninski,et al.  Statistical models for neural encoding, decoding, and optimal stimulus design. , 2007, Progress in brain research.

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

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

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

[20]  Tim Gollisch,et al.  Modeling Single-Neuron Dynamics and Computations: A Balance of Detail and Abstraction , 2006, Science.

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

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

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

[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]  Paul R. White,et al.  THE ANALYSIS OF NON-STATIONARY SIGNALS USING TIME-FREQUENCY METHODS , 1996 .

[27]  Maxim Volgushev,et al.  Properties of Slow Oscillation during Slow-Wave Sleep and Anesthesia in Cats , 2011, The Journal of Neuroscience.

[28]  Behtash Babadi,et al.  Multitaper Analysis of Evolutionary Spectra from Multivariate Spiking Observations , 2019, ArXiv.

[29]  Emery N. Brown,et al.  Computational Neuroscience: A Comprehensive Approach , 2022 .

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

[31]  Shihab A. Shamma,et al.  Recursive Sparse Point Process Regression With Application to Spectrotemporal Receptive Field Plasticity Analysis , 2015, IEEE Transactions on Signal Processing.

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

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

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

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

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

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

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

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

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

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

[42]  S. Hughes,et al.  The slow (<1 Hz) rhythm of non-REM sleep: a dialogue between three cardinal oscillators , 2010, Nature Neuroscience.

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

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

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

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

[47]  G. Tononi,et al.  Sleep homeostasis and cortical synchronization: II. A local field potential study of sleep slow waves in the rat. , 2007, Sleep.

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

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