多元时间序列分析中一个重要概念是引导与滞后关系 为此， 用互相关阵来衡量时间序列之间的线性关系的强度 \(k\)元弱平稳列\(\boldsymbol r_t\)嘚滞后\(l\)的互协方差阵定义为 \[\begin{aligned} \Gamma_l =

时间序列之间的线性相依性的分类

两个分量\(r_{it}\)和\(r_{js}\)的线性相依关系可以分类为如下情况：

则\(r_{it}\)和\(r_{js}\)之间存在相互的反馈關系， 互为引导和滞后

多元时间序列的一些计算使用蔡瑞胸(R.S. Tsay)教授的MTS扩展包。

考虑IBM股票与标普500指数从1926年1月到2008年12月的月对数收益率 共996个观測值。单位为百分之一

IBM的月对数收益率时间序列图：

标普500的月对数收益率时间序列图：

两个序列同时刻的散点圖：

同时刻具有明显的正相关，样本相关系数：

同时刻的相关系数为0.65显著。

IBM对滞后一步的标普500：

标普500对滞后一步的IBM：

用MTS包的MTSplot()函数绘制时间序列图：

计算两個序列的单样本简单样本统计量：

从超额峰度(Kurtosis)来看两个序列有重尾分布 而标普的重尾分布更明显。

用MTS包的ccm()函数显示简约的互相关阵序列囷图形 R.S Tsay教授的MTS包的ccm()函数的源代码见后面， 略作修改

随后给出了\(l=1,2,\dots\)的CCM(即\(\hat{\boldsymbol\rho}_l\))的简化表示。 两个序列的同步的相关性较强 相关系数为0.645； IBM对数收益率的序列自相关（矩阵的左上角元素）都不显著； 标普对数收益率的序列自相关（矩阵的右下角元素）在滞后1,3,5,9等位置显著； IBM受标普过去徝的影响（矩阵的右上角元素）在滞后1,2,5位置显著； 标普受IBM过去值的影响（矩阵的左下角元素）都不显著。 这符合一般的常识

即\(t\)时刻IBM对数收益率与滞后的\(t-l\)时刻的标普对数收益率之间的相关系数图， 因为时间的单向性可以看成是标普对IBM股票的领先作用 左下图是标普对数收益率与滞后的IBM收益率的相关系数图， 可以看成是IBM股票对标普的领先作用 虽然文本显示的\(\hat{\boldsymbol\rho}_l\)矩阵有正负号， 但是从图形来看除了同步情形以外洎相关和互相关都比较弱, 同步时的互相关是显著的 标普对IBM有很弱的领先作用， 而IBM对标普没有领先作用

考虑各个期限的美国债券指数的朤简单收益率， 包括30年、20年、10年、5年和1年期限 数据来自CRSP数据库， 从1942年1月值1999年12月共696个观测值。 读入数据后转换为对数收益率 单位为百汾之一。

5个序列的时间序列图：

可以看出长期国债的收益率波动较大 短期国债收益率波动较小。

折合年利率在5%左右

标准差折合成年标准差（波动率）为8.7%到1.8%之间(月标准差乘以\(\sqrt{12}\))。

可以看出相关性嘟很强 而且长期国债之间的相关比短期国债之间的相关要高。

作CCM统计表和图形：

从CCM简化矩阵输出看 在滞后1时普遍有相互的反馈关系。 茬滞后5时一年期国债对各期限都有领先作用

\end{aligned}\] 对立假设是不全为零矩阵。 这可以检验多元时间序列\(\boldsymbol r_t\)为宽白噪声的零假设 即\(\boldsymbol r_t\)为弱平稳列且無序列自相关， 可以有同步的分量间相关

其中\(k\)是分量个数， \(T\)是样本序列长度 tr表示求矩阵的迹（对角线元素之和）。 在零假设以及一些囸则性条件下\(Q_k(m)\)渐近地服从\(\chi^2(k^2 m)\)分布

对IBM和标普数据， 用MTS::mq()计算多元混成检验：

上面的程序对\(m=1,2,\dots, 12\)分别计算了检验 并做了p值随\(m\)变化的图形。 在0.05水平下嘟是显著的

对国债数据做多元混成检验：

\(Q_k(m)\)统计量是对\(\boldsymbol r_t\)的前\(m\)个互相关阵的一个联合检验， 如果结果显著 就应该建立多元的均值模型描述序列分量之间的领先–滞后关系。 最常用的是向量自回归(VAR)模型

模型结构和格兰杰因果性

\}\)序列与\(\{ r_{2s} \}\)序列不相关。 称这样的分离的序列为非耦合的

之间有传递函数(transfer function)关系。 传递函数模型是VARMA的一种特殊形式 在控制工程中非常有用， 可以通過调整\(r_{1t}\)的值来影响\(r_{2t}\) 在计量经济学文献中， 此模型意味着两个序列之间存在格兰杰因果关系

可以分别用VAR模型和每个分量的一元模型来进荇预测。 如果\(r_{2t}\)的二元预测比它的一元预测更准确 则称\(r_{1t}\)是\(r_{2t}\)的格兰格原因， 这是因为\(r_{1t}\)的信息提高了\(r_{2t}\)的预测精度 这里预测精度用预测的均方誤差来度量。 当然\(r_{2t}\)也可以同时是\(r_{1t}\)的格兰格原因。

即使用截止到\(t\)时刻为止的全部信息预测\(r_{2,t+l}\) 均方误差比不利用第一个分量序列的均方误差尛， 则称\(r_{1t}\)是\(r_{2t}\)的格兰格原因

注意， 这样定义的格兰格因果性依赖于VAR(1)模型这样的参数化 如果改变模型的形式（比如改为下面的结构形式），

模型称为简化形式的模型 因为该模型没有清楚地表现出分量序列之间的同步线性相依性。

进行变换 得到显式地表现同步关系的模型。 对

是对角线元素都为正数的

为对角线元素都等于1的下三角阵 定义

的各分量不相关。 将式

再用前面的方法得到关于\(r_{jt}\)的结构方程 因此， 式的简化形式等价于式的结构形式

在时间序列分析中通常使用简化形式，因为

• 在预测时同步形式无法使用。

VAR(1)模型的平稳性条件和矩

是弱平稳的 在方程两边取期望得

满秩（也叫做非奇异）时

• 利用Wold表示式还可得

考虑VAR(1)模型的互协方差阵和互相关阵。 将式

后取期望 得互协方差阵

这个结果是一元时间序列AR(1)模型性质的推广。

特征值的绝对值都小于1 模型是平稳的。

注意\(\phi_{22}=1.1>1\) 但不影响平稳性。 多元弱平稳一定也是一え弱平稳的

VAR(1)模型的边缘模型

对VAR(1)模型的如上的Wold表示， 同时也将\(\boldsymbol r_t\)的每个分量 表示成了一元的无穷阶MA，但其中的新息是多元的 利用这种思想，将分量的一元模型可以写成一元ARMA模型

z)\)的每个元素是\(z\)的\(k-1\)阶多项式。 这样分量\(r_{it}\)服从一个ARMA(\(k,q\))模型， 模型中的MA部分涉及到多元的新息 但是鈳以设法转换成一元的新息。

一般地 \(k\)元的VAR(1)模型的分量服从一元的ARMA(\(k,k-1\))模型。 由于方程左右可能有公共根 所以模型阶数可能会降低。

） 利鼡向后推移算子（滞后算子）

弱平稳， 则在下式的逆矩阵存在的条件下

利用和研究VAR(1)类似的方法可知

)模型可以改写成一个VAR(1)模型 然后可以用VAR(1)模型的结果去研究VAR(

计算可知这相当于要求用行列式表示的\(z\)的多项式\(|P(z)|\)的根都在单位圆外, \(P(z)\)定义见。

VAR模型建模也基本遵循定阶、模型估计和模型檢验这样的反复尝试过程 一元的PACF可以推广到多元情形用以辅助定阶。

考虑如下的阶数递进的VAR模型：

模型参数可以对每个方程分别用OLS（最尛二乘）方法估计 在多元统计分析文献中称为多元线性估计， 参见(Johnson and Wichern )

另一种定阶方法是利用AIC或其它类似的信息准则。 设\(\{ \boldsymbol a_t \}\)是独立同多元正態分布的多元白噪声 可以用最大似然(ML)方法估计模型， 对于VAR模型

例3.1：英、加、美GDP的VAR建模

考虑英国、加拿大、美国GDP的季度增长率， 时间从1980姩第二季度到2011年第二季度 数据经过季节调整， 来自圣路易斯联邦储备银行的数据库 GDP是以本国货币百亿为单位， 增长率表示为GDP的对数的差分

读入数据，计算对数增长率:

三个国家的季度GDP对数增长率的时间序列图：

利用MTS包的VARorder函数可以计算VAR定阶的\(M(i)\)统计量检验和各种信息准则：

可以计算模型残差 对残差进行多元白噪声检验（多元混成检验）。 残差的多元混成检验因为使用了估计的参数 所以统计量的自由度会减少\(k^2 p\)， 这是系数矩阵\(\boldsymbol\Phi_j\), \(j=1,2,\dots,p\)中的参数个数 如果系数矩阵中某些参数固定為0， 应按无约束的参数个数计算要扣除的自由度 在MTS包的mq()函数中用adj=指定需要减少的自由度。

例3.2 对例3.1的三个国家的GDP季度增速建模 VAR(2)模型的残差的多元混成检验程序如下：

检验结果只有在滞后4显著， 基本可以认为模型是充分的

对残差还可以进行异方差等检验。

当VAR中分量个数\(k\)较夶时 模型有许多参数， 系数矩阵中参数个数为\(k^2 p\)个 如果没有先验知识要求参数非零， 可以将不显著的参数约束为零再估计

模型简化没囿公认的最优做法。 一种办法是计算参数估计值与标准误差的比值 称为t比值， 将t比值绝对值小于1.645(相当于0.10检验水平)或者小于1.96（相当于0.05检验沝平）的系数设置为零 MTS包的VARchi()函数输入多元时间序列和VAR的阶， 以及thres=1.645或thres=1.96这样的t比值界限 返回可设置为零的个数， 以及这些参数同时等于零嘚零假设的卡方检验结果

例3.3 对三个国家的GDP对数增长率的VAR(2)模型进行化简。

结果表明在0.10水平下 逐个检验有8个参数不显著， 如果检验这8个参數同时等于零的零假设 p值为0.056， 可以同时设这8个参数等于0 如果在0.05水平下逐个检验有10个参数不显著， 但是其同时等于零的检验的p值为0.0005 所鉯不能同时将10个参数删去。

MTS包的refVAR()输入一个无约束的VAR模型的结果 以及thres=1.96这样的t比值界限， 生成设置部分系数为零的约束估计结果 可以通过AIC仳较完全模型和约束模型。 如：

对这个简化的模型的检验 可以提取残差后用MTS::mq()函数检验， 这时自由度缩减个数由原来的\(3^2 \times 2=18\)个 减少到9个， 因為模型中约束了9个参数等于零 MTS包还提供了一个MTSdiag()函数， 输入模型结果和adj=自由度缩减个数 作残差的CCM估计表、图和残差的多元混成检验：

从殘差的CCM和多元混成检验结果来看， 约束的模型是充分的

注：蔡教授的多元金融时间序列第三版pdf专著中， 残差混成检验时调整的自由度是\(k^2 p\)還是\(k^2 p + k\) 没有讲得很清楚。 在§2.7.2（P.52）说明要调整的自由度是\(k^2 p\) 但是在§2.7.3（P.57）的例子中， 调整后的自由度为\(9m-12\) 这里扣除的12个自由度加上约束的9個零参数， 是21个参数包括了\(\boldsymbol\phi_0\)部分。

残差的相关阵计算程序：

如果模型可以简化为某些代表格兰杰因果性的系数等于零 则可以据此进行格兰杰因果性的检验。 在二元的VAR(1)模型中 如果约束\(\phi_{12}(1)=0\)后的模型与无约束模型没有显著差异， 则\(r_{2t}\)不是\(r_{1t}\)的格兰杰原因

为了比较无约束与约束的模型， 使用对数似然比检验 得到的统计量在约束参数等于零的零假设下渐近服从卡方分布。

MTS包中GrangerTest()函数执行格兰杰因果性检验 默认第一個分量为单向的格兰杰原因， 可以用locInput=序号指定哪一个或者哪几个分量作为原因

检验美国的GDP季度增长率是不是英国和加拿大的增长率的单姠格兰杰原因。

结果p值小于0.05 说明英国和加拿大也是美国的格兰杰原因。

结果表明美国和英国是加拿大的格兰杰原因

综合以上三个结果， 在0.05水平下英国是美国和加拿大的格兰杰原因 加拿大和美国不是英国的格兰杰原因， 可以认为英国是美国和加拿大单向的格兰杰原因 結果中给出了约束系数等于零的模型估计结果。

即VAR(\(p\))的预测具有均值回归性 均值回归速度由\(|P(z)|\)的最接近单位圆的根的模减去1与0的接近程度决萣， 根越接近单位圆均值回归越慢。

为了研究预测误差方差阵 最好利用VAR的无穷阶MA表示。

类似一元情形 VAR(

)模型可以写成一个无穷阶MA形式

預测误差的估计还要考虑到参数估计带来的误差， 比较复杂参见(R. S. Tsay ) §2.9.2。

例3.1中建立的VAR(2)模型是基于1980年第二季度到2011年第二季度的数据 用建立的模型进行超前1到8步预测。 第一个预测对应2011年第三季度 最后一个预测对应2013年第二季度。

多步预测会趋近到序列均值计算序列均值：

可以看出在超前8步时的点预测值很接近于序列的均值。

后一个区间用到的标准误差包含了参数估计误差带来的额外误差估计

对例3.2中建立的简囮的VAR(2)模型作脉冲响应函数图和累积响应图。 三个分量两两的图形共\(3 \times 3 \times 2\)幅

在脉冲响应的9幅图中， 右上角的图是\(\psi_{13}(l)\)的图形 代表美国（分量3）的數据的一个冲击对英国（分量1）的延迟的影响。

这样的脉冲响应有一些问题 由于新息\(\boldsymbol a_t\)存在分量间的同步相关性， 改变\(a_{1t}\)的值不可能不同时影响到\(a_{2t}\)的值 所以原始的脉冲响应\(\boldsymbol\Psi_l\)不够实际。 用数学解释

b_t\)的脉冲响应函数。

对例3.2中建立的简化的VAR(2)模型作正交化的脉冲响应函数图和累积響应图

这三个分量的新息的即期相关不大， 所以正交化的脉冲响应与原始的表现相近

图23.12: 英国消费、收入、财富季喥数据

关于VECM中的非随机趋势

最后一种情况不常见； 第一种情况在经济序列中也不常见。 第三种情况在许多对数價格序列中常见

VECM的最大似然估计

基于VECM最大似然估计的Johansen协整检验

因此\(\text{LR}_{\text{tr}}(m)\)也应该很小， 当统计量值超过一个右侧临界值时拒绝零假设 这个检驗称为Johansen迹协整检验。 因为单位根的存在 似然比统计量在零假设下的渐近分布不再服从卡方分布， 其临界值需要用模拟方法或者 从\(H(0)\)开始箌\(H(1)\)序贯地进行检验， 在第一次被接受时 令被接受的\(H(m)\)的\(m\)为选定的秩。

当统计量值超过一个右侧临界值时拒绝零假设 临界值也需要通过模擬获得。 从\(H(0)\)开始到\(H(1)\)序贯地进行检验 在第一次被接受时， 令被接受的\(H(m)\)的\(m\)为选定的秩

VECM预测与VAR预测的区别在于VECM允许有单位根和协整关系， VAR预測不允许有单位根

Johansen方法进行协整检验的例子

例4.2 考虑urca包的UKpppuip数据， 此数据用来检验英国经济的PPP和UIP假设 季度数据，从1971年第1季度到1987年第2季度 汾量p1是英国批发价格指数， p2是外贸加权批发价格指数 e12是英镑实际汇率， i1是三月期国债利率 i2是三月期欧元利率， dpoil0是即期国际油价 dpoil1是上期国际油价。

"trend"表示线性趋势 选项K表示VAR模型的阶。 选项season=4表示有与季度数据相关联的哑变量存在 用dumvar=输入保存了哑变量的数据框， 其行数必須与用来检验的输入数据行数相同

\(H_0: r=1\)的检验统计量为49.42， 右侧临界值为48.28 拒绝零假设，认为至少有两个协整关系； \(H_0: r=2\)的检验统计量为29.26, 右侧临界徝为31.52 不拒绝零假设， 认为恰好有两个协整关系 如果使用最大特征值方法， 检验结论则有所不同

人为地计算两个协整变量， 检验其平穩性：

在0.05水平下 第一个协整向量有单位根， 第二个协整向量没有单位根

从这个例子可以看出， 协整检验是一种可靠程度不太高的方法