A qualitative comparison of some diffusion models for neural activity via stochastic ordering

A number of diffusion processes have been proposed as a continuous analog of Stein's model for the subthreshold membrane potential of a neuron. Interspike intervals are then described as the first-passage-time of the corresponding diffusion model through a suitable threshold. Various biological considerations suggest the use of more sophisticated models in lieu of the Ornstein-Uhlenbeck model. However, the advantages of the additional complexity are not always clear. Comparisons among different models generally use numerical methods in specific examples without a general sensitivity analysis on the role of the model parameters. Here, we compare the distribution of interspike intervals from different models using the method of stochastic ordering. The qualitative comparison of the role of each parameter extends the results obtained from numerical simulations. One result on neurons with high positive net excitation is that the reversal potential models considered do not greatly differ from the Ornstein-Uhlenbeck model. For neurons with increased inhibition, the models give greater differences among the interspike interval distributions. In particular, when the mean trajectories are matched, the Feller model gives shorter times than the Ornstein-Uhlenbeck model but longer times than our double reversal potential model.