The above forecast formulas and functions for univariate time series can be generalized directly to the multivariate case. Here the forecast
is a vector, and the mean square errors are now error covariance matrices and (
7.8) becomes
. The function
newforecast we wrote earlier for updating the forecast should be modified to incorporate the matrix nature of the
and
weights.