Parameters of stochastic diffusion processes estimated from observations of first-hitting times: application to the leaky integrate-and-fire neuronal model.

A theoretical model has to stand the test against the real world to be of any practical use. The first step is to identify parameters in the model estimated from experimental data. In many applications where renewal point data are available, models of first-hitting times of underlying diffusion processes arise. Despite the seemingly simplicity of the model, the problem of how to estimate parameters of the underlying stochastic process has resisted solution. The few attempts have either been unreliable, difficult to implement, or only valid in subsets of the relevant parameter space. Here we present an estimation method that overcomes these difficulties, is computationally easy and fast to implement, and also works surprisingly well on small data sets. The method is illustrated on simulated and experimental data. Two common neuronal models--the Ornstein-Uhlenbeck and Feller models--are investigated.