Implications for multi-neuron recording applications Accuracy of spike sorter is crucial to retaining information in single-unit responses. Error tolerance: to maintain 80% information - < 2% to maintain 50% information - < 13-14%
were added at a percent error rate of
% error, ε
Figure 6) Misclassification error. Exponential curve (solid) fit to data to model İ(s,ŝ) (top) and coding fraction (bottom) as a function of total percent spike sorting error. Dashed lines demarcate where error causes information rate to fall by half its information rate (top) and where coding fraction crosses 0 (bottom).
% error, ε
Figure 7) Comparing effect of error on information rate (bottom) when average firing rate (top) is kept constant and when it is changing linearly.
% error, ε
Figure 5)Effect of false positive (x’s) and false negative (dots) misdetection error on information rate (top)coding fraction (middle), and average firing rate (bottom).
Figure 3) The optimal linear estimate of the stimulus was computed by applying the Wiener-Hopf filter to ρ(t) . From the error between s and ŝ, the signal-to-noise ratio (SNR) was computed and İ(s, ρ) was estimated by İ(s, ŝ). Coding fraction was also calculated from the mean square error and standard deviation of the stimulus,
according to .
Figure 4) Stimulus signal s(t) (top solid) is input into IAF neuron which outputs a spike train. The spike train is contaminated by random sorting error. Output with no error is shown in bottom panel. The spike train is then input into an optimal reconstruction filter to yield stimulus estimate ŝ(t) (top dashed).
Deborah S. Won, David Y. Chong, and Patrick D. Wolf
Department of Biomedical Engineering, Duke University, Durham, NC 27708
EFFECTS OF SPIKE SORTING ERROR ON INFORMATION CONTENT IN MULTI-NEURON RECORDINGS
Figure 1) The current prototype for multi-neuron acquisition and processing systems involves sorting the multi-unit signal on each channel into single-unit signals before most of the data is discarded and only the onset times of action potentials from each unit are retained. The goal is to retain as much of the relevant information as possible in the original spike trains. The information may be about the stimulus that elicited the response or the behavior that was to be elicited by the neural activity.. Ideally, the spike output from the neurons are reproduced at the output of the acquisition system, and the information transmitted to the computational decoders remains unchanged. However, spike sorting does not exactly replicate the spike trains due to the error inevitably associated with spike sorting.
The aim of this investigation is to quantify the effects of incorrectly discriminated spikes on encoded information so as to set tolerable limits on sorting error and glean insight into the necessity of sorting. To address how sorting error affects the ability to decode the multi-neuron response, spike trains were simulated with random detection error and with classification error and changes in the Shannon information content were observed.
 Y. Tock and G. Inbar, Modern Techniques in Neuroscience Research, Springer-Verlag, 1999.
 F. Rieke, D.A. Bodnar, and W. Bialek, Proceedings of the Royal Society London B (1995) 262: 259-265.
 R. Wessel, C. Koch, and F. Gabbiani, Journal of Neurophysiology (1996) 75: 2280-2293.
 S. Seung, http://hebb.mit.edu/courses/9.29/lectures/ lecture02.html.
 K. D. Harris, D. A. Henze, J. Csicsvari, H. Hirase, and G. Buzsaki, Journal of Neurophysiology (2000) 84: 401-414.
 B. Wheeler and W. Heetderks, IEEE Trans on Biomedical Engineering (1982).
Methods To compute the information transmission rate İ(s, ρ) and coding fraction γ, we utilize a signal-reconstruction paradigm commonly used in information-theoretic analysis -. Average values were found for results from 10 unique input signals.
Information rate is fit well by exponential curve for percent errors > roughly 2% (R2 > .9). For smaller percent errors, the decay rate is steeper. At percent errors greater than 17% misdetection and 19% misclassification, the coding fraction γ drops below 0. Since γ provides a measure of how much of the stimulus variance is accounted for by the estimate variance, the error at which γbecomes less than 0 is the point beyond which spikes act as completely random occurrences, uncorrelated with the stimulus input. Information per spike also changed similarly to information rate starting at 0.8-0.9 bits/spike and decreasing to an asymptotic value of 0.1-0.2 bits/spike. Figures display results for one trial.
Funding provided by National Science Foundation and DARPA contract #N66-001-02-C-8022.