A multiple filter test for the detection of rate changes in renewal processes with varying variance. (English) Zbl 1454.62365
Summary: Nonstationarity of the event rate is a persistent problem in modeling time series of events, such as neuronal spike trains. Motivated by a variety of patterns in neurophysiological spike train recordings, we define a general class of renewal processes. This class is used to test the null hypothesis of stationary rate versus a wide alternative of renewal processes with finitely many rate changes (change points). Our test extends ideas from the filtered derivative approach by using multiple moving windows simultaneously. To adjust the rejection threshold of the test, we use a Gaussian process, which emerges as the limit of the filtered derivative process. We also develop a multiple filter algorithm, which can be used when the null hypothesis is rejected in order to estimate the number and location of change points. We analyze the benefits of multiple filtering and its increased detection probability as compared to a single window approach. Application to spike trains recorded from dopamine midbrain neurons in anesthetized mice illustrates the relevance of the proposed techniques as preprocessing steps for methods that assume rate stationarity. In over 70% of all analyzed spike trains classified as rate nonstationary, different change points were detected by different window sizes.
MSC:
| 62P10 | Applications of statistics to biology and medical sciences; meta analysis |
| 62M10 | Time series, auto-correlation, regression, etc. in statistics (GARCH) |
Keywords:
stochastic processes; renewal processes; change point detection; nonstationary rate; multiple filters; multiple time scalesReferences:
| [1] | Antoch, J. and Hušková, M. (1994). Procedures for the detection of multiple changes in series of independent observations. In Asymptotic Statistics ( Prague , 1993). Contrib. Statist. 3-20. Physica, Heidelberg. · doi:10.1007/978-3-642-57984-4_1 |
| [2] | Basseville, M. and Nikiforov, I. V. (1993). Detection of Abrupt Changes : Theory and Application . Prentice Hall, Englewood Cliffs, NJ. |
| [3] | Benda, J. and Herz, A. V. M. (2003). A universal model for spike-frequency adaptation. Neural Comput. 15 2523-2564. · Zbl 1052.92010 · doi:10.1162/089976603322385063 |
| [4] | Bertrand, P. (2000). A local method for estimating change points: The “hat-function.” Statistics 34 215-235. · Zbl 0955.62031 · doi:10.1080/02331880008802714 |
| [5] | Bertrand, P. R., Fhima, M. and Guillin, A. (2011). Off-line detection of multiple change points by the filtered derivative with \(p\)-value method. Sequential Anal. 30 172-207. · Zbl 1215.62076 · doi:10.1080/07474946.2011.563710 |
| [6] | Billingsley, P. (1999). Convergence of Probability Measures , 2nd ed. Wiley, New York. · Zbl 0944.60003 |
| [7] | Brodsky, B. E. and Darkhovsky, B. S. (1993). Nonparametric Methods in Change-Point Problems . Kluwer Academic, Dordrecht. · Zbl 1274.62512 |
| [8] | Brody, C. D. (1999). Correlations without synchrony. Neural Comput. 11 1537-1551. |
| [9] | Brown, E. N., Kass, R. E. and Mitra, P. P. (2004). Multiple neural spike train data analysis: State-of-the-art and future challenges. Nat. Neurosci. 7 456-461. |
| [10] | Csörgő, M. and Horváth, L. (1987). Asymptotic distributions of pontograms. Math. Proc. Cambridge Philos. Soc. 101 131-139. · Zbl 0645.60040 · doi:10.1017/S0305004100066470 |
| [11] | Csörgő, M. and Horváth, L. (1997). Limit Theorems in Change-Point Analysis . Wiley, Chichester. · Zbl 0884.62023 |
| [12] | Dreyer, J. K., Herrik, K. F., Berg, R. W. and Hounsgaard, J. D. (2010). Influence of phasic and tonic dopamine release on receptor activation. J. Neurosci. 30 14273-14283. |
| [13] | Farkhooi, F., Strube-Bloss, M. and Nawrot, M. P. (2009). Serial correlation in neural spike trains: Experimental evidence, stochastic modelling, and single neuron variability. Phys. Rev. E 79 021905. |
| [14] | Gonon, F. G. (1988). Nonlinear relationship between impulse flow and dopamine released by rat midbrain dopaminergic neurons as studied by in vivo electrochemistry. Neuroscience 24 19-28. |
| [15] | Gourévitch, B. and Eggermont, J. J. (2007). A nonparametric approach for detection of bursts in spike trains. J. Neurosci. Methods 160 349-358. |
| [16] | Grün, S., Diesmann, M. and Aertsen, A. (2002). “Unitary Events” in multiple single-neuron activity. II. Non-Stationary data. Neural Comput. 14 81-119. · Zbl 1005.92006 · doi:10.1162/089976602753284464 |
| [17] | Grün, S., Riehle, A. and Diesmann, M. (2003). Effect of cross-trial nonstationarity on joint-spike events. Biol. Cybernet. 88 335-351. · Zbl 1075.92013 · doi:10.1007/s00422-002-0386-2 |
| [18] | Grün, S. and Rotter, S., eds. (2010). Analysis of Parallel Spike Trains. Springer Series in Computational Neuroscience 7 . Springer, New York. |
| [19] | Howe, M. W., Tierney, P. L., Sandberg, S. G., Phillips, P. E. M. and Graybiel, A. M. (2013). Prolonged dopamine signalling in striatum signals proximity and value of distant rewards. Nature 500 575-579. |
| [20] | Hušková, M. and Slabý, A. (2001). Permutation tests for multiple changes. Kybernetika ( Prague ) 37 605-622. · Zbl 1264.62038 |
| [21] | Johnson, D. H. (1996). Point process models of single-neuron discharges. J. Comput. Neurosci. 3 275-299. |
| [22] | Jones, M. C., Marron, J. S. and Sheather, S. J. (1996). A brief survey of bandwidth selection for density estimation. J. Amer. Statist. Assoc. 91 401-407. · Zbl 0873.62040 · doi:10.2307/2291420 |
| [23] | Kass, R. E., Ventura, V. and Brown, E. N. (2005). Statistical issues in the analysis of neuronal data. J. Neurophysiol. 94 8-25. |
| [24] | Kendall, D. G. and Kendall, W. S. (1980). Alignments in two-dimensional random sets of points. Adv. in Appl. Probab. 12 380-424. · Zbl 0425.60009 · doi:10.2307/1426603 |
| [25] | Legéndy, C. R. and Salcman, M. (1985). Bursts and recurrences of bursts in the spike trains of spontaneously active striate cortex neurons. J. Neurophysiol. 53 926-939. |
| [26] | Messer, M., Kirchner, M., Schiemann, J., Roeper, J., Neininger, R. and Schneider, G. (2014). Supplement to “A multiple filter test for the detection of rate changes in renewal processes with varying variance.” . · Zbl 1454.62365 |
| [27] | Mileykovskiy, B. and Morales, M. (2011). Duration of inhibition of ventral tegmental area dopamine neurons encodes a level of conditioned fear. J. Neurosci. 31 7471-7476. |
| [28] | Muhsal, B. (2013). Change-point methods for multivariate autoregressive models and multiple structural breaks in the mean. Dissertation, available at . |
| [29] | Nawrot, M. P., Aertsen, A. and Rotter, S. (1999). Single-trial estimation of neuronal firing rates: From single-neuron spike trains to population activity. J. Neurosci. Methods 94 81-92. |
| [30] | Nawrot, M. P., Boucsein, C., Rodriguez Molina, V., Riehla, A., Aertsen, A. and Rotter, S. (2008). Measurement of variability dynamics in cortical spike trains. Journal of Neuroscience Methods 169 347-390. |
| [31] | Perkel, D. H., Gerstein, G. L. and Moore, G. P. (1967a). Neuronal spike trains and stochastic point processes. I. The single spike train. Biophys. J. 7 391-417. |
| [32] | Perkel, D. H., Gerstein, G. L. and Moore, G. P. (1967b). Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains. Biophys. J. 7 419-440. |
| [33] | Petrov, V. V. (1995). Limit Theorems of Probability Theory : Sequences of Independent Random Variables . Clarendon, Oxford. · Zbl 0826.60001 |
| [34] | Pollard, D. (1984). Convergence of Stochastic Processes . Springer, New York. · Zbl 0544.60045 |
| [35] | Redgrave, P., Rodriguez, M., Smith, Y., Rodriguez-Oroz, M. C., Lehericy, S., Bergman, H., Agid, Y., DeLong, M. R. and Obeso, J. A. (2010). Goal-directed and habitual control in the basal ganglia: Implications for Parkinson’s disease. Nat. Rev. , Neurosci. 11 760-772. |
| [36] | Rieke, F., Warland, D., de Ruyter van Steveninck, R. and Bialek, W. (1999). Spikes : Exploring the Neural Code . MIT Press, Cambridge, MA. · Zbl 0912.92004 |
| [37] | Schiemann, J., Klose, V., Schlaudraff, F., Bingmer, M., Seino, S., Magill, P. J., Schneider, G., Liss, B. and Roeper, J. (2012). K-ATP channels control in vivo burst firing of dopamine neurons in the medial substantia nigra and novelty-induced behavior. Nat. Neurosci. 15 1272-1280. |
| [38] | Schneider, G. (2008). Messages of oscillatory correlograms: A spike train model. Neural Comput. 20 1211-1238. · Zbl 1138.92308 · doi:10.1162/neco.2007.12-06-424 |
| [39] | Shimazaki, H. and Shinomoto, S. (2007). A method for selecting the bin size of a time histogram. Neural Comput. 19 1503-1527. · Zbl 1115.92014 · doi:10.1162/neco.2007.19.6.1503 |
| [40] | Staude, B., Rotter, S. and Grün, S. (2008). Can spike coordination be differentiated from rate covariation? Neural Comput. 20 1973-1999. · Zbl 1140.92008 · doi:10.1162/neco.2008.06-07-550 |
| [41] | Steinebach, J. and Eastwood, V. R. (1995). On extreme value asymptotics for increments of renewal processes. J. Statist. Plann. Inference 45 301-312. · Zbl 0831.60060 · doi:10.1016/0378-3758(94)00079-4 |
| [42] | Steinebach, J. and Zhang, H. Q. (1993). On a weighted embedding for pontograms. Stochastic Process. Appl. 47 183-195. · Zbl 0779.60031 · doi:10.1016/0304-4149(93)90013-T |
| [43] | Tokdar, S., Xi, P., Kelly, R. C. and Kass, R. E. (2010). Detection of bursts in extracellular spike trains using hidden semi-Markov point process models. J. Comput. Neurosci. 29 203-212. |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
