Measuring real-time synchronization in both spike trains and continuous time series
© Kreuz et al; licensee BioMed Central Ltd. 2011
Published: 18 July 2011
A wide variety of approaches to quantify the dissimilarity between two spike trains has been proposed. The Victor-Purpura metric  evaluates the cost needed to transform one spike train into the other using only certain elementary steps. Another metric proposed by van Rossum  measures the Euclidean distance between the two spike trains after convolution of the spikes with an exponential function. These methods, like many others, involve one parameter that sets the time-scale. In contrast, two more recent approaches, the ISI- and the SPIKE-distance, are parameter free and time-scale adaptive [3–5]. While the ISI-distance relies on the relative length of interspike intervals, the SPIKE-distance is sensitive to spike coincidences. Both measures are easy to visualize in a time-resolved manner and can be applied to more than two spike trains, either as averages over pairwise distances or as truly multivariate measures .
The SPIKE-distance is calculated from instantaneous values of spike train dissimilarity for which at each time moment not only the last previous spike but also the first following spike is taken into account. This non-causal dependence on ‘future’ spiking does not allow for a real-time calculation. Furthermore, the SPIKE-distance is designed to be applied to sequences of discrete events and does not foresee an application to continuous data (such as the electroencephalogram, EEG).
Here we address both issues and present two extensions of the SPIKE-distance. In the first extension the original method is modified such that instantaneous value of dissimilarity for two or more spike trains can be calculated in real-time. In the second extension the SPIKE-distance is transformed into a measure of dissimilarity which can be applied to two or more channels of continuous data. Both variants can also be combined such that the instantaneous synchronization among continuous time series is monitored in real-time.
DC has been supported by grant 2010FI-B2 00079 of the ’Generalitat de Catalunya’ and European Social Funds, and RGA by the Ramón y Cajal program.
- Victor JD, Purpura KP: 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 Comput. 2001, 13: 751-763. 10.1162/089976601300014321.View ArticlePubMedGoogle Scholar
- Kreuz T, Haas JS, Morelli A, Abarbanel HDI, Politi A: Measuring spike train synchrony. J Neurosci Methods. 2007, 165: 151-161. 10.1016/j.jneumeth.2007.05.031.View ArticlePubMedGoogle Scholar
- Kreuz T, Chicharro D, Andrzejak RG, Haas JS, Abarbanel HDI: Measuring multiple spike train synchrony. J Neurosci Methods. 2009, 183: 287-299. 10.1016/j.jneumeth.2009.06.039.View ArticlePubMedGoogle Scholar
- Kreuz T, Chicharro D, Greschner M, Andrzejak RG: Time-resolved and time-scale adaptive measures of spike train synchrony. J Neurosci Methods. 2011, 195: 92-106. 10.1016/j.jneumeth.2010.11.020.View ArticlePubMedGoogle Scholar
- . The Matlab source code for calculating and visualizing the ISI- and the SPIKE-distance as well as information about their implementation can be found under http://www.fi.isc.cnr.it/users/thomas.kreuz/sourcecode.html
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.