Infinite Systems of Interacting Chains with Memory of Variable Length—A Stochastic Model for Biological Neural Nets

We consider a new class of non Markovian processes with a countable number of interacting components. At each time unit, each component can take two values, indicating if it has a spike or not at this precise moment. The system evolves as follows. For each component, the probability of having a spike at the next time unit depends on the entire time evolution of the system after the last spike time of the component. This class of systems extends in a non trivial way both the interacting particle systems, which are Markovian (Spitzer in Adv. Math. 5:246–290, 1970) and the stochastic chains with memory of variable length which have finite state space (Rissanen in IEEE Trans. Inf. Theory 29(5):656–664, 1983). These features make it suitable to describe the time evolution of biological neural systems. We construct a stationary version of the process by using a probabilistic tool which is a Kalikow-type decomposition either in random environment or in space-time. This construction implies uniqueness of the stationary process. Finally we consider the case where the interactions between components are given by a critical directed Erdös-Rényi-type random graph with a large but finite number of components. In this framework we obtain an explicit upper-bound for the correlation between successive inter-spike intervals which is compatible with previous empirical findings.

[1]  T. E. Harris,et al.  On chains of infinite orde , 1955 .

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

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

[4]  P. Brémaud,et al.  STABILITY OF NONLINEAR HAWKES PROCESSES , 1996 .

[5]  F. Spitzer.  Interaction of Markov processes , 1970 .

[6]  H. Berbee,et al.  Chains with infinite connections: Uniqueness and Markov representation , 1987 .

[7]  Anders Öberg,et al.  Square summability of variations of g-functions and uniqueness of g-measures , 2003 .

[8]  JORMA RISSANEN,et al.  A universal data compression system , 1983, IEEE Trans. Inf. Theory.

[9]  Robert Fortet,et al.  Sur des chaînes à liaisons complètes , 1937 .

[10]  P. Ferrari,et al.  Perfect simulation for interacting point processes, loss networks and Ising models , 1999, math/9911162.

[11]  J. Goldberg,et al.  RESPONSE OF NEURONS OF THE SUPERIOR OLIVARY COMPLEX OF THE CAT TO ACOUSTIC STIMULI OF LONG DURATION. , 1964, Journal of neurophysiology.

[12]  Decay of Correlations for Non-Hölderian Dynamics. A Coupling Approach , 1998, math/9806132.

[13]  Nancy L. Garcia,et al.  Loss Network Representation of Peierls Contours , 1998, math/9806131.

[14]  Wulfram Gerstner,et al.  Spiking Neuron Models , 2002 .

[15]  Béla Bollobás,et al.  Random Graphs: Notation , 2001 .

[16]  Chains with Complete Connections and One-Dimensional Gibbs Measures , 2003, math/0305025.

[17]  B Cessac,et al.  A discrete time neural network model with spiking neurons: II: Dynamics with noise , 2010, Journal of mathematical biology.

[18]  D. Brillinger.  Maximum likelihood analysis of spike trains of interacting nerve cells , 2004, Biological Cybernetics.

[19]  A. Hawkes.  Point Spectra of Some Mutually Exciting Point Processes , 1971 .

[20]  Sandro Gallo,et al.  Chains with unbounded variable length memory: perfect simulation and a visible regeneration scheme , 2009, Advances in Applied Probability.

[21]  N. Garcia,et al.  Kalikow-type decomposition for multicolor infinite range particle systems , 2010, 1008.2740.

[22]  F. A. Galves.  Coupling, renewal and perfect simulation of chains of infinite order , 2001 .

[23]  J. Rasmussen,et al.  Perfect simulation of Hawkes processes , 2005, Advances in Applied Probability.

[24]  Pablo A. Ferrari,et al.  Processes with long memory: Regenerative construction and perfect simulation , 2002 .

[25]  Béla Bollobás,et al.  Random Graphs , 1985 .

[26]  Alejandro Maass,et al.  Cesàro mean distribution of group automata starting from measures with summable decay , 2000, Ergodic Theory and Dynamical Systems.