Statistical properties of superimposed stationary spike trains

The Poisson process is an often employed model for the activity of neuronal populations. It is known, though, that superpositions of realistic, non- Poisson spike trains are not in general Poisson processes, not even for large numbers of superimposed processes. Here we construct superimposed spike trains from intracellular in vivo recordings from rat neocortex neurons and compare their statistics to specific point process models. The constructed superimposed spike trains reveal strong deviations from the Poisson model. We find that superpositions of model spike trains that take the effective refractoriness of the neurons into account yield a much better description. A minimal model of this kind is the Poisson process with dead-time (PPD). For this process, and for superpositions thereof, we obtain analytical expressions for some second-order statistical quantities—like the count variability, inter-spike interval (ISI) variability and ISI correlations—and demonstrate the match with the in vivo data. We conclude that effective refractoriness is the key property that shapes the statistical properties of the superposition spike trains. We present new, efficient algorithms to generate superpositions of PPDs and of gamma processes that can be used to provide more realistic background input in simulations of networks of spiking neurons. Using these generators, we show in simulations that neurons which receive superimposed spike trains as input are highly sensitive for the statistical effects induced by neuronal refractoriness.

[1]  Moritz Helias,et al.  Equilibrium and response properties of the integrate-and-fi re neuron in discrete time , 2022 .

[2]  Johannes Schemmel,et al.  Spike-Frequency Adapting Neural Ensembles: Beyond Mean Adaptation and Renewal Theories , 2007, Neural Computation.

[3]  Robert E. Kass,et al.  A Spike-Train Probability Model , 2001, Neural Computation.

[4]  Benjamin Lindner,et al.  Superposition of many independent spike trains is generally not a Poisson process. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  Farzad Farkhooi,et al.  Adaptation reduces variability of the neuronal population code. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Bernard C. Picinbono Output Dead-Time in Point Processes , 2009, Commun. Stat. Simul. Comput..

[7]  A. Reyes,et al.  Relation between single neuron and population spiking statistics and effects on network activity. , 2006, Physical review letters.

[8]  Wulfram Gerstner,et al.  Spiking Neuron Models: An Introduction , 2002 .

[9]  D H Johnson,et al.  The transmission of signals by auditory-nerve fiber discharge patterns. , 1983, The Journal of the Acoustical Society of America.

[10]  A. Aertsen,et al.  Beyond the Cortical Column: Abundance and Physiology of Horizontal Connections Imply a Strong Role for Inputs from the Surround , 2011, Front. Neurosci..

[11]  Carl van Vreeswijk,et al.  Temporal Correlations in Stochastic Networks of Spiking Neurons , 2002, Neural Computation.

[12]  Wulfram Gerstner,et al.  SPIKING NEURON MODELS Single Neurons , Populations , Plasticity , 2002 .

[13]  A. R. Kamat DISTRIBUTION THEORY OF TWO ESTIMATES FOR STANDARD DEVIATION BASED ON SECOND VARIATE DIFFERENCES , 1954 .

[14]  P. Dayan,et al.  Supporting Online Material Materials and Methods Som Text Figs. S1 to S9 References the Asynchronous State in Cortical Circuits , 2022 .

[15]  Jaime de la Rocha,et al.  Supplementary Information for the article ‘ Correlation between neural spike trains increases with firing rate ’ , 2007 .

[16]  M. Nawrot,et al.  Serial correlation in neural spike trains: experimental evidence, stochastic modeling, and single neuron variability. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Benjamin Lindner,et al.  Interspike interval statistics of neurons driven by colored noise. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[18]  S. Ostojic Interspike interval distributions of spiking neurons driven by fluctuating inputs. , 2011, Journal of neurophysiology.

[19]  Sonja Grün,et al.  Analysis of Parallel Spike Trains , 2010 .

[20]  A. Holden Models of the stochastic activity of neurones , 1976 .

[21]  Sonja Grün,et al.  Serial interval statistics of spontaneous activity in cortical neurons in vivo and in vitro , 2007, Neurocomputing.

[22]  S. Levin Lectu re Notes in Biomathematics , 1983 .

[23]  G L GERSTEIN,et al.  An approach to the quantitative analysis of electrophysiological data from single neurons. , 1960, Biophysical journal.

[24]  J. Assad,et al.  Beyond Poisson: Increased Spike-Time Regularity across Primate Parietal Cortex , 2009, Neuron.

[25]  C. Vreeswijk Stochastic Models of Spike Trains , 2010 .

[26]  Robert Rosenbaum,et al.  Mechanisms That Modulate the Transfer of Spiking Correlations , 2011, Neural Computation.

[27]  H. Barlow,et al.  MAINTAINED ACTIVITY IN THE CAT'S RETINA IN LIGHT AND DARKNESS , 1957, The Journal of general physiology.

[28]  David R. Cox,et al.  The statistical analysis of series of events , 1966 .

[29]  Wilson Truccolo,et al.  Stochastic Models for Multivariate Neural Point Processes: Collective Dynamics and Neural Decoding , 2010 .

[30]  Nicolas Brunel,et al.  Dynamics of Networks of Excitatory and Inhibitory Neurons in Response to Time-Dependent Inputs , 2011, Front. Comput. Neurosci..

[31]  Walter L. Smith,et al.  ON THE SUPERPOSITION OF RENEWAL PROCESSES , 1954 .

[32]  Stefan Rotter,et al.  Dependence of Neuronal Correlations on Filter Characteristics and Marginal Spike Train Statistics , 2008, Neural Computation.

[33]  Benjamin Lindner,et al.  How Noisy Adaptation of Neurons Shapes Interspike Interval Histograms and Correlations , 2010, PLoS Comput. Biol..

[34]  R. Douglas,et al.  A Quantitative Map of the Circuit of Cat Primary Visual Cortex , 2004, The Journal of Neuroscience.

[35]  Gert Cauwenberghs,et al.  Neuromorphic Silicon Neuron Circuits , 2011, Front. Neurosci.

[36]  Néstor Parga,et al.  Response of integrate-and-fire neurons to noisy inputs filtered by synapses with arbitrary timescales: firing rate and correlations , 2009, Neural Computation.

[37]  J. Donoghue,et al.  Collective dynamics in human and monkey sensorimotor cortex: predicting single neuron spikes , 2009, Nature Neuroscience.

[38]  Eero P. Simoncelli,et al.  Spatio-temporal correlations and visual signalling in a complete neuronal population , 2008, Nature.

[39]  Don H. Johnson,et al.  Point process models of single-neuron discharges , 1996, Journal of Computational Neuroscience.

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

[41]  Daniel P. Heyman,et al.  Stochastic models in operations research , 1982 .

[42]  Stefan Rotter,et al.  Measurement of variability dynamics in cortical spike trains , 2008, Journal of Neuroscience Methods.

[43]  Shigeru Shinomoto,et al.  Differences in Spiking Patterns Among Cortical Neurons , 2003, Neural Computation.

[44]  Henry C. Tuckwell,et al.  Introduction to theoretical neurobiology , 1988 .

[45]  H. Sompolinsky,et al.  Chaos in Neuronal Networks with Balanced Excitatory and Inhibitory Activity , 1996, Science.

[46]  I. Miller Probability, Random Variables, and Stochastic Processes , 1966 .

[47]  Olivier D. Faugeras,et al.  A Constructive Mean-Field Analysis of Multi-Population Neural Networks with Random Synaptic Weights and Stochastic Inputs , 2008, Front. Comput. Neurosci..

[48]  D. Cox,et al.  The statistical analysis of series of events , 1966 .

[49]  Nicolas Brunel,et al.  Dynamics of Sparsely Connected Networks of Excitatory and Inhibitory Spiking Neurons , 2000, Journal of Computational Neuroscience.

[50]  N. Brunel,et al.  Firing frequency of leaky intergrate-and-fire neurons with synaptic current dynamics. , 1998, Journal of theoretical biology.

[51]  Mounir Boukadoum,et al.  Mechanisms Gating the Flow of Information in the Cortex: What They Might Look Like and What Their Uses may be , 2010, Front. Comput. Neurosci..

[52]  Michael J. Berry,et al.  Refractoriness and Neural Precision , 1997, The Journal of Neuroscience.

[53]  Moritz Helias,et al.  Instantaneous Non-Linear Processing by Pulse-Coupled Threshold Units , 2010, PLoS Comput. Biol..

[54]  Moritz Helias,et al.  Nonequilibrium dynamics of stochastic point processes with refractoriness. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[55]  Moritz Helias,et al.  Finite Post Synaptic Potentials Cause a Fast Neuronal Response , 2010, Front. Neurosci..

[56]  Eric Li,et al.  A nonstationary Poisson point process describes the sequence of action potentials over long time scales in lateral-superior-olive auditory neurons , 2004, Biological Cybernetics.