Computes residuals of a Vector ARMA model.
varmares(w, phi, theta)¶
w (NxK matrix) – time series.
phi ((K*P)xK matrix) – AR coefficient matrices.
theta ((K*Q)xK matrix) – MA coefficient matrices.
res (NxK matrix) –
residuals. If the calculation fails res is set to missing value with error code:
Reason for Failure
\(M < 1\)
\(N < 1\)
\(P < 0\)
\(Q < 0\)
\(P = 0\) and \(Q = 0\)
floating point work space too small
integer work space too small
AR parameters too close to stationarity boundary
model not stationary
model not invertible
\(I+M'H'HM\) not positive definite
varmares() is adapted from code developed by Jose Alberto Mauricio of the
Universidad Complutense de Madrid. It was published as Algorithm AS311
in Applied Statistics. Also described in “Exact Maximum Likelihood
Estimation of Stationary Vector ARMA Models,” JASA, 90:282-264.