Fast and reliable methods for extracting functional connectivity in large populations

TY - ABST

T1 - Fast and reliable methods for extracting functional connectivity in large populations

AU - Roudi, Yasser

AU - Tyrcha, Joanna

AU - Hertz, John

N1 - Udgivelsesdato: 13 July 2009 Volumne: 10

PY - 2009

Y1 - 2009

N2 - (Dansk abstrakt findes ikke)The simplest model for describing multi-neuron spike statistics is the pairwise Ising model. To start, one divides the spike trains into small time bins and to each neuron i and each time bin t assigns a binary variables si(t)= -1 if neuron i has not emitted any spikes in that time bin, and 1 if it has emitted one spike or more. One then can construct an Ising model, P(s )=Z-1exp{h.s+sJs} for the spike patterns with the same means and pair correlations as the data, using Boltzmann learning, which is in principle exact.  The elements Jij , of the matrix J can be considered to be functional couplings. However, Boltzmann learning is prohibitively time-consuming for large networks. Here, we compare the results from five fast approximate methods for finding the couplings with those from Boltzmann learning.    We used data from a simulated network of spiking neurons operating in a balanced state of asynchronous firing with a mean rate of ~10 Hz for excitatory neurons. Employing a bin size of 10 ms, we performed Boltzmann learning to fit Ising models for populations of size N up to 200 excitatory neurons chosen randomly from the 800 in the simulated network.  We studied the following methods:  A) a naive mean-field approximation, for which J is equal to minus the inverse of the covariance matrix. B) an independent-pair approximation, C) a low rate, small-population approximation (the low-rate limit of (B), which is valid generally in the limit of small Nrt, where r is the average rate (spikes/time bin) and t is the bin width,  D) inversion of the TAP equations from spin-glass theory, and  E) a weak-correlation approximation proposed recently by Sessak and Monasson.  We quantified the quality of these approximations, as functions of N, by computing the RMS error and R2, treating the Boltzmann couplings as the true ones. We found that while all the approximations are good for small N, the TAP, Sessak-Monasson, and, in particular, their average outperform the others by a relatively large margin for N.  Thus, these methods offer a useful tool for fast analysis of multineuron spike data.

AB - (Dansk abstrakt findes ikke)The simplest model for describing multi-neuron spike statistics is the pairwise Ising model. To start, one divides the spike trains into small time bins and to each neuron i and each time bin t assigns a binary variables si(t)= -1 if neuron i has not emitted any spikes in that time bin, and 1 if it has emitted one spike or more. One then can construct an Ising model, P(s )=Z-1exp{h.s+sJs} for the spike patterns with the same means and pair correlations as the data, using Boltzmann learning, which is in principle exact.  The elements Jij , of the matrix J can be considered to be functional couplings. However, Boltzmann learning is prohibitively time-consuming for large networks. Here, we compare the results from five fast approximate methods for finding the couplings with those from Boltzmann learning.    We used data from a simulated network of spiking neurons operating in a balanced state of asynchronous firing with a mean rate of ~10 Hz for excitatory neurons. Employing a bin size of 10 ms, we performed Boltzmann learning to fit Ising models for populations of size N up to 200 excitatory neurons chosen randomly from the 800 in the simulated network.  We studied the following methods:  A) a naive mean-field approximation, for which J is equal to minus the inverse of the covariance matrix. B) an independent-pair approximation, C) a low rate, small-population approximation (the low-rate limit of (B), which is valid generally in the limit of small Nrt, where r is the average rate (spikes/time bin) and t is the bin width,  D) inversion of the TAP equations from spin-glass theory, and  E) a weak-correlation approximation proposed recently by Sessak and Monasson.  We quantified the quality of these approximations, as functions of N, by computing the RMS error and R2, treating the Boltzmann couplings as the true ones. We found that while all the approximations are good for small N, the TAP, Sessak-Monasson, and, in particular, their average outperform the others by a relatively large margin for N.  Thus, these methods offer a useful tool for fast analysis of multineuron spike data.

U2 - 10.1186/1471-2202-10-S1-O9

DO - 10.1186/1471-2202-10-S1-O9

M3 - Conference abstract in journal

SN - 1471-2202

VL - 10, suppl 1

SP - 09

JO - BMC Neuroscience

JF - BMC Neuroscience

T2 - 18th Annual Computational Neuroscience Meeting CNS*2009

Y2 - 18 July 2009 through 23 July 2009

ER -