Monthly
288 pp. per issue
6 x 9, illustrated
ISSN
0899-7667
E-ISSN
1530-888X
2014 Impact factor:
2.21

Neural Computation

May 2010, Vol. 22, No. 5, Pages 1245-1271
(doi: 10.1162/neco.2009.07-08-823)
© 2009 Massachusetts Institute of Technology
Decoding Poisson Spike Trains by Gaussian Filtering
Article PDF (460.86 KB)
Abstract

The temporal waveform of neural activity is commonly estimated by low-pass filtering spike train data through convolution with a gaussian kernel. However, the criteria for selecting the gaussian width σ are not well understood. Given an ensemble of Poisson spike trains generated by an instantaneous firing rate function λ(t), the problem was to recover an optimal estimate of λ(t) by gaussian filtering. We provide equations describing the optimal value of σ using an error minimization criterion and examine how the optimal σ varies within a parameter space defining the statistics of inhomogeneous Poisson spike trains. The process was studied both analytically and through simulations. The rate functions λ(t) were randomly generated, with the three parameters defining spike statistics being the mean of λ(t), the variance of λ(t), and the exponent α of the Fourier amplitude spectrum 1/fα of λ(t). The value of σopt followed a power law as a function of the pooled mean interspike interval I, σopt = aIb, where a was inversely related to the coefficient of variation CV of λ(t), and b was inversely related to the Fourier spectrum exponent α. Besides applications for data analysis, optimal recovery of an analog signal waveform λ(t) from spike trains may also be useful in understanding neural signal processing in vivo.