高铁梅计量经济学课件

1第九章向量自回归和误差修正模型传统的经济计量方法是以经济理论为基础来描述变量关系的模型。但是，经济理论通常并不足以对变量之间的动态联系提供一个严密的说明，而且内生变量既可以出现在方程的左端又可以出现在方程的右端使得估计和推断变得更加复杂。为了解决这些问题而出现了一种用非结构性方法来建立各个变量之间关系的模型。本章所要介绍的向量自回归模型(vectorautoregression，VAR)和向量误差修正模型(vectorerrorcorrectionmodel，VEC)就是非结构化的多方程模型。2向量自回归(VAR)是基于数据的统计性质建立模型，VAR模型把系统中每一个内生变量作为系统中所有内生变量的滞后值的函数来构造模型，从而将单变量自回归模型推广到由多元时间序列变量组成的“向量”自回归模型。VAR模型是处理多个相关经济指标的分析与预测最容易操作的模型之一，并且在一定的条件下，多元MA和ARMA模型也可转化成VAR模型，因此近年来VAR模型受到越来越多的经济工作者的重视。§9.1向量自回归理论3VAR(p)模型的数学表达式是(9.1.5)其中：yt是k维内生变量向量，p是滞后阶数，样本个数为T。kk维矩阵A1，…，Ap是要被估计的系数矩阵。t是k维扰动向量，它们相互之间可以同期相关，但不与自己的滞后值相关及不与等式右边的变量相关，假设是t的协方差矩阵，是一个(kk)的正定矩阵。9.1.1VAR模型的一般表示Tt,,2,1tptpttεyAyAy114如果行列式det[A(L)]的根都在单位圆外，则式(9.1.5)满足稳定性条件，可以将其表示为无穷阶的向量动平均(VMA(∞))形式(9.1.6)其中ttLεCy)(1)()(LLAC2210)(LLLCCCCkIC05对VAR模型的估计可以通过最小二乘法来进行，假如对矩阵不施加限制性条件，由最小二乘法可得矩阵的估计量为(9.1.7)其中：当VAR的参数估计出来之后，由于A(L)C(L)=Ik，所以也可以得到相应的VMA(∞)模型的参数估计。ttTεεΣˆˆ1ˆptpttttyAyAyAyεˆˆˆˆ22116由于仅仅有内生变量的滞后值出现在等式的右边，所以不存在同期相关性问题，用普通最小二乘法(OLS)能得到VAR简化式模型的一致且有效的估计量。即使扰动向量t有同期相关，OLS仍然是有效的，因为所有的方程有相同的回归量，其与广义最小二乘法(GLS)是等价的。注意，由于任何序列相关都可以通过增加更多的yt的滞后而被消除（absorbed），所以扰动项序列不相关的假设并不要求非常严格。7例9.1我国货币政策效应实证分析的VAR模型为了研究货币供应量和利率的变动对经济波动的长期影响和短期影响及其贡献度，根据我国1995年1季度～2004年4季度的季度数据，设居民消费价格指数为P(1990年=100)、居民消费价格指数变动率为PR(P/P-1-1)*100)、实际GDP的对数，ln(GDP/P)为ln(gdp)、实际M1的对数，ln(M1/P)为ln(m1)和实际利率rr(一年期贷款利率R-PR）。利用VAR(3)模型对ln(gdp)，ln(m1)和rr，3个变量之间的关系进行实证研究，其中实际GDP和实际M1以对数的形式出现在模型中，而实际利率没有取对数。8EViews软件中VAR模型的建立和估计1．建立VAR模型为了创建一个VAR对象，应选择Quick/EstimateVAR…或者选择Objects/Newobject/VAR或者在命令窗口中键入var。便会出现下图的对话框(以例9.1为例)：9可以在对话框内添入相应的信息：(1)选择模型类型（VARType）：(2)在EstimationSample编辑框中设置样本区间(3)输入滞后信息在LagIntervalsforEndogenous编辑框中输入滞后信息，表明哪些滞后变量应该被包括在每个等式的右端。这一信息应该成对输入：每一对数字描述一个滞后区间。例如，滞后对14表示用系统中所有内生变量的1阶到4阶滞后变量作为等式右端的变量。102．VAR估计的输出VAR对象的设定框填写完毕，单击OK按纽，EViews将会在VAR对象窗口显示如下估计结果：11表中的每一列对应VAR模型中一个内生变量的方程。对方程右端每一个变量，EViews会给出系数估计值、估计系数的标准差(圆括号中)及t-统计量(方括号中)。例如，在D(logGDPTC_P)的方程中RR_TC(-1)的系数是0.000354。同时，有两类回归统计量出现在VAR对象估计输出的底部：12输出的第一部分显示的是每个方程的标准OLS回归统计量。根据各自的残差分别计算每个方程的结果，并显示在对应的列中。输出的第二部分显示的是VAR模型的回归统计量。13例9.1结果如下：3个方程调整的拟合优度分别为：可以利用这个模型进行预测及下一步的分析。)ln()1ln(928.0038.000035.0068.0029.10048.091.618.23865.10118.00145.064.1)ln()1ln(111ttttttGDPMRRGDPMRR)ln()1ln(752.0138.00029.0329.0562.00086.084.652.1626.1222tttGDPMRRtttttteeeGDPMRR321333)ln()1ln(208.0169.00017.0215.0093.00047.003.1951.15298.0697.0,746.0,986.02212GDPMRRRR14同时，为了检验扰动项之间是否存在同期相关关系，可用残差的同期相关矩阵来描述。用ei表示第i个方程的残差，i=1，2，3。其结果如表9.1所示。表9.1残差的同期相关矩阵e1e2e3e11-0.23-0.504e2-0.2310.274e3-0.5040.274115从表中可以看到实际利率rr、实际M1的ln(m1)方程和实际GDP的ln(gdp)方程的残差项之间存在的同期相关系数比较高，进一步表明实际利率、实际货币供给量(M1)和实际GDP之间存在着同期的影响关系，尽管得到的估计量是一致估计量，但是在本例中却无法刻画它们之间的这种同期影响关系。169.1.2结构VAR模型(SVAR)在式(9.1.1)或式(9.1.3)中，可以看出，VAR模型并没有给出变量之间当期相关关系的确切形式，即在模型的右端不含有当期的内生变量，而这些当期相关关系隐藏在误差项的相关结构之中，是无法解释的，所以将式(9.1.1)和式(9.1.3)称为VAR模型的简化形式。本节要介绍的结构VAR模型(StructuralVAR，SVAR)，实际是指VAR模型的结构式，即在模型中包含变量之间的当期关系。171．两变量的SVAR模型为了明确变量间的当期关系，首先来研究两变量的VAR模型结构式和简化式之间的转化关系。如含有两个变量(k=2)、滞后一阶(p=1)的VAR模型结构式可以表示为下式ztttttxtttttuzxxbbzuzxzbbx12212121201121111210(9.1.8)Tt,,2,118在模型(9.1.8)中假设：（1）变量过程xt和zt均是平稳随机过程；（2）随机误差uxt和uzt是白噪声序列，不失一般性，假设方差x2=z2=1；（3）随机误差uxt和uzt之间不相关，cov(uxt,uzt)=0。式(9.1.8)一般称为一阶结构向量自回归模型(SVAR(1))。19它是一种结构式经济模型，引入了变量之间的作用与反馈作用，其中系数b12表示变量zt的单位变化对变量xt的即时作用，21表示xt-1的单位变化对zt的滞后影响。虽然uxt和uzt是单纯出现在xt和zt中的随机冲击，但如果b210，则作用在xt上的随机冲击uxt通过对xt的影响，能够即时传到变量zt上，这是一种间接的即时影响；同样，如果b120，则作用在zt上的随机冲击uzt也可以对xt产生间接的即时影响。冲击的交互影响体现了变量作用的双向和反馈关系。ztttttxtttttuzxxbbzuzxzbbx12212121201121111210202．多变量的SVAR模型下面考虑k个变量的情形，p阶结构向量自回归模型SVAR(p)为tptptttuyΓyΓyΓyB22110(9.1.13)其中：,,111212211120kkkkbbbbbbBpiikkikikikiiikiii,,2,1,)()(2)(1)(2)(22)(21)(1)(12)(11Γkttttuuu21u21可以将式(9.1.13)写成滞后算子形式kttttELIuuuyB)'(,)((9.1.14)其中：B(L)=B01L2L2…pLp，B(L)是滞后算子L的kk的参数矩阵，B0Ik。需要注意的是，本书讨论的SVAR模型，B0矩阵均是主对角线元素为1的矩阵。如果B0是一个下三角矩阵，则SVAR模型称为递归的SVAR模型。22不失一般性，在式(9.1.14)假定结构式误差项(结构冲击)ut的方差-协方差矩阵标准化为单位矩阵Ik。同样，如果矩阵多项式B(L)可逆，可以表示出SVAR的无穷阶的VMA(∞)形式其中：ttLuDy)((9.1.15)1)()(LLBD2210)(LLLDDDD100BD23式(9.1.15)通常称为经济模型的最终表达式，因为其中所有内生变量都表示为外生变量的分布滞后形式。而且外生变量的结构冲击ut是不可直接观测得到，需要通过yt各元素的响应才可观测到。可以通过估计式(9.1.5)，转变简化式的误差项得到结构冲击ut。从式(9.1.6)和式(9.1.15)，可以得到ttLLuDεC)()((9.1.16)24上式对于任意的t都是成立的，称为典型的SVAR模型。由于C0=Ik，可得式(9.1.17)两端平方取期望，可得所以我们可以通过对D0施加约束来识别SVAR模型。ttεuD0(9.1.17)ΣDD00(9.1.18)259.2结构VAR(SVAR)模型的识别条件前面已经提到，在VAR简化式中变量间的当期关系没有直接给出，而是隐藏在误差项的相关关系的结构中。自Sims的研究开始，VAR模型在很多研究领域取得了成功，在一些研究课题中，VAR模型取代了传统的联立方程模型，被证实为实用且有效的统计方法。然而，VAR模型存在参数过多的问题，如式(9.1.1)中，一共有k(kp+d)个参数，只有所含经济变量较少的VAR模型才可以通过OLS和极大似然估计得到满意的估计结果。26为了解决这一参数过多的问题，计量经济学家们提出了许多方法。这些方法的出发点都是通过对参数空间施加约束条件从而减少所估计的参数。SVAR模型就是这些方法中较为成功的一种。9.2.1VAR模型的识别条件在经济模型的结构式和简化式之间进行转化时，经常遇到模型的识别性问题，即能否从简化式参数估计得到相应的结构式参数。27对于k元p阶简化VAR模型利用极大似然方法，需要估计的参数个数为tptpttεyAyAy11(9.2.1)222kkpk(9.2.2)而对于相应的k元p阶的SVAR模型来说，需要估计的参数个数为(9.2.4)tptpttuyΓyΓyB110(9.2.3)22kpk28要想得到结构式模型惟一的估计参数，要求识别的阶条件和