Detection and localization of multiple rate changes in Poisson spike trains
BMC Neuroscience volume 12, Article number: P268 (2011)
In statistical spike train analysis, stochastic point process models usually assume stationarity, in particular that the underlying spike train shows a constant firing rate (e.g. ). However, such models can lead to misinterpretation of the associated tests if the assumption of rate stationarity is not met (e.g. ). Therefore, the analysis of nonstationary data requires that rate changes can be located as precisely as possible. However, present statistical methods focus on rejecting the null hypothesis of stationarity without explicitly locating the change point(s) (e.g. ).
We propose a test for stationarity of a given spike train that can also be used to estimate the change points in the firing rate. Assuming a Poisson process with piecewise constant firing rate, we propose a Step-Filter-Test (SFT) which can work simultaneously in different time scales, accounting for the high variety of firing patterns in experimental spike trains. Formally, we compare the numbers N1=N1(t,h) and N2=N2(t,h) of spikes in the time intervals (t-h,t] and (h,t+h]. By varying t within a fine time lattice and simultaneously varying the interval length h, we obtain a multivariate statistic D(h,t):=(N1-N2)/√(N1+N2), for which we prove asymptotic multivariate normality under homogeneity. From this a practical, graphical device to spot changes of the firing rate is constructed.
Our graphical representation of D(h,t) (Figure 1A) visualizes the changes in the firing rate. For the statistical test, a threshold K is chosen such that under homogeneity, |D(h,t)|<K holds for all investigated h and t with probability 0.95. This threshold can indicate potential change points in order to estimate the inhomogeneous rate profile (Figure 1B). The SFT is applied to a sample data set of spontaneous single unit activity recorded from the substantia nigra of anesthetized mice. In this data set, multiple rate changes are identified which agree closely with visual inspection. In contrast to approaches choosing one fixed kernel width , our method has advantages in the flexibility of h.
Brown EN, Kass RE, Mitra P: Multiple neural spike train analysis: state-of-the-art and future challenges. Nat Neurosci. 2004, 7 (5): 456-61. 10.1038/nn1228.
Brody CD: Correlations without synchrony. Neural Comp. 1999, 11: 1537-51. 10.1162/089976699300016133.
Gourévitch B, Eggermont JJ: A simple indicator of nonstationarity of firing rate in spike trains. J Neurosci Methods. 2007, 163 (1): 181-7. 10.1016/j.jneumeth.2007.02.021.
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-2.
We thank Brooks Ferebee for stimulating discussions. This work was supported by the LOEWE-Schwerpunkt “Neuronale Koordination Forschungsschwerpunkt Frankfurt” (MM, JR), by the BMBF Project Bernstein Fokus: Neurotechnologie Frankfurt, FKZ 01GQ0841 (MB) and by the Deutsche Forschungsgemeinschaft SFB 815 (JR).
About this article
Cite this article
Tillmann, M., Messer, M., Bingmer, M. et al. Detection and localization of multiple rate changes in Poisson spike trains. BMC Neurosci 12 (Suppl 1), P268 (2011). https://doi.org/10.1186/1471-2202-12-S1-P268