Skip to Main content Skip to Navigation
Journal articles

Computational aspects of Bayesian spectral density estimation

Abstract : Gaussian time-series models are often specified through their spectral density. Such models present several computational challenges, in particular because of the nonsparse nature of the covariance matrix. We derive a fast approximation of the likelihood for such models. We propose to sample from the approximate posterior (i.e., the prior times the approximate likelihood), and then to recover the exact posterior through importance sampling. We show that the variance of the importance sampling weights vanishes as the sample size goes to infinity. We explain why the approximate posterior may typically be multimodal, and we derive a Sequential Monte Carlo sampler based on an annealing sequence to sample from that target distribution. Performance of the overall approach is evaluated on simulated and real datasets. In addition, for one real-world dataset, we provide some numerical evidence that a Bayesian approach to semiparametric estimation of spectral density may provide more reasonable results than its frequentist counterparts. The article comes with supplementary materials, available online, that contain an Appendix with a proof of our main Theorem, a Python package that implements the proposed procedure, and the Ethernet dataset.
Complete list of metadatas

https://hal-hec.archives-ouvertes.fr/hal-01026131
Contributor : Antoine Haldemann <>
Submitted on : Saturday, July 19, 2014 - 10:31:36 PM
Last modification on : Thursday, March 5, 2020 - 5:50:51 PM

Links full text

Identifiers

Collections

Citation

Nicolas Chopin, J. Rousseau, Brunero Liseo. Computational aspects of Bayesian spectral density estimation. Journal of Computational and Graphical Statistics, Taylor & Francis, 2013, 22 (3), pp.533-557. ⟨10.1080/10618600.2013.785293⟩. ⟨hal-01026131⟩

Share

Metrics

Record views

357