An application of reversible-jump Markov chain Monte Carlo to spike classification of multi-unit extracellular recordings.
- 16 January 2003
- journal article
- research article
- Published by Taylor & Francis Ltd in Network: Computation in Neural Systems
- Vol. 14 (1), 61-82
- https://doi.org/10.1088/0954-898x/14/1/304
Abstract
Multi-electrode recordings in neural tissue contain the action potential waveforms of many closely spaced neurons. While we can observe the action potential waveforms, we cannot observe which neuron is the source for which waveform nor how many source neurons are being recorded. Current spike-sorting algorithms solve this problem by assuming a fixed number of source neurons and assigning the action potentials given this fixed number. We model the spike waveforms as an anisotropic Gaussian mixture model and present, as an alternative, a reversible-jump Markov chain Monte Carlo (MCMC) algorithm to simultaneously estimate the number of source neurons and to assign each action potential to a source. We derive this MCMC algorithm and illustrate its application using simulated three-dimensional data and real four-dimensional feature vectors extracted from tetrode recordings of rat entorhinal cortex neurons. In the analysis of the simulated data our algorithm finds the correct number of mixture components (sources) and classifies the action potential waveforms with minimal error. In the analysis of real data, our algorithm identifies clusters closely resembling those previously identified by a user-dependent graphical clustering procedure. Our findings suggest that a reversible-jump MCMC algorithm could offer a new strategy for designing automated spike-sorting algorithms.Keywords
This publication has 15 references indexed in Scilit:
- Particle filters for state-space models with the presence of unknown static parametersIEEE Transactions on Signal Processing, 2002
- A Global Geometric Framework for Nonlinear Dimensionality ReductionScience, 2000
- Trajectory Encoding in the Hippocampus and Entorhinal CortexNeuron, 2000
- Multi-channel spike detection and sorting using an array processing techniqueNeurocomputing, 1999
- Reversible jump Markov chain Monte Carlo computation and Bayesian model determinationBiometrika, 1995
- Bayesian Density Estimation and Inference Using MixturesJournal of the American Statistical Association, 1995
- Bayesian Modeling and Classification of Neural SignalsNeural Computation, 1994
- A Bayesian method for classification and discriminationThe Canadian Journal of Statistics / La Revue Canadienne de Statistique, 1992
- Illustration of Bayesian Inference in Normal Data Models Using Gibbs SamplingJournal of the American Statistical Association, 1990
- Monte Carlo sampling methods using Markov chains and their applicationsBiometrika, 1970