×

Time series decomposition into oscillation components and phase estimation. (English) Zbl 1418.62339

Summary: Many time series are naturally considered as a superposition of several oscillation components. For example, electroencephalogram (EEG) time series include oscillation components such as alpha, beta, and gamma. We propose a method for decomposing time series into such oscillation components using state-space models. Based on the concept of random frequency modulation, Gaussian linear state-space models for oscillation components are developed. In this model, the frequency of an oscillator fluctuates by noise. Time series decomposition is accomplished by this model like the Bayesian seasonal adjustment method. Since the model parameters are estimated from data by the empirical Bayes’ method, the amplitudes and the frequencies of oscillation components are determined in a data-driven manner. Also, the appropriate number of oscillation components is determined with the Akaike information criterion (AIC). In this way, the proposed method provides a natural decomposition of the given time series into oscillation components. In neuroscience, the phase of neural time series plays an important role in neural information processing. The proposed method can be used to estimate the phase of each oscillation component and has several advantages over a conventional method based on the Hilbert transform. Thus, the proposed method enables an investigation of the phase dynamics of time series. Numerical results show that the proposed method succeeds in extracting intermittent oscillations like ripples and detecting the phase reset phenomena. We apply the proposed method to real data from various fields such as astronomy, ecology, tidology, and neuroscience.

MSC:

62M10 Time series, auto-correlation, regression, etc. in statistics (GARCH)
62P10 Applications of statistics to biology and medical sciences; meta analysis

Software:

astsa; CircStat
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Akaike, H. (1980). Likelihood and the Bayes procedure. In J. M. Bernardo, M. H. DeGroot, D. V. Lindley, & A. F. M. Smith (Eds.), Bayesian statistics (pp. 1-13). Durhom, NC: International Society of Bayesian Statistics. · Zbl 0471.62033
[2] Akaike, H., & Ishiguro, M. (1980). BAYSEA: A Bayesian seasonal adjustment program. Tokyo: Institute of Statistical Mathematics.
[3] Ba, D., Babadi, B., Purdon, P. L., & Brown, E. N. (2014). Robust spectrotemporal decomposition by iteratively reweighted least squares. Proceedings of the National Academy of Sciences, 111, E5336-E5345. ,
[4] Berens, P. (2009). CircStat: A Matlab toolbox for circular statistics. Journal of Statistical Software, 31, 1-21. ,
[5] Brockwell, P. J., & Davis, R. A. (2009). Time series: Theory and methods. New York: Springer. · Zbl 1169.62074
[6] Bruns, A. (2004). Fourier-, Hilbert- and wavelet-based signal analysis: Are they really different approaches? Journal of Neuroscience Methods, 137, 321-332. ,
[7] Buzsáki, G. (2011). Rhythms of the brain. Oxford: Oxford University Press. · Zbl 1204.92017
[8] Campbell, M. J., & Walker, A. M. (1977). A survey of statistical work on the Mackenzie River Series of annual Canadian lynx trappings for the years 1821-1934 and a new analysis. Journal of the Royal Statistical Society A, 140, 411-431. ,
[9] Cohen, M. X. (2014). Analyzing neural time series data. Cambridge MA: MIT Press.
[10] Hanslmayr, S., Klimesch, W., Sauseng, P., Gruber, W., Doppelmayr, M., Freunberger, … Birbaumer, N. (2007). Alpha phase reset contributes to the generation of ERPs. Cerebral Cortex, 17, 1-8. ,
[11] Harvey, M. X. (1989). Forecasting, structural time series models and the Kalman filter. Cambridge: Cambridge University Press.
[12] Harvey, M. X., & Trimbur, T. M. (2003). General model-based filters for extracting cycles and trends in economic time Series. Review of Economics and Statistics, 85, 244-255. ,
[13] Huerta, G., & West, M. (1999a). Priors and component structures in autoregressive time series models. Journal of the Royal Statistical Society B, 61, 881-899. , · Zbl 0940.62079
[14] Huerta, G., & West, M. (1999b). Bayesian inference on periodicities and component spectral structure in time series. Journal of Time Series Analysis, 20, 401-416. , · Zbl 0958.62090
[15] Kitagawa, G. (2010). Introduction to time series modeling. London: Chapman & Hall. , · Zbl 1204.62159
[16] Kitagawa, G., & Gersch, W. (1984). A smoothness priors-state space modeling of time series with trend and seasonality. Journal of the American Statistical Association, 79, 378-389.
[17] Le Van Quyen, M., Foucher, J., Lachaux, J. P., Rodriguez, E., Lutz, A., Martinerie, J., & Varela, F. J. (2001). Comparison of Hilbert transform and wavelet methods for the analysis of neuronal synchrony. Journal of Neuroscience Methods, 111, 83-98. ,
[18] Lopour, B. A., Tavassoli, A., Fried, I., & Ringach, D. L. (2013). Coding of information in the phase of local field potentials within human medial temporal lobe. Neuron, 79, 594-606. ,
[19] Makeig, S., Westerfield, M., Jung, T. P., Enghoff, S., Townsend, J., Courchesne, E., & Sejnowski, T. J. (2002). Dynamic brain sources of visual evoked responses. Science, 295, 690-694. ,
[20] Mizuseki, K., Sirota, A., Pastalkova, E., & Buzsaki, G. (2009a). Multi-unit recordings from the rat hippocampus made during open field foraging. http://dx.doi.org/10.6080/K0Z60KZ9
[21] Mizuseki, K., Sirota, A., Pastalkova, E., & Buzsaki, G. (2009b). Theta oscillations provide temporal windows for local circuit computation in the entorhinal-hippocampal loop. Neuron, 64, 267-280. ,
[22] Nguyen, D. P., Wilsona, M. A., Brown, E. N., & Barbieri, R. (2009). Measuring instantaneous frequency of local field potential oscillations using the Kalman smoother. Journal of Neuroscience Methods, 184, 365-374. ,
[23] O’Keefe, J., & Recce, M. L. (1993). Phase relationship between hippocampal place units and the EEG theta rhythm. Hippocampus, 3, 317-330. ,
[24] Quinn, B. G., & Hannan, E. J. (2013). The estimation and tracking of frequency. Cambridge: Cambridge University Press. · Zbl 1256.62054
[25] Shumway, R. H., & Stoffer, D. S. (1982). An approach to time series smoothing and forecasting using the EM algorithm. Journal of Time Series Analysis, 3, 253-264. , · Zbl 0502.62085
[26] Siapas, A. G., Lubenov, E. V., & Wilson, M. A. (2005). Prefrontal phase locking to hippocampal theta oscillations. Neuron, 46, 141-151. ,
[27] Tamura, Y., Sato, T., Ooe, M., & Ishiguro, M. (1991). A procedure for tidal analysis with a Bayesian information criterion. Geophysical Journal International, 104, 507-516. ,
[28] Tong, H. (1977). Some comments on the Canadian lynx data. Journal of the Royal Statistical Society A, 140, 432-436. ,
[29] Tong, H. (1990). Nonlinear time series: A dynamical system approach. Oxford: Oxford University Press. · Zbl 0716.62085
[30] West, M. (1997). Time series decomposition. Biometrika, 84, 489-494. , · Zbl 0882.62088
[31] West, M., Prado, R., & Krystal, A. D. (1999). Evaluation and comparison of EEG traces: Latent structure in nonstationary time series. Journal of the American Statistical Association, 94, 1083-1095. ,
[32] Whittle, P. (1953). Estimation and information in stationary time series. Arkiv for Matematik, 2, 423-434. , · Zbl 0053.41003
[33] Wiener, N. (1966). Nonlinear problems in random theory. Cambridge, MA: MIT Press.
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.