  The inverse Ising model is used in computational neuroscience to infer
probability distributions of the synchronous activity of large neuronal
populations. This method allows for finding the Boltzmann distribution with
single neuron biases and pairwise interactions that maximizes the entropy and
reproduces the empirical statistics of the recorded neuronal activity. Here we
apply this strategy to large populations of retinal output neurons (ganglion
cells) of different types, stimulated by multiple visual stimuli with their own
statistics. The activity of retinal output neurons is driven by both the inputs
from upstream neurons and the recurrent connections. We first apply the
standard inverse Ising model approach, and show that it accounts well for the
system's collective behavior when the input visual stimulus has short-ranged
spatial correlations, but fails for long-ranged ones. This happens because
stimuli with long-ranged spatial correlations synchronize the activity of
neurons over long distances. This effect cannot be accounted for by pairwise
interactions, and so by the pairwise Ising model. To solve this issue, we apply
a previously proposed framework that includes a temporal dependence in the
single neurons biases to model how neurons are driven in time by the stimulus.
Thanks to this addition, the stimulus effects are taken into account by the
biases, and the pairwise interactions allow for characterizing the network
effect in the population activity and for reproducing the structure of the
recurrent functional connections in the retinal architecture. We found that the
retinal architecture splits into weakly interacting subpopulations composed of
strongly interacting neurons. Overall, this temporal framework fixes the
problems of the standard, static, inverse Ising model and accounts for the
system's collective behavior, for stimuli with either short or long-range
correlations.
