Truncation thresholds: a pair of spike detection thresholds computed using truncated probability distributions

We describe a method for computing a pair of spike detection thresholds, called `truncation thresholds', using truncated probability distributions, for extracellular recordings. In existing methods the threshold is usually set to a multiple of an estimate of the standard deviation of the noise in the recording, with the multiplication factor being chosen between 3 and 5 according to the researcher's preferences. Our method has the following advantages over these methods. First, because the standard deviation is usually estimated from the entire recording, which includes the spikes, it increases with ring rate. By contrast, truncation thresholds decrease in absolute value with increasingring rate, thereby capturing more of the signal. Second, the parameters of the selected noise distribution are estimated more accurately by maximum likelihood tting of the truncated distribution to the data delimited by the truncation thresholds. Third, the computation of the truncation thresholds is completely data-driven. It does not involve a user- de ned multiplication factor. Fourth, methods that use a threshold that is proportional to the estimated standard deviation of the noise assume that the noise distribution is symmetrical around the mean. By contrast, truncation thresholds are not linked to each other by an assumption of symmetry about some axis. Fifth, existing methods do not verify that subthreshold data obey a noise distribution. Truncation thresholds, however, are de ned by the fact that the distribution of the data they delimit is statistically indistinguishable, according to the Kolmogorov{Smirnov test, from a selected distribution, truncated at those thresholds. Application of the method is illustrated using recordings from cortical area M1 in awake behaving rats, as well as in simulated recordings. Source code and executables of a software suite that computes the truncation thresholds are provided for the case when the noise distribution is modeled as truncated normal.

PDF

___

[1] Chorev E, Epsztein J, Houweling AR, Lee AK, Brecht M. Electrophysiological recordings from behaving animals{ going beyond spikes. Curr Opin Neurobiol 2009; 19: 513-519.
[2] Obien ME, Deligkaris K, Bullmann T, Bakkum DJ, Frey U. Revealing neuronal function through microelectrode array recordings. Front Neurosci 2015; 8: 423.
[3] Buzsaki G. Large-scale recording of neuronal ensembles. Nat Neurosci 2004; 7: 446-451.
[4] Engel AK, Moll CK, Fried I, Ojemann GA. Invasive recordings from the human brain: clinical insights and beyond. Nat Rev Neurosci 2005; 6: 35-47.
[5] Rowland NC, Breshears J, Chang EF. Neurosurgery and the dawning age of Brain-Machine Interfaces. Surg Neurol Int 2013; 4 (Suppl. 1): S11-S14.
[6] Lewicki MS. A review of methods for spike sorting: the detection and classi cation of neural action potentials. Network 1998; 9: R53-R78.
[7] Quiroga RQ, Nadasdy Z, Ben-Shaul Y. Unsupervised spike detection and sorting with wavelets and superparamag- netic clustering. Neural Comput 2004; 16: 1661-1687.
[8] Vargas-Irwin C, Donoghue JP. Automated spike sorting using density grid contour clustering and subtractive waveform decomposition. J Neurosci Methods 2007; 164: 1-18.
[9] Takekawa T, Isomura Y, Fukai T. Accurate spike sorting for multi-unit recordings. Eur J Neurosci 2010; 31: 263-272.
[10] Jackel D, Frey U, Fiscella M, Franke F, Hierlemann A. Applicability of independent component analysis on high- density microelectrode array recordings. J Neurophysiol 2012; 108: 334-348.
[11] Todorova S, Sadtler P, Batista A, Chase S, Ventura V. To sort or not to sort: the impact of spike-sorting on neural decoding performance. J Neural Eng 2014; 11: 056005.
[12] Okatan M, Kocaturk, M. Action potential detection in extracellular recordings using truncation thresholds". In: Medical Technologies National Conference (TIPTEKNO); 15{18 October 2015; Bodrum, Turkey. New York, NY, USA: IEEE. pp. 1-4 (in Turkish with English abstract).
[13] Okatan M, Kocaturk, M. Firing rate dependence of truncation thresholds". In: Medical Technologies National Conference (TIPTEKNO); 15{18 October 2015; Bodrum, Turkey. New York, NY, USA: IEEE. pp. 1-4 (in Turkish with English abstract).
[14] Thakur PH, Lu H, Hsiao SS, Johnson KO. Automated optimal detection and classi cation of neural action potentials in extra-cellular recordings. J Neurosci Methods 2007; 162: 364-376.
[15] Yang Z, Liu W, Keshtkaran MR, Zhou Y, Xu J, Pikov V, Guan C, Lian Y. A new EC-PC threshold estimation method for in vivo neural spike detection. J Neural Eng 2012; 9: 046017.
[16] Press WH, Flannery BB, Teukolsky SA, Vetterling WT. Numerical Recipes in C: The Art of Scienti c Computing. 2nd ed. Cambridge, UK: Cambridge University Press, 1992.
[17] Kocaturk M, Gulcur HO, Canbeyli R. Toward building hybrid biological/in silico neural networks for motor neuroprosthetic control. Frontiers in Neurorobotics 2015; 9: 8.
[18] Donoho D, Johnstone IM. Ideal spatial adaptation by wavelet shrinkage. Biometrika 1994; 81: 425-455.
[19] Uzun R, Ozer M, Perc M. Can scale-freeness offset delayed signal detection in neuronal networks? Europhys Lett 2014; 105: 60002.
[20] Perc M, Marhl M. Ampli cation of information transfer in excitable systems that reside in a steady state near a bifurcation point to complex oscillatory behaviour. Phys Rev E 2005; 71: 026229.