时间序列模型diqizhangppt课件

1、X12季节调整方法 X12对X11方法进展改良：1 扩展贸易日和节假日影响的调理功能，添加季节、趋势循环和不规那么要素分解模型的选择功能；2 新的季节调整结果稳定性诊断功能；3 添加X12-ARIMA模型的建模和模型选择功能。 例子：1 2021年国庆房地产销售？2 农产品销售？3 月饼的销售？4 羊肉和狗肉销售？ 第七专题时间序列模型二本讲要点：一、季节调整、分解与平滑二、单整与协整三、VAR和误差修正模型1234一 X-N季节调整方法 1954年，美国商务部国势普查局BCDCC在美国国家经济研讨局NBER战前研讨的挪动平均比法The Ratio-Moving Average Method的

2、根底上，开发了关于季节调整的最初的电子计算机程序，开场大规模地对时间序列进展季节调整。以后，季节调整方法不断改良，每次改良都以X再加上序号表示。1960年，发表X-3方法，该方法和以前的程序相比，特异项的替代方法和季节要素的计算方法略有不同。1961年，发表X-10方法。该方法思索到了根据不规那么变动和季节变动的相对大小来选择计算季节要素的挪动平均项数。1965年，发表X-11方法，成为一种常用的季节调整方法。 课件讨论教师：周靖祥单位：湘潭大学商学院Email1:zjx2021sinaEmail2:zhoujingx126 CourseZJX202107一、经济时间序列的季节调整、分解与平滑

3、 内容安排：一 X-N季节调整方法 二X11方法三趋势分解 四指数平滑 中国的春节、国庆节等节假日对经济时间序列也会产生影响。例如，春节影响可以添加当周或前一周商品的零售额。在X12方法中，贸易日和节假日影响可以从不规那么要素中同时估计得到。在X12方法中，可以对不规那么要素建立ARIMAX模型，包括贸易日和节假日影响的回归变量，而且还可以指明奇特值的影响，并在估计其他回归影响的同时消除。留意：E-views中的节假日调整只针对美国，不能运用于其他国家中国：圣诞节和万圣节。2 节假日影响的调整利用E-views 软件对一个月度或季度时间序列进展季节调整的操作方法。在E-views任务环境中，翻

4、开一个月度或季度时间序列的任务文件，双击需进展数据处置的序列名，进入存放时间序列的任务表中，在序列窗口的工具栏中单击Proc按钮将显示菜单：3 季节调整相关操作 E-views软件1 Census X12方法 E-views进展季节调整时将执行以下步骤： 【1】给出一个被调整序列的阐明文件和数据文件； 【2】利用给定的信息执行X12程序； 【3】前往一个输出文件，将调整后的结果存在E-views任务文件中。 X12的E-views接口菜单只是一个简短的描画，E-views还提供一些菜单不能实现的接口功能，更普通的命令接口程序。 调用X12季节调整过程，在序列窗口选择Procs/Seasonal

5、 Adjustment / Census X12，翻开一个对话框： 什么是季节调整？季节性变动的缘由：气候、社会制度及风俗习惯定义：月度和季度数据都含有季节变动要素，由于季节要素的影响呵斥年度周期性变化，遮盖或混淆经济开展中其他客观变化规律，经济分析时必需去掉季节动摇的影响，将季节要素从原序列中剔除，这个过程就是“季节调整 Seasonal Adjustment季节性调整方法：1X-N季节调整方法【1】 X11方法X11 Method指定季节调整分解的方式：乘法；加法；伪加法此方式必需伴随ARIMA阐明；对数加法。留意乘法；伪加法和对数加法不允许有零和负数。【2】 季节滤波Seasonal F

6、ilter当估计季节因子时，允许选择季节挪动平均滤波，缺省是X12自动确定。 留意：假设序列短于20年，X12不允许指定315的季节滤波。2季节调整选择Seasonal Adjustment Option【4】 存调整后的分量序列名Component Series to save X12将被调整的序列名作为缺省列在Base name框中，可以改动序列名。多项选择钮中选择要保管的季节调整后分量序列，X12将加上相应的后缀存在任务文件中： 最终的季节调整后序列SA； 最终的季节因子SF； 最终的趋势循环序列TC； 最终的不规那么要素分量IR； 季节/贸易日因子D16； 假日/贸易日因子D18；【3

7、】趋势滤波Trend Filter (Henderson) 当估计趋势循环分量时，允许指定亨德松挪动平均的项数，可以输入大于1和小于等于101的奇数，缺省是由X12自动选择。3ARIMA选择ARIMA Option 点击ARIMA Option标签，可出现以下对话框： X12允许在季节调整前对被调整序列建立一个适宜的ARIMA模型。【1】 数据转换Data Transformation 在配备一个适宜的ARMA模型之前允许转换序列： * 缺省是不转换； * Auto选择是根据计算出来的AIC准那么自动确定是不做转换还是进展对数转换； * Logistic选择将序列 y 转换为 log(y/(1

8、-y)，序列的值被定义在0和1之间； * Box-Cox power选择要求提供一个参数 ，做以下转换：由每天经济活动的总和组成的月度时间序列受该月各周的影响，这种影响称为贸易日影响或周任务日影响。例如，对于零售业在每周的星期一至星期五的销售额比该周的星期六、星期日要少得多。北京周一周五商场不拥堵？因此，在某月假设多出的星期天数是一周的前五天，那么该月份销售额将较低；假设多出的星期天数是一周的星期六、星期日，那么该月份销售额将较高。又如，在流量序列中平均每天的影响将产生“月长度影响。由于在每年中二月份的长度是不一样的，所以这种影响不能够完全被季节要素接受。二月份残留的影响被称为润年影响28和2

9、9的差别。1 贸易日和节假日影响【2】ARIMA阐明(ARIMA Spec) 允许在2种不同的方法中选择ARIMA模型。 Specify in-line 选择 要求提供ARIMA模型阶数的阐明p d q(P D Q) p 非季节的AR阶数 d 非季节的差分阶数 q 非季节的MA阶数 P 季节AR阶数 D 季节差分阶数 Q 季节MA阶数 【3】回归因子选择Regressors 允许在ARIMA模型中指定一些外生回归因子，利用多项选择钮可选择常数项，或季节虚拟变量，事先定义的回归因子可以捕捉贸易日和节假日的影响。 可以在进展季节调整和利用ARIMA模型得到用于季节调整的向前/向后预测值之前，先去掉

10、确定性的影响例如节假日和贸易日影响。选择Adjustment Option能否进展这项调整？，确定在那一个步骤里调整：在ARIMA步骤，还是X-11步骤？4贸易日和节假日影响【1】Trading Day Effects消除贸易日影响有2种选择，依赖于序列是流量序列还是存量序列诸如存货。对于流量序列还有2种选择，是对周任务日影响进展调整还是对仅对周日-周末影响进展调整。存量序列仅对月度序列进展调整，需给出被观测序列的月天数。【2】Holiday effects 仅对流量序列做节假日调整。对每一个节日，必需提供一个数，是到这个节日之前影响的继续天数。 Easter 复活节 Labor 美国、加拿大

11、的劳工节，九月第一个星期一 Thanksgiving 感恩节在美国为11月第4个星期4；加拿大为10月第2个星期1 Christmas 圣诞节 留意这些节日只针对美国，不能运用于其他国家。提供的各种诊断： 季节要素的稳定性分析Stability Analysis of Seasonals Sliding spans 挪动间距 检验被调整序列在固定大小的挪动样本上的变化； Historical revisions 历史修正检验被调整序列添加一个新观测值，即添加一个样本时的变化。 其他诊断Other Diagnostics 还可以选择显示各种诊断输出。5 诊断Diagnostics二X11方法 X

12、-11法是美国商务部规范的季节调整方法乘法模型、加法模型，序列可被分解为趋势项与季节项的乘积或和。乘法模型只适用于序列值都为正的情形。 假设在季节调整对话框中选择X-11选项，调整后的序列及因子序列会被自动存入E-views任务文件中，在过程的结尾X-11简要的输出及错误信息也会在序列窗口中显示。 关于调整后的序列的名字。E-views在原序列名后加SA，但也可以改动调整后的序列名，这将被存储在任务文件中。 需求留意，季节调整的观测值的个数是有限制的。X-11只作用于含季节数据的序列，需求至少4整年的数据，最多能调整20年的月度数据及30年的季度数据。 TramoTime Series Reg

13、ression with ARIMA Noise, Missing Observation, and Outliers是对具有缺失观测值，ARIMA误差、几种外部影响的回归模型完成估计、预测和插值的程序。SeatsSignal Extraction in ARIMA Time Series)是基于ARIMA模型的将可观测时间序列分解为不可观测分量的程序。这两个程序是有Victor Gomez 和Agustin Maravall 开发的。中选择Pross/Seasonal Adjustment/Tramo Seats 时， E-views执行外部程序，将数据输给外部程序，然后将结果前往E-vie

14、ws。 三 Tramo/Seats方法 三趋势分解 如何将趋势和循环要素进展分解？测定长期趋势有多种方法，比较常用的方法有回归分析方法、挪动平均法、阶段平均法(phase average，PA方法)、HP滤波方法和频谱滤波方法frequency (band-pass) filer， BP滤波。1、Hodrick-PrescottHP滤波 在宏观经济分析中，经常需求分解序列组成成分中的长期趋势，Hodrick-Prescott滤波是被广泛运用的一种方法。该方法在Hodrick and Prescott(1980) 分析战后美国经济周期的论文中初次运用。设Yt是包含趋势成分和动摇成分的经济时间序列

15、，YtT是其中含有的趋势成分， YtC是其中含有的动摇成分。那么 计算HP滤波就是从Yt中将YtT 分别出来 。 普通地，时间序列Yt中的不可观测部分趋势YtT常被定义为下面最小化问题的解： 其中：c(L)是延迟算子多项式： 结合两式，HP滤波的问题就归结为使下面损失函数最小，即： 最小化问题用c(L)YtT2 来调整趋势的变化，并随着 的增大而增大。这里存在一个权衡问题，要在趋势要素对实践序列的跟踪程度和趋势光滑度之间作一个选择。 = 0 时，满足最小化问题的趋势等于序列Yt； 添加时，估计趋势中的变化总数相对于序列中的变化减少，即 越大，估计趋势越光滑； 趋于无穷大时，估计趋势将接近线性函

16、数。普通阅历地， 的取值如下： HP滤波处置不赖于经济周期峰和谷确实定。把经济周期看成宏观经济动摇对某些缓慢变动途径的偏离，这种途径在期间内单调地增长，所以称之为趋势。运用Hodrick-Prescott滤波来平滑序列，选择Procs/Hodrick Prescott Filter出现下面的HP滤波对话框：对平滑后的序列给一个变量名， E-views有默许值。然后给定平滑参数的值，年度数据取100，季度和月度数据分别取1600和14400。不允许填入非整数的数据。点击OK后，E-views与原序列一同显示处置后的序列。留意只需包括在当前任务文件样本区间内的数据才被处置，平滑后序列区间外的数据都

17、为NA。 例：利用HP滤波方法求潜在产出和产出缺口 设Yt为我国的季度GDP目的，利用季节调整方法将GDP中的季节要素和不规那么要素去掉，得到GDP_TC序列。潜在产出Y*,即趋势利用HP滤波计算出来的YtT来替代，GDP的循环要素Yt序列由式下式计算：GDP的循环要素YtC序列实践上就是围绕趋势线上下的动摇，称为GDP缺口序列。它是一个绝对量的产出缺口。也可以用相对量表示产出缺口，本例用Gapt来表示相对产出缺口，可由下式计算得到：四指数平滑 指数平滑是可调整预测的简一方法。当只需少数观测值时这种方法是有效的。与运用固定系数的回归预测模型不同，指数平滑法的预测用过去的预测误差进展调整。要用指

18、数平滑法预测，选择Procs/Exponential Smoothing 显示如下对话框：【1】平滑方法：在5种方法中选择一种方法。【2】平滑参数：可以指定平滑参数也可以让E-views估计。要估计参数，在填充区内输入字母e，E-views估计使误差平方和最小的参数值。假设估计参数值趋于1，这阐明序列趋于随机游走，最近的值对估计未来值最有用。要指定参数值，在填充区内输入参数值，一切参数值在0-1之间，假设他输入的参数值超出这一区间， E-views将会估计这个参数。 【3】平滑后的序列名：可以为平滑后的序列指定一个名字，E-views在原序列后加SM指定平滑后的序列名，也可以改动。【4】估计样

19、本：必需指定预测的样本区间不论能否选择估计参数。缺省值是当前任务文件的样本区间。E-views将从样本区间末尾开场计算预测值。【5】季节循环：可以改动每年季节数缺省值为每年12个月、4个季度。这个选项允许预测不规那么间距的数据，在空白处输入循环数。 1 单指数平滑一个参数 这种单指数平滑方法适用于序列值在一个常数均值上下随机动摇的情况，无趋势及季节要素。yt 平滑后的序列 计算式如下： , , t = 2, 3, , T其中: , 为平滑因子。 越小， 越平缓，反复迭代可得到： 缘由：由此可知为什么这种方法叫指数平滑，y 的预测值是 y 过去值的加权平均，而权数被定义为以时间为指数的方式。 2

20、 双指数平滑一个参数 这种方法是将单指数平滑进展两次运用一样的参数。适用于有线性趋势的序列。序列y的双指数平滑以递归方式定义为： 其中: 0 1, St 是单指数平滑后的序列，Dt 是双指数平滑序列。留意双指数平滑是阻尼因子为 0 1 的单指数平滑方法。双指数平滑的预测如下： 最后一个表达式阐明双指数平滑的预测有线性趋势，截距为 2ST DT ，斜率为 (ST DT )/(1 )， T 是估计样本的期末值。 3 Holt-Winters 无季节趋势两个参数 这种方法适用于具有线性时间趋势无季节变差的情形。这种方法与双指数平滑法一样以线性趋势无季节成分进展预测。双指数平滑法只用了一个参数，这种方

21、法用两个参数。yt 平滑后的序列 由下式给出： 其中: a 表示截距；b表示斜率，即趋势。这两个参数由如下递归式定义：其中: k 0 , , 在0-1之间，为阻尼因子。这是一种有两个参数的指数平滑法。 预测值计算如下： 这些预测值具有线性趋势，截距为 aT ，斜率为 bT ， T 是估计样本的期末值。 4 Holt-Winter加法模型三个参数 该方法适用于具有线性时间趋势和加法模型的季节变差。yt 平滑后的序列 由下式给出：其中：at 表示截距，bt 表示斜率， at + bt k 表示趋势，St 为加法模型的季节因子，s 表示季节周期长度，月度数据 s =12，季度数据 s = 4。需求用

22、简单的方法给出季节因子的第一年初值，以及截距和斜率的初值。这三个系数由下面的递归式定义：其中：k 0， 在01之间，为阻尼因子。预测值由下式计算：其中：ST+k-s用样本数据最后一年的季节因子，T 是估计样本的期末值。 5 Holt-winters乘法模型三个参数 这种方法适用于序列具有线性趋势和乘法季节变化。yt 的平滑序列 由下式给出：其中：at 表示截距，bt 表示斜率， at + bt k 表示趋势，St 为乘法模型的季节因子，s 表示季节周期长度，月度数据 s =12，季度数据 s = 4。需求用简单的方法给出季节因子的第一年的初值，以及截距和斜率的初值。这三个系数定义如下：其中：k

23、 0， 在01之间，为阻尼因子。预测值由下式计算 其中：ST+k-s 用样本数据最后一年的季节因子，T 是估计样本的期末值。例高铁梅，2006：指数平滑方法运用 本例利用指数平滑方法对我国上证收盘指数时间范围：1991年1月-2003年3月的月度时间序列 (sh_s) 进展拟合和预测。采用五种平滑模型对1991年1月-2002年9月的数据做指数平滑，并利用预测公式得到2002年10月-2003年3月半年的预测值。二、单整与协整内容安排： 一随机过程和平稳性原理 二平稳性检验的详细方法 三非平稳性数据的处置 四协整实际和协整检验 五误差修正模型一随机过程和平稳性原理普通称依赖于参数时间t的随机变

24、量集合 为随机过程。例，假设样本察看值y1，y2,yt是来自无穷随机变量序列y-2, y-1,y0 ,y1 ,y2 的一部分，那么这个无穷随机序列称为随机过程。 1 随机过程随机过程中有一特殊情况叫白噪音，其定义如下：假设随机过程服从的分布不随时间改动，且：对一切t 对一切t 那么，这一随机过程称为白噪声。 假设一个随机过程的均值和方差在时间过程上都是常数，并且在任何两时期的协方差值仅依赖于该两时期间的间隔或滞后，而不依赖于计算这个协方差的实践时间，就称它其为平稳的。2 平稳性原理平稳随机过程的性质：均值 对一切t方差 对一切t协方差 对一切t其中 即滞后k的协方差或自(身)协方差， 是 和

25、，也就是相隔k期的两值之间的协方差。 数据整理搜集和加工Unit Root 检验Vector Autoregression Theory (VAR)CointegrationError Correction ModelsImpulse Response Functions Variance DecompositionGranger Causality检验非平稳平稳将一个随机游走变量即非平稳数据对另一个随机游走变量进展回归能够导致荒唐的结果，传统的显著性检验将告知我们变量之间的关系是不存在的。有时候时间序列的高度相关仅仅是由于二者同时随时间有向上或向下变动的趋势，并没有真正的联络。这种情况就称为

26、“伪回归Spurious Regression。3 伪回归二平稳性检验的详细方法单位根检验的根本原理： David Dickey和Wayne Fuller的单位根检验unit root test即迪基富勒DF检验，是在对数据进展平稳性检验中比较经常用到的一种方法。1 单位根检验 DF检验的根本思想：从思索如下模型开场：* 其中 即前面提到的白噪音零均值、恒定方差、非自相关的随机误差项。由式*我们可以得到： * * * 依次将式(*)(*)、(*)代入相邻的上式，并整理，可得： &根据 值的不同，可以分三种情况思索：1假设 1，那么当T时， 0，即对序列的冲击将随着时间的推移其影响逐渐减弱，此时

27、序列是稳定的。 2假设 1，那么当T时， ，即对序列的冲击随着时间的推移其影响反而是逐渐增大的，很显然，此时序列是不稳定的。3 假设 =1，那么当T时， =1，即对序列的冲击随着时间的推移其影响是不变的，很显然，序列也是不稳定的。 对于式*，DF检验相当于对其系数的显著性检验，所建立的零假设是：H0 ： 假设回绝零假设，那么称Yt没有单位根，此时Yt是平稳的；假设不能回绝零假设，我们就说Yt具有单位根，此时Yt被称为随机游走序列random walk series是不稳定的。 方程*也可以表达成： 其中 = - ， 是一阶差分运算因子。此时的零假设变为：H0： =0。留意到假设不能回绝H0，那

28、么 = 是一个平稳序列，即 一阶差分后是一个平稳序列，此时我们称一阶单整过程integrated of order 1序列，记为I (1)。 I (1)过程在金融、经济时间序列数据中是最普遍的，而I (0)那么表示平稳时间序列。从实际与运用的角度，DF检验的检验模型有如下的三个： * 其中t是时间或趋势变量，在每一种方式中，建立的零假设都是：H0： 或H0： ，即存在一单位根。上述差别在于能否包含有常数截距和趋势项。假设误差项是自相关的，就把 * 修正如下：上式中添加了 的滞后项，建立在此根底上的DF检验又被称为增广的DF检验augmented Dickey-Fuller，简记ADF。ADF检

29、验统计量和DF统计量有同样的渐近分布，运用一样的临界值。首先，看如何判别检验模型能否应该包含常数项和时间趋势项。阅历做法是：调查数据图形。其次，如何判别滞后项数m。在实证中，常用的方法有两种： 2 ADF模型检验1渐进t检验。该种方法需求选择一个较大的m值，然后用t检验确定系数能否显著，假设是显著的，那么选择滞后项数为m；假设不显著，那么减少m直到对应的系数值是显著的2信息准那么。常用的信息准那么有AIC信息准那么、SC信息准那么，普通而言，选择给出了最小信息准那么值的m值普通是经过差分处置来消除数据的不平稳性。即对时间序列进展差分，然后对差分序列进展回归。对于金融数据做一阶差分后，即由总量数

30、据变为增长率，普通会平稳。数据处置的两难：但往往会丧失总量数据的长期信息，而这些信息对分析问题来说又是必要的。这就是通常我们所说的时间序列检验的两难问题。三非平稳性数据的处置 四协整实际和协整检验 思想：一些经济变量能够是非平稳的，但是它们的线性组合却有能够是平稳的，这种平稳的线性组合被称为协整方程，并且可被解释为变量之间的长期稳定的平衡关系。自变量和因变量之间存在协整关系，就是因变量能被自变量的线性组合所解释，两者之间存在平衡关系，因变量不能被自变量的所解释的部分构成一个残差序列，这个残差序列应该是平稳的。 格兰杰因果关系实际及格兰杰因果关系检验 格兰杰因果关系实际的思想是，x能否引起了y，

31、主要看如今的y 可以在多大的程度上被过去的x解释，即参与x的滞后值能否使解释程度提高。 格兰杰因果关系检验的原理 ：H0：1 特征根迹检验(trace检验) 由于r个最大特征根可得到r个协整向量，而对于其他k r个非协整组合来说，r+1，k应该为0，于是可得到原假设、备选假设为 相应的检验统计量为r称为特征根迹统计量。1 Johansen协整检验 2当 1 不显著时，接受H10 ，阐明只需1个协整向量，依次进展下去，直到接受Hr0，阐明存在r个协整向量。这r个协整向量就是对应于最大的r个特征根的经过正规化的特征向量。依次检验这一系列统计量的显著性： 1当 0 不显著时(即 0 值小于某一显著性

32、程度下的Johansen分布临界值)，接受H00 (r = 0)，阐明有k个单位根，0个协整向量(即不存在协整关系)。当 0 显著时(即 0 值大于某一显著性程度下的Johansen分布临界值)，回绝H00 ，那么阐明至少有一个协整向量，必需接着检验 1 的显著性。 根据右边假设检验，大于临界值回绝原假设。继续检验的过程可归纳为如下的序贯过程：1 临界值，回绝H10 ，阐明至少有2个协整向量； r 临界值，接受Hr0，阐明只需r个协整向量。 2最大特征值检验 对于Johansen协整检验，另外一个类似的检验方法是 ：检验统计量是基于最大特征值的，其方式为 其中 r 称为最大特征根统计量，简记为

33、-max统计量。 检验从下往上进展，首先检验0 ，假设 0 临界值，回绝H00 ，至少有1个协整向量。接受H00 (r = 0)，阐明最大特征根为0，无协整向量，否那么接受H01，至少有1个协整向量；假设 1 显著，回绝H10，接受至少有2个协整向量的备择假设H11；依次进展下去，直到接受Hr0，共有r个协整向量。 3协整方程的方式 与单变量时间序列能够出现均值非零、包含确定性趋势或随机趋势一样，协整方程也可以包含截距和确定性趋势。假设方程能够会出现如下情况Johansen，1995：【1】VAR模型 没有确定趋势，协整方程没有截距： 【2】VAR模型没有确定趋势，协整方程有截距项 ：【3】V

34、AR模型有确定性线性趋势，但协整方程只需截距： 【4】VAR模型和协整方程都有线性趋势，协整方程的线性趋势表示为 ： 【5】VAR模型有二次趋势，协整方程仅有线性趋势： 其中 是k ( kr )阶矩阵，它被称为 的正交互余矩阵(orthogonal complement) ，即 。 与 有关的项是协整关系的外部确定项，当确定项同时出如今协整关系的内部和外部时， 的分解不是独一可识别的。Johansen(1995)指出可将属于误差修正项内的那部分外生项正交地投影于 空间上，所以 是 的0空间，即 。 留意细节：【1】 Johansen协整检验的临界值对k =10的序列都是有效的。而且临界值依赖于

35、趋势假设，对于包含其他确定性回归量的模型能够是不适宜。例如，VAR模型中假设包含转移(变化)虚拟变量，能够使程度系列yt 产生一个不延续的线性趋势。【2】 迹统计量和最大特征值统计量的结论能够产生冲突。对这样的情况，建议检验估计得到的协整向量，并将选择建立在协整关系的解释才干上。 为了实现协整检验，从VAR对象或Group(组)对象的工具栏中选择View/Cointegration Test 即可。协整检验仅对知非平稳的序列有效，所以需求首先对VAR模型中每一个序列进展单位根检验。E-views软件中协整检验实现的实际根底是Johansen (1991, 1995a)协整实际。在Cointeg

36、ration Test Specification的对话框以下图中将提供关于检验的详细信息： 2 协整检验在E-views软件中的实现3协整检验的设定【1】 确定性趋势的阐明 序列也许会有非零均值，或与随机趋势一样有确定趋势。类似地，协整方程也能够会有截距和确定趋势，关于协整的LR检验统计量的渐近分布不再是通常的2分布，它的分布依赖于与确定趋势有关的假设。因此，为了完成这个检验，需求提供关于根本数据的趋势假设。E-views在Deterministic Trend assumption of test对话框中，对前面的5种能够方式提供了检验。假设不能确定用哪一个趋势假设，可以选择Summary

37、 of all 5 trend assumption第6个选择协助确定趋势假设的选择。这个选项在5种趋势假设的每一个下面都标明协整关系的个数，可以看到趋势假设检验结果的敏感性。 【2】 外生变量对话框还允许指定包含于VAR模型中的附加的外生变量Xt。常数和线性趋势不应被列在该编辑框中，由于它们在5个Trend Specification选项中得到了指定。假设确实包含外生变量，该当认识到EViews算出的临界值并没有思索这些变量。 【3】滞后区间加+1 该当用一对数字确定协整检验的滞后区间。需求留意的是：滞后设定是指在辅助回归中的一阶差分的滞后项，不是指原序列。例如，假设在编辑栏中键入“1 2，

38、协整检验用yt对yt-1，yt-2和其他指定的外生变量作回归，此时与原序列yt有关的最大的滞后阶数是3。对于一个滞后阶数为1的协整检验，在编辑框中应键入“0 0。 【4】协整检验结果的解释* 协整关系的数量 输出结果的第一部分给出了协整关系的数量，并以两种检验统计量的方式显示：第一种检验结果是所谓的迹统计量，列在第一个表格中；第二种检验结果是最大特征值统计量，列在第二个表格中。对于每一个检验结果：第一列显示了在原假设成立条件下的协整关系数；第二列是 矩阵按由大到小排序的特征值；第三列是迹检验统计量或最大特征值统计量；第四列是在5%显著性程度下的临界值；最后一列是根据MacKinnon-Haug

39、-Michelis (1999) 提出的临界值所得到的P值。* 协整关系 输出的第二部分给出协整关系 和调整参数 的估计。假设不强加一些恣意的正规化条件，协整向量 是不可识别的。在第一块中报告了基于正规化 其中S11在Johansen(1995a)中作出了定义的 和 的估计结果。留意：在Unrestricted Cointegrating Coefficients下 的输出结果：第一行是第一个协整向量，第二行是第二个协整向量，以此类推。其他的部分是在每一个能够的协整关系数(r =0，1，k -1)正规化后的估计输出结果。一个可选择的正规化方法是：在系统中，前r个变量作为其他k r个变量的函数。

40、近似的规范误差在可识别参数的圆括号内输出。五误差修正模型 两个经济变量存在协整关系，即两者之间有长期平衡关系。但是，在短期内也许会出现失衡。因此，可把协整方程中的误差项看作“平衡误差，并利用误差项把经济变量之间的短期行为和它的长期值联络起来。假设变量间存在协整关系，即阐明这些变量间存在着长期稳定的关系，而这种长期稳定的关系是在短期动态过程的不断调整下得以维持。 Engle和Granger将协整与误差修正模型结合起来，建立了向量误差修正模型。只需变量之间存在协整关系，可以由自回归分布滞后模型导出误差修正模型。而在VAR模型中的每个方程都是一个自回归分布滞后模型，因此，可以以为VEC模型是含有协整

41、约束的VAR模型，多运用于具有协整关系的非平稳时间序列建模。 1 向量误差修正模型(VEC) 其中的每一个方程都是一个误差修正模型。ecmt-1 = yt-1是误差修正项，反映变量之间的长期平衡关系，系数向量 反映变量之间的平衡关系偏离长期平衡形状时，将其调整到平衡形状的调整速度。一切作为解释变量的差分项的系数反映各变量的短期动摇对作为被解释变量的短期变化的影响，可以剔除其中统计不显著的滞后差分项。 其中每个方程的误差项 i (i =1，2，k) 都具有平稳性。一个协整体系由多种表示方式，用误差修正模型表示是当前处置这种问题的普遍方法，即： 假设yt 所包含的k个I (1)过程存在协整关系，不

42、包含外生变量可改写为 ： 思索一个两变量(y1，y2)的包含误差修正项、但没有滞后差分项的VEC模型。误差修正项是： 那么VEC模型为： 其中： ，写成一方程方式为： 其中，系数1，2 代表调整速度。在这个简单的模型中，等式右端独一的变量是误差修正项。在长期平衡中，这一项为0。然而，假设y1，y2 在上一期偏离了长期平衡，那么误差修正项非零，1和2会将其向平衡形状调整。 由于序列y1t，y2t的不同特征，模型可以指定成不同的方式：以下是简单的VEC模型，可以构造构造VEC模型，还可以思索VEC模型的Granger因果检验、脉冲呼应函数和方差分解。关于VAR模型和VEC模型更多的讨论，见Davi

43、dson和Mackinnon1993及汉密尔顿1999。 【1】 假设两个内生变量y1和y2不含趋势项，并且协整方程有截距，那么VEC模型有如下方式 ： 【2】 假设在序列中有线性趋势，那么VEC模型有如下方式： 【3】 协整方程中能够有趋势项 t，其方式为： 【4】 假设序列中存在着隐含的二次趋势项 t，等价于VEC模型的括号外也存在线性趋势项，其方式为：1如何估计VEC模型 由于VEC模型的表达式仅仅适用于协整序列，应先运转Johansen协整检验，并确定协整关系数，需求提供协整信息作为VEC对象定义的一部分。 2 VEC)模型在E-views软件中的实现向量在VAR对象设定框中，从VAR

44、 Type中选择Vector Error Correction项。在VAR Specification栏中，除特殊情况外应该提供与无约束的VAR模型一样的信息：【1】常数或线性趋势项不应包括在Exogenous Series的编辑框中，应定义在Cointegration栏中。【2】在VEC模型中滞后间隔的阐明指一阶差分的滞后。例如，滞后阐明“1 1将包括VEC模型右侧的变量的一阶差分项的滞后，即VEC模型是两阶滞后约束的VAR模型 。为估计没有一阶差分项的VEC模型，指定滞后的方式为：“0 0。 2在Cointegration栏中定义VEC模型常数和趋势时，必需从5个趋势假设阐明中选择一个，也

45、必需在适当的编辑框中填入协整关系的个数，应该是一个小于VEC模型中内生变量个数的正数。 3假设想强加约束于协整关系或(和)调整参数，用Restrictions栏以下图。留意：假设没在VAR Specification栏中单击Vector Error Correction项，这一栏将是灰色的。 含义：在有两个协整方程的情况，约束第三个变量外生于协整方程，两个协整方程第一个变量的系数均为1。单击OK按纽即可估计VEC模型，分两步完成：第一步，从Johansen所用的协整检验估计协整关系；第二步，用所估计的协整关系构造误差修正项，并估计包括误差修正项作为回归量的一阶差分方式的VAR模型。 VEC模型

46、估计的输出包括两部分。第一部分，假设不强加约束，EViews将会用系统默许的能可以识别一切的协整关系的正规化方法。系统默许的正规化表述为：将VEC模型中前r个变量作为剩余k r个变量的函数，其中r表示协整关系数，k是VEC模型中内生变量的个数。第二部分输出是在第一步之后以误差修正项作为回归量的一阶差分的VAR模型。误差修正项以CointEq1，CointEq2，表示方式输出。VEC模型输出结果的底部，有系统的两个对数似然值。第一个值标有Log Likelihood(d.f.adjusted)，其计算用自在度修正的残差协方差矩阵，这是无约束的VAR模型的对数似然值。标有Log Likelihoo

47、d的值是以没有修正自在度的残差协方差矩阵计算的，与协整检验所输出的值是可比较的。 4 VEC系数的获得 对于VEC模型，系数的估计保管在三个不同的二维数组中：A，B和C。A包含调整参数；B包含协整向量；C包含短期参数一阶差方项滞后的系数。 【1】A的第一个目的是VEC的方程序号，第二个目的是协整方程的序号。例如，A2，1表示：VEC的第二个方程中的第一个协整方程的调整系数。 【2】B的第一个目的是协整方程序号，第二个目的是协整方程的变量序号。例如，B2，1表示：第二个协整方程中第一个变量的系数。 【3】 C的第一个目的是VEC的方程序号，第二个目的是VEC中一阶差分回归量的变量序号。例如，C(

48、2 , 1)表示：VEC第二个方程中第一个一阶差分回归量的系数。 在VEC模型的名字后面加一个点号和系数元素，就可以获得这些系数，如： var01.a(2,1) var01.b(2,1) var01.c(2,1) 察看A , B和C的每一个元素和被估计系数的对应关系，从VAR的工具栏中选择 View/Representations 即可。三、 向量自回归和误差修正模型传统的经济计量方法是以经济实际为根底来描画变量关系的模型。但是，经济实际通常并缺乏以对变量之间的动态联络提供一个严密的阐明，而且内生变量既可以出如今方程的左端又可以出如今方程的右端使得估计和推断变得更加复杂。为理处理这些问题而出现

49、了一种用非构造性方法来建立各个变量之间关系的模型。 向量自回归(VAR)是基于数据的统计性质建立模型，VAR模型把系统中每一个内生变量作为系统中一切内生变量的滞后值的函数来构造模型，从而将单变量自回归模型推行到由多元时间序列变量组成的“向量自回归模型。VAR模型是处置多个相关经济目的的分析与预测最容易操作的模型之一，并且在一定的条件下，多元MA和ARMA模型也可转化成VAR模型。1 向量自回归实际 VAR(p) 模型的数学表达式是 其中：yt 是 k 维内生变量向量，Xt 是d 维外生变量向量，p是滞后阶数，样本个数为T 。kk维矩阵A1，Ap和kd维矩阵B是要被估计的系数矩阵。t是k维扰动向量，它们相互之间可以同期相关，但不与本人的滞后值相关及不与等式右边的变量相关，假设 是t的协方差矩阵，是一个(kk)的正定矩阵。式(9.1.1)可以用矩阵表示为 1VAR模型的普通表示： 即含有k个时间序列变量