Premium
ESTIMATION OF MULTIVARIATE TIME SERIES
Author(s) -
Shea B. L.
Publication year - 1987
Publication title -
journal of time series analysis
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.576
H-Index - 54
eISSN - 1467-9892
pISSN - 0143-9782
DOI - 10.1111/j.1467-9892.1987.tb00423.x
Subject(s) - mathematics , univariate , autoregressive–moving average model , kalman filter , series (stratigraphy) , recursion (computer science) , multivariate statistics , covariance matrix , algorithm , covariance , autoregressive model , estimation theory , moving average , statistics , paleontology , biology
. The algorithm proposed here is a multivariate generalization of a procedure discussed by Pearlman (1980) for calculating the exact likelihood of a univariate ARMA model. Ansley and Kohn (1983) have shown how the Kalman filter can be used to calculate the exact likelihood function when not all the observations are known. In Shea (1983) it is shown that this algorithm is much quicker than that of Ansley and Kohn (1983) for all ARMA models except an ARMA (2, 1) and a couple of low‐order AR processes and therefore when we have no missing observations this algorithm should be used instead. The Fortran subroutine G13DCF in the NAG (1987) Library fits a vector ARMA model using an adaptation of this algorithm. Experience in the use of this routine suggests that having reasonably good initial estimates of the ARMA parameter matrices, and in particular the residual error covariance matrix, can not only substantially reduce the computing time but more important improve the convergence properties of the minimization procedure. We therefore propose a method of calculating initial estimates of the ARMA parameters which involves using a generalization of the concept of inverse cross covariances from the univariate to the multivariate case. Finally theory is put into practice with the fitting of a bivariate model to a couple of real‐life time series.