Take the general forecasting problem to be:

Yt+h = E[Yt-1 | Xt-1]

where Y is the predicted outcome of interest, E is the expectation operator, X is the information set, and t is the time index. Then, let the process be state-dependent:

Yt+h = E[Yt-1 | Xt-1, st-1]

where s is a discrete, that is, mutually exclusive, state of the world. For example, in an economic context, one may consider expansions and recessions as two different states of the world, and in the US, these labels are formally declared by the NBER.

The sovereign package provides the tools and methods to solves this state-dependent forecasting problem. First, using either parametric or non-parametric machine learning techniques, users may sort, identify, and classify discrete states of the world. Second, using vector auto-regressions or direct projections, a user can create both single- and multi-state time series forecasting models. Third, users may analyze their models and forecasts through impulse response functions, forecast error variance decompositions, and historical error decompositions.

Available tools and techniques may be reviewed under the Tools tab. While vignettes and other extended documentation of the sovereign package’s capabilities may be found under the Workflow tab.


One may install sovereign through the package’s GitHub page with


or through CRAN via



If you should have questions, concerns, or wish to collaborate, please contact Tyler J. Pike