An improved method for the estimation of firing rate dynamics using an optimal digital filter

In most neural systems, neurons communicate by means of sequences of action potentials or ‘spikes’. Information encoded by spike trains is often quantified in terms of the firing rate which emphasizes the frequency of occurrence of action potentials rather than their exact timing. Common methods for estimating firing rates include the rate histogram, the reciprocal interspike interval, and the spike density function. In this study, we demonstrate the limitations of these aforementioned techniques and propose a simple yet more robust alternative. By convolving the spike train with an optimally designed Kaiser window, we show that more robust estimates of firing rate are obtained for both low and high-frequency inputs. We illustrate our approach by considering spike trains generated by simulated as well as experimental data obtained from single-unit recordings of first-order sensory neurons in the vestibular system. Improvements were seen in the prevention of aliasing, phase and amplitude distortion, as well as in the noise reduction for sinusoidal and more complex input profiles. We review the generality of the approach, and show that it can be adapted to describe neurons with sensory or motor responses that are characterized by marked nonlinearities. We conclude that our method permits more robust estimates of neural dynamics than conventional techniques across all stimulus conditions.