Joint modeling of observed inter-arrival times and waveform data with multiple hidden states for neural spike-sorting