时间序列分析|VAR向量自回归

1 基本概述

1.1 概念

向量自回归模型（vector autoregressive model，简称VAR模型）是非结构性方程组模型，用于估计多个变量之间的动态关系。向量自回归模型把系统中每一个内生变量作为系统中所有内生变量的滞后值的函数来构造模型，从而实现了将单变量自回归模型推广到由多元时间序列变量组成的“向量”自回归模型。比如说存在一个系统，系统内有多个变量，VAR 模型分别将每一个变量作为因变量 Y，而系统内所有变量的滞后值作为自变量来建立方程。这样的话，系统内具有多少个变量就能够建立多少个方程式，这些式子就能描述多个变量之间的动态关系。

VAR模型常用于预测相互联系的时间序列系统以及分析随机扰动对变量系统的动态影响，主要应用于宏观经济学。是处理多个相关经济指标的分析与预测中最容易操作的模型之一。

由于向量自回归模型把每个内生变量作为系统中所有内生变量滞后值的函数来构造模型，从而避开了结构建模方法中需要对系统每个内生变量关于所有内生变量滞后值的建模问题。

1.2 内生变量与外生变量

1.21 内生变量

• 内生变量是具有某种概率分布的随机变量，它的参数是联立方程系统估计的元素，是由模型系统决定的，同时也对模型系统产生影响。
• 内生变量–般都是明确经济意义变量。
• 一般情况下，内生变量与随机项相关
• 在联立方程模型中，内生变量既作为被解释变量，又可以在不同的方程中作为解释变量。

1.22 外生变量

• 外生变量一般是确定性变量，或者是具有临界概率分布的随机变量，其参数不是模型系统研究的元素。
• 外生变量影响系统，但本身不受系统的影响。
• 外生变量一般是经济变量、政策变量、虚拟变量。
• 一般情况下，外生变量与随机项不相关。

注意：一个变量是内生变量还是外生变量，由经济理论和经济意义决定，不是从数学形式决定。

1.3 基本步骤

1.31 平稳性检验

VAR模型建立之前需要对各时间序列变量进行平稳性检验（ADF检验）。若各时间序列均是平稳序列，则可建立 VAR 模型；否则得到的向量自回归模型是伪回归。若是各数据不满足平稳性，但通过了协整检验，也可建立向量自回归模型。

1.32 确定滞后阶数

如果滞后期太少，误差项的自相关会很严重，并导致参数的非一致性估计。在VAR模型中适当加大k值（增加滞后变量个数），可以消除误差项中存在的自相关；但从另一方面看，k值又不宜过大。k值过大会导致自由度减小，直接影响模型参数估计量的有效性。

有两种方法可以做滞后阶数检验：

1. 分析各种准则，最后确定滞后阶数，AIC准则、SC准则、HQ准则、LogL准则、最终预测误差（FPE），选择k值的原则是在增加k值的过程中使相应的值达到最小。
2. 分析似然比（LR），这种方法不会出现第一种方法的无效结果。

1.33 构建VAR模型并对参数进行估计

1.34 模型稳定性检验

VAR模型稳定的充分与必要条件是所有特征值都要在单位圆以内（在以横轴为实数轴，纵轴为虚数轴的坐标体系中，以原点为圆心半径为1的圆称为单位圆），模型可以进一步做脉冲响应分析和方差分解。

1.35 脉冲响应

它描述的是VAR模型中的一个内生变量（冲击变量）的冲击给另一个内生变量（受冲击变量）所带来的影响。

1.36 方差分解

VAR的方差分解能够给出随机新息的相对重要性信息。方差分解是分析预测残差的标准差由不同冲击影响的比例，也是对应内生变量对标准差的贡献比例。

1.4 格兰杰因果检验

在有些情况下，时间序列分析也会出现伪相关问题，也就是可以计算出较大的相关系数的变量实际上并不相关。

针对此问题，格兰杰因果检验由此而生。格兰杰因果检验用于检验时间序列之间是否存在相关关系，它是能否建立脉冲函数的前提。

在VAR模型中，格兰杰检验的因果关系不是通常所说的因果关系（并非真正汉语意义上的“因果关系”），而是说先发生的事情对后发生的事情有一定的影响，或者说某个变量是否可以用来提高对其他相关变量的预测能力。所以，格兰杰因果关系的实质是一种“预测”关系。其实质是考量一个变量的滞后量能否加入到其他变量的公式中。当一个变量确实受到其他变量的滞后量影响时，可以称这两个变量具有格兰杰因果关系。

格兰杰因果检采取以下方式验证是否是真正的相关关系：

（1）估计当前的Y值被Y本身滞后期取值所能解释的程度；

（2）检验加入X的滞后期后，Y的被解释程度是否提高；

（3）如果满足条件（2），则X是Y的格兰杰成因，此时X的滞后期系数具有统计显著性。

由于格兰杰因果关系检验是在向量自回归模型的基础上进行的，因此向量自回归模型本身的合理性对格兰杰因果关系检验的结果也是非常重要的。例如，向量自回归模型本身应当具有恰当的滞后期。

几个问题：

1.如果三个变量原序列和一阶差分后序列都各有一个不平稳，但二阶差分后都平稳，请问如果要建立var模型，使用原序列还是二阶差分后序列。
（1）如果同阶单整且协整就可以建立var模型，这三个变量之间存在协整关系，所以可以用原序列建立var模型。
（2）还有一种说法是，不建立var模型，协整后直接建立vec误差修正模型。

首先，二阶差分没有经济意义。其次，若原序列有两个都平稳，而剩余的这个一阶差分平稳，就可以用两原序列加一个一阶差分的数据进行回归检验。再次，若没出现上述情况，就使用二阶数据进行协整检验，再建立误差修正模型即可。

2.为什么二阶差分不能用VAR？

因为二阶差分后无法解释其经济含义，比如说，GDP不平稳但是GDP增长率平稳，这个OK，有经济学意义，但是利率是不平稳的，利率的差分什么意义呢？

3.什么叫同阶单整呢？为什么同阶单整才可以建立VAR模型呢？

（1）单整序列
我们在第5篇笔记的第3节里对指数的对数收益率序列做了一次差分。实际我们可以在差分序列上再做差分。

当某一个不平稳的时间序列做了n次差分后第一次变成平稳序列了，那么我们就称这个时间序列为n阶单整。

根据重要统计量（之前介绍过的均值、方差、自相关系数）的可加性，我们有以下两个结论：

平稳序列的线性组合也是平稳的；
平稳序列和一阶单整序列的线性组合是一阶单整的；
那么问题来了，两个一阶单整序列的线性组合还是一阶单整的吗？答案是不一定。

（2）协整关系
但是有一种情况我们可能用两个一阶单整的序列构造出平稳序列（平稳序列的一阶差分一定也是平稳的）。

对于两个走势很像的时间序列，比如，格力和美的的股价走势，它们的对数收益率序列可以看成股价序列的一阶差分。

从前面的笔记可以假设对数收益率序列是平稳的。所以格力和美的的股价走势可以说是一阶单整的。

因为格力和美的的股价走势强相关，那么可以想象它们股价的差值序列似乎应该是平稳的。如果是平稳的，我们就说格力和美的的股价走势序列具有协整关系（Cointegration）。

更为严格的定义如下。

描述协整关系的模型和AR模型很像，都是线性模型。区别在于AR模型描述的是一个序列的不同时间点的关系；协整模型描述的是两个序列的相同时间点的关系。

协整关系通常用在配对交易中，用来判断两个走势相近的标的价格差值序列是否平稳，进而判断这两个标的走势是否平稳。

在实际应用中，我们通常会对标的价格取对数，以减小量纲的影响。

（3）模型解释与协整检验
我们针对两个序列寻找它们的协整关系，我们在它们之间建立了下面的关系：

​

由上式，我们已经假设这两个序列之间是线性关系。究竟具不具有线性关系，实际上需要验证的。所以当我们拟合出模型后，我们首先应该看看模型是否有不错的解释效果（实际上很多人都忽略了这一点）；然后我们再进行协整的检验。

验证线性模型的解释效果有很多指标，SSE（和方差）、MSE（均方误差）、RMSE（均方根误差）、MAE（平均绝对误差）、R-Squared（R平方检验，为R平方，修正R平方等）。

因为本篇笔记还是以协整为核心，所以上述方法不详述了。后面碰到实际问题再说。

那么协整性怎么检验呢？一般有两个方法，约翰森（Johansen）检验法与恩格尔-格兰杰（Engel-Granger）检验法。

约翰森（Johansen）检验法针对含有滞后项的情况，如下：

本篇笔记考虑更为简单的恩格尔-格兰杰（Engel-Granger）检验法。对于下述关系式进行两步检验：

4.VAR模型里要求所有变量都要同阶差分后平稳吗？
国内目前还是要求所有变量都平稳，然后才能做VAR，这里不是要求都统一次数差分平稳，而只是要求平稳，至于差分几次，就看实际情况了。关键问题其实并不在这儿，问题在经济意义。比如说，GDP不平稳但是GDP增长率平稳，这个OK，有经济学意义，但是利率是不平稳的，利率的差分什么意义呢？

多说一句，之所以最开始提了个国内，就是因为国外现在会直接用VAR，不考虑平稳性问题，这样能最大限度保留数据的信息。比如，Christiano et al.(2005)。在国内这么干可是会desk reject的。Stiglitz(2017)对宏观经济学家处理时间序列的方式进行了抨击，讽刺宏观经济学家只能处理平稳数据，丧失了数据信息。不过Christiano et al.(2018)对此进行了反驳，有兴趣可以找来看看。