Measuring neuronal identification quality in ensemble recordings (isoitools) (Neymotin et al. 2011)

 Download zip file 
Help downloading and running models
"... Here we describe information theoretic measures of action potential waveform isolation applicable to any dataset, that have an intuitive, universal interpretation, and that are not dependent on the methods or choice of parameters for single unit isolation, and that have been validated using a dataset."
1 . Neymotin SA, Lytton WW, Olypher AO, Fenton AA (2011) Measuring the quality of neuronal identification in ensemble recordings Journal of Neuroscience 31(45):16398-16409 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type:
Brain Region(s)/Organism:
Cell Type(s):
Gap Junctions:
Simulation Environment: C or C++ program;
Model Concept(s): Action Potentials; Methods; Extracellular Fields;
Implementer(s): Neymotin, Sam [samn at];
isoitools - command-line tools for calculating cluster quality.
The software is intended for use with clusters formed from of feature
vectors of neural extracellular field recordings.

These programs were tested/developed on LINUX systems, but can be
compiled to run on Microsoft Windows or Mac OS. To compile, you will
need a standard C++ compiler and the make utility. Compile the C++
files from the command line with the make command. That will produce
executable files. The intermediate .o files can be removed with make
clean. make install will copy the binary files to the top-level
directory (isoitools).

layout of files:

  src - source code directory, along with documentation in readme
  files in subdirectories

    isoi - source code for program that calculates Isolation
    Information cluster quality measures

    isorat - source code for program that calculates Isolation
    Distance and LRatio cluster quality measures

    wave2f - source code for program that

  wavex.txt - short recording of extracellular waveforms using a
              tetrode in rat CA1. The format of the file allows using
              it with wave2f. The first line of the file is the number
              of microwires used for recording (4, since a
              tetrode). The next line is the number of samples on each
              channel (64). The following line specifies the sampling
              frequency (32 KHZ). The remaining lines are the waveform
              data for each spike (4 channels at 64 samples per
              channel, so 256 numbers per line).

To build the software, change to the src directory and use the make
command. Then run the make install command. This will copy the
binaries for the 3 programs above into the top-level directory. The
usage instructions and further information on each of the programs is
located in the subdirectories as readme.txt files.

wavex.txt is a sample data file that can be used with these programs.

The following instructions demonstrate how to use the programs:

 1. use wave2f to convert the waveforms into feature vectors & save
 results to file:
  wave2f wavex.txt wavefeat.txt
 2. use isoi to obtain the Isolation Information quality measures &
 save results to file:
  isoi wavefeat.txt wave_IsoI.txt
 3. use isorat to obtain the Isolation Distance and LRatio quality
 measures & save results to file:
  isorat wavefeat.txt wave_IsoDLRat.txt

Please consult the individual readme.txt files for more information on
the individual programs (located in the src subdirectories).

For more information contact Sam Neymotin ( samn at neurosim dot
 downstate dot edu ) or Andre Fenton ( afenton at nyu dot edu ).


 The methods used in these program are described in an article at The
  Journal of Neuroscience, Measuring the quality of neuronal
  identification in ensemble recordings by Neymotin SA, Lytton WW,
  Olypher AO, Fenton AA (2011).

<license info>
Copyright (C) 2003-2011 Sam Neymotin & BioSignal Group

This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or (at
your option) any later version.

This program is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of
General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program.  If not, see <>.
</license info>

Neymotin SA, Lytton WW, Olypher AO, Fenton AA (2011) Measuring the quality of neuronal identification in ensemble recordings Journal of Neuroscience 31(45):16398-16409[PubMed]

References and models cited by this paper

References and models that cite this paper

Bentley JL (1975) Multidimensional binary search trees used for associative searching Communications Of ACM 18:509-517

Bryant JL, Roy S, Heck DH (2009) A technique for stereotaxic recordings of neuronal activity in awake, head-restrained mice. J Neurosci Methods 178:75-9

Bulmer MG (1979) Principles of Statistics

Buzsaki G (2004) Large-scale recording of neuronal ensembles. Nat Neurosci 7:446-51 [PubMed]

Cormen TH, Leiserson CE, Rivest RL, Stein C (2001) Introduction to Algorithms

Csicsvari J, Henze DA, Jamieson B, Harris KD, Sirota A, Bartho P, Wise KD, Buzsaki G (2003) Massively parallel recording of unit and local field potentials with silicon-based electrodes. J Neurophysiol 90:1314-23 [Journal] [PubMed]

Dickey AS, Suminski A, Amit Y, Hatsopoulos NG (2009) Single-unit stability using chronically implanted multielectrode arrays. J Neurophysiol 102:1331-9 [PubMed]

Emondi AA, Rebrik SP, Kurgansky AV, Miller KD (2004) Tracking neurons recorded from tetrodes across time. J Neurosci Methods 135:95-105 [PubMed]

Fenton AA, Kao HY, Neymotin SA, Olypher A, Vayntrub Y, Lytton WW, Ludvig N (2008) Unmasking the CA1 ensemble place code by exposures to small and large environments: more place cells and multiple, irregularly arranged, and expanded place fields in the larger space. J Neurosci 28:11250-62 [PubMed]

Gray CM, Maldonado PE, Wilson M, McNaughton B (1995) Tetrodes markedly improve the reliability and yield of multiple single-unit isolation from multi-unit recordings in cat striate cortex. J Neurosci Methods 63:43-54 [PubMed]

Hampson RE, Simeral JD, Deadwyler SA (2001) What ensemble recordings reveal about functional hippocampal cell encoding. Prog Brain Res 130:345-57 [PubMed]

Harris KD, Csicsvari J, Hirase H, Dragoi G, Buzsaki G (2003) Organization of cell assemblies in the hippocampus. Nature 424:552-6 [PubMed]

Harris KD, Henze DA, Csicsvari J, Hirase H, Buzsaki G (2000) Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements. J Neurophysiol 84:401-14 [Journal] [PubMed]

Henze DA, Borhegyi Z, Csicsvari J, Mamiya A, Harris KD, ki G (2000) Intracellular features predicted by extracellular recordings in the hippocampus in vivo. J Neurophysiol 84:390-400 [Journal] [PubMed]

Hill DN, Mehta SB, Kleinfeld D (2011) Quality metrics to accompany spike sorting of extracellular signals. J Neurosci 31:8699-705 [PubMed]

Johnson DH, Sinanovic S (2001) Symmetrizing the Kullback-Leibler distance IEEE Trans. Info. Theory

Kullback S, Leibler RA (1951) On information and sufficiency Ann Math Stat 22:79-86

Lewicki MS (1998) A review of methods for spike sorting: the detection and classification of neural action potentials. Network 9:R53-78 [PubMed]

Lytton WW (2008) Computer modelling of epilepsy. Nat Rev Neurosci 9:626-37 [Journal] [PubMed]

Mahalanobis PC (1936) On the generalized distance in statistics Proceedings of the National Institute of Science, India

O'Keefe J, Recce ML (1993) Phase relationship between hippocampal place units and the EEG theta rhythm. Hippocampus 3:317-30 [PubMed]

Olypher AV, Klement D, Fenton AA (2006) Cognitive disorganization in hippocampus: a physiological model of the disorganization in psychosis. J Neurosci 26:158-68 [PubMed]

Park E, Dvorak D, Fenton AA (2011) Ensemble place codes in hippocampus: CA1, CA3, and dentate gyrus place cells have multiple place fields in large environments. PLoS One 6:e22349-705 [PubMed]

Paxinos G, Watson M (2007) The rat brain in stereotaxic coordinates

Press WH, Flannery BP, Teukolsky SA, Vetterling WT (1992) Numerical Recipes in C: The Art of Scientific Computing

Quirk MC, Wilson MA (1999) Interaction between spike waveform classification and temporal sequence detection. J Neurosci Methods 94:41-52 [PubMed]

Quiroga RQ, Nadasdy Z, Ben-Shaul Y (2004) Unsupervised spike detection and sorting with wavelets and superparamagnetic clustering. Neural Comput 16:1661-87 [PubMed]

Schmitzer-Torbert N, Jackson J, Henze D, Harris K, Redish AD (2005) Quantitative measures of cluster quality for use in extracellular recordings. Neuroscience 131:1-11 [PubMed]

Schmitzer-Torbert N, Redish AD (2004) Neuronal activity in the rodent dorsal striatum in sequential navigation: separation of spatial and reward responses on the multiple T task. J Neurophysiol 91:2259-72 [PubMed]

Shoham S, Fellows MR, Normann RA (2003) Robust, automatic spike sorting using mixtures of multivariate t-distributions. J Neurosci Methods 127:111-22 [PubMed]

Thompson LT, Best PJ (1990) Long-term stability of the place-field activity of single units recorded from the dorsal hippocampus of freely behaving rats. Brain Res 509:299-308 [PubMed]

Tolias AS, Ecker AS, Siapas AG, Hoenselaar A, Keliris GA, Logothetis NK (2007) Recording chronically from the same neurons in awake, behaving primates. J Neurophysiol 98:3780-90

Wang Q, Kulkarni SR, Verdu S (2006) A Nearest-Neighbor Approach to Estimating Divergence between Continuous Random Vectors ISIT 2006 Seattle, USA

Wheeler BC (1999) Automatic Discrimination of Single Units chapter in Methods for Neural Ensemble Recordings

Wilson MA, McNaughton BL (1993) Dynamics of the hippocampal ensemble code for space. Science 261:1055-8 [PubMed]

Wood F, Black MJ, Vargas-Irwin C, Fellows M, Donoghue JP (2004) On the variability of manual spike sorting. IEEE Trans Biomed Eng 51:912-8

(36 refs)