Optimal filtering for spike sorting of multi-site electrode recordings.

We derive an optimal linear filter, to reduce the distortions of the peak amplitudes of action potentials in extracellular multitrode recordings, which are due to background activity and overlapping spikes. This filter is being learned very efficiently from the raw recordings in an unsupervised manner and responds to the average waveform with an impulse of minimal width. The average waveform does not have to be known in advance, but is learned together with the optimal filter. The peak amplitude of a filtered waveform is a more reliable estimate for the amplitude of an action potential than the peak of the biphasic waveform and can improve the accuracy of the event detection and clustering procedures. We demonstrate a spike-sorting application, in which events are detected using the Mahalanobis distance in the N-dimensional space of filtered recordings as a distance measure, and the event amplitudes of the filtered recordings are clustered to assign events to individual units. This method is fast and robust, and we show its performance by applying it to real tetrode recordings of spontaneous activity in the visual cortex of an anaesthetized cat and to realistic artificial data derived therefrom.

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

[2]  Bruce L. McNaughton,et al.  The stereotrode: A new technique for simultaneous isolation of several single units in the central nervous system from multiple unit records , 1983, Journal of Neuroscience Methods.

[3]  J.P. Miller,et al.  Optimal discrimination and classification of neuronal action potential waveforms from multiunit, multichannel recordings using software-based linear filters , 1994, IEEE Transactions on Biomedical Engineering.

[4]  B. McNaughton,et al.  Tetrodes markedly improve the reliability and yield of multiple single-unit isolation from multi-unit recordings in cat striate cortex , 1995, Journal of Neuroscience Methods.

[5]  Terrence J. Sejnowski,et al.  An Information-Maximization Approach to Blind Separation and Blind Deconvolution , 1995, Neural Computation.

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

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

[8]  Aapo Hyvärinen,et al.  A Fast Fixed-Point Algorithm for Independent Component Analysis , 1997, Neural Computation.

[9]  Maneesh Sahani,et al.  Tetrodes for monkeys , 1997 .

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

[11]  Jean-Francois Cardoso,et al.  Blind signal separation: statistical principles , 1998, Proc. IEEE.

[12]  Sonja Grün,et al.  Detecting unitary events without discretization of time , 1999, Journal of Neuroscience Methods.

[13]  P. Maher,et al.  Handbook of Matrices , 1999, The Mathematical Gazette.

[14]  Richard A. Andersen,et al.  Latent variable models for neural data analysis , 1999 .

[15]  Maneesh Sahani,et al.  Simultaneous paired intracellular and tetrode recordings for evaluating the performance of spike sorting algorithms , 1999, Neurocomputing.

[16]  Kenneth D. Miller,et al.  Cross-channel correlations in tetrode recordings: implications for spike-sorting , 1999, Neurocomputing.

[17]  Terrence J. Sejnowski,et al.  Blind source separation of more sources than mixtures using overcomplete representations , 1999, IEEE Signal Processing Letters.

[18]  Terrence J. Sejnowski,et al.  Learning Overcomplete Representations , 2000, Neural Computation.

[19]  Christophe Andrieu,et al.  Bayesian deconvolution of noisy filtered point processes , 2001, IEEE Trans. Signal Process..

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

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

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

[23]  Christophe Pouzat,et al.  Improved spike-sorting by modeling firing statistics and burst-dependent spike amplitude attenuation: a Markov chain Monte Carlo approach. , 2004, Journal of neurophysiology.