The spontaneous-rate histogram of the auditory nerve can be explained by only two or three spontaneous rates and...

Estimates of the spontaneous discharge rate (SR) of auditory-nerve (AN) fibers are often based on measurements of the average rate over a long (e.g., 30 s) interval. These measurements are important because SR is apparently correlated with other AN properties, such as threshold to acoustic stimuli, shape of rate-level function, recovery from prior stimulation, and certain anatomical characteristics. Furthermore, histograms of SR estimates from large numbers of fibers suggest that they can be divided into two (i.e., low and high) or three (i.e., low, medium, and high) SR classes. Yet, even "simple" statistical estimates, such as average rate, can behave surprisingly poorly for processes with long-range dependence (LRD), which has been found in the spontaneous activity of AN fibers. In particular, LRD greatly increases the variability of estimates of mean discharge rate. We investigated the implications of this effect of LRD for our understanding of the SRs of AN fibers. The fractional-Gaussian-noise-driven Poisson process (fGnDP) was originally developed to model the LRD action-potential trains of AN fibers. Using rate estimates computed from this model, we were able to reproduce the shape of published histograms of SR using only three fixed SR values. Moreover, by using a Poisson-equivalent integrate-and-fire (IF) model in place of the inhomogeneous Poisson process in the fGnDP model, we were able to reproduce SR histograms using only two fixed SR values. These results suggest that AN fibers may have only two or three possible values for their long-term average spontaneous discharge rates. In other words, all "high-SR" neurons may actually have the same underlying SR. Furthermore, both "low-SR" and "medium-SR" neurons may have a single "true" SR value, or these two classes may have two different "true" SR values. Furthermore, the Poisson-equivalent IF model may prove useful in other applications involving the modeling of trains of action potentials.