Semiparametric Bayesian Inference for Time Series with Mixed Spectra

A Bayesian analysis is presented of a time series which is the sum of a stationary component with a smooth spectral density and a deterministic component consisting of a linear combination of a trend and periodic terms. The periodic terms may have known or unknown frequencies. The advantage of our approach is that different features of the data such as the regression parameters, the spectral density, unknown frequencies, and missing observations are combined in a hierarchical Bayesian framework and estimated simultaneously. A Bayesian test to detect the presence of deterministic components in the data is also constructed. By using an asymptotic approximation to the likelihood, the computation is carried out efficiently using Markov chain Monte Carlo in O(Mn) operations, where n is the sample size and M and is the number of iterations. We show empirically that our approach works well on real and simulated examples