The mbsts package: Multivariate Bayesian Structural Time Series Models in R

Structural time series models are state space models for time series data. They are constructed in terms of components each of which has a direct interpretation. For example, one may consider a decomposition in which a series can be seen as the sum of trend and regression components. The multivariate Bayesian structural time series (MBSTS) model (Qiu et al., 2018) is a generalized version of many structural time series models and is constructed as the sum of a trend component, a seasonal component, a cycle component, a regression component, and an error term, where each component provides an independent and additional effect. Users have flexibility in choosing these components and are free to construct their specific forms, for example adding on a regression component with predictors generated through data mining as that in (Jammalamadaka et al., 2019). The MBSTS model uses the Bayes selection technique via Markov chain Monte Carlo (MCMC) methods to select among a set of contemporary predictors, thus one does not need to commit to a fixed set of predictors. Specifically, the variable selection technique uses a "spike and slab" approach, through which a different set of predictors can be selected in each MCMC iteration. Then important predictors will be selected according to their overall frequency of numbers being selected over the total number of MCMC iterations. The multivariate structure and the Bayesian framework allow the model to take advantage of the association structure among target series.