“…1: (1) detection of temporal segments of the voltage trace that are likely to contain spikes, (2) estimation of a set of features for each segment, and (3) classification of the segments according to these features. A variety of methods exist for solving each step (e.g., (1) thresholding based on absolute value (Obeid and Wolf, 2004), squared values (Rutishauser et al, 2006), Teager energy (Choi et al, 2006), or other nonlinear operators (Rebrik et al, 1999), (2) features such as peak-to-peak width/amplitude, projections onto principal components (Lewicki, 1998), or wavelet coefficients (Quiroga et al, 2004;Kwon and Oweiss, 2011), and (3) classification methods such as K-means (Lewicki, 1998), mixture models (Sahani, 1999;Shoham et al, 2003), or superparamagnetic methods (Quiroga et al, 2004)). …”