Pu-Ming Zhang | Search PubMed

PubMed2006

A robust method for spike sorting with automatic overlap decomposition.

Spike sorting is the mandatory first step in analyzing multiunit recording signals for studying information processing mechanisms within the nervous system. Extracellular recordings usually contain overlapped spikes produced by a number of neurons adjacent to the electrode, together with unknown background noise, which in turn induce some difficulties in neural signal identification. In this paper, we propose a robust method to deal with these problems, which employs an automatic overlap decomposition technique based on the relaxation algorithm that requires simple fast Fourier transforms. The performance of the presented system was tested at various signal-to-noise ratio levels based on synthetic data that were generated from real recordings.

Action Potentials↗

PubMed2004

Spike sorting based on automatic template reconstruction with a partial solution to the overlapping problem.

A new method for spike sorting is proposed which partly solves the overlapping problem. Principal component analysis and subtractive clustering techniques are used to estimate the number of neurons contributing to multi-unit recording. Spike templates (i.e. waveforms) are reconstructed according to the clustering results. A template-matching procedure is then performed. Firstly all temporally displaced templates are compared with the spike event to find the best-fitting template that yields the minimum residue variance. If the residue passes the chi(2)-test, the matching procedure stops and the spike event is classified as the best-fitting template. Otherwise the spike event may be an overlapping waveform. The procedure is then repeated with all possible combinations of two templates, three templates, etc. Once one combination is found, which yields the minimum residue variance among the combinations of the same number of component templates and makes the residue pass the chi(2)-test, the matching procedure stops. It is unnecessary to check the remaining combinations of more templates. Consequently, the computational effort is reduced and the over-fitting problem can be partly avoided. A simulated spike train was used to assess the performance of the proposed method, which was also applied to a real recording of chicken retina ganglion cells.

Action Potentials↗