第八专题时间序列模型PPT课件

1、一、结构VAR模型(SVAR) 内容安排： 第1页/共105页VAR模型并没有给出变量之间当期相关关系的确切形式，即在模型的右端不含有内生变量，而这些当期相关关系隐藏在随机误差项中，无法被观察到。模型中的误差项t是不可观测的，可以被看作是不可解释的随机扰动。结构VAR模型(Structural VAR，SVAR)，实际是指VAR模型的结构式，即在模型中包含变量之间的当期关系。 第2页/共105页含有两个变量(k=2)、滞后一阶(p=1)的VAR模型结构式可以表示为下式： 10121111212021211221ttttxtttttztxbb zxzuzbb xxzu称为一阶结构向量自回归模型（

2、SVAR（1）。 结构式经济模型，引入变量之间的作用与反馈作用，系数 b12表示变量zt的单位变化对变量xt的即时作用，21表示xt-1的单位变化对zt的滞后影响。冲击的交互影响体现了变量作用的双向和反馈关系。 第3页/共105页假设：【1】变量过程xt和zt均是平稳随机过程；【 2 】 随 机 误 差 ux t和 uz t 是 白 噪 声 序 列 ， 方差 ；【3】随机误差uxt 和uzt 之间不相关， 。 122zx0),cov(ztxtuu第4页/共105页 为导出VAR模型的简化式方程，将上述模型表示为矩阵形式： 该模型可以简单地表示为：ztxtttttuuzxbbzxbb112221

3、12112010211211tttuyyB1100第5页/共105页假设B0可逆，可导出简化式方程为 ：其中：tttuByBBy101110010ttyAA11020100100aaBA222112111101aaaaBAtttt2110uB第6页/共105页 从而可以看到，简化式扰动项t是结构式扰动项ut的线性组合，因此代表一种复合冲击。因为uxt 和uzt是不相关的白噪声序列，则可以断定上述1t和 2t 也是白噪声序列，并且均值和方差为：2211222122112222121111111)var(, 0)(, 0)(bbbbbbtsEEzxttst221122122211222122222

4、2111)var(, 0)(, 0)(bbbbbbtsEExzttst第7页/共105页同期的1t和 2 t之间的协方差为： 可以看出当b12 0或b21 0时，VAR模型简化式中的扰动项不再像结构式中那样不相关。当b12 = b21 = 0时，即变量之间没有即时影响，上述协方差为0，相当于对B0矩阵施加约束。 22112122122112212221212111)(),cov(bbbbbbbbEzxtttt第8页/共105页考虑k个变量的情形，p阶结构向量自回归模型SVAR(p)为： tptptttuyyyyB22110其中： 111212211120kkkkbbbbbbBpiikkikik

5、ikiiikiii,2, 1)()(2)(1)(2)(22)(21)(1)(12)(11第9页/共105页可以将上写成滞后算子形式： kttttELIuuuyB) (,)( 其中： ，B(L)是滞后算子L的 kk 的参数矩阵，B0 Ik。如果B0 是一个下三角矩阵，则SVAR模型称为递归的SVAR模型。 ppLLLLBB2210)(第10页/共105页假定结构式误差项(结构冲击) ut 的方差-协方差矩阵标准化为单位矩阵Ik。同样，如果矩阵多项式B(L)可逆，可以表示出SVAR的无穷阶的VMA()形式： 其中： ttL uDy)(1)()(LLBD2210)(LLLDDDD100 BD第11页

6、/共105页自Sims的研究开始，VAR模型开始取代了传统的联立方程模型，被证实为实用且有效的统计方法。VAR模型存在参数过多的问题，只有所含变量较少的VAR模型才可以通过OLS和极大似然估计得到“满意”的估计结果。 为了解决参数过多的问题，通过对参数空间施加约束条件从而减少所估计的参数，SVAR模型就是其中的一种。 第12页/共105页模型的识别性问题，即能否从简化式参数估计得到相应的结构式参数。对于k元p阶简化VAR模型： 利用极大似然方法，需要估计的参数个数为：tptpttyAyAy11222kkpk对于相应的k元p阶的SVAR模型： 来说，需要估计参数个数为： tptpttuyyyB1

7、1022kpk第13页/共105页 要想得到结构式模型惟一的估计参数，要求识别的阶条件和秩条件，即简化式的未知参数不比结构式的未知参数多。因此，如果不对结构式参数加以限制，将出现模型不可识别的问题。 对于k元p阶SVAR模型，需要对结构式施加的限制条件个数为：k(k-1)/2。这些约束条件可以是同期(短期)的，也可以是长期的。 第14页/共105页【1】短期约束 短期约束通常直接施加在矩阵D0 上，表示经济变量对结构冲击的同期响应，常见的可识别约束是简单的0约束排除方法。 通过Cholesky-分解建立递归形式的短期约束： Sims提出使D0 矩阵的上三角为0的约束方法，这是一个简单的对协方差

8、矩阵 的Cholesky-分解。Cholesky-分解的基本思想： 第15页/共105页 对于任意实对称正定矩阵 ，存在惟一一个主对角线元素为1的下三角形矩阵G和惟一一个主对角线元素为正的对角矩阵Q使得： （*） 利用这一矩阵G可以构造一个k维向量ut ，构造方法为 ut =G-1t，设 ： GGQ)(ttE第16页/共105页 则：QGGGQGGGGGuut111111)()(tttEE由于Q是对角矩阵，可得ut 的元素互不相关，其（j, j）元素是ujt 的方差。令Q 1/2表示其（j, j）元素为u jt 的标准差的矩阵。上式(*)可写为： PPGQGQ2/12/1 其中：P=GQ1/2

9、是一个下三角矩阵，上式被称为Cholesky (乔利斯基)分解。 第17页/共105页 Sims施加约束的基本过程是： 由于 是正定矩阵，所以可得到Cholesky因子P，即 。而且，当给定矩阵 时，Cholesky因子P是惟一确定的。 对于VAR模型两边都乘以P 1，得到：PPttLyA)(),(0ktVWN第18页/共105页tttLLuyAyAP)()(*1 其中： 。由于： ttPu1kt tttttEEEIPPPPPPuut111111)()()()(ut 是协方差为单位矩阵的白噪声向量，即： ),(kktVWNI0u第19页/共105页依据经济理论假设的短期约束 一般短期约束的施加

10、不必是下三角形式的。SVAR模型中的同期表示矩阵B0 是D0 的逆，即 ，因此可以通过对B0 施加限制条件实现短期约束。 例如：对于税收(y1t)、政府支出(y2t)和产出(y3t)的三变量SVAR模型来说 ，由于模型中包含3个内生变量，则k(k-1)/2= 3，因此需要对模型施加3个约束条件，才能识别出结构冲击。 100 BD第20页/共105页【2】 长期约束 关于长期约束的概念最早是由Blanchard 和 Quah在1989年提出的，是为了识别模型供给冲击对产出的长期影响。施加在结构VMA()模型的系数矩阵Di (i=1，2，)上的约束通常称为长期约束。最常见的长期约束的形式是对 的第

11、i行第j列元素施加约束，典型的是0约束形式，表示第i个变量对第j个变量的累积乘数影响为0。0iiD第21页/共105页SVAR模型根据其建模特点，主要分3种类型：K-型，C-型和AB-型，其中AB-型是最通常的类型，而K-型、C-型都可视为是AB-型的特殊形式。【1】 K-型 假定K是一个(kk)的可逆矩阵，K矩阵左如下形式的VAR模型： ( )ttLAy第22页/共105页其中： 扰动项 t 转变为正交扰动项ut (协方差矩阵是一个单位阵)，向量yt中各元素间的当期相关关系是由可逆矩阵K来决定的。假定知道t 的方差-协方差矩阵的真实形式： ( )ttLKAyKttuK ktE0u)(kttE

12、Iuu)(ttuK kttttIEE)()(uuKK第23页/共105页从而有： 这意味着对矩阵施加了k(k+1)/2个非线性的限制，K中剩下k(k1)/2个自由参数，还须给出k(k1)/2个短期约束条件。kIKK【2】 C-型型 假定C是一个(kk)的可逆矩阵，对于VAR模型： ( )ttLAy第24页/共105页如果满足下列条件： 则称上述模型为C-型SVAR模型。 在这一模型中，ut是相互独立的扰动，而t是独立正交的扰动项ut的线性组合。与K-型模型所不同的是：在这个模型中，内生变量之间没有同期关系，每个变量对正交扰动项的响应是通过矩阵C模拟的。 ttCu ktE0u)(kttEIuu)

13、(第25页/共105页 由 ，可以得到 。假定 的形式已知， 意味着对C矩阵施加了k(k+1)/2个非线性的限制性条件，C中剩下k(k1)/2个自由参数。如果C矩阵是下三角矩阵，则C矩阵就相当于Cholesky-分解的P矩阵。 )()(CuCuttttEECCCC【3】 AB-型 假定A、B是(kk)的可逆矩阵，A矩阵左乘上式形式的VAR模型，则得： 第26页/共105页如果A、B满足下列条件：则称上述模型为AB-型SVAR模型。 注意到AB-模型是最典型的SVAR模型，可以涵盖K-模型和C-模型。如果AB-模型中的A矩阵为单位矩阵，则AB-模型就转化为C-模型。如果AB-模型中的B矩阵为单位

14、矩阵，则此AB-模型为K-模型。 ttyLAAA)(ttBuA ktE0u)(kttEIuu)(第27页/共105页由：得到：如果 的形式已知，则 是对矩阵A、B的参数施加了k(k+1)/2个非线性限制条件，剩下2k2 k (k+1)/2个自由参数。 )()(BuBuAAttttEEBBAABBAA第28页/共105页在VAR窗口中选择：Procs /Estimate Structural Factorization 即可。 在E-views中SVAR模型采用AB-型： 其中et，ut是k维向量，et是可观测到的（或简化式的）残差，相当于前文的t，而ut 是不可观测的结构新息(结构式残差)。A

15、、B是待估计的k k矩阵。简化式残差et的协方差矩阵为 ttBuAe （五）在E-views中估计SVAR模型 第29页/共105页第30页/共105页用文本形式表示的短期约束： 对于更一般的约束，可用文本形式指定可识别的约束。在文本形式中，以一系列的方程表示关系：Aet = But ，并用特殊的记号识别et和ut向量中的每一个元素。A、B矩阵中被估计的元素必须是系数向量中被指定的元素。例如：对于有3个变量的VAR模型，约束A矩阵为主对角线是1的下三角矩阵，B矩阵是一对角矩阵。在这些约束条件下，Aet = But 的关系式可以写为下面的形式。 第31页/共105页tttttttttubeaea

16、eubeaeube333232131322212121111 为了以文本形式指定这些约束，从VAR对象窗口选择Procs/Estimate Structure Factorization，并单击Text按钮，在编辑框中，应键入下面的方程： 3*)6(2*)5(1*)4(32*)3(1*)2(21*) 1 (1ucececeuceceuce第32页/共105页特殊的关键符“e1”， “e2”， “e3”分别代表et向量中的第一、第二、第三个元素，而“ u1”， “ u2”， “ u3”分别代表ut 向量中的第一、第二、第三个元素。在这个例子中，A、B矩阵中的未知元素以系数向量c中的元素来代替。并

17、且对A、B矩阵的约束不必是下三角形式，可以依据具体的经济理论来建立约束。第33页/共105页一旦建立了模板矩阵，在VAR对象窗口的菜单中选择Procs/Estimate Structural Factorization，在SVAR Option对话框中，选择Matrix和Long-run Pattern按钮，并在相应的的编辑框中键入模版矩阵的名字。 用文本形式表示的长期约束： 在VAR对象窗口的菜单中选择Procs/Estimate Structural Factorization，并击活Text按钮，在编辑框中键入下面的形式： 第34页/共105页 lr2(u1)=0 zero LR res

18、ponse of 2nd variable to 1st shock 在撇号后面的内容是注释。这个约束以特殊的关键字“1r #”开始，数字代表受约束的响应变量；在圆括号内，必须指定脉冲关键字 u和扰动项序号，在其后紧跟等号和响应值（通常是0）。需注意：当需列出多个长期约束时，不要混淆短期与长期约束。第35页/共105页A、B矩阵的估计： 根据上述所描述的任何一种形式的可识别约束，单击SVAR Options对话框的OK按钮，可以估计A、B矩阵。为了使用脉冲响应和方差分解的结构选项，必须先估计这两个矩阵。 假定扰动项是多元正态的，EViews使用极大似然估计法估计A、B矩阵。使用不受限制的参数代

19、替受限制的参数计算似然值。对数似然值通过得分方法最大化，在这儿梯度和期望信息矩阵使用解析法计算。 第36页/共105页最优化控制(Optimization Control) 最优化过程控制的选项在SVAR Options对话框的Optimization Control栏下提供。可以指定初始值、迭代的最大数和收敛标准。 估计的输出 一旦估计收敛，E-views会在VAR对象窗口中显示估计的结果，包括：估计值、标准误差和被估计无约束参数的Z统计量及对数似然的最大值。基于被估计的信息矩阵的逆（Hessian的负的期望值）所估计的标准误差在最后的估计中计算。第37页/共105页 VAR模型中一个重要的

20、问题就是滞后阶数的确定。在选择滞后阶数p时，一方面想使滞后数足够大，以便能完整反映所构造模型的动态特征。但是另一方面，滞后数越大，需要估计的参数也就越多，模型的自由度就减少。所以通常进行选择时，需要综合考虑，既要有足够数目的滞后项，又要有足够数目的自由度。事实上，这是VAR模型的一个缺陷，在实际中常常会发现，将不得不限制滞后项的数目，使它少于反映模型动态特征性所应有的理想数目。二、滞后阶数p的确定 第38页/共105页确定滞后阶数的LR(似然比)检验 ：)(|ln|)ln(221kmTLRjj LR (Likelihood Ratio) 检验方法，从最大的滞后数开始，检验原假设：在滞后数为j时

21、，系数矩阵Aj的元素均为0；备择假设为：系数矩阵Aj中至少有一个元素显著不为0。2 (Wald)统计量如下： 其中m是可选择的其中一个方程中的参数个数：m =d+ kj，d是外生变量的个数，k是内生变量个数， 和 分别表示滞后阶数为(j 1)和 j 的VAR模型的残差协方差矩阵的估计。 1jj第39页/共105页 从最大滞后数开始，比较LR统计量和5%水平下的临界值，如果LR 时，拒绝原假设，表示统计量显著，此时表示增加滞后值能够显著增大极大似然的估计值；否则，接收原假设。每次减少一个滞后数，直到拒绝原假设。205. 0AIC信息准则和信息准则和SC准则准则 实际研究中，大家比较常用的方法还有

22、AIC信息准则和SC信息准则，其计算方法可由下式给出： 第40页/共105页其中在VAR模型中n = k(d + pk)是被估计的参数总数，k是内生变量个数，T是样本长度，d是外生变量的个数，p是滞后阶数，l是由下式确定的：TnTl22AICTTnTlln2SCln22ln12TTkl第41页/共105页在E-views软件中滞后阶数p的确定： 完成VAR模型的估计，选择View/Lag Structure/Lag Length Criteria，需要指定较大的滞后阶数，表中将显示出直至最大滞后数的各种信息标准（如果在VAR模型中没有外生变量，滞后从1开始，否则从0开始）。表中用“*”表示从每

23、一列标准中选的滞后数。在47列中，是在标准值最小的情况下所选的滞后数。第42页/共105页在E-views软件中关于残差的各种检验 【1】相关图相关图(Correlogram) VAR模型在指定的滞后数的条件下得到的残差的交叉相关图（样本自相关）。交叉相关图3种形式示：有两种表格形式，一种是以变量来显示（Tabulate by Variable），另一种是以滞后阶数来显示(Tabulate by Lag)。曲线图（Graph）显示交叉相关图矩阵形式。点线代表滞后的相关系数加减两倍的渐近标准误差的曲线图 。第43页/共105页【2】混合的自相关检验混合的自相关检验 计算与指定阶数所产生的残差序列

24、相关的多变量Box-Pierce/Ljung-Box Q统计量。 同时计算出Q统计量和调整后的Q统计量（即：小样本修正）。在原假设是滞后h期残差不存在序列相关的条件下，两个统计量都近似的服从自由度为k2 (h p)的2 统计量，其中p为VAR模型的滞后阶数。 第44页/共105页【3 3】自相关自相关LM检验检验 计算与直到指定阶数所产生的残差序列相关的多变量LM检验统计量。滞后h阶数的检验统计量是通过残差t 关于原始右侧回归量和滞后残差 t-h的辅助回归运算得到的，这里t-h 缺少的前h个值被赋予0。参考Johansen (1995)LM统计量的计算公式。在原假设是滞后h期没有序列相关的条件

25、下，LM统计量渐近地服从自由度为k2的2 分布。 第45页/共105页【4】 正态性检验 这是J-B残差正态检验在多变量情形下的扩展，这种检验主要是比较残差的第三、第四阶残差矩与来自正态分布的那些矩。 【5 5】White异方差检验异方差检验 回归检验是通过残差序列对每一个回归量及回归量交叉项乘积的回归来实现的，并检验回归的显著性。第46页/共105页 在实际应用中，由于VAR模型是一种非理论性的模型，因此在分析VAR模型时，往往不分析一个变量的变化对另一个变量的影响如何，而是分析当一个误差项发生变化，或者说模型受到某种冲击时对系统的动态影响，这种分析方法称为脉冲响应函数方法(impulse

26、response function，IRF)。三、脉冲响应与方差分解 第47页/共105页 用时间序列模型来分析影响关系的一种思路，是考虑扰动项的影响是如何传播到各变量的。两变量的VAR 模型来说明脉冲响应函数的基本思想。 （一）脉冲响应函数的基本思想 ttttttttttttzdzdxcxczzbzbxaxax222112211122112211其中，ai，bi，ci，di是参数， 是扰动项，假定是具有下面这样性质的白噪声向量： ),(21ttt第48页/共105页 假定上述系统从0期开始活动，且设x-1 = x-2 = z-1 = z-2= 0，又设于第0期给定扰动项10 =1，20 =0

27、，并且其后均为0，即 1t =2t =0(t ，)，称此为第0期给x以脉冲，下面讨论xt 与zt 的响应，t = 0时： stttsttttt,0)(,00)()var(,0)(2221第49页/共105页0, 100zx将其结果代入上式，当t = 1时：1111,czax再把此结果代入式，当t=2时 ：,112212cbaax112112cdcacz继续这样计算下去，设求得结果为称为由x的脉冲引起的x的响应函数。同样求得x的脉冲引起的z的响应函数： ,43210 xxxxx当t = 0时：,43210zzzzz第50页/共105页将上述讨论推广到多变量的VAR(p)模型： （二）VAR模型的

28、脉冲响应函数 ttpptLLLLCCIAAIy)()(22111 VMA()表达式的系数可按下面的方式给出，由于VAR的系数矩阵A和VMA的系数矩阵C必须满足下面关系： 第51页/共105页kpPLLLLICCIAAI)(2211kLLII221其中：1 = 2 = = 0。关于q的条件递归定义MA系数： PqPqqqCACACACACACAC2211211211kpqkpqpqpq0CIC，令，令若00第52页/共105页考虑VMA()的表达式： yt的第i个变量yit可以写成：其中k是变量个数。 tktLLCCIy)(221)(3)3(2)2(1)1()0(1jtijjtijjtijjti

29、jkjitccccy第53页/共105页 仅考虑两个变量的情形： ， q =1 , 2 , 3 ,， i , j = 1 , 2 现在假定在基期给 y1 一个单位的脉冲，即： )()(qijqcC2221)2(22)2(21)2(12)2(111211)1(22)1(21)1(12)1(1121)0(22)0(21)0(12)0(1121ttttttttccccccccccccyy第54页/共105页2 1 0 1 2 3 4 5 t其他, 00, 11tt， 210,02tt则由 y1的脉冲引起的y2的响应函数为：)2(2122)1(2121)0(2120,2, 1, 0cytcytcyt第

30、55页/共105页 因此，一般地，由yj的脉冲引起的yi的响应函数可以求出如下： ,)4()3()2()1()0(ijijijijijccccc 且由yj的脉冲引起的yi的累积(accumulate)响应函数可表示为 ：0)(qqijc第56页/共105页Cq的第i行、第j列元素还可以表示为 ： jtqtiqijyc,)(作为q的函数，它描述了在时期t，其他变量和早期变量不变的情况下yi,t+q对yjt的一个冲击的反应(对应于经济学中的乘数效应)，把它称作脉冲响应函数。也可以用矩阵的形式表示为： tqtqyC第57页/共105页 脉冲响应函数在E-views软件中的实现： 为了得到脉冲响应函数

31、，先建立一个VAR模型，然后在VAR工具栏中选择View/Impulse Response或者在工具栏选择Impulse，并得到下面的对话框，有两个菜单：Display 和 Impulse Definition。第58页/共105页Display菜单提供下列选项： 【1】显示形式（Display Format） 选择以图或表来显示结果。如果选择Combined Graphs 则Response Standard Error选项是灰色，不显示标准误差。而且应注意：输出表的格式是按响应变量的顺序显示，而不是按脉冲变量的顺序。 【2】显示信息（Display Information） 输入产生冲击的

32、变量（Impulses）和希望观察其脉冲响应的变量（Responses）。可以输入内生变量的名称，也可以输入变量的对应的序数。如果想显示累计的响应，则需要单击Accumulate Response选项。对于稳定的VAR模型，脉冲响应函数应趋向于0，且累计响应应趋向于某些非0常数。 第59页/共105页【3】脉冲响应标准差（Response Standard Error） 提供计算脉冲响应标准误差的选项。解析的或Monte Carlo标准误差对一些Impulse选项和误差修正模型（VEC）一般不一定有效。若选择了Monte Carlo，还需在下面的编辑框确定合适的迭代次数。 如果选择表的格式，被

33、估计的标准误差将在响应函数值下面的括号内显示。如果选择以多图来显示结果，曲线图将包括关于脉冲相应的正负（+/-）两个标准偏离带。在Combined Graphs中将不显示标准误差偏离带。 第60页/共105页Impulse Definition菜单提供了转换脉冲的选项： 【1】 Residual-One Unit 设置脉冲为残差的一个单位的冲击。这个选项忽略了VAR模型残差的单位度量和相关性，所以不需要转换矩阵的选择。这个选项所产生的响应函数是VAR模型相对应VMA()模型的系数。 【2】 Residual-One Std.Dev 设置脉冲为残差的一个标准偏差的冲击。这个选项忽略了VAR模型残

34、差的相关性。 第61页/共105页【3】 Cholesky 用残差协方差矩阵的Cholesky 因子的逆来正交化脉冲。这个选项为VAR模型的变量强加一个次序，并将所有影响变量的公共因素归结到在VAR模型中第一次出现的变量上。注意：如果改变变量的次序，将会明显地改变响应结果。可以在Cholesky Ordering 的编辑框中重新定义VAR模型中变量的次序。 第62页/共105页【5】结构分解（Structural Decomposition） 用结构因子分解矩阵估计的正交转换矩阵。如果没有先估计一个结构因子分解矩阵，或者没有对模型施加约束，这个选项不能用。 【4】广义脉冲（Gneralized

35、 Impluses） 描述Pesaran和Shin(1998)构建的不依赖于VAR模型中变量次序的正交的残差矩阵。应用按上面的Cholesky顺序计算的第j个变量的Cholesky因子得到第j个变量的扰动项的广义脉冲响应。 第63页/共105页【2】用户指定（User Specified） 这个选项允许用户定义脉冲。建立一个包含脉冲的矩阵（或向量），并在编辑框中输入矩阵的名字。如果VAR模型中有k个内生变量，则脉冲矩阵必须是k行和1列或k列的矩阵，每一列代表一个脉冲向量。 例：一个有k（= 3）个变量的VAR模型，希望同步对第一个变量有一个正的一个单位的冲击，给第二个变量一个负的一个单位的冲击

36、，建立一个31的脉冲矩阵SHOCK. 值分别为：1，1，0。在编辑框中键入矩阵的名字:SHOCK。 第64页/共105页 脉冲响应函数描述的是VAR模型中的一个内生变量的冲击给其他内生变量所带来的影响。而方差分解(variance decomposition)是通过分析每一个结构冲击对内生变量变化(通常用方差来度量)的贡献度，进一步评价不同结构冲击的重要性。因此，方差分解给出对VAR模型中的变量产生影响的每个随机扰动的相对重要性的信息。其基本思想如下所述。 （三）方差分解 第65页/共105页 脉冲响应函数是随着时间的推移，观察模型中的各变量对于冲击是如何反应的，然而对于只是要简单地说明变量间

37、的影响关系又稍稍过细了一些。因此，Sims（1980）依据VMA()表示，提出方差分解方法，定量地但是相当粗糙地把握变量间的影响关系。其思路如下：根据式： 可知各个括号中的内容是第j个扰动项j从无限过去到现在时点对yi影响的总和。求方差，假定j无序列相关，则：)(3)3(2)2(1)1()0(1jtijjtijjtijjtijkjitccccy第66页/共105页这是把第j个扰动项对第i个变量从无限过去到现在时点的影响，用方差加以评价的结果。此处还假定扰动项向量的协方差矩阵 是对角矩阵，则yi的方差是上述方差的k项简单和： jjqijqjtijjtijjtijccccE2)(022)2(1)1

38、()0()()()()var(2)(01jjqjjqkjitcy第67页/共105页 yi的方差可以分解成k种不相关的影响，因此为测定各个扰动项相对yi的方差有多大程度的贡献，定义了如下尺度： )()()var()()(2)(012)(02)(0jjqijqkjjjqijqitjjqijqijccycRVC即相对方差贡献率(relative variance contribution，RVC)是根据第j个变量基于冲击的方差对yi的方差的相对贡献度来观测第j个变量对第i个变量的影响。 第68页/共105页实际上，不可能用直到s = 的项和来评价。如果模型满足平稳性条件，则随着q的增大呈几何级数性

39、的衰减，所以只需取有限的s项。VAR(p)模型的前s期的预测误差是：可得近似的相对方差贡献率(RVC)： 001122110ICCCCCststtt)()()(2)(1012)(10jjqijsqkjjjqijsqijccsRVC第69页/共105页其中RVCji (s)具有如下的性质： 【1】【2】如果RVCji (s)大时，意味着第j个变量对第i个变量的影响大，相反地，RVCji (s)小时，可以认为第j个变量对第i个变量的影响小。 kjisRVCij， 21,1)(01)(1sRVCkjij第70页/共105页方差分解在Eviews软件中的实现 VAR的工具栏中选View/Varianc

40、e decomposition项。注意，因为非正交的因子分解所产生的分解不具有较好的性质，所以所选的因子分解仅限于正交的因子分解。 第71页/共105页 与脉冲响应函数一样，如果改变VAR模型中变量的顺序，基于Cholesky 因子的方差分解能有明显的改变。例如，排在第一个变量的第一期分解完全依赖于它自己的扰动项。 只有像在SVAR模型中那样估计一个结构因子分解矩阵时，基于结构正交化的因子分解才是有效的。注意：如果SVAR模型是恰好可识别的，那么预测的标准误差将等同于Cholesky因子分解的标准误差。对于过度识别的SVAR模型，为了保持更有效的性质，所预测的标准误差可能不同于Cholesky

41、因子分解的标准误差。 第72页/共105页四、AR系列动态模型内容安排： 第73页/共105页（一）ADLM模型ADLM（autoregressive distributed lag Model）称为自回归分布滞后模型。模型的特点： 相比于标准的协整检验，不论变量是否同为过程，或 同 为 过 程 ， 既 不 需 要 变 量 同 阶 单 整 ， 都 可 以 用ADLM模型来检验变量之间的长期关系。第74页/共105页模型结构模型结构典型的 模型的结构如下： 其中：L是滞后算子(lag operator，定义如下：1( ,)( ,)ktiiitttiL P yL q xwu12( ,.)kADLM

42、 p q qq212( ,)1.ppL PLLL 212( ,)1.iiqiiiiiqL qLLL 1ttL yy第75页/共105页建模方法建模方法 ADLM建模的方法包括两个阶段： 【1】建立与该ADLM模型相对应的误差修正模型（ECM），并计算出ECM模型中的F统计量，判断变量间是否存在长期稳定的关系。 【2】运用ARDL模型，估计变量之间长期关系的系数。第76页/共105页所谓ARIMA模型，是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对其滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA模型根据原序列是否平稳以及回归中所含部分的不同，具体又包括：移动平均过程

43、（MA）、自 回 归 过 程 （ A R ） 、 自 回 归 移 动 平 均 过 程（ARMA）以及ARIMA过程。（二）ARIMA模型第77页/共105页（1）MA过程q阶的移动平均（MA）过程可用下式表示： 其中u为常数项， 为白噪音过程 引入滞后算子L，原式可以写成: 或者 其中 1- 12- 2-.ttttqtqYueq eq eq eqititti=1Y = u+L+ttY = u+(L)2q12q(L )= 1+L +L.L t第78页/共105页MA（q）过程的特征 【1】 【2】 【3】自协方差 当kq时 0 当kq时，ACF(j)=0，此现象为截尾，是MA(q)过程的一个特征

44、。222j1j+12j+2qq-j12q1 j=0ACF(j)(.) (1.)0 jq 第86页/共105页（4）ARIMA模型的估计矩估计：利用样本自协方差函数和样本自相关函数，对模型的参数作估计。极大似然估计：包括无条件极大似然估计、条件极大似然估计、精确似然估计等方法。非线性估计：利用迭代搜索思想。 最小二乘估计：对于不包含MA部分的ARMA模型（即AR模型），可以利用普通最小二乘法对参数进行估计。第87页/共105页（5）ARMA模型的预测以平稳的AR(2)过程为例：其中 为零均值白噪音过程由模型的平稳性，我们有： t1t-12t-2tY = c+Y+Y+ utut+ 11t2t-1t

45、+ 1Y= c+Y +Y+ ut+ 21t+ 12tt+ 2Y= c+Y+Y + u第88页/共105页（三）ARCH模型Autoregressive Conditional Heteroscedasticity ,简称ARCH）模型,反映随机过程的一种特殊特性：即方差随时间变化而变化，且具有波动性。ARCH模型已广泛地应用于金融领域。第89页/共105页ARCH模型:由均值方程和条件方差方程给出： 表示t-1时刻所有可得信息的集合， 为条件方差用极大似然估计法对方程进行估计。tttyx222101122var( |).tttttp t phaaaa1tth第90页/共105页（四）GARCH

46、模型 一般的GARCH(p,q)模型如下表示 可用极大似然估计法估计。 GARCH(p,q)的推广tttyx2201111.ttpt ptqt qhaaahh第91页/共105页ARCH效应的识别ARCH LM Test：【1】ARCH效应的识别通常是对于残差项中是否存在自回归条件异方差现象的拉格朗日乘数检验（Lagrange multiplier test, Engle 1982）【2】ARCH LM Test 中的统计量F统计量Obs*R-squared 统计量第92页/共105页ARCH效应的识别：效应的识别： 原假设 ： 辅助回归方程： 其中e是残差，ARCH LM Test 对最小二

47、乘法，两阶段最小二乘法，非线性最小二乘法都适用。0H01.0q222201122.tttqt qteeeev第93页/共105页（G）ARCH模型的估计方法（最大似然估计模型的估计方法（最大似然估计法）：法）： 对于AR(1)- GARCH (1,1)模型： 可以通过最大化下述对数函数来估计模型参数：1201111ttttttyyhh 2111()11log(2 )log( )222TTttttttyyTLhh 第94页/共105页X-ARCH模型的类型：模型的类型：【1】GARCH-M模型【2】TARCH模型【3】EGARCH模型可以分析：【1】股市（汇市）收益波动【2】使用TARCH和EGARCH模型度量股市（汇市）收益波动非对称性【3】股市（汇市）波动溢出效应的研究第95页/共105页利用状态空间形式表示动态系统主要有两个优点： 第一，状态空间模型将不可观测的变量(状态变量)并入可观测模型并与其一起得到估计结果； 第二，状态空间模型是利用强有效的递归算法卡尔曼滤波来估计的