Unbiased estimation of precise temporal correlations between spike trains

A key issue in systems neuroscience is the contribution of precise temporal inter-neuronal interactions to information processing in the brain, and the main analytical tool used for studying pair-wise interactions is the cross-correlation histogram (CCH). Although simple to generate, a CCH is influenced by multiple factors in addition to precise temporal correlations between two spike trains, thus complicating its interpretation. A Monte-Carlo-based technique, the jittering method, has been suggested to isolate the contribution of precise temporal interactions to neural information processing. Here, we show that jittering spike trains is equivalent to convolving the CCH derived from the original trains with a finite window and using a Poisson distribution to estimate probabilities. Both procedures over-fit the original spike trains and therefore the resulting statistical tests are biased and have low power. We devise an alternative method, based on convolving the CCH with a partially hollowed window, and illustrate its utility using artificial and real spike trains. The modified convolution method is unbiased, has high power, and is computationally fast. We recommend caution in the use of the jittering method and in the interpretation of results based on it, and suggest using the modified convolution method for detecting precise temporal correlations between spike trains.

[1]  Stuart N. Baker,et al.  An Accurate Measure of the Instantaneous Discharge Probability, with Application to Unitary Joint-Event Analysis , 2000, Neural Computation.

[2]  A. Aertsen,et al.  On the significance of correlations among neuronal spike trains , 2004, Biological Cybernetics.

[3]  G. Buzsáki,et al.  Behavior-dependent short-term assembly dynamics in the medial prefrontal cortex , 2008, Nature Neuroscience.

[4]  Sonja Grün,et al.  Detecting unitary events without discretization of time , 1999, Journal of Neuroscience Methods.

[5]  E. S. Pearson,et al.  On questions raised by the combination of tests based on discontinuous distributions. , 1950, Biometrika.

[6]  George L Gerstein,et al.  Searching for significance in spatio-temporal firing patterns. , 2004, Acta neurobiologiae experimentalis.

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

[8]  Sonja Grün,et al.  Unitary Events in Multiple Single-Neuron Spiking Activity: I. Detection and Significance , 2002, Neural Computation.

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

[10]  D. Perkel,et al.  Nerve-impulse patterns: A quantitative display technique for three neurons , 1975, Brain Research.

[11]  Jerzy Neyman,et al.  Joint Statistical Papers , 1967 .

[12]  M. Abeles,et al.  Detecting precise firing sequences in experimental data , 2001, Journal of Neuroscience Methods.

[13]  DiesmannMarkus,et al.  Unitary events in multiple single-neuron spiking activity , 2002 .

[14]  Jerzy Neyman,et al.  The testing of statistical hypotheses in relation to probabilities a priori , 1933, Mathematical Proceedings of the Cambridge Philosophical Society.

[15]  Amir Globerson,et al.  Correlations between Groups of Premotor Neurons Carry Information about Prehension , 2008, The Journal of Neuroscience.

[16]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. I. The single spike train. , 1967, Biophysical journal.

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

[18]  Asohan Amarasingham,et al.  At what time scale does the nervous system operate? , 2003, Neurocomputing.

[19]  Sonja Grün,et al.  Detecting synfire chain activity using massively parallel spike train recording. , 2008, Journal of neurophysiology.

[20]  R N Lemon,et al.  Synchronization in monkey motor cortex during a precision grip task. I. Task-dependent modulation in single-unit synchrony. , 2001, Journal of neurophysiology.

[21]  W. Singer,et al.  Synchronization of neuronal responses in primary visual cortex of monkeys viewing natural images. , 2008, Journal of neurophysiology.

[22]  Carlos D. Brody,et al.  Correlations Without Synchrony , 1999, Neural Computation.

[23]  M. Abeles.  Local Cortical Circuits: An Electrophysiological Study , 1982 .

[24]  M. Abeles.  Quantification, smoothing, and confidence limits for single-units' histograms , 1982, Journal of Neuroscience Methods.

[25]  Hagai Bergman,et al.  Trial to trial variability in either stimulus or action causes apparent correlation and synchrony in neuronal activity , 2001, Journal of Neuroscience Methods.

[26]  George L. Gerstein,et al.  Improvements to the Sensitivity of Gravitational Clustering for Multiple Neuron Recordings , 2000, Neural Computation.

[27]  A. Aertsen,et al.  Spike synchronization and rate modulation differentially involved in motor cortical function. , 1997, Science.

[28]  E. Vaadia,et al.  Spatiotemporal structure of cortical activity: properties and behavioral relevance. , 1998, Journal of neurophysiology.

[29]  Akira Date,et al.  On the temporal resolution of neural activity , 1998 .

[30]  A. Tamhane,et al.  Multiple Comparison Procedures , 1989 .

[31]  Z. Nadasdy,et al.  Neurons of the cerebral cortex exhibit precise interspike timing in correspondence to behavior. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[32]  A. Tamhane,et al.  Multiple Comparison Procedures , 2009 .

[33]  M. Abeles.  The Quantification and Graphic Display of Correlations Among Three Spike Trains , 1983, IEEE Transactions on Biomedical Engineering.

[34]  Sonja Grün,et al.  Bounds of the Ability to Destroy Precise Coincidences by Spike Dithering , 2007, BVAI.

[35]  E. Stark,et al.  Encoding of reach and grasp by single neurons in premotor cortex is independent of recording site. , 2007, Journal of neurophysiology.