向量自回归模型讲义

第8章VAR模型与协整1980年Sims提出向量自回归模型(vectorautoregressivemodel)。这种模型采用多方程联立的形式，它不以经济理论为基础，在模型的每一个方程中，内生变量对模型的全部内生变量的滞后值进行回归，从而估计全部内生变量的动态关系。8.1 向量自回归(VAR)模型定义8.1.1 模型定义VAR模型是自回归模型的联立形式，所以称向量自回归模型。假设y*y2t之间存在关系，如果分别建立两个自回归模型yi,t=f(yi,t-i,yi,t-2,)y2,t=f(y2,t-i,y2,t-2,)则无法捕捉两个变量之间的关系。如果采用联立的形式，就可以建立起两个变量之间的关系。VAR模型的结构与两个参数有关。一个是所含变量个数N，一个是最大滞后阶数k。以两个变量yit，y2t滞后i期的VAR模型为例，y1,t=C1+11.1y1,t-1+12.1y2,t-1+U1ty2,t=C2+21.1y1,t-1+22.1y2,t-1+U2t(8.1)IID (0,2), Cov(U1 t, U2 t) = 0。写其中U1t,U2t成矩阵形式是,设,则，那么,y1ty2tC1C2yty2t11.121.1C1, c = c2 ,12.122.1Yt = c +1 Yt-1 + Ut含有示如下:Yt = c +y1,t 1y2,t 111.121.1N个变量滞后k期的1 Yt-1 +2 Yt-2 + +Ut IID (0,)其中，Yt = (y1, ty2, t C = (C1C2 cn)yN, t)11.j21. j12.j22.j1N.j2N.jN1.jN2.jU1tU2t12.122.1(8.2)U1t,Ut = U2t(8.3)VAR模型表k Yt-k + Ut,(8.4),j = 1,2, k,NN.jUt=(U1tU2,tUNt),Yt为N1阶时间序列列向量。C为N1阶常数项列向量。i,，k均为NN阶参数矩阵，utIID(0,)是N1阶随机误差列向量，其中每一个元素都是非自相关的，但这些元素，即不同方程对应的随机误差项之间可能存在相关。因VAR模型中每个方程的右侧只含有内生变量的滞后项，他们与ut是渐近不相关的，所以可以用OLS法依次估计每一个方程，得到的参数估计量都具有一致性。估计VAR的EViews4.1操作：打开工作文件，点击Quick键,选EstimateVAR功能。作相应选项后，即可得到VAR的表格式输出方式。在VAR模型估计结果窗口点击View选representation功能可得到VAR的代数式输出结果。8.1.2 VAR模型的特点是：（ 1） 不以严格的经济理论为依据。在建模过程中只需明确两件事：共有哪些变量是相互有关系的，把有关系的变量包括在VAR模型中；确定滞后期ko使模型能反映出变量间相互影响的绝大部分。（ 2） VAR模型对参数不施加零约束。（对无显着性的参数估计值并不从模型中剔除，不分析回归参数的经济意义。）（ 3） VAR模型的解释变量中不包括任何当期变量，所有与联立方程模型有关的问题在VAR模型中都不存在（主要是参数估计量的非一致性问题）。（ 4） VAR模型的另一个特点是有相当多的参数需要估计。比如一个VAR模型含有三个变量，最大滞后期k=3，则有kN2=332=27个参数需要估计。当样本容量较小时，多数参数的估计量误差较大。（5）无约束VAR模型的应用之一是预测。由于在VAR模型中每个方程的右侧都不含有当期变量，这种模型用于样本外一期预测的优点是不必对解释变量在预测期内的取值做任何预测。（6）用VAR模型做样本外近期预测非常准确。做样本外长期预测时，则只能预测出变动的趋势，而对短期波动预测不理想。西姆斯（Sims）认为VAR模型中的全部变量都是内生变量。近年来也有学者认为具有单向因果关系的变量，也可以作为外生变量加入VAR模型。附录：（file:B8c1）VAR模型静态预测的EViews操作：点击Procs选MakeModel功能。点击Solve。在出现的对话框的Solutionoption（求解选择）中选择Staticsolution（静态解）。VAR模型动态预测的EViews操作：点击Procs选MakeModel功能（工作文件中如果已经有Model，则直接双击Model）。点击Solve。在出现的对话框的Solutionoption（求解选择）中选择Dynamicsolution（静态解）。注意：Model窗口中的第一行，“ASSIGNALLF”表示模拟结果保存在原序列名后加F的新序列中，以免原序列中的数据被覆盖掉。静态预测的效果非常好。动态预测的表现是前若干期预测值很接近真值，以后则只能准确预测变化的总趋势，而对动态的变化特征预测效果较差。综上所述，用VAR做样本外动态预测1，2期则预测效果肯定是非常好的。8.2 VAR模型稳定的条件VAR模型稳定的充分与必要条件是(1见(8.(3) 式)的所有特征值都要在单位圆以内(在以横轴为实数轴，纵轴为虚数轴的坐标体系中，以原点为圆心，半径为1的圆称为单位圆)，或特征值的模都要小于1。1 先回顾单方程情形。以AR(2)过程yt=1yt-1+2yt-2+ut(8.11)为例。改写为(1-1L-2L2)yt=(L)yt=ut(8.12)yt稳定的条件是(L)=0的根必须在单位圆以外。2 对于VAR模型，也用特征方程判别稳定性。以(8.3)式，Yt=c+iYt-i+ut,为例，改写为(I-1L)Yt=c+ut(8.13)保持VAR模型稳定的条件是|I-1L|=0的根都在单位圆以外。|I-iL|=0在此称作相反的特征方程(reversecharacteristicfunction)。(第2章称特征方程)例8.1以二变量(N=2)，k=1的VAR模型y1t5/81/2yi,tiu1ty2t=1/45/8y2,ti+u2t(8.14)5/81/2其中1=1/45/8为例分析稳定性。相反的特征方程是II.1L|=10(5/8)L(1/2)L=I11I01(1/4)L(5/8)L=(1-(5/8)L)2-1/8L2=(1-0.978L)(1-0.27L)=0(8.15)求解得L1=1/0.978=1.022,L2=1/0,27=3.690因为L1,L2都大于1,所以对应的VAR模型是稳定的。3 .VAR模型稳定的另一种判别条件是，特征方程|1-I|=0的根都在单位圆以内。特征方程|1-I|=0的根就是1的特征值。例8.2仍以VAR模型(8.14)为例，特征方程表达如下：5/8 1/25/81/21/4 5/801/45/8(5/8-)2-1/8=(5/8-)2二.1/8)2=(0.978-)(0.271-)=0(8.16)得1=0.9786,2=0.2714。1,2是特征方程|1-I|=0的根，是参数矩阵1的特征值。因为1=0.978,2=0.271,者B小于1,该VAR模型是稳定的。注意：(1)因为L1=1/0.978=1/1,L2=1/0.27=1/2,所以特征方程与相反的特征方程的根互为倒数，L=1/。(2)在单方程模型中，通常用相反的特征方程(L)=0的根描述模型的稳定性，即单变量过程稳定的条件是(相反的)特征方程(L)=0的根都要在单位圆以外；而在VAR模型中通常用特征方程|1-I|=0的根描述模型的稳定性。VAR模型稳定的条件是，特征方程|1-I|=0的根都要在单位圆以内，或相反的特征方程|I-L1|=0的根都要在单位圆以4.对于k1的k阶VAR模型可以通过友矩阵变换(companionform),改写成1阶分块矩阵的VAR模型形式。然后利用其特征方程的根判别稳定性。具体变换过程如下。给出k阶VAR模型，Yt=c+1Yt-1+2Yt-2+kYt-k+ut(8.17)再配上如下等式，Yt-1=Yt-1Yt-2=Yt-2Yt-k+1=Yt-k+1把以上k个等式写成分块矩阵形式，Ytc口1口 2Yt 10I0Yt 2=0+0IMYt k 1 NK 10 NK 100k1nkYt1ut00Yt2000Yt3+0I0 NK NKYtk NK 10 NK 1(8.18)其中每一个元素都表示一个向量或矩阵。令Yt=(Yt-1Yt-2Yt-k+i)nk1C=(c000)nk1口i口2n；k1iikI00000A=0I000NKNKUt = (ut0 0) NK 1上式可写为Yt=C+AYt-1+Ut（8.19）注意，用友矩阵变换的矩阵（向量）用正黑体字母表示。k阶VAR模型用友矩阵表示成了1阶分块矩阵的VAR模型。例如，2变量2阶VAR模型的友矩阵变换形式是YtYt 1+ I1Yt 1 utYt 2 + 0(8.20)其中等式的每一个元素（项）都表示一个41阶向量或44阶矩阵。例如，2变量3阶VAR模型的友矩阵变换形式是(8.21)Ytc123Yt1utY-=0+I00Yt2+0Yt200I0Yt30其中等式的每一个元素（项）都表示一个61阶向量或66阶矩阵。VAR模型的稳定性要求A的全部特征值，即特征方程|A-I|=0的全部根必须在单位圆以内或者相反的特征方程|I-LA|=0的全部根必须在单位圆以外。注意：特征方程中的A是NkNk阶的。特征方程中的I也是NkNk阶的。以2阶VAR模型的友矩阵变换为例,I1L2ILI=I-1L-2L2=0(8.22)的全部根必须在单位圆以外。以3阶VAR模型的友矩阵变换为例,I00|I-AL|=0I0I00l00I0I0IiL2L3LILI00ILI=|I-1L-2L2-3L3|=0(8.23)的全部根必须在单位圆以外。因此，对于k阶VAR模型的友矩阵变换形式，特征方程是，|I-1L-2L2-kLk|=0(8.24)附录：求VAR模型特征根的EViews4.1操作：在VAR模型估计结果窗口点击View选LagStructrure,ARRootsTable功能，即可得到VAR模型的全部特征根。若选LagStructrure,ARRootsGraph功能，即可得到单位圆曲线以及VAR模型全部特征根的位置图。VARStabilityConditionCheckRootsofCharacteristicPolynornialEndogenousriables:LNGPLNCPLNIPE?（口genQusvariables；CLagspecification:13Date:05/16/05Time:22:23RaolModulus1.654971.0154970.466948-0.673264I118193440.466948+0.673264i0,8W344口,0057360.805736-0.220604-0.488752i0.536232-0.220604+0.4S3752i0.5362320,440194C4401940.158B44-022935210.27B9870J68844+0.22935210,278987Warning:Allestonerootoulsidstheunitcircle.VARd口日5notsatisfythestabilitycoridition.InverseRootsofARCharacteristicPolynomialVARStabilityConditionCheckRootsofCharacteristicPolynomialEndogenousvariablesPHOQHONHOExogenousYariribleEi：CLagspecificalion:18Date：05/1M5Time：22:19RootModulus0. S41175-0.491343i 9741630.841175+0.491343) 9741630.966912 9669120.9446240.9446240.414853 + 079351210 8954180.4148B3-0793512! 3954180768707 +Q.J32410i,062064-0766787-0.43241010.9020540.115255+0.8173291 B2541501115256-0,8173291 32S415-0.4E2S31 -0.662775i 308211-0.462531 +0.68277510.908211-0.146717-0.7B5193i0.7987930.146717+0.78513317937S30.692431 - 0.337520i0.7934940.692431 + 0.387520 i0.7934940.7556B3 - 0.2272750.7091210L7556B3 + D,227275i7B91210764318 +Q10Bie?i771657-0764318-0.106167!,7716570 0905B2 十 0.58723110 5941760L0905B2-0.537231 i0.594176.0.419544 +0.3756991,563176-0.419544-0.376699!0.563176Norootliasoutsidetheunitcircle.VARsatisfiesIhestabilitycondilion.Inverse Roots of AR Characteristic Polynomial8.3VAR模型的稳定性(stability)特征现在讨论VAR模型的稳定性特征。稳定性是指当把一个脉动冲击施加在VAR模型中某一个方程的新息(innovation)过程上时，随着时间的推移，这个冲击会逐渐地消失。如果是不消失，则系统是不稳定的。下面分析一阶VAR模型Yt = c+ 1 Yt-1 + u t(8.29)t=1时，有Y1=c+1Y0+u1(8.30)t=2时，采用迭代方式计算，Y2 = c +1 Y1 + u2 = c += (I + 1) c +12 Y0 +当 t = 3 时，进一步迭代，1 (c +1 Y0 + u1) + u21 u1 + u2(8.31)Y3=c+1Y2+u3=c+1(I+1)c+12Y0+1u1+u2+u3=(I+1+12)c+13Y0+12u1+1u2+u3(8.32)对于t期，按上述形式推导t1iYt=(I+i+12+it-1)c+itYo+i产ut-i(8.33)由上式可知，10=I。通过上述变换，把Yt表示成了漂移项向量、初始值向量Yo和新息向量ut的函数。可见系统是否稳定可以通过观察漂移项向量c、初始值向量Yo和新息向量ut经受冲击后的表现。假定模型是稳定的，将有如下3个结论。( 1) 假设t=1时，对c施加一个单位的冲击，那么到t期的影响是(I+1+12+1t-1)当t时，此影响是一个有限值，(I-1)-1。( 2) 假设在初始值Y0上施加一个单位的冲击。到t期的影响是it。随着t,1t0,影响消失(因为对于平稳的VAR模型，1中的元素小于1，所以随着t，取t次方后，1t0)。t1i(3)从i0n，ut-i项可以看出，白噪声中的冲t1击离t期越远，影响力就越小。nii=(I-1)-1,i0t1称作长期乘子矩阵，是对i0n1iut-i求期望得到的。对单一方程的分析知道，含有单位根的自回归过程对新息中的脉动冲击有长久的记忆能力。同理，含有单位根的VAR模型也是非平稳过程。当新息中存在脉动冲击时，VAR模型中内生变量的响应不会随时间的推移而消失。平稳变量构成的一定是稳定(stability)的模型，但稳定的模型不一定由平稳变量构成。也可能由非平稳(nonstationary)变量(存在协整关系)构成。8.4VAR模型滞后期k的选择建立VAR模型除了要满足平稳性条件外，还应该正确确定滞后期k。如果滞后期太少，误差项的自相关会很严重，并导致参数的非一致性估计。正如在第4章介绍ADF检验的原理一样，在VAR模型中适当加大k值(增加滞后变量个数)，可以消除误差项中存在的自相关。但从另一方面看，k值又不宜过大。k值过大会导致自由度减小，直接影响模型参数估计量的有效性。下面介绍几种选择k值的方法。1)用LR统计量选择k值。LR(似然比)统计量定义为，22LR=-2(logL(k)-logL(k+1)(N)(8.34)其中logL(k)和logL(k+1)分别是VAR(k)和VAR(k+1)模型的极大似然估计值。k表示VAR模型中滞后变量的最大滞后期。LR统计量渐近服从2(N2)分布。显然当VAR模型滞后期的增加不会给极大似然函数值带来显着性增大时，即LR统计量的值小于临界值时，新增加的滞后变量对VAR模型毫无意义。应该注意，当样本容量与被估参数个数相比不够充分大时，LR的有限样本分布与LR渐近分布存在很大差异。2)用赤池(Akaike)信息准则(AIC)选择k值。cTiUt2处AIC=logT+T(8.34)其中ut表示残差，T表示样本容量，k表示最大滞后期。选择k值的原则是在增加k值的过程中使AIC的值达到最小。EViews3.0的计算公式是10gL2kAIC=-2T+73)用施瓦茨(Schwartz)准则(SC)选择k值。：WklogTSC=logt-t+(8.35)其中ut表示残差，T表示样本容量，k表示最大滞后期。选择最佳k值的原则是在增加k值的过程中使SC值达到最小。EViews3.0的计算公式是logLklogTsc=-2-y+例8.3以第8章案例为例，k=1、2、3、4时的logL、AkaikeAIC和SchwarzSC的值见下表。10gL-2 (log L(k) - logL(k+1)Akaike AICSchwarz SCVAR(1) VAR(2) VAR(3) VAR(4)184.6198.9200.0207.828.6-7.84-7.362.215.6-8.27-7.41-8.09-6.852(9)16.9-8.23-6.6LacjL)ogLLRFPEAICSCHQ022 93003MA7 75E-05-0.951430-0.327311-0 9059351175.5519276 06258.34E-0a-77631677.291710闱中IF2195 793233,67880*4.93 日00*8.321857*-7.453032*-8003404*3199 96546.409457G.31E-08-8 093591-6.852398-7,B38G45VARLagOrderSlettiarCnlariaEndogenousvariables:LNGPLNCPLNIPExogenous.enables:CDate:05f25ffiTimeZ3:52Sampler19631907Includedobser/ation&-42+indicatesImgorderselectedbythecriterionLR:sequentialmodifiedLRteststatistic(eachtestat5%level)FPE:FinalpredictionerrorAIC:AkaikeinformationcriterionSC:SchwarzirifarmationcriterionHQ:Hannan-Quinninformalioricriterion建立滞后2期的VAR模型是可以的附录：考察VAR模型最大滞后期的EViews4.1操作：在VAR模型估计结果窗口点击View选LagStructrure,LagLengyhCriteria功能，即可得到5个评价统计量的值。8.5VAR模型的脉冲响应函数由于VAR模型参数的OLS估计量只具有一致性，单个参数估计值的经济解释是很困难的。要想对一个VAR模型做出分析，通常是观察系统的脉冲响应函数(1)脉冲响应函数。脉冲响应函数描述一个内生变量对误差冲击的反应。具体地说，它描述的是在随机误差项上施加一个标准差大小的冲击后对内生变量的当期值和未来值所带来的影响。对于如下VAR模型，yi,t表示GDP,y2,t表示货币供应量，yi,t=ci+ii.iyi,t-i+i2.iy2,t-i+uit、y2,t=C2+2i.iyi,t-i+22.iy2,t-i+U2t(8.36)在模型(8.36)中，如果误差uit和u2t不相关，就很容易解释。uit是yi,t的误差项；u2t是y2,t的误差项。脉冲响应函数衡量当期u1t和u2t一个标准差的冲击分别对GDP和货币存量的当前值和未来值的影响。对于每一个VAR模型都可以表示成为一个无限阶的向量MA(8)过程。具体方法是对于任何一个VAR(k)模型都可以通过友矩阵变换改写成一个VAR(1)模型(见8.1.2节)。Yt=A1Yt-1+Ut(I-LA1)Yt=UtYt=(I-LA1)-1Ut=Ut+AiUt-i+Ai2Ut-2+AisUt-s+这是一个无限阶的向量MA(oo)过程。或写成，Yt+s=Ut+s+AiUt+s-i+Ai2Ut+s-2+AisUt+全部的移动平均参数矩阵用改用j,(j=i,-s)表示，Yt+s=Ut+s+iUt+s-i+2Ut+s-2+sUt+(8.37)其中i=Ai,2=Ai2,，s=Ais,显然，由(8.37)式有下式成立，_Ytss=Uts中第i行第j列元素表示的是，令其它误差项在任何时期都不变的条件下，当第j个变量yjt对应的误差项ujt在t期受到一个单位的冲击后，对第i个内生变量yit在t+s期造成的影响。把s中第i行第j列元素看作是滞后期s的函数y，t sujt,s = 1,2, 3,称作脉冲响应函数(impulse-responsefunction),脉冲响应函数描述了其它变量在t期以及以前各期保持不变的前提下，yi,t+s对uj,t时一次冲击的响应过程。对脉冲响应函数的解释出现困难源于实际中各方程对应的误差项从来都不是完全非相关的。当误差项相关时，它们有一个共同的组成部分，不能被任何特定的变量所识别。为处理这一问题，常引入一个变换矩阵M与ut相乘，vt=Mut(0,)从而把ut的方差协方差矩阵变换为一个对角矩阵。现在有多种方法。其中一种变换方法称作乔利斯基（Cholesky）分解法，从而使误差项正交。原误差项相关的部分归于VAR系统中的第一个变量的随机扰动项。在上面的例子里，u1t和u2t的共同部分完全归于Ult,因为Ult在U2t、*、八之前。虽然乔利斯基分解被广泛应用，但是对于共同部分的归属来说，它还是一种很随意的方法。所以方程顺序的改变将会影响到脉冲响应函数。因此在解释脉冲响应函数时应小心。注意：对于Ut中的每一个误差项，内生变量都对应着一个脉冲响应函数。这样，一个含有4个内生变量的VAR将有16个脉冲响应函数。附录：VAR模型残差序列及其方差、协方差矩阵的求法。点击VAR窗口中的Procs键，选MakeResiduals（生成残差）功能，工作文件中就会生成以residOI,resid02,为编号的残差序列（残差序列的顺序与VAR模型估计对话框中输入的变量顺序相一致），并打开残差序列数据组窗口。在这个残差序列数据组窗口中点击View键，选择Covariances功能，即可得到残差序列的方差、协方差矩阵。选择Correlation功能，即可得到残差序列的相关系数矩阵。CovarianceMatrixRESID01RESID02RESID03RESID012585272-.056293*7.166736RESID02QQ5G29315240536.196976RESID03厅1667366.1909753E7537ECorrelationMatrixRESID01RESIDQ2RE5ID03RESID011.DOOQOO-0.002836-0.232497RESID02-0.002B3G1.0000000.261936RESID03-0.232497,2616351.000000附录：脉冲响应的EViews操作(file:VAR01)点击VAR窗口中的Impulse键。在随后弹出的对话框中做出各项选择后点击OK键。例8.4美国民用燃油价格、生产量、储量的脉冲呼1应图。ResponseNFHOt?OneS.D.ImovitiansRuponseQHOtoOneS.D,InnovationsRupon01NHOtoOrw轧D.InriQVitNnf2i&IO1214-lg1B20222424AglID12I*Ifi192022242g白IO12，*IS恒2D2M图1-ofthreeinnovations图2图3图1表示美国民用燃油价格（PHO）分别对燃油价格（PHO）、生产量（QHO）、储量（NHO）3个方程相应新息过程一个标准差冲击的响应。图2表示美国燃油生产量（QHO）分别对燃油价格（PHO）、生产量（QHO）、储量（NHO）3个方程相应新息过程一个标准差冲击的响应。图3表示美国燃油储量（NHO）分别对燃油价格（PHO）、生产量（QHO）、储量（NHO）3个方程相应新息过程一个标准差冲击的响应。8.6格兰杰非因果性检验VAR模型还可用来检验一个变量与另一个变量是否存在因果关系。经济计量学中格兰杰（Granger）非因果性定义如下：格兰杰非因果性：如果由yt和xt滞后值所决定的yt的条件分布与仅由yt滞后值所决定的条件分布相同，即（ytyt-i,xt-1,）=（ytyt-1,）,（8.38）则称xt-1对yt存在格兰杰非因果性。格兰杰非因果性的另一种表述是其它条件不变，若加上xt的滞后变量后对yt的预测精度不存在显着性改善，则称xt-1对yt存在格兰杰非因果性关系。为简便，通常总是把xt-1对yt存在非因果关系表述为xt(去掉下标-1)对yt存在非因果关系(严格讲，这种表述是不正确的)。在实际中，除了使用格兰杰非因果性概念外，也使用“格兰杰因果性”概念。顾名思义，这个概念首先由格兰杰(Granger1969)提出。西姆斯(Sims1972)也提出因果性定义。这两个定义是一致的。根据以上定义，xt对yt是否存在因果关系的检验可通过检验VAR模型以yt为被解释变量的方程中是否可以把xt的全部滞后变量剔除掉而完成。比如VAR模型中以yt为被解释变量的方程表示如下：kkyt=i1iyti+i1ixti+u1t(8.39)如有必要，常数项，趋势项，季节虚拟变量等都可以包括在上式中。则检验xt对yt存在格兰杰非因果性的零假设是Ho：1=2=k=0显然如果(8.39)式中的xt的滞后变量的回归参数估计值全部不存在显着性，则上述假设不能被拒绝。换句话说，如果xt的任何一个滞后变量的回归参数的估计值存在显着性，则结论应是Xt对yt存在格兰杰因果关系。上述检验可用F统计量完成。(SSESSEu)kF=SSEu(TkN)其中SSEr表示施加约束(零假设成立)后的残差平方和。SSEu表示不施加约束条件下的残差平方和。k表示最大滞后期。N表示VAR模型中所含当期变量个数，本例中N=2,T表示样本容量。在零假设成立条件下，F统计量近似服从F(k,t-kn)分布。用样本计算的F值如果落在临界值以内，接受原假设，即xt对yt不存在格兰杰因果关系。例8.5:(file:stock)以661天(1999.1.4-2001.10.5)的上海(SH)和深圳(SZ)股票收盘价格综合指数为例，7 0 06 0 05 0 04 0 03 0 01 0 02 0 03 0 04 0 05 0 06 0 02 0 0 015 0 010 0 0滞后10期的Granger因果性检验结果如下:（当概率小于0.05时，表示推翻原假设）Paiiv/isoGrangerCausalityTestsDato:06/1B/04Timo:19:56Sample:1661Lags；:10NullHypothesis：ObsF-StatisticProbabilitySHd。日普notGrangerCsuseSZSZdoetnotGraingerCauseSH氏11.3637523,43950.193160.00000上表中概率定义为,P(F1.36)=0.19316图示如下:图7P（F23.44）=0.00000因为F值（1.36）落在原假设接受域，所以原假设“上海股票价格综合指数对深圳股票价格综合指数不存在Granger因果关系”被接受。因为F值（23.44）落在原假设拒绝域，所以原假设“深圳股票价格综合指数对上海股票价格综合指数不存在Granger因果关系”被推翻。附录：格兰杰因果关系检验EViews操作方法是，打开数剧组窗口，点View键，选GrangerCausility。在打开的对话窗口中填上滞后期（上面的结果取滞后期为10）,点击OK键。用滞后5,10,15,20,25期的检验式分别检验，结果见下表：k=5k=10k=15k=20k=25Ho：SHdoesnotGranger1.081.361.211.291.40接受HoCauseSZHo:SZdoesnotGranger43.923.415.912.610.3拒名HoCauseSH结论都是上海股票价格综合指数不是深圳股票价格综合指数变化的原因，但深圳股票价格综合指数是上海股票价格综合指数变化的原因。注意：（1）滞后期k的选取是任意的。实质上是一个判断性问题。一般来说要试检验若干个不同滞后期k的格兰杰因果关系检验，且结论相同时，才可以最终下结论。（2）当做xt是否为导致yt变化的格兰杰原因检验时，如果zt也是yt变化的格兰杰原因，且zt又与xt相关，这时在xt是否为导致yt变化的格兰杰因果关系检验式的右端应加入zt的滞后项（实际上是3个变量VAR模型中的一个方程）。（3）不存在协整关系的非平稳变量之间不能进行格兰杰因果关系检验。8.7VAR模型与协整如果VAR模型Yt=1Yt-1+2Yt-1+kYt-k+ut,utIID（0,）的内生变量都含有单位根，那么可以用这些变量的一阶差分序列建立一个平稳的VAR模型。Yt=1*Yt-i+2*Yt-2+Tk*Yt-k+ut*然而，当这些变量存在协整关系时，采用差分的方法构造VAR模型虽然是平稳的，但不是最好的选择。如果YtI(1)，且非平稳变量间存在协整关系。那么由这些非平稳变量组成的线性组合则是平稳的。建立单纯的差分VAR模型将丢失重要的非均衡误差信息。因为变量间的协整关系给出了变量间的长期关系。同时用这种非均衡误差以及变量的差分变量同样可以构造平稳的VAR模型。从而得到一类重要的模型，这就是向量误差修正模型。对于k阶VAR模型，Yt=iYt-i+2Yt-2+kYt-k+ut,其向量误差修正模型(VEC)的表达式是Yt=(1+2+k-I)Yt-i-(2+3+k)Yt-i-(3+k)Yt-2-kYt-(k-1)+ut(8.45)k=i-I=1+2+k-Iiikj=-ij1i,j=1,2,k-1)Yt=Yt-1+1Yt-1+2Yt-2+k-1Yt-(k-1)+ut(8.47)称为压缩矩阵(impactmatrix，影响矩阵)。是全部参数矩阵的和减一个单位阵。为多项式矩阵，其中每一个元素都是一个多项式。运算规则于一般矩阵相同。滞后期的延长不影响对协整向量个数的分析。根据Granger定理，向量误差修正模型(VEC)的表达式是A?(L)(1-L)Yt=Yt-1+d(L)ut(8.48)其中A?(L)是多项式矩阵A(L)分离出因子(1-L)后降低一阶的多项式矩阵，d(L)是由滞后算子表示的多项式矩阵。上式与(8.47)式完全相同。其中A?(L)(1-L)Yt=A?(L)Yt=Yt-1Yt-1-2Yt-2-k-1Yt-(k-1)d(L)ut=ut在这里d(L)退化为单位列向量。若YtCI(1,1)，比较(8.47)和(8.48)式必然有其中是协整矩阵，是调整系数矩阵。和都是Nr阶矩阵。表示有r个协整向量，1,2，r,存在r个协整关系。因为YtI(1),所以YtI(0)。从VAR模型变换为模型向量误差修正模型称为协整变换(cointegratingtransformation)。压缩矩阵决定模型VAR模型中是否存在，以及以什么规模存在协整关系。若YtCI(1,1)，YtI(0)，所以除了Yt-1，模型中各项都是平稳的。而对于Yt-1有如下三种可能。1)当Yt的分量不存在协整关系，的特征根为零，=0。2)若rank()=N(满秩)，保证Yt-1平稳的唯一一种可能是YtI(0)。3)当YtI(1),若保证Yt-i平稳，只有一种可能，即Yt的分量存在协整关系。YtI(0)rank()=rN,可以分解=和都是Nr阶矩阵。其中是协整矩阵，是调整系数矩阵。表示有r个协整向量，2,，存在r个协整关系。假定YtI(1)具有一般性。如果某个变量的单整阶数高于1，可通过差分取其相应单整阶数为1的序列加入模型。上式也可以加入位移项与趋势项。若=成立，且存在r个协整关系，则Yt-1的一般表达式是Yt-1 =111r212rN1NrNr111r212rN1Nr11 ( 11 y1,t 1.N1 ( 11 y1,t 1.8.8 单位根与Yt-1111Nr1rN11 y1,t 1.r1 y1,t 1.NryN,t 1)1r ( r1yN,t 1)Nr ( r1y1,t 1y2,t 1r N yN,t 1 N 1.1NyN,t 1.rN yN,t 1 r 11,t 1. rN yN,t 1)1,t 1. rN yN ,t 1 ) N 1(oVAR模型中存在单位根与压缩矩阵降秩的关系。以k=1的VAR模型为例，Yt=1Yt-1+ut（8.56）它的VEC表达式是Yt=Yt-1+ut。如果VAR模型中存在单位根，一定降秩。8.9VAR模型中协整向量的估计与检验8.9.1 VAR模型中协整向量的估计（略）8.9.2 VAR模型中协整向量的检验检验存在r个协整向量（存在r个协整关系）,即（N-r）个非协整向量，或者（N-r）个单位根，可以表达为特征方程|-I|=0相应（N-r）个特征值，r+1,，N,为零。几个等价检验:1）存在r个协整关系应特征方程|-I|=0存在r个非零特征值（根）2）存在r个协整关系冬VAR模型中存在（N-r）个单位根3）存在r个协整关系应特征方程|-I|=0相应（N-r）个特征值r+1,，N,为零。协整检验是一个连续检验过程。Yt=Yt-i+iYt-i+k-1Yt-（k-1）+Dt+ut（1） 首先从检验r=0开始。意即在VAR模型中不存在协整向量（含有N个单位根）。如果r=0不能被拒绝，说明N个变量间不存在协整关系。检验到此终止。不能建立VEC模型。如果r=0被拒绝，则应继续进行下面的检验。（2） r1。意即在VAR模型中存在1个协整向量（含有N-1个单位根）。如果r1不能被拒绝，检验到此终止。如果r1被拒绝，则应进一步作如下检验。（3） rNT。意即在VAR模型中存在NT个协整向量（含有1个单位根）。如果rN-1不能被拒绝，检验到此终止。如果rN-1被拒绝，说明r=N。在检验过程中，比如rr*-1已经被拒绝，但rr*不能被拒绝，则结论是VAR模型（8.61）中存在r*个协整向量。（4） 4）协整检验过程中的每一步检验都属于右单端检验。8.3.3VEC模型中确定项的处理1常数项的处理VEC模型中常数项的位置可分3种情形讨论。位置不同，相应的协整检验用表也不同。（ 1）常数项完全属于协整空间。（ 2）常数项的一部分进入协整空间，一部分属于数据空间（VAR的常数项）。（ 3）常数项只进入数据空间（VAR的常数项），不进入协整空间。2趋势项的处理对时间趋势项t进入协整空间，表示变量协整关系中也存在线性趋势。进入数据空间（VAR的常数项）。表示原变量中存在二次方的时间趋势项，或差分变量中存在一次方的时间趋势项，8.9.4协整检验用表根据和t所在位置不同，检验协整关系的LR统计量的分布也不同。检验时应选择相应的临界值表。附表1给出了5种模型条件下所对应的临界值。附表1VAR模型协整检验临界值表(迹统计量)单位根个数模型类型N-r0.100.050.0112.863.846.51模型(1)210.4712.5316.31=0,=0321.6324.3129.75协整空间中无常数项、无趋势项。436.5839.8945.58数据空间中无均值、无趋势项。555.4459.4666.52678.3682.4990.457104.77109.99119.808135.24141.20152.329169.45175.77187.3110206.05212.67226.4011248.45255.27269.8117.529.2412.97模型(2)217.8519.9624.6010,2=0,=0332.0034.9141.07协整空间中有常数项、无趋势项。449.6553.1260.16数据空间中无均值、无趋势项。571.8676.0784.45697.18102.14111.017126.58131.70143.098159.48165.58177.209196.37202.92215.7410236.54244.15257.6811282.45291.40307.6412.693.766.65模型(3)213.3315.4120.0410,20,=0326.7929.6835.65协整空间中有常数项、无趋势项。443.9547.2154.46数据空间中有线性趋势、无二次趋势项。564.8468.5276.07689.4894.15103.187118.50124.24133.578150.53156.00168.369186.39192.89204.9510225.85233.13247.1811269.96277.71293.44110.4912.2516.26模型(4)222.7625.3230.4510,20,10,2=0339.0642.4448.45协整空间中有常数项、有线性趋势项。459.1462.9970.05数据空间中有线性趋势、无二次趋势项。583.2087.3196.586110.42114.90124.757141.01146.76158.498176.67182.82196.089215.17222.21234.4110256.72263.42279.0711303.13310.81327.4512.573.746.40模型（5）216.0618.1723.4610,20,10,20331.4234.5540.49协整空间中有常数项、有线性趋势项。450.7454.6461.24数据空间中有一次趋势项。573.4077.7485.786100.14104.94114.367130.84136.61146.998164.34170.80182.519201.95208.97222.4610244.12250.84263.9411288.08295.99312.58注：1.模型(1)-(5)分别摘自Osterwald-Lenum(1992)表0,表1*,表1,表2*,表2。2.表示检验水平，N表示VAR模型中变量个数，r表示协整向量个数。案例分析1:关于中国GDP、宏观消费与基本建设投资的VEC模型分析。（293-302）用EViews估计VAR、VEC模型。1.建立VAR模型对任何一组有关系的经济变量都可以直接建立VAR模型。最大滞后期k的选择可以依据LR检验、赤池准则、Schwartz准则。建立VAR模型的EViews步骤是（1）点击Quick键，选EstimateVAR功能，得如下对话框：EViews4.1View的子菜单VARS-tabililyC6iidilionCheckolsofCharacteristicPolynomialEndogenousvariables:LhlGPLNCPLNIPExogenousvariables:CLagspecificalion:13Date:11/UM4Tine:16:0SRootModulus1.0154971.015497,4GE94e-0,6732G4i0.9193440.466948+0,673264108193J40805736U.HUb/JbW.220604-0.48976210.5362320.220604+0.489752i0.5362320.4401940.4401940.150344-0.229352102739870.190044+0,2293521027rB987Warning: At least one root outside the unit circle VAR does not satisfy thia stability condition.InveiseRootsofARChducteristicPolTm&midl451.005-O.o-，口由-A.0-1.5-1.5H.Ci-口?0.00.511.01.5问题：（1）非平稳经济变量之间可以建立VAR模型吗？若不存在协整关系不可以；若存在协整关系，在滞后项充分多的前提下可以建立VAR模型。这相当于每个方程都是AEG协整回归式。2.检验变量间是否存在协整关系。从工作文件中选中变量，打开数据组窗口，点击View键，选CointegrationTest功能，得如下对话框：其中有5种选择。协整空间无常数项、无时间趋势项；协整空间有常数项、无时间趋势项，数据空间无常数项；协整空间有常数项、无时间趋势项；协整空间有常数项、有时间趋势项，数据空间无时间趋势项；协整空间有常数项、有时间趋势项，数据空间有时间趋势项。上述5种情形总览。根据变量的实际情况作出选择JohansenComtegratj.onTes