Skip to content


  • Poster presentation
  • Open Access

Using spike train distances to identify the most discriminative neuronal subpopulation

BMC Neuroscience201314 (Suppl 1) :P35

  • Published:


  • Firing Rate
  • Individual Neuron
  • Discrimination Performance
  • Spike Timing
  • Pairwise Similarity

During the last decade spike train distances [14] have become an essential means to characterize neural coding in a wide range of neurophysiological contexts. In a typical setup different stimuli are presented repeatedly and spike train distances are used to carry out a pairwise similarity analysis in order to evaluate whether responses to the same stimulus exhibit smaller distances than responses to different stimuli. With the increasing availability of multi-neuron recordings these kinds of analyses can now be performed not just for individual neurons but rather for simultaneously recorded neuronal populations.

Recently, the bivariate Victor-Purpura and the van Rossum distances [1, 2] have been extended to quantify dissimilarities between multi-unit responses. These new population measures [5, 6] have been designed such that they can estimate the discrimination performance of either the population as a whole ("summed population") or of individual neurons ("labelled line") or of interpolations between these two extremes. However, the question that these approaches fail to answer is the following: In cases the two extremes fail, which subpopulation within the larger population discriminates the presented stimuli best?

In this study we thus follow a complimentary approach and present an algorithm which addresses exactly this question. The brute-force approach of calculating the pairwise distance matrices and the stimulus discrimination performance for every possible neuronal subpopulation is not feasible even for moderate numbers of neurons. Instead, our algorithm relies on an iterative scheme which considerably restricts the number of subpopulations for which the pairwise distance matrices and the stimulus discrimination performance actually have to be calculated.

The spike train distance that we use to evaluate whether the responses elicited by different stimuli can be distinguished is the SPIKE-distance [4, 7]. Whereas distances like the Victor-Purpura or the van Rossum spike train distance rely on a time-scale parameter, the SPIKE-distance is parameter-free and time-scale independent. This allows for easy comparability of results obtained for vastly different firing rates (which depend on the size of the pooled subpopulation). In contrast to the ISI-distance [3] it is sensitive to spike timing which is an important property that will be needed in some of the setups used.



TK and NB acknowledge funding support from the European Commission through the Marie Curie Initial Training Network 'Neural Engineering Transformative Technologies (NETT)', project 289146. TK would like to thank the Italian Ministry of Foreign Affairs for support regarding the activity of the Joint Italian-Israeli Laboratory on Neuroscience.

Authors’ Affiliations

Institute for Complex Systems, CNR, Sesto Fiorentino, Italy


  1. 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
  2. van Rossum MCW: A novel spike distance. Neural Computation. 2001, 13: 751-763. 10.1162/089976601300014321.View ArticlePubMedGoogle Scholar
  3. 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
  4. Kreuz T, Chicharro D, Houghton C, Andrzejak RG, Mormann F: Monitoring spike train synchrony. J Neurophysiol. 2013, Doi: 10.1152/jn.00873.2012Google Scholar
  5. 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.View ArticlePubMedGoogle Scholar
  6. Houghton C, Sen K: A new multineuron spike train metric. Neural Computation. 2008, 20: 1495-1511. 10.1162/neco.2007.10-06-350.PubMed CentralView ArticlePubMedGoogle Scholar
  7. The Matlab source code for calculating and visualizing the SPIKE-distance as well as information about the implementation can be found under. []


© Kreuz and Bozanic; licensee BioMed Central Ltd. 2013

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.