A Unified Framework of Third Order Time and Frequency Domain Analysis for Neural Spike Trains

Third order time and frequency analysis has exhibited great potential for correlation analysis of multi-sensor datasets, but is usually presented as separate time domain and frequency domain approaches. A combined framework of both frequency domain and time domain has rarely been used. This paper proposes a non-parametric third order time and frequency domain framework which used two dimensional Fourier transforms to bridge the gap between time domain and frequency domain. A unified framework offers flexibility and efficiency to apply to data. In this paper we study neural spike train data treated as stochastic point processes. In time domain direct analysis, third order cumulant densities of spike trains are applied, which need all first-, secondand third order product densities to be calculated before constructing the third order cumulant density, which brings additional challenges. The novelty in this study is that a new framework is proposed which can offer an alternative approach without calculating lower order quantities and can reveal nonlinear relationship between neural recordings. The results show that the present framework provides a novel non-parametric method to estimate both time and frequency domain measurements which is applicable to neural spike trains.

[1]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains. , 1967, Biophysical journal.

[2]  Olaf Sporns,et al.  Complex network measures of brain connectivity: Uses and interpretations , 2010, NeuroImage.

[3]  K. S. Lii,et al.  Cross-Bispectrum Computation and Variance Estimation , 1981, TOMS.

[4]  M.R. Raghuveer,et al.  Bispectrum estimation: A digital signal processing framework , 1987, Proceedings of the IEEE.

[5]  D. Brillinger The spectral analysis of stationary interval functions , 1972 .

[6]  Gary R. Wilson,et al.  Time delay estimation using the cross bispectrum , 1990, IEEE Trans. Signal Process..

[7]  M K Habib,et al.  Dynamics of neuronal firing correlation: modulation of "effective connectivity". , 1989, Journal of neurophysiology.

[8]  Stefan Rotter,et al.  Higher-Order Statistics of Input Ensembles and the Response of Simple Model Neurons , 2003, Neural Computation.

[9]  Ifije E. Ohiorhenuan,et al.  Sparse coding and high-order correlations in fine-scale cortical networks , 2010, Nature.

[10]  David M. Halliday Spike-Train Analysis for Neural Systems , 2005 .

[11]  Shan Yu,et al.  Higher-Order Interactions Characterized in Cortical Activity , 2011, The Journal of Neuroscience.

[12]  Stefan Rotter,et al.  Interplay between Graph Topology and Correlations of Third Order in Spiking Neuronal Networks , 2016, PLoS Comput. Biol..

[13]  J. R. Rosenberg,et al.  The Fourier approach to the identification of functional coupling between neuronal spike trains. , 1989, Progress in biophysics and molecular biology.

[14]  C. Granger Investigating causal relations by econometric models and cross-spectral methods , 1969 .

[15]  Jerry M. Mendel,et al.  Tutorial on higher-order statistics (spectra) in signal processing and system theory: theoretical results and some applications , 1991, Proc. IEEE.

[16]  A M Amjad,et al.  A framework for the analysis of mixed time series/point process data--theory and application to the study of physiological tremor, single motor unit discharges and electromyograms. , 1995, Progress in biophysics and molecular biology.

[17]  R. Mason,et al.  Inhibition of endocannabinoid metabolism attenuates enhanced hippocampal neuronal activity induced by kainic acid , 2008, Synapse.

[18]  B. Conway,et al.  Detection of weak synaptic interactions between single Ia afferent and motor‐unit spike trains in the decerebrate cat. , 1993, The Journal of physiology.