L33 Some Specific Multivariate Time Series Models
Linear Filtering¶
- \(X_{t}\) is a "r" dimensional input series
- \(Y_{t}\) is a "k" dimensional output series
Multivariate linear (Time invariant) filter relating \(X_{t}\) and \(Y_{t}\)
where, \(\Psi_{j}\) are \(k\times r\) matrices for all \(j\).
If \(\Psi_{j}=0\) for all \(j \lt 0\)
$$
Y_{t} = \sum_{j=0}^\infty \Psi_{j} X_{t-j}
$$
\(Y_{t}\) is expressible only in terms of present and past values of the input series \(X_{t}\).
- Let \(\lVert A \rVert^{2}\) denote the "norm" of \(A\) = \(tr(A'A)\)
- Linear filter is "stable" \(\impliedby\) \(\sum_{j=0}^\infty \lVert \Psi_{j} \rVert \lt \infty\)
Theorem
If
- Linear filter is stable
- Input random vectors \(\{ X_{t} \}\) have finite second moments
Then
- Filtering representations exists uniquely and converges to mean square
If the linear filter is stable and \(X_{t}\) is stationary with cross covariances \(\Gamma_{x}(l)\), then \(Y_{t}\) is stationary with cross covariances,
Wold (MA) representation of the Series¶
- "infinite MA representation of a stationary vector process"
Let \(\{ Y_{t} \}\) be a multivariate stationary process with mean \(\mu\)
$$
Y_{t} - \mu = \sum_{j=0}^\infty \Psi_{j}e_{t-j}
$$
$$
Y_{t} = \mu + \sum_{j=0}^\infty \Psi(B)e_{t}
$$
Where \(\Psi(B)\) is a \(k \times k\) matrix of backward shift operator. (= \(\sum_{j=0}^\infty \Psi_{j}B^j\))
Vector Autoregressive Moving Average (VARMA)¶
\(VARMA(p,q)\) process
where,
- \(\Phi_{p}(B)= \Phi_{0} - \Phi_{1}(B) - \dots \Phi_{p}B^p\)
- and \(\Phi_{i}\) here are matrices and not single quantities (parameters)
- \(\Theta_{q}(B)= \Theta_{0} - \Theta_{1}(B) - \dots \Theta_{q}B^q\)
- and \(\Theta_{i}\) here are matrices and not single quantities (parameters)
- \(q=0\) \(\implies\) \(VAR(p)\)
- \(p=0\) \(\implies\) \(VMA(q)\)
A process is stationary if we can represent this as a convergent vector moving average process of infinite order.
$$
Y_{t} = \mu + e_{t} + \sum_{j=1}^\infty \Psi_{j}e_{t-j} = \mu + \Phi(B)e_{t}
$$
By taking \(Y_{t}-\mu\) on one side we get
And so,
$$
\Phi(B) = 1 + \sum_{j=1}^\infty \Phi_{j}B^j = \Phi_{p}(B)^{-1}\Theta_{q}(B)
$$
- This is not VARMA because it doesn't have any \(Y_{t-i}\) terms on the RHS. It is just VMA.
- If we can write this as VMA, it is stationary.
- where \(\Phi^*(z) =\text{adj}(\Phi(z))\)
Thus the process is
The process is stationary if \(\{ \det(\Phi(z))^{-1} \}\) is convergent for \(\lvert z \rvert \lt 1\).
VMA (order \(q\))¶
where
Vector MA(1)¶
Assume, \(\mu = 0\)
By recursive substitution,
A Vector Infinite autoregressive process can be rewritten as VMA(1)