Abstract
This paper develops a nonparametric method for estimating the spectral density of multivariate stationary time series using basis expansion. A likelihood-based approach is used to fit the model through the minimization of a penalized Whittle negative log-likelihood. Then, a Newton-type algorithm is developed for the computation. In this method, we smooth the Cholesky factors of the multivariate spectral density matrix in a way that the reconstructed estimate based on the smoothed Cholesky components is consistent and positive-definite. In a simulation study, we have illustrated and compared our proposed method with other competitive approaches. Finally, we apply our approach to two real-world problems, Electroencephalogram signals analysis, \(El\ Ni\tilde{n}o\) Cycle.
Similar content being viewed by others
References
Akaike H (1974) A new look at the statistical model identification. IEEE Trans Autom Control 19(6):716–723
Berkowitz J, Diebold FX (1998) Bootstrapping multivariate spectra. Rev Econ Stat 80(4):664–666
Brillinger DR (2001) Time series: data analysis and theory. SIAM, Philadelphia
Brockwell PJ, Davis RA (1991) Stationary time series, time series: theory and methods. Springer, Berlin, pp 1–41
Calder M, Davis RA (1997) Introduction to whittle (1953) the analysis of multiple stationary time series, breakthroughs in statistics. Springer, Berlin, pp 141–169
Capon J (1979) Maximum-likelihood spectral estimation. Nonlinear Methods Spectral Anal 155–179
Dai M, Guo W (2004) Multivariate spectral analysis using Cholesky decomposition. Biometrika 91(3):629–643
Davies RB (1973) Asymptotic inference in stationary Gaussian time-series. Adv Appl Probab 5(3):469–497
Dunsmuir W (1979) A central limit theorem for parameter estimation in stationary vector time series and its application to models for a signal observed with noise. Ann Stat 490–506
Eilers PH, Marx BD et al (1996) Flexible smoothing with b-splines and penalties. Stat Sci 11(2):89–121
Fan J, Kreutzberger E (1998) Automatic local smoothing for spectral density estimation. Scand J Stat 25(2):359–369
Güler I, Kiymik MK, Akin M, Alkan A (2001) AR spectral analysis of EEG signals by using maximum likelihood estimation. Comput Biol Med 31(6):441–450
Hansen M, Huang J, Kooperberg C, Stone C, Truong Y (2001) Statistical modeling with spline functions methodology and theory
Krafty RT, Collinge WO (2013) Penalized multivariate whittle likelihood for power spectrum estimation. Biometrika 100(2):447–458
Maadooliat M, Huang JZ, Hu J (2015) Integrating data transformation in principal components analysis. J Comput Graph Stat 24(1):84–103
Maadooliat M, Sun Y, Chen T (2018) Nonparametric collective spectral density estimation with an application to clustering the brain signals. Stat Med 37(30):4789–4806
Maturana-Russel P, Meyer R (2021) Bayesian spectral density estimation using p-splines with quantile-based knot placement. Comput Stat 36:2055–2077
Müller B, Reinhardt J, Strickland MT (1995) Neural networks: an introduction. Springer, Berlin
Pawitan Y (1996) Automatic estimation of the cross-spectrum of a bivariate time series. Biometrika 83(2):419–432
Pawitan Y, O’sullivan F (1994) Nonparametric spectral density estimation using penalized whittle likelihood. J Am Stat Assoc 89(426):600–610
Priestley MB (1981) Spectral analysis and time series: probability and mathematical statistics. Number 04; QA280, P7
Quiroga RQ, Kraskov A, Kreuz T, Grassberger P (2002) Performance of different synchronization measures in real data: a case study on electroencephalographic signals. Phys Rev E 65(4):041903
Rosen O, Stoffer DS (2007) Automatic estimation of multivariate spectra via smoothing splines. Biometrika 94(2):335–345
Shumway R, Stoffer D (2011) ARIMA models’. Time series analysis and its applications. Springer, New York
Shumway RH, Stoffer DS, Stoffer DS (2000) Time series analysis and its applications, vol 3. Springer, New York
van Luijtelaar G (1997) The WAG/Rij rat model of absence epilepsy: ten years of research: a compilation of papers. Nijmegen University Press
Whittle P (1953) Estimation and information in stationary time series. Arkiv för matematik 2(5):423–434
Whittle P (1961) Gaussian estimation in stationary time seris. Bull. Int. Stat. Inst. 39:105–129
Zhang S (2022) Automatic estimation of spatial spectra via smoothing splines. Comput Stat 37(2):565–590
Acknowledgements
The authors express sincere gratitude to the anonymous referees for their valuable and insightful comments, which significantly contributed to the enhancement of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Nezampour, S., Nematollahi, A., Krafty, R.T. et al. A new approach to nonparametric estimation of multivariate spectral density function using basis expansion. Comput Stat 39, 3625–3641 (2024). https://doi.org/10.1007/s00180-023-01451-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00180-023-01451-4