- Poster presentation
- Open Access
Network inference from non-stationary spike trains
BMC Neuroscience volume 12, Article number: P150 (2011)
Current approaches to the problem of inferring network connectivity from spike data [1, 2] make a stationarity assumption, which is often not valid. Here we describe a method for inferring both the connectivity of a network in the presence of nonstationarity state and the time-dependent external drive that causes it.
Consider an experiment in which the network is subjected repeatedly to a potentially unknown external input (such as would arise from sensory stimulation). We assume the spikes to be binned in time and represented by a binary array: S i (t,r) = ±1, according to whether neuron i fires or not in time bin t of repetition r of the measurement. We fit these data to the simplest kind of binary stochastic model: At time step t of repetition r, each formal neuron receives a net input, H i (t,r) = h i (t) + ∑ j J ij S j (t,r), and it takes the value +1 at the next step with a probability given by a logistic sigmoidal function 1/[1+exp(-2H i (t,r))] of H i (t,r). Maximizing the likelihood of the data leads to learning rules
for the model parameters -- the couplings J ij and external inputs h i (t). For weak coupling and/or densely connected networks, we have developed faster alternative algorithms . These are based on expanding the learning rules around mean-field and TAP  equations for m i (t) = ‹S i (t,r)› r . (TAP equations are a generalization of the usual mean-field equations for highly connected random networks.)
We have applied this method, as well as conventional ones assuming stationarity, to data sets from (1) the stochastic model itself, (2) a realistic computational model of a small cortical network, and (3) salamander retina under visual stimulation. In all three cases, we find that if we perform the reconstruction assuming stationarity, the algorithms effectively invent fictitious couplings to explain stimulus-induced correlations: The couplings in the network are systematically overestimated. The nonstationary treatment outlined above enables us, for sufficient data, to find both the correct (weaker) couplings and to extract the time-dependence of the external input. To illustrate this, figure 1 shows the J ij s found using the nonstationary algorithm plotted against those found using the stationary one, based on spike trains of 40 salamander retinal neurons stimulated by 120 repetitions of a 26.5-second clip from a film.. The mean J ij is reduced, from 0.0471 to -0.0028, and the large positive J ij s found assuming stationarity are reduced by a facto of 2-3 when nonstationarity is taken into account.
Schneidman E, Berry M, Segev R, Bialek W: Weak pairwise correlations imply strongly correlated networks states in a neural population. Nature. 2006, 440: 1007-1012. 10.1038/nature04701.
Roudi Y, Tyrcha J, Hertz J: Ising model for neural data: Model quality and approximate methods for extracting functional connectivity. Phys Rev E. 2009, 79: 051915-10.1103/PhysRevE.79.051915.
Roudi Y, Hertz J: Mean field theory for nonequilibrium network reconstruction. Phys Rev Lett. 2011, 106: 048702-10.1103/PhysRevLett.106.048702.
Thouless DJ, Anderson PW, Palmer RG: Solution of ‘solvable model of a spin glass’. Phil Mag. 1977, 35: 593-601. 10.1080/14786437708235992.
We thank Michael Berry for providing the salamander retinal data.