Volume 10 Supplement 1
Measuring spike train synchrony between neuronal populations
© Kreuz et al; licensee BioMed Central Ltd. 2009
Published: 13 July 2009
With the increasing availability of multi-unit recordings the focus of attention starts to shift from bivariate methods towards methods that provide the possibility to study patterns of activity across many neurons. Measures of multi-neuron spike train synchrony are becoming indispensable tools for addressing issues such as network synchronization, spike timing reliability and neuronal coding. However, many multi-neuron synchrony measures are extensions of bivariate measures. Two of the most prominent bivariate approaches are the spike train metrics by Victor-Purpura and van Rossum [1, 2]. The former evaluates the cost needed to transform one spike train into the other using only certain elementary steps , while the latter measures the Euclidean distance between the two spike trains after convolution of the spikes with an exponential function . Both methods involve one parameter that sets the time scale. In contrast, a more recent bivariate approach, the ISI-distance , is time scale independent and self-adaptive. Another essential difference is that the ISI-distance relies on the relative length of interspike intervals (ISI) and not on the timing of spikes. Finally, this method also allows the visualization of the relative firing pattern in a time-resolved manner.
Recently, the Victor-Purpura and the van Rossum distances have been extended to quantify dissimilarities between multi-unit responses [4, 5]. To calculate the multi-unit Victor-Purpura metric, simultaneous spikes are labeled by the neuron that fired them, but this label can be changed at an additional cost which is determined by a second parameter. By varying this population parameter the metric is shifted from a "labeled line" (LL) code metric in which the distance is defined as the sum of the distances of the single neurons, to a "summed population" (SP) code metric in which the spike trains are superimposed before the distance is calculated . In the extended van Rossum metric, the spike trains of each population are located in a space of vector fields (with a different unit vector assigned to each neuron). In this case interpolation between the LL and the SP coding is achieved by varying the angle (a second parameter) between unit vectors .
Here we present an analogous extension for the ISI-distance  that also interpolates between the LL and the SP codes. This multi-neuron ISI-distance inherits all the basic properties of the bivariate ISI-distance described above; in particular, it is also time scale independent and thus, in contrast to the other two multi-neuron metrics, depends on one population parameter only. In this study we compare all three multi-neuron distances using both controlled simulations and real data. We stress the advantages of our extension with respect to visualization, computational cost and applicability to larger numbers of spike trains with higher numbers of spikes.
TK has been supported by the Marie Curie Individual Intra-European Fellowship "DEAN," project No 011434, DC by the I.U.E. Department of the Generalitat of Catalunya and the European Social Fund, and RGA by the Ramón y Cajal program.
- Victor J, Purpura K: Nature and precision of temporal coding in visual cortex: A metric-space analysis. J Neurophysiol. 1996, 76: 1310-1326.PubMedGoogle Scholar
- van Rossum MCW: A novel spike distance. Neural Computation. 2001, 13: 751-763. 10.1162/089976601300014321.PubMedView ArticleGoogle Scholar
- Kreuz T, Haas JS, Morelli A, Abarbanel HDI, Politi A: Measuring spike train synchronization. J Neurosci Methods. 2007, 165: 151-161. 10.1016/j.jneumeth.2007.05.031.PubMedView ArticleGoogle Scholar
- Aronov D, Reich DS, Mechler F, Victor JD: Neural coding of spatial phase in V1 of the macaque monkey. J Neurophysiol. 2003, 89: 3304-3327. 10.1152/jn.00826.2002.PubMedView ArticleGoogle Scholar
- Houghton C, Sen K: A new multineuron spike train metric. Neural Computation. 2008, 20: 1495-1511. 10.1162/neco.2007.10-06-350.PubMed CentralPubMedView ArticleGoogle Scholar
- The Matlab source code for calculating and visualizing all ISI-distances as well as information about their implementation can be found under. [http://inls.ucsd.edu/~kreuz/Source-Code/Spike-Sync.html]
This article is published under license to BioMed Central Ltd.