VARMAX (Vector Autoregressive Moving Average model with eXogenous variables) extends the ARMA/ARIMA model in two ways:
- To work with time series with multiple response variables (vector time series).
- To work with exogenous variables, or variables that are independent of the other variables in the system.
The model includes both the dynamic relationship between the multiple response variables and the relationship between the dependent and independent variables.
Nonseasonal VARMAX Model
This formula represents a nonseasonal VARMAX model:
| Variable | Description |
|---|---|
| Yt | Time series made stationary. First term is autoregressive component, second term is exogenous component, third term is moving average component, the fourth term (C) is vector of constants, and fifth term (Et) is vector of residual errors. |
| Yt | Vector of n response variables. |
| Xt | Vector of m exogenous variables. |
| p | Number of previous periods of endogenous variables included in model. |
| q | Number of previous periods included in moving average. |
| b | Number of previous periods of exogenous variables included. |
| Φi | n * n matrix of autoregressive parameters. |
| Bi | n * m matrix of exogenous variable parameters. |
| Θi | n * n matrix of moving average parameters. |
| Et | Difference between actual and predicted value of Yt, (Yt - Ŷt). |
Seasonal VARMAX Model
This formula represents a seasonal VARMAX model:
(1 - Φ1Back - … - ΦpBackp)(1 - Φ1Backm - … - ΦspBackm*sp)(1 - Back)d(1 - Backm)sd(Yt) =
C + (1 + Θ1Back + … + ΘqBackq) (1 + Θ1Backm + … + ΘsqBackm*sq) Et +
(B1 + B2Back + … + Bb-1Backb) (1 - Back)d (1 - Backm)sd Xt
Variables are as in the nonseasonal VARMAX Model and the following table.
| Variable | Description |
|---|---|
| Back | Back-shift operator, that is: Back(yt) = yt-1 |
| m | Number of periods in each season. |
| d | Number of differencing steps performed to stationarize time series. |
| sp | Seasonal parameter corresponding to p. |
| sq | Seasonal parameter corresponding to q. |
| sd | Seasonal parameter corresponding to d. |
Difference between VARMAX and R Function VARMA
The results of the VARMAX function and the VARMA function in the R package MTS might not match exactly, because VARMA applies a boundary to the numerical optimization and VARMAX does not. Also, the default optimization method (the Method syntax element) is different for VARMA and VARMAX.