A nonparametric Bayesian alternative to spike sorting

The analysis of extra-cellular neural recordings typically begins with careful spike sorting and all analysis of the data then rests on the correctness of the resulting spike trains. In many situations this is unproblematic as experimental and spike sorting procedures often focus on well isolated units. There is evidence in the literature, however, that errors in spike sorting can occur even with carefully collected and selected data. Additionally, chronically implanted electrodes and arrays with fixed electrodes cannot be easily adjusted to provide well isolated units. In these situations, multiple units may be recorded and the assignment of waveforms to units may be ambiguous. At the same time, analysis of such data may be both scientifically important and clinically relevant. In this paper we address this issue using a novel probabilistic model that accounts for several important sources of uncertainty and error in spike sorting. In lieu of sorting neural data to produce a single best spike train, we estimate a probabilistic model of spike trains given the observed data. We show how such a distribution over spike sortings can support standard neuroscientific questions while providing a representation of uncertainty in the analysis. As a representative illustration of the approach, we analyzed primary motor cortical tuning with respect to hand movement in data recorded with a chronic multi-electrode array in non-human primates. We found that the probabilistic analysis generally agrees with human sorters but suggests the presence of tuned units not detected by humans.

[1]  G. Radons,et al.  Analysis, classification, and coding of multielectrode spike trains with hidden Markov models , 2004, Biological Cybernetics.

[2]  Thomas L. Griffiths,et al.  Particle Filtering for Nonparametric Bayesian Matrix Factorization , 2006, NIPS.

[3]  Paul Fearnhead,et al.  Particle filters for mixture models with an unknown number of components , 2004, Stat. Comput..

[4]  S. Roweis,et al.  Time-Varying Topic Models using Dependent Dirichlet Processes , 2005 .

[5]  Carl E. Rasmussen,et al.  The Infinite Gaussian Mixture Model , 1999, NIPS.

[6]  L. Frank,et al.  An application of reversible-jump Markov chain Monte Carlo to spike classification of multi-unit extracellular recordings. , 2003, Network.

[7]  J. Pitman Combinatorial Stochastic Processes , 2006 .

[8]  Eric R. Ziecel Aspects of Uncertainty , 1995 .

[9]  J. Csicsvari,et al.  Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements. , 2000, Journal of neurophysiology.

[10]  Michael A. West,et al.  Hierarchical priors and mixture models, with applications in regression and density estimation , 2006 .

[11]  Radford M. Neal Connectionist Learning of Belief Networks , 1992, Artif. Intell..

[12]  Carl E. Rasmussen,et al.  Modelling Spikes with Mixtures of Factor Analysers , 2004, DAGM-Symposium.

[13]  Matthew Fellows,et al.  On the variability of manual spike sorting , 2004, IEEE Transactions on Biomedical Engineering.

[14]  Matthew Fellows,et al.  Robustness of neuroprosthetic decoding algorithms , 2003, Biological Cybernetics.

[15]  Michael J. Black,et al.  Automatic spike sorting for neural decoding , 2004, The 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[16]  R. Segev,et al.  A method for spike sorting and detection based on wavelet packets and Shannon's mutual information , 2002, Journal of Neuroscience Methods.

[17]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[18]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[19]  S. MacEachern,et al.  A semiparametric Bayesian model for randomised block designs , 1996 .

[20]  Partha P. Mitra,et al.  Automatic sorting of multiple unit neuronal signals in the presence of anisotropic and non-Gaussian variability , 1996, Journal of Neuroscience Methods.

[21]  Shy Shoham,et al.  Robust, automatic spike sorting using mixtures of multivariate t-distributions , 2003, Journal of Neuroscience Methods.

[22]  Wei Wu,et al.  Bayesian Population Decoding of Motor Cortical Activity Using a Kalman Filter , 2006, Neural Computation.

[23]  Richard A. Andersen,et al.  On the Separation of Signals from Neighboring Cells in Tetrode Recordings , 1997, NIPS.

[24]  S. MacEachern,et al.  Estimating mixture of dirichlet process models , 1998 .

[25]  Mati Joshua,et al.  Quantifying the isolation quality of extracellularly recorded action potentials , 2007, Journal of Neuroscience Methods.

[26]  M S Lewicki,et al.  A review of methods for spike sorting: the detection and classification of neural action potentials. , 1998, Network.

[27]  Thomas L. Griffiths,et al.  A more rational model of categorization , 2006 .

[28]  Michael I. Jordan,et al.  An Introduction to Variational Methods for Graphical Models , 1999, Machine-mediated learning.

[29]  A P Georgopoulos,et al.  On the relations between the direction of two-dimensional arm movements and cell discharge in primate motor cortex , 1982, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[30]  M. Escobar,et al.  Markov Chain Sampling Methods for Dirichlet Process Mixture Models , 2000 .

[31]  Adrian E. Raftery,et al.  Bayesian Regularization for Normal Mixture Estimation and Model-Based Clustering , 2007, J. Classif..

[32]  Francesco Ventriglia,et al.  Synaptic fusion pore structure and AMPA receptor activation according to Brownian simulation of glutamate diffusion , 2003, Biological Cybernetics.

[33]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[34]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[35]  J. E. Griffin,et al.  Order-Based Dependent Dirichlet Processes , 2006 .

[36]  Timothy J. Robinson,et al.  Sequential Monte Carlo Methods in Practice , 2003 .

[37]  Thomas L. Griffiths,et al.  Infinite latent feature models and the Indian buffet process , 2005, NIPS.

[38]  Yoshio Sakurai,et al.  Automatic sorting for multi-neuronal activity recorded with tetrodes in the presence of overlapping spikes. , 2003, Journal of neurophysiology.

[39]  P. Fearnhead,et al.  On‐line inference for hidden Markov models via particle filters , 2003 .