Setting Adaptive Spike Detection Threshold for Smoothed TEO Based on Robust Statistics Theory

We propose a novel approach aimed at adaptively setting the threshold of the smoothed Teager energy operator (STEO) detector to be used in extracellular recording of neural signals. In this proposed approach, to set the adaptive threshold of the STEO detector, we derive the relationship between the low-order statistics of its input signal and the ones of its output signal. This relationship is determined with only the background noise component assumed to be present at the input. Robust statistics theory techniques were used to achieve an unbiased estimation of these low-order statistics of the background noise component directly from the neural input signal. In this paper, the emphasis is made on extracellular neural recordings. However, the proposed method can be used in the analysis of different biomedical signals where spikes are important for diagnostic (e.g., ECG, EEG, etc.). We validated the efficacy of the proposed method using synthetic neural signals constructed from real neural recordings signals. Four different sets of extracellular recordings from four distinct neural sources have been exploited to that purpose. The first dataset is recorded from an adult male monkey using the Utath 10×10 microelectrode array implemented in the prefrontal cortex, the second one was obtained from the visual cortex of a rat using a stainless-steel-tipped microelectrode, the third dataset came from recording in a human medial lobe using intracranial electrode, and finally, the fourth one was extracted from recordings in a macaque parietal cortex using a single tetrode. Simulation results show that our approach is effective and robust, and outperforms state-of-the-art adaptive detection methods in its category (i.e., efficient and simple, and do not require a priori knowledge about neural spike waveforms shapes).

[1]  Peter J. Huber,et al.  Robust Statistics , 2005, Wiley Series in Probability and Statistics.

[2]  Steven Kay,et al.  Fundamentals Of Statistical Signal Processing , 2001 .

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

[4]  Ramanathan Gnanadesikan,et al.  Methods for statistical data analysis of multivariate observations , 1977, A Wiley publication in applied statistics.

[5]  M. Genton,et al.  Highly Robust Estimation of Dispersion Matrices , 2001 .

[6]  L. Isserlis ON A FORMULA FOR THE PRODUCT-MOMENT COEFFICIENT OF ANY ORDER OF A NORMAL FREQUENCY DISTRIBUTION IN ANY NUMBER OF VARIABLES , 1918 .

[7]  J. Gower,et al.  Methods for statistical data analysis of multivariate observations , 1977, A Wiley publication in applied statistics.

[8]  R.R. Harrison,et al.  A Low-Power Integrated Circuit for a Wireless 100-Electrode Neural Recording System , 2006, IEEE Journal of Solid-State Circuits.

[9]  I. Miller Probability, Random Variables, and Stochastic Processes , 1966 .

[10]  R. Quian Quiroga,et al.  Unsupervised Spike Detection and Sorting with Wavelets and Superparamagnetic Clustering , 2004, Neural Computation.

[11]  P. Rousseeuw,et al.  Alternatives to the Median Absolute Deviation , 1993 .

[12]  M. Sawan,et al.  An Ultra Low-Power CMOS Automatic Action Potential Detector , 2009, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

[13]  R. Northrop Signals and Systems Analysis In Biomedical Engineering , 2003 .

[14]  Min Dai,et al.  A Method for Solving Overlap Problem in Spike Sorting Based on Genetic Algorithm , 2009, 2009 2nd International Congress on Image and Signal Processing.

[15]  Shih-Tseng Lee,et al.  Detection of neuronal spikes using an adaptive threshold based on the max–min spread sorting method , 2008, Journal of Neuroscience Methods.

[16]  Patrick D. Wolf,et al.  Evaluation of spike-detection algorithms fora brain-machine interface application , 2004, IEEE Transactions on Biomedical Engineering.

[17]  S. Mukhopadhyay,et al.  A new interpretation of nonlinear energy operator and its efficacy in spike detection , 1998, IEEE Transactions on Biomedical Engineering.

[18]  I. Bankman,et al.  Optimal detection, classification, and superposition resolution in neural waveform recordings , 1993, IEEE Transactions on Biomedical Engineering.

[19]  I. Obeid Comparison of Spike Detectors based on Simultaneous Intracellular and Extracellular Recordings , 2007, 2007 3rd International IEEE/EMBS Conference on Neural Engineering.

[20]  Patrick D. Wolf,et al.  Optimizing the automatic selection of spike detection thresholds using a multiple of the noise level , 2009, Medical & Biological Engineering & Computing.

[21]  Eivind Kvedalen Signal processing using the Teager Energy Operator and other nonlinear operators , 2003 .

[22]  Petros Maragos,et al.  AM-FM energy detection and separation in noise using multiband energy operators , 1993, IEEE Trans. Signal Process..

[23]  Don H. Johnson,et al.  Fundamental Detection and Estimation Limits in Spike Sorting , 2007, 2007 IEEE International Conference on Acoustics, Speech and Signal Processing - ICASSP '07.

[24]  Bruno O. Shubert,et al.  Random variables and stochastic processes , 1979 .

[25]  Peter J. Rousseeuw,et al.  Time-Efficient Algorithms for Two Highly Robust Estimators of Scale , 1992 .

[26]  Hanzhang Lu,et al.  Automated optimal detection and classification of neural action potentials in extra-cellular recordings , 2007, Journal of Neuroscience Methods.

[27]  Taejeong Kim,et al.  A new action potential detector using the MTEO and its effects on spike sorting systems at low signal-to-noise ratios , 2006, IEEE Transactions on Biomedical Engineering.

[28]  Monson H. Hayes,et al.  Statistical Digital Signal Processing and Modeling , 1996 .

[29]  Joel W. Burdick,et al.  Spike detection using the continuous wavelet transform , 2005, IEEE Transactions on Biomedical Engineering.

[30]  Ramanathan Gnanadesikan Methods for Statistical Data Analysis of Multivariate Observations: Gnanadesikan/Methods , 1997 .