Home Estimation of neural firing rate: the wavelet density estimation approach
Article
Licensed
Unlicensed Requires Authentication

Estimation of neural firing rate: the wavelet density estimation approach

  • Abed Khorasani and Mohammad Reza Daliri EMAIL logo
Published/Copyright: August 8, 2013

Abstract

The computation of neural firing rates based on spike sequences has been introduced as a useful tool for extraction of an animal’s behavior. Different estimating methods of such neural firing rates have been developed by neuroscientists, and among these methods, time histogram and kernel estimators have been used more than other approaches. In this paper, the problem in the estimation of firing rates using wavelet density estimators has been considered. The results of simulation study in estimation of underlying rates based on spike sequences sampled from two different variable firing rates show that the proposed wavelet density method provides better and more accurate estimation of firing rates with smooth results compared to two other classical approaches. Furthermore, the performance of a different family of wavelet density estimators in the estimation of the underlying firing rate of biological data have been compared with results of both time histogram and kernel estimators. All in all, the results show that the proposed method can be useful in the estimation of firing rate of neural spike trains.


Corresponding author: Mohammad Reza Daliri, Biomedical Engineering Department, Faculty of Electrical Engineering, Iran University of Science and Technology (IUST), Narmak, 16846-13114 Tehran, Iran, Phone: +98-21-73225738, Fax: +98-2173225777

References

[1] Abeles M. Quantification, smoothing, and confidence-limits for single-units histograms. J Neurosci Methods 1982; 5: 317–325.10.1016/0165-0270(82)90002-4Search in Google Scholar

[2] Bair W, Cavanaugh JR, Smith MA, Movshon JA. The timing of response onset and offset in macaque visual neurons. J Neurosci 2002; 22: 3189–3205.10.1523/JNEUROSCI.22-08-03189.2002Search in Google Scholar

[3] Breiman L, Meisel W, Purcell E. Variable kernel estimates of multivariate densities. Technometrics 1977; 19: 135–144.10.1080/00401706.1977.10489521Search in Google Scholar

[4] Brown SL, Joseph J, Stopfer M. Encoding a temporally structured stimulus with a temporally structured neural representation. Nat Neurosci 2005; 8: 1568–1576.10.1038/nn1559Search in Google Scholar

[5] Brown EN, Kass RE, Mitra PP. Multiple neural spike train data analysis: state-of-the-art and future challenges. Nat Neurosci 2004; 7: 456–461.10.1038/nn1228Search in Google Scholar

[6] Cencov NN. Evaluation of an unknown distribution density from observations. Soviet Math Dokl 1962; 3: 1559–1562.Search in Google Scholar

[7] Dayan P, Abbott L. Theoretical neuroscience: computational and mathematical modeling of neural systems. Cambridge, MA: MIT Press 2001.Search in Google Scholar

[8] DiMatteo I, Genovese C, Kass R. Bayesian curve-fitting with free-knot splines. Biometrika 2001; 88: 1055–1071.10.1093/biomet/88.4.1055Search in Google Scholar

[9] Donoho DL, Kerkyacharian IM, Picard D. Density estimation by wavelet thresholding. Ann Stat 1996; 24: 508–539.10.1214/aos/1032894451Search in Google Scholar

[10] Gerstein G, Kiang NYS. An approach to the quantitative analysis of electrophysiological data from single neurons. Biophys J 1960; 1: 15–28.10.1016/S0006-3495(60)86872-5Search in Google Scholar

[11] Gerstein GL, Perkel DH. Simultaneously recorded trains of action potentials—analysis and functional interpretation. Science 1969; 3: 828–830.Search in Google Scholar

[12] Habbema J, Remme J. Variable kernel density estimation in discriminant analysis. Compstat. Proceedings in Computational Statistics. Vienna: Physica Verlag 1978.Search in Google Scholar

[13] Levi R, Varona R, Arshavsky YI, Rabinovich MI, Selverston AI. The role of sensory network dynamics in generating a motor program. J Neurosci 2005; 25: 9807–9815.10.1523/JNEUROSCI.2249-05.2005Search in Google Scholar

[14] Nawrot M, Aertsen A, Rotter S. Single-trial estimation of neuronal firing rates: from single-neuron spike trains to population activity. J Neurosci Methods 1999; 94: 81–92.10.1016/S0165-0270(99)00127-2Search in Google Scholar

[15] Parzen E. Estimation of a probability density-function and mode. Ann Math Stat 1962; 33: 1065–1076.10.1214/aoms/1177704472Search in Google Scholar

[16] Peter A, Rangarajan A. Maximum likelihood wavelet density estimation with applications to image and shape matching. IEEE Trans Image Process 2008; 17: 458–468.10.1109/TIP.2008.918038Search in Google Scholar

[17] Richmond BJ, Optican LM, Spitzer H. Temporal encoding of two-dimensional patterns by single units in primate primary visual cortex. I. Stimulus-response relations. J Neurophysiol 1990; 64: 351–369.10.1152/jn.1990.64.2.351Search in Google Scholar

[18] Sanderson A. Adaptive filtering of neuronal spike train data. IEEE Trans Biomed Eng 1980; 27: 271–274.10.1109/TBME.1980.326633Search in Google Scholar

[19] Shimazaki H, Shinomoto S. A method for selecting the bin size of a time histogram. Neural Comput 2007; 19: 1503–1527.10.1162/neco.2007.19.6.1503Search in Google Scholar

[20] Shimazaki H, Shinomoto S. Kernel bandwidth optimization in spike rate estimation. J Comput Neurosci 2010; 29: 171–182.10.1007/s10827-009-0180-4Search in Google Scholar

[21] Vannucci M, Shinomoto S. Nonparametric density estimation using wavelets. ISDS, Duke University 1998: 26–95.Search in Google Scholar

[22] Wiener M, Richmond B. Model based decoding of spike trains. Biosystems 2002; 67: 295–300.10.1016/S0303-2647(02)00087-4Search in Google Scholar

Received: 2013-1-8
Accepted: 2013-7-24
Published Online: 2013-08-08
Published in Print: 2013-08-01

©2013 by Walter de Gruyter Berlin Boston

Downloaded on 12.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/bmt-2013-0060/html
Scroll to top button