with the EM method and the Kalman filter.
You may obtain the list of {time, estimated rate, estimated regularity}.
BayesRR.zip
Function 'BayesRR' returns the list of {time, estimated rate, estimated regularity}.
- Estimating instantaneous irregularity of neuronal firing
T. Shimokawa and S. Shinomoto, Neural Computation (2009) 21:1931-1951.- A characterization of the time-rescaled gamma process as a model for spike trains
T. Shimokawa, S. Koyama, and S. Shinomoto, Journal of Computational Neuroscience (2010) DOI: 10.1007/s10827-009-0194-y.
If you have any questions, or have suggestions for improving the programs, please contact Shigeru Shinomoto or Takeaki Shimokawa.