R语言建立VAR模型分析联合内生变量的动态关系

VAR模型是向量自回归模型的简称，是基于数据的统计性质建立的一种常用的计量经济模型，它把系统中每一个内生变量作为系统中所有内生变量的滞后值的函数来构造模型，从而将单变量自回归模型推广到由多元时间序列变量组成的“向量”自回归模型。本实验运用 R 语言来建立两变量的向量自回归模型，首先是检验两变量序列的平稳性，然后进行协整检验，确定待拟合模型的滞后阶数，再拟合VAR模型和做脉冲响应分析，最终对拟合的VAR模型进行分析于预测。通过本实验学会用VAR模型处理多个相关经济指标的分析与预测。

本课程难度为较难，属于中级级别课程，适合具有 R 语言基础并且有一定经管背景的用户，熟悉 R 语言基础知识并加深巩固，学会用 R 语言建立计量经济模型分析经济数据。

打开Xfce 终端，下载实验所需的数据，并启动 R ：

以下实验原理介绍来自 MBA 智库百科：向量自回归模型

向量自回归模型简称VAR模型，是一种常用的计量经济模型，1980年由克里斯托弗·西姆斯（Christopher Sims）提出。VAR模型是用模型中所有当期变量对所有变量的若干滞后变量进行回归。VAR模型用来估计联合内生变量的动态关系，而不带有任何事先约束条件。它是AR模型的推广，此模型目前已得到广泛应用。

向量自回归(VAR)是基于数据的统计性质建立模型，VAR模型把系统中每一个内生变量作为系统中所有内生变量的滞后值的函数来构造模型，从而将单变量自回归模型推广到由多元时间序列变量组成的“向量”自回归模型。VAR模型是处理多个相关经济指标的分析与预测最容易操作的模型之一，并且在一定的条件下，多元MA和ARMA模型也可转化成VAR模型，因此近年来VAR模型受到越来越多的经济工作者的重视。

本实验利用 2005 年 8 月汇率改革后消费者价格指数（CPI）和人民币名义有效汇率（NEER）的月度数据建立两变量 VAR 模型来分析人民币汇率变动与CPI之间的动态关系。

可以看到数据长度为 137，即包含了 137 个月度 CPI 和 NEER 数据。

下面我们绘制两个变量的时序图，观察消费者价格指数和名义汇率的波动情况。

可以看到消费者价格指数（CPI）的波动幅度较大，而名义有效汇率波动较小但明显呈现一个波动上升的趋势，两个变量序列大致是不平稳的。

在拟合 VAR 模型之前，需要对变量进行平稳性检验，如果要拟合的内生变量都是平稳的或者同阶单整的才可进行 VAR 模型的拟合。但在平稳性检验之前需要先项两变量取对数，以消除时间序列的异方差的影响。

判断是否平稳主要看详细拟合结果的最后两个部分，即：

1 是检验统计量的值，2 是对应的显著性水平下检验统计量的临界值。

单位根检验的原假设是序列存在单位根。

对于 lncpi 和 lnneer 的单位根检验结果:

由于两个变量都存在单位根，对序列差分后检验序列的平稳性。

由于检验的详细结果太长，截取需要的部分一阶差分平稳性检验结果：

lncpi 的一阶差分平稳性检验

由于 dlncpi 的检验统计量的值为 -8.1895，在 1%、5%、10% 的显著性水平上都小于对应的临界值，因此拒绝存在单位根的原假设，即 lncpi 是一阶差分平稳的。

lnneer 的一阶差分平稳性检验

由于 dlnneer 的检验统计量的值为 -7.0409，在 1%、5%、10% 的显著性水平上都小于对应的临界值，因此拒绝存在单位根的原假设，即 lncpi 是一阶差分平稳的。

我们再来看一阶差分后的时间序列图：

可以看到差分后的两个变量序列没有明显的波动聚集或则上升下降的趋势，比较平稳。

两个内生变量都是一阶差分平稳的，即都是一阶单整的，是不平稳的时间序列，不能做 Granger 因果检验，只能做协整检验。

协整检验主要针对非平稳的单个序列，但它们的线性组合可能是平稳的。几个变量之间可能存在的一种长期均衡关系进行检验，表现为存在某个协整方程。

由于所有变量都是一阶单整的，是非平稳时间序列，因此各变量之间可能存在协整关系，如果要对所选择的内生变量进行VAR模型的构建，需要进行协整检验，以判断各个变量之间是否存在长期稳定的协整关系，处理各变量之间的是否存在伪回归问题。在本实验我们运用E-G两步法协整检验，需要用到 lmtest 包中的 dwtest 函数。

由于拟合的模型变量和截距项的 t 检验在 5% 的显著性水平上都是显著的，且模型的 F 检验在5%的显著性水平上也是显著的，因此拟合的线性模型是合理的。

协整回归方程为：lncpi=4.80296-0.03666 lnneer+εt

由于 Durbin-Watson test 检验的 p-value 的值几乎为 0，该数值小于 0.05 ，说明在 5% 的显著性水平上残差序列不独立，具有自相关性。

检验残差序列的平稳性：

根据残差序列的平稳性检验结果，在 5% 的显著性水平上拒绝残差序列存在单位根的原假设，即残差序列是平稳的，说明 CPI 和名义有效汇率两个序列之间存在协整关系，意味着我国的消费者价格指数和名义有效汇率之间具有长期均衡关系，增长或者减少具有协同效应。

。

E-G两步法：第二步：误差修正模型的建立

协整回归方程：

误差修正项的系数为负，符合误差修正机制，反映了上一期偏离长期均衡的数量将在下一期得到反向修正，这也符合之前证明的协整关系。

在拟合 VAR 模型之前还需要确定拟合几阶 VAR 模型，也就是确定滞后阶数。确定滞后阶数需要用到 vars 包中的 VARselect 函数。

根据结果，不同的信息准则有不同的滞后阶数，选择 2 阶或者 4 阶都是可以的，一般来说选择在相同条件下更加简洁的模型，因此选择 2 阶滞后。

在确定好最优滞后阶数以后我们就可以拟合模型：

根据模型拟合结果，所得的 VAR(2) 的模型方程（保留两位小数）为：

也可以用 coef 函数来查看模型估计的简明结果。

会拟合出一致的 VAR(2) 方程（系数保留两位小数）。

在查看完拟合结果的图形之后，我们来绘制拟合结果的脉冲响应图，需要用到 var 包中的 irf 函数。然后通过 plot 函数画出图形。

根据图形 lncpi 自身以及 lnneer 的波动对 lncpi 有正向的冲击。 lncpi 对自身的影响没有滞后期，并且自身波动的影响随着时期的增加会越来越小。lnneer 波动对 lncpi 的影响在第一期以前是逐渐减少的，但在第一期减少为 0 之后随着时间的增加影响越来越大。

lncpi 波动对 lnneer 有一个负向的冲击，并且随着时期数的增加负向的影响会越来越大。lnneer 波动对自身有一个正向的冲击，这个冲击从开始先增加，在第二期达到最大值以后又逐渐减少。

模型预测通常需要用到 predict 函数，具体用法如下：

这样我们就得到了两变量（lncpi、lnneer）的 VAR(2)模型的滞后 10 期预测结果：

通过本次实验深入了解向量自回归模型，学会运用 R 语言来进行 VAR模型建模。学习如何检验一个序列的平稳性，以及如何运用协整检验来分析非平稳时间序列的关系。学会运用 var 包中的包括 VARselect() 函数、VAR() 函数、coef() 函数以及 irf() 函数各种函数来进行 VAR 模型的建模，选择最优滞后阶数，拟合 VAR 模型，对模型结果进行脉冲响应分析以及模型的预测。希望同过实验能够独立的运用 R 语言进行 VAR 模型建模于分析，巩固 R 的同时加深自己的金融知识。