Unsupervised Spike Detection and Sorting with Wavelets and Superparamagnetic Clustering

This study introduces a new method for detecting and sorting spikes from multiunit recordings. The method combines the wave let transform, which localizes distinctive spike features, with super paramagnetic clustering, which allows automatic classification of the data without assumptions such as low variance or gaussian distributions. Moreover, an improved method for setting amplitude thresholds for spike detection is proposed. We describe several criteria for implementation that render the algorithm unsupervised and fast. The algorithm is compared to other conventional methods using several simulated data sets whose characteristics closely resemble those of in vivo recordings. For these data sets, we found that the proposed algorithm outperformed conventional methods.

[1]  M. Abeles,et al.  Multispike train analysis , 1977, Proceedings of the IEEE.

[2]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[3]  U. Wolff.  Comparison Between Cluster Monte Carlo Algorithms in the Ising Model , 1989 .

[4]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[5]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[6]  Amara Lynn Graps,et al.  An introduction to wavelets , 1995 .

[7]  V. Samar,et al.  Multiresolution Analysis of Event-Related Potentials by Wavelet Decomposition , 1995, Brain and Cognition.

[8]  D. Kleinfeld,et al.  Variability of extracellular spike waveforms of cortical neurons. , 1996, Journal of neurophysiology.

[9]  Eytan Domany,et al.  Superparamagnetic Clustering of Data , 1996 .

[10]  Eytan Domany,et al.  Data Clustering Using a Model Granular Magnet , 1997, Neural Computation.

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

[12]  Eytan Domany,et al.  Superparamagnetic clustering of data — The definitive solution of an ill-posed problem , 1999 .

[13]  Erkki Oja,et al.  Independent component analysis: algorithms and applications , 2000, Neural Networks.

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

[15]  E Hulata,et al.  Detection and sorting of neural spikes using wavelet packets. , 2000, Physical review letters.

[16]  J. Letelier,et al.  Spike sorting based on discrete wavelet transform coefficients , 2000, Journal of Neuroscience Methods.

[17]  R Quian Quiroga,et al.  Wavelet Transform in the analysis of the frequency composition of evoked potentials. , 2001, Brain research. Brain research protocols.

[18]  Gilles Laurent,et al.  Using noise signature to optimize spike-sorting and to assess neuronal classification quality , 2002, Journal of Neuroscience Methods.

[19]  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.

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .

[21]  R Quian Quiroga,et al.  Performance of different synchronization measures in real data: a case study on electroencephalographic signals. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  R. Quian Quiroga,et al.  Single-trial event-related potentials with wavelet denoising , 2003, Clinical Neurophysiology.