《随机时间序列分析模型课件.ppt》由会员分享，可在线阅读，更多相关《随机时间序列分析模型课件.ppt（81页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、随机时间序列分析模型第1页，此课件共81页哦经典计量经济学模型与时间序列模型经典计量经济学模型与时间序列模型确定性时间序列模型与随机性时间序列确定性时间序列模型与随机性时间序列模型模型第2页，此课件共81页哦一、时间序列模型的基本概念及其适用性一、时间序列模型的基本概念及其适用性第3页，此课件共81页哦（一）时间序列模型的基本概念（一）时间序列模型的基本概念 随随机机时时间间序序列列模模型型（time series modeling）是指仅用它的过去值及随机扰动项所建立起来的模型，其一般形式为 Xt=F(Xt-1,Xt-2,t)建立具体的时间序列模型，需解决如下三个问题建立具体的时间序列模型，
2、需解决如下三个问题：(1)模型的具体形式模型的具体形式 (2)时序变量的滞后期时序变量的滞后期 (3)随机扰动项的结构随机扰动项的结构 例如，取线性方程、一期滞后以及白噪声随机扰动项（t=t），模型将是一个1阶自回归过程阶自回归过程AR(1)：Xt=Xt-1+t这里，t特指一白噪声一白噪声。第4页，此课件共81页哦 一般的p阶自回归过程阶自回归过程AR(p)是 Xt=1Xt-1+2Xt-2+pXt-p+t (*)(1)如果随机扰动项是一个白噪声(t=t)，则称(*)式为一纯纯AR(p)过程（过程（pure AR(p)process），记为 Xt=1Xt-1+2Xt-2+pXt-p+t (2)如
3、果t不是一个白噪声，通常认为它是一个q阶的移移动平均（动平均（moving average）过程）过程MA(q)：t=t-1t-1-2t-2-qt-q 该式给出了一个纯纯MA(q)过过程程（pure MA(p)process）。第5页，此课件共81页哦 将纯AR(p)与纯MA(q)结合，得到一个一般的自回归移动平均自回归移动平均（autoregressive moving average）过程）过程ARMA（p,q）：Xt=1Xt-1+2Xt-2+pXt-p+t-1t-1-2t-2-qt-q 该式表明：该式表明：（1）一一个个随随机机时时间间序序列列可可以以通通过过一一个个自自回回归归移移动动
4、平平均均过过程程生生成成，即该序列可以由其自身的过去或滞后值以及随机扰动项来解释。（2）如如果果该该序序列列是是平平稳稳的的，即它的行为并不会随着时间的推移而变化，那么就可以通过该序列过去的行为来预测未来。那么就可以通过该序列过去的行为来预测未来。这也正是随机时间序列分析模型的优势所在。第6页，此课件共81页哦 经典回归模型的问题：经典回归模型的问题：迄迄今今为为止止，对一个时间序列Xt的变动进行解释或预测，是通过某个单方程回归模型或联立方程回归模型进行的，由于它们以因果关系为基础，且具有一定的模型结构，因此也常称为结结构构式式模模型型（structural model）。然然而而，如果Xt波
5、动的主要原因可能是我们无法解释的因素，如气候、消费者偏好的变化等，则利用结构式模型来解释Xt的变动就比较困难或不可能，因为要取得相应的量化数据，并建立令人满意的回归模型是很困难的。有有时时，即使能估计出一个较为满意的因果关系回归方程，但由于对某些解释变量未来值的预测本身就非常困难，甚至比预测被解释变量的未来值更困难，这时因果关系的回归模型及其预测技术就不适用了。（二）时间序列分析模型的适用性（二）时间序列分析模型的适用性第7页，此课件共81页哦 例例如如，时时间间序序列列过过去去是是否否有有明明显显的的增增长长趋趋势势，如果增长趋势在过去的行为中占主导地位，能否认为它也会在未来的行为里占主导地
6、位呢？或者时时间间序序列列显显示示出出循循环环周周期期性性行行为为，我们能否利用过去的这种行为来外推它的未来走向？随随机机时时间间序序列列分分析析模模型型，就就是是要要通通过过序序列列过过去去的的变变化化特特征征来预测未来的变化趋势来预测未来的变化趋势。使用时间序列分析模型的另一个原因在于使用时间序列分析模型的另一个原因在于：如果经济理论正确地阐释了现实经济结构，则这一结构可以写成类似于ARMA(p,q)式的时间序列分析模型的形式。在这些情况下，采用另一条预测途径在这些情况下，采用另一条预测途径：通过时间序列的历史数据，得出关于其过去行为的有关结论，进而对时间序列未来行为进行推断。第8页，此课
7、件共81页哦 例如，例如，对于如下最简单的宏观经济模型：这里，Ct、It、Yt分别表示消费、投资与国民收入。Ct与与Yt作作为为内内生生变变量量，它它们们的的运运动动是是由由作作为为外外生生变变量的投资量的投资It的运动及随机扰动项的运动及随机扰动项 t的变化决定的。的变化决定的。第9页，此课件共81页哦上述模型可作变形如下：两个方程等式右边除去第一项外的剩余部分可看成一个综合性的随机扰动项，其特征依赖于投资项It的行为。如如果果It是是一一个个白白噪噪声声，则消费序列Ct就成为一个1阶阶自自回回归归过过程程AR(1)，而收入序列Yt就成为一个(1,1)阶的自回归移动平均过程阶的自回归移动平均
8、过程ARMA(1,1)。第10页，此课件共81页哦二、随机时间序列模型的平稳性条件二、随机时间序列模型的平稳性条件第11页，此课件共81页哦 自回归移动平均模型（ARMA）是随机时间序列分析模型的普遍形式，自回归模型（AR）和移动平均模型（MA）是它的特殊情况。关于这几类模型的研究，是时时间间序序列列分分析析的的重重点点内内容容：主主要要包包括括模型的平稳性分析模型的平稳性分析、模型的识别模型的识别和和模型的估计模型的估计。（一）一）AR(p)AR(p)模型的平稳性条件模型的平稳性条件 随随机机时时间间序序列列模模型型的的平平稳稳性性，可可通通过过它它所所生生成成的的随随机机时时间间序序列列的
9、的平稳性来判断平稳性来判断。如如果果一个p阶自回归模型AR(p)生成的时间序列是平稳的，就说该AR(p)模型是平稳的，否则否则，就说该AR(p)模型是非平稳的。第12页，此课件共81页哦考虑p阶自回归模型AR(p)Xt=1Xt-1+2Xt-2+pXt-p+t (*)引入滞后算子（滞后算子（lag operator）L：LXt=Xt-1,L2Xt=Xt-2,LpXt=Xt-p(*)式变换为 (1-1L-2L2-pLp)Xt=t 记(L)=(1-1L-2L2-pLp),则称多项式方程 (z)=(1-1z-2z2-pzp)=0为AR(p)的特征方程特征方程(characteristic equati
10、on)(characteristic equation)。可以证明，可以证明，如果该特征方程的所有根在单位圆外如果该特征方程的所有根在单位圆外（根的模大于（根的模大于1 1），则），则AR(p)AR(p)模型是平稳的。模型是平稳的。第13页，此课件共81页哦 例例5.1 AR(1)模型的平稳性条件。对1阶自回归模型AR(1)方程两边平方再求数学期望，得到Xt的方差由于Xt仅与t相关，因此，E(Xt-1t)=0。如果该模型稳定，则有E(Xt2)=E(Xt-12)，从而上式可变换为：在稳定条件下，该方差是一非负的常数，从而有|1。第14页，此课件共81页哦 而AR(1)的特征方程的根为 z=1/A
11、R(1)稳定，即|1，意味着特征根大于1。例例5.2 AR(2)模型的平稳性。对AR(2)模型 方程两边同乘以Xt，再取期望得：第15页，此课件共81页哦又由于于是 同样地，由原式还可得到于是方差为 第16页，此课件共81页哦由平稳性的定义，该方差必须是一不变的正数，于是有 1+21,2-11,|2|1这就是AR(2)的平稳性条件的平稳性条件，或称为平稳域平稳域。它是一顶点分别为（-2,-1），（2,-1），（0,1）的三角形。第17页，此课件共81页哦对应的特征方程1-1z-2z2=0 的两个根z1、z2满足：z1z2=-1/2 ,z1+z2=-1/2 AR(2)模型解出1，2由AR(2)的
12、平稳性，|2|=1/|z1|z2|1，有于是|z2|1。由 2-1 1可推出同样的结果。第18页，此课件共81页哦 对高阶自回模型对高阶自回模型AR(p)来说来说，多数情况下没有必要直接计算其特征方程的特征根，但有一些有用的规则一些有用的规则可用来检验高阶自回归模型的稳定性可用来检验高阶自回归模型的稳定性：(1)AR(p)模型稳定的必要条件是模型稳定的必要条件是：1+2+p1 (2)(2)由于i(i=1,2,p)可正可负，AR(p)模型稳模型稳定的充分条件是：定的充分条件是：|1|+|2|+|p|1 第19页，此课件共81页哦 对于移动平均模型MR(q)：Xt=t-1t-1-2t-2-qt-q
13、 其中t是一个白噪声，于是 （二）（二）MA(q)模型的平稳性模型的平稳性 当滞后期大于q时，Xt的自协方差系数为0。因此:有限阶移动平均模型总是平稳的有限阶移动平均模型总是平稳的。第20页，此课件共81页哦 由于ARMA(p,q)模型是AR(p)模型与MA(q)模型的组合：Xt=1Xt-1+2Xt-2+pXt-p+t-1t-1-2t-2-qt-q （三）（三）ARMA(p,q)模型的平稳性模型的平稳性 而而MA(q)模型总是平稳的，因此模型总是平稳的，因此ARMA(p,q)模型的平稳性取决模型的平稳性取决于于AR(p)部分的平稳性。部分的平稳性。当当AR(p)部分平稳时，则该部分平稳时，则该
14、ARMA(p,q)模型是平稳的，否则，模型是平稳的，否则，不是平稳的。不是平稳的。第21页，此课件共81页哦 最后最后 （1 1）一个平稳的时间序列总可以找到生成它的平稳的随机过程或）一个平稳的时间序列总可以找到生成它的平稳的随机过程或模型；模型；（2 2）一个非平稳的随机时间序列通常可以通过差分的方法将它变换为）一个非平稳的随机时间序列通常可以通过差分的方法将它变换为平稳的，对差分后平稳的时间序列也可找出对应的平稳随机过程或模型。平稳的，对差分后平稳的时间序列也可找出对应的平稳随机过程或模型。因此，因此，如果我们将一个非平稳时间序列通过如果我们将一个非平稳时间序列通过d d次差分，将它变次差
15、分，将它变为平稳的，然后用一个平稳的为平稳的，然后用一个平稳的ARMA(p,q)ARMA(p,q)模型作为它的生成模型，则模型作为它的生成模型，则我们就说该原始时间序列是一个我们就说该原始时间序列是一个自回归单整移动平均自回归单整移动平均（autoregressive integrated moving averageautoregressive integrated moving average）时间序列，记为）时间序列，记为ARIMA(p,d,q)ARIMA(p,d,q)。例如，例如，一个一个ARIMA(2,1,2)ARIMA(2,1,2)时间序列在它成为平稳序列之前先得差分时间序列在它成
16、为平稳序列之前先得差分一次，然后用一个一次，然后用一个ARMA(2,2)ARMA(2,2)模型作为它的生成模型的。模型作为它的生成模型的。当然，当然，一个一个ARIMA(p,0,0)ARIMA(p,0,0)过程表示了一个纯过程表示了一个纯AR(p)AR(p)平稳过程；一平稳过程；一个个ARIMA(0,0,q)ARIMA(0,0,q)表示一个纯表示一个纯MA(q)MA(q)平稳过程。平稳过程。第22页，此课件共81页哦三、随机时间序列模型的识别三、随机时间序列模型的识别第23页，此课件共81页哦 所所谓谓随随机机时时间间序序列列模模型型的的识识别别，就是对于一个平稳的随机时间序列，找出生成它的合
17、适的随机过程或模型，即判断该时间序列是遵循一纯AR过程、还是遵循一纯MA过程或ARMA过程。所所使使用用的的工工具具主要是时间序列的自自相相关关函函数数（autocorrelation function，ACF）及偏偏自自相相关关函函数数（partial autocorrelation function，PACF）。第24页，此课件共81页哦（一）一）AR(p)AR(p)过程过程 1.1.自相关函数自相关函数ACFACF 1阶自回归模型阶自回归模型AR(1)Xt=Xt-1+t 的k阶滞后自协方差自协方差为：=1,2,因此，AR(1)模型的自相关函数自相关函数为=1,2,由由AR(1)的稳定性知
18、的稳定性知|1，因此，因此，k k时，呈指数形衰减，时，呈指数形衰减，直到零直到零。这种现象称为拖尾拖尾或称AR(1)有无穷记忆有无穷记忆（infinite memory）。注意注意，0时，呈振荡衰减状。第25页，此课件共81页哦 Xt=1Xt-1+2Xt-2+t该模型的方差0以及滞后1期与2期的自协方差1,2分别为阶自回归模型阶自回归模型AR(2)类似地,可写出一般的一般的k期滞后自协方差期滞后自协方差：(K=2,3,)于是,AR(2)的k 阶自相关函数阶自相关函数为：(K=2,3,)其中:1=1/(1-2),0=1如果如果AR(2)AR(2)稳定，则由稳定，则由 1 1+2 211知知|k
19、 k|衰减趋于零，呈拖尾状。衰减趋于零，呈拖尾状。至于衰减的形式，要看至于衰减的形式，要看AR(2)AR(2)特征根的实虚性，特征根的实虚性，若为实根，则呈若为实根，则呈单调或振荡型衰减，若为虚根，则呈正弦波型衰减。单调或振荡型衰减，若为虚根，则呈正弦波型衰减。第26页，此课件共81页哦一般地，p阶自回归模型阶自回归模型AR(p)Xt=1Xt-1+2Xt-2+pXt-p+tk期滞后协方差为:从而有自相关函数:可见，无论无论k k有多大，有多大，k k的计算均与其到的计算均与其到p p阶滞后的自相阶滞后的自相关函数有关关函数有关，因此呈拖尾状呈拖尾状。如果如果AR(p)AR(p)是稳定的，则是稳
20、定的，则|k k|递减且趋于零递减且趋于零。第27页，此课件共81页哦 其中：1/zi是AR(p)特征方程(z)=0的特征根，由AR(p)平稳的条件知，|zi|p，Xt与Xt-k间的偏自偏自相关系数相关系数为零。AR(p)的一个主要特征是的一个主要特征是:kp时，时，k*=Corr(Xt,Xt-k)=0 即即 k*在在p以后是截尾的。以后是截尾的。一随机时间序列的识别原则：一随机时间序列的识别原则：若若XtXt的的偏偏自自相相关关函函数数在在p p以以后后截截尾尾，即即kp时时，k*=0=0，而而它它的的自自相关函数相关函数 k是拖尾的，则此序列是自回归是拖尾的，则此序列是自回归AR(p)AR
21、(p)序列。序列。第30页，此课件共81页哦 在实际识别时，由于样本偏自相关函数rk*是总体偏自相关函数k*的一个估计，由于样本的随机性，当kp时，rk*不会全为0，而是在0的上下波动。但可以证明，当kp时，rk*服从如下渐近正态分布:rk*N(0,1/n)式中n表示样本容量。因此，如果计算的rk*满足 需指出的是需指出的是，我们就有95.5%的把握判断原时间序列在p之后截尾。第31页，此课件共81页哦 对MA(1)过程（二）MA(q)MA(q)过程过程 可容易地写出它的自协方差系数自协方差系数：于是，MA(1)过程的自相关函数自相关函数为：可见，当当k1时，时，k k0，即，即Xt与与Xt-
22、k不相关，不相关，MA(1)MA(1)自相关自相关函数是截尾的。函数是截尾的。第32页，此课件共81页哦 MA(1)过程可以等价地写成过程可以等价地写成 t t关于无穷序列关于无穷序列X Xt t，X Xt-1t-1，的线性组的线性组合的形式：合的形式：或（*）(*)是一个AR()过程，它的偏自相关函数非截尾但却趋于零，因此MA(1)MA(1)的偏自相关函数是非截尾但却趋于零的。的偏自相关函数是非截尾但却趋于零的。注意注意:(*)式只有当|1时才有意义，否则意味着距Xt越远的X值，对Xt的影响越大，显然不符合常理。因此，我们把把|1|q时，Xt与与Xt-k不相关，即存在截尾现象，因此，当当kq
23、时，时，k k=0是是MA(q)的一个特征的一个特征。于是：可可以以根根据据自自相相关关系系数数是是否否从从某某一一点点开开始始一一直直为为0 0来来判判断断MA(q)MA(q)模型的阶。模型的阶。第34页，此课件共81页哦 与MA(1)相仿，可以验证MA(q)过程的偏自相关函数是非截尾但趋于零的。MA(q)模型的识别规则：模型的识别规则：若随机序列的自相关函数截尾，即自若随机序列的自相关函数截尾，即自q q以后，以后，k k=0=0（kqkq）；而它的偏自相关函数是拖尾的，则此序列）；而它的偏自相关函数是拖尾的，则此序列是滑动平均是滑动平均MA(q)MA(q)序列。序列。同样需要注意的是同样
24、需要注意的是：在实际识别时，由于样本自相关函数rk是总体自相关函数k的一个估计，由于样本的随机性，当kq时，rk不会全为0，而是在0的上下波动。但可以证明，当kq时，rk服从如下渐近正态分布:rkN(0,1/n)式中n表示样本容量。因此，如果计算的如果计算的r rk k满足：满足：我们就有就有95.5%95.5%的把握判断原时间序列在的把握判断原时间序列在q q之后截尾之后截尾。第35页，此课件共81页哦 ARMA(p,q)的自相关函数的自相关函数，可以看作MA(q)的自相关函数和AR(p)的自相关函数的混合物。当当p=0时，它具有截尾性质时，它具有截尾性质；当当q=0时，它具有拖尾性质；时，
25、它具有拖尾性质；当当p、q都不为都不为0时，它具有拖尾性质时，它具有拖尾性质 从识别上看，通常：从识别上看，通常：ARMA(p，q)过程的偏自相关函数（过程的偏自相关函数（PACF）可能在可能在p阶滞后前有阶滞后前有几项明显的尖柱（几项明显的尖柱（spikes），但从），但从p阶滞后项开始逐渐趋向于零；阶滞后项开始逐渐趋向于零；而而它的自相关函数（它的自相关函数（ACF）则是在则是在q阶滞后前有几项明显的尖柱，阶滞后前有几项明显的尖柱，从从q阶滞后项开始逐渐趋向于零。阶滞后项开始逐渐趋向于零。（三）（三）ARMA(p,q)ARMA(p,q)过程过程 第36页，此课件共81页哦 表表5.1 AR
26、MA(p,q)模型的模型的ACF与与PACF理论模式理论模式 模型 ACF PACF 白噪声 0=kr 0*=kr AR(p)衰减趋于零（几何型或振荡型）P 阶后截尾：0*=kr，kp MA(q)q阶后截尾：，0=kr，kq 衰减趋于零（几何型或振荡型）ARMA(p,q)q阶后衰减趋于零（几何型或振荡型）p阶后衰减趋于零（几何型或振荡型）第37页，此课件共81页哦 图图4.2 ARMA(p,q)模型的模型的ACF与与PACF理论模式理论模式 ACF PACF 模型模型1：tttXXe+=-17.00.00.20.40.60.812345678ACF10.00.20.40.60.81234567
27、8PACF1第38页，此课件共81页哦第39页，此课件共81页哦第40页，此课件共81页哦四、随机时间序列模型的估计四、随机时间序列模型的估计第41页，此课件共81页哦 AR(p)、MA(q)、ARMA(p,q)模型的估计方法较多，大大体体上上分分为为3类：类：（1）最小二乘估计；）最小二乘估计；（2）矩估计；）矩估计；（3）利用自相关函数的直接估计）利用自相关函数的直接估计。下面有选择地加以介绍。结构阶数模型识别确定估计参数第42页，此课件共81页哦（一）（一）AR(p)AR(p)模型的模型的Yule WalkerYule Walker方程估计方程估计 在AR(p)模型的识别中，曾得到 利用
28、k=-k，得到如下方程组：此方程组被称为此方程组被称为Yule Walker方程组方程组。该方程组建立了该方程组建立了AR(p)AR(p)模型的模型参数模型的模型参数 1 1,2 2,p p与自相关函数与自相关函数 1 1,2 2,p p的关系，的关系，第43页，此课件共81页哦 利用实际时间序列提供的信息，利用实际时间序列提供的信息，首先首先求得自相关函数的估计值求得自相关函数的估计值 然后然后利用利用Yule Walker方程组，求解模型参数的估计值方程组，求解模型参数的估计值由于 于是 从而可得 2 2的估计值的估计值 在具体计算时，可用样本自相关函数rk替代。第44页，此课件共81页哦
29、（二）（二）MA(q)MA(q)模型的矩估计模型的矩估计 将MA(q)模型的自协方差函数中的各个量用估计量代替，得到：首先首先求得自协方差函数的估计值，(*)是一个包含(q+1)个待估参数(*)的非线性方程组，可以用直接法直接法或迭代法迭代法求解。常用的迭代方法有常用的迭代方法有线性迭代法线性迭代法和和Newton-Raphsan迭代法迭代法。第45页，此课件共81页哦 1.MA(1)1.MA(1)模型的直接算法模型的直接算法 对于MA(1)模型，（*）式相应地写成于是 或有于是有解 由于参数估计有两组解，可根据可逆性条件|1|1的MA(q)模型，一般用迭代算法估计参数：由（*）式得 第一步第
30、一步，给出的一组初值，比如代入（*）式，计算出第一次迭代值（*）第47页，此课件共81页哦 第二步第二步，将第一次迭代值代入（*）式，计算出第二次迭代值 按此反复迭代下去，直到第m步的迭代值与第m-1步的迭代值相差不大时（满足一定的精度），便停止迭代，并用第m步的迭代结果作为（*）的近似解。第48页，此课件共81页哦（三）（三）ARMA(p,q)ARMA(p,q)模型的矩估计模型的矩估计 在 ARMA(p,q)中 共 有(p+q+1)个 待 估 参 数1,2,p与1,2,q以及2，其估计量计算步骤及公式如下：第一步第一步，估计1,2,p 是总体自相关函数的估计值，可用样本自相关函数rk代替。第
31、49页，此课件共81页哦 第二步，第二步，改写模型，求1,2,q以及2的估计值 将模型 改写为：令 于是(*)可以写成：(*)构成一个MA模型。按照估计MA模型参数的方法，可以得到1,2,q以及2的估计值。第50页，此课件共81页哦（四）（四）AR(p)AR(p)的最小二乘估计的最小二乘估计 假设模型AR(p)的参数估计值已经得到，即有 残差的平方和为：(*)根据最小二乘原理，所要求的参数估计值是下列方程组的解：即 j=1,2,p (*)解该方程组，就可得到待估参数的估计值。第51页，此课件共81页哦 为了与AR(p)模型的Yule Walker方程估计进行比较，将(*)改写成：j=1,2,p
32、由自协方差函数的定义，并用自协方差函数的估计值 代入，上式表示的方程组即为：或 j=1,2,pj=1,2,p第52页，此课件共81页哦解该方程组，得到：即为参数的最小二乘估计。Yule Walker方程组的解比较发现，当n足够大时，二者是相似的。2的估计值为：第53页，此课件共81页哦 需要说明的是，需要说明的是，在上述模型的平稳性、识别与估计的讨论中，在上述模型的平稳性、识别与估计的讨论中，ARMA(p,q)模型中均未包含常数项。模型中均未包含常数项。如果包含常数项，该常数项并不影响模型的原有性质如果包含常数项，该常数项并不影响模型的原有性质，因为通过适当的变形，可将包含常数项的模型转换为不
33、含常数项的模型。下面以一般的ARMA(p,q)模型为例说明。对含有常数项的模型 方程两边同减/(1-1-p)，则可得到 其中第54页，此课件共81页哦五、模型的检验五、模型的检验第55页，此课件共81页哦 由于ARMA(p,q)模型的识别与估计是在假设随机扰动项是一白噪声的基础上进行的，因此，如果估计的模型确认正如果估计的模型确认正确的话，残差应代表一白噪声序列确的话，残差应代表一白噪声序列。如果通过所估计的模型计算的样本残差不代表一白噪声，如果通过所估计的模型计算的样本残差不代表一白噪声，则说明模型的识别与估计有误，需重新识别与估计。则说明模型的识别与估计有误，需重新识别与估计。在实际检验时
34、，主要检验残差序列是否存在自相关在实际检验时，主要检验残差序列是否存在自相关。（一）残差项的白噪声检验（一）残差项的白噪声检验 可用可用QLB的统计量进行的统计量进行 2检验检验：在给定显著性水平下，可计算不同滞后期的QLB值，通过与 2分布表中的相应临界值比较，来检验是否拒绝残差序列为白噪声的假设。若大于相应临界值，则应拒绝所估计的模型，需重新识别与估计。第56页，此课件共81页哦（二）（二）AICAIC与与SBCSBC模型选择标准模型选择标准 另外一个遇到的问题是，在实际识别ARMA(p,q)模型时，需多次反复偿试，有可能存在不止一组（p,q）值都能通过识别检验。显然，增加增加p与与q的阶
35、数，可增加拟合优度的阶数，可增加拟合优度，但却同时降低了自由度但却同时降低了自由度。因此，对可能的适当的模型，存在着模型的对可能的适当的模型，存在着模型的“简洁性简洁性”与模型与模型的拟合优度的权衡选择问题。的拟合优度的权衡选择问题。第57页，此课件共81页哦 其中，n为待估参数个数（p+q+可能存在的常数项），T为可使用的观测值，RSS为残差平方和（Residual sum of squares）。在选择可能的模型时，在选择可能的模型时，AIC与与SBC越小越好越小越好 显显然然，如如果果添添加加的的滞滞后后项项没没有有解解释释能能力力，则则对对RSSRSS值值的的减减小小没没有有多多大大帮
36、帮助助，却却增增加加待待估估参参数数的的个个数数，因因此此使使得得AICAIC或或SBCSBC的的值值增增加。加。需需注注意意的的是是：在不同模型间进行比较时，必须选取相同的时间段。常用的模型选择的判别标准有：常用的模型选择的判别标准有：赤池信息法赤池信息法（Akaike information criterion，简记为简记为AIC）与施瓦兹贝叶斯法施瓦兹贝叶斯法（Schwartz Bayesian criterion，简记为简记为SBC）：第58页，此课件共81页哦 由第一节知：中国支出法GDP是非平稳的，但它的一阶差分是平稳的，即支出法GDP是I(1)时间序列。可以对经过一阶差分后的GD
37、P建立适当的ARMA(p,q)模型。记GDP经一阶差分后的新序列为GDPD1，该新序列的样本自相关函数图与偏自相关函数图如下：例例4.3 中国支出法中国支出法GDP的的ARMA(p,q)模型估计。模型估计。第59页，此课件共81页哦 图形：图形：样本自相关函数图形呈正弦线型衰减波，而偏自相关函数图形则在滞后两期后迅速趋于0。因此可初步判断该序列满足可初步判断该序列满足2 2阶自回归过程阶自回归过程AR(2)AR(2)。自相关函数自相关函数与偏自相关函数偏自相关函数的函数值：函数值：相关函数具有明显的拖尾性；偏自相关函数值在k2以后，可认为：可认为：偏自相关函数是截尾的。再次验证了一阶差分后的偏
38、自相关函数是截尾的。再次验证了一阶差分后的GDPGDP满足满足AR(2)AR(2)随机过程。随机过程。表表5.2 中国中国GDP一阶差分序列的样本自相关函数与偏自相关函数一阶差分序列的样本自相关函数与偏自相关函数kkr*krkkr*krkkr*kr10.8590.8597-0.034-0.25213-0.361-0.08620.622-0.4418-0.1120.01214-0.3630.07630.378-0.0659-0.1750.0415-0.3080.04340.1910.06610-0.228-0.11716-0.216-0.02250.0870.07711-0.282-0.1921
39、7-0.128-0.04860.036-0.05112-0.32-0.0218-0.059-0.002第60页，此课件共81页哦设序列GDPD1的模型形式为 有如下Yule Walker 方程：解为：用用OLSOLS法回归的结果为：法回归的结果为：（7.91)(-3.60)r2=0.8469 R2=0.8385 DW=1.15第61页，此课件共81页哦 有时，在用回归法时，也可加入常数项有时，在用回归法时，也可加入常数项。本例中加入常数项的回归为：（1.99）（7.74）（-3.58）r2=0.8758 R2=0.8612 DW.=1.22 第62页，此课件共81页哦模型检验模型检验 下表列出
40、三模型的残差项的自相关系数及QLB检验值。模型1与模型3的残差项接近于一白噪声，但模型2存在4阶滞后相关问题，Q统计量的检验也得出模型2拒绝所有自相关系数为零的假设。因此:模型模型1 1与与3 3可作为描述中国支出法可作为描述中国支出法GDPGDP一阶差分序列的随机生成过程。一阶差分序列的随机生成过程。表表 5.3 模型残差项的自相关系数及模型残差项的自相关系数及Q检验值检验值 模型1 模型2 模型3 K Resid-ACF Q Resid-ACF Q Resid-ACF Q 1 0.382 3.3846 0.258 1.5377 0.257 1.5263 2 0.014 3.3893-0.1
41、39 2.0077-0.040 1.5646 3-0.132 3.8427-0.246 3.5677-0.059 1.6554 4-0.341 7.0391-0.529 11.267-0.328 4.6210 5-0.170 7.8910-0.300 13.908-0.151 5.2864 6 0.253 9.9097 0.271 16.207 0.345 9.0331 7 0.144 10.613 0.158 17.051 0.155 9.8458 8 0.057 10.730 0.116 17.541 0.076 10.059 9-0.019 10.745 0.097 17.914 0.0
42、11 10.064 10-0.146 11.685-0.036 17.969-0.123 10.728 11-0.233 14.329-0.136 18.878-0.230 13.319 12-0.049 14.461 0.064 19.104-0.012 13.328 第63页，此课件共81页哦用建立的用建立的AR(2)模型对中国支出法模型对中国支出法GDP进行外推预测。进行外推预测。模型模型1可作如下展开：于是，当已知t-1、t-2、t-3期的GDP时，就可对第t期的GDP作出外推预测。模型模型3的预测式与此相类似，只不过多出一项常数项。对对20012001年中国支出法年中国支出法GDPG
43、DP的预测结果（亿元）的预测结果（亿元）预测值预测值 实际值实际值 误差误差 模型模型1 95469 95933 -0.48%1 95469 95933 -0.48%模型模型3 97160 95933 1.28%3 97160 95933 1.28%第64页，此课件共81页哦 由由于于中国人均居民消费（CPC）与 人 均 国 内 生 产 总 值（GDPPC）这两时间序列是非平稳的，因此不宜直接建立它们的因果关系回归方程。但但 它它 们们 都都 是是I(2)I(2)时时 间间 序序 列列，因 此 可 以 建 立 它 们 的ARIMA(p,d,q)模型。下面只建立下面只建立中国人均居民消费（中国人
44、均居民消费（CPCCPC）的随机时间序列模型。的随机时间序列模型。中国人均居民消费（CPC）经过二次差分后的新序列记为CPCD2，其自相关函数、偏自相关函数及Q统计量的值列于下表：例例4.4 中国人均居民消费的中国人均居民消费的ARMA(p,q)ARMA(p,q)模型模型第65页，此课件共81页哦 在5%的显著性水平下，通过Q统计量容易验证该序列本身就接近于一白噪声，因此可考虑采用零阶考虑采用零阶MA(0)MA(0)模型模型:由于k=2时，|r2|=|-0.29|因此,也可考虑采用下面的MA模型：表表5.45.4 CPCD2CPCD2 序列的自相关函数、偏自相关函数与序列的自相关函数、偏自相关
45、函数与 Q Q 统计量值统计量值 k ACF PACF Q k ACF PACF Q 1 0.125 0.125 0.269 7 0.196 0.014 6.286 2-0.294-0.314 1.882 8-0.218-0.335 8.067 3-0.034 0.060 1.906 9-0.010 0.024 8.072 4-0.213-0.350 2.919 10 0.102-0.147 8.650 5-0.258-0.193 4.576 11-0.071 0.001 9.025 6 0.131 0.017 5.057 12 0.006-0.119 9.029 第66页，此课件共81页哦
46、当然，还可观察到自相关函数在滞后4、5、8时有大于0.2的函数值，因此，可考虑在模型中增加MA(4)、MA(5)、MA(8)。不同模型不同模型的回归结果列于表的回归结果列于表4.54.5。可以看出可以看出:在纯在纯MAMA模型中，模型模型中，模型4 4具有较好的性质，但由于具有较好的性质，但由于MA(5)MA(5)的的t t检验偏小，因此可选取模型检验偏小，因此可选取模型3 3。表表 5.55.5 中国居民人均消费水平的中国居民人均消费水平的ARMAARMA模型模型 模型 a MA(2)MA(4)MA(5)MA(8)AR(1)R2 SSR AIC 1 24.57 0 93137.4 8.94
47、2 32.4-0.89 0.42 53699.9 8.54 (3.62)(-7.43)3 14.07-0.72-1.71 0.7 28128.8 8.03 (8.75)(-3.07)(-5.08)4 11.73-1.09-1.99-1.3 0.82 17480.8 7.7 (17.81)(-3.38)(-4.61)(-1.58)5 11.79-1.07-1.91-1.25-0.34 0.81 17402.7 7.84 (14.93)(-3.10)(-2.56)(-1.42)(-0.15)6 14.95-0.66-1.27 -1.99 0.75 22924.2 7.97 (5.16)(-2.14
48、)(-1.77)(-1.29)7 214.25-2.53-2.45 -6.52 1.39 0.99 8943.7 7.06 (63.83)(-2.25)(-2.53)(-2.23)(98.26)第67页，此课件共81页哦 最后，给出通过模型最后，给出通过模型3 3的外推预测。的外推预测。模型3的展开式为：即 由于t表示预测期的随机扰动项，它未知，可假设为0，于是t期的预测式为:为模型3中滞后2期与滞后4期的相应残差项的估计值。第68页，此课件共81页哦 表5.6列出了采用模型3对中国居民人均居民消费水平的2期外推预测。为了对照，表中也同时列出了采用2.10的模型的预测结果。表表 5.6 5.6
49、 中国居民人均消费水平中国居民人均消费水平2 2 期外推预测比较（单位：元）期外推预测比较（单位：元）实际值 ARMA模型 因果关系模型 预测值 相对误差（%）预测值 相对误差（%）1997 2834 3048 7.6 2822-0.4 1998 2972 3407 14.6 2977 0.2 第69页，此课件共81页哦例例例例5.5 5.5 利用自相关和偏自相关系数利用自相关和偏自相关系数利用自相关和偏自相关系数利用自相关和偏自相关系数 识别识别识别识别ARMAARMA模型并检验序列相关性模型并检验序列相关性模型并检验序列相关性模型并检验序列相关性 考考察察美美国国1947年年第第1季季度度
50、1995年年第第1季季度度的的消消费费和和GDP之之间间的的关关系系，发发现现残残差差存存在在序序列列相相关关，并并且且例例4.4用用AR(3)模模型型修修正正了了残残差差的的序序列列相相关关性性。但但是是，并并没没有有说说明明是是通通过过怎怎样样的的方方法法来来判判断断残残差差服服从从一一个个AR(3)模模型型。这这个个例例子子将将借借助助Q-统统计计量量、自自相相关关系系数数和和偏偏自自相相关关系系数数图图，说说明明如如何何判判断断模模型型的的阶阶数数。回回归归方方程：程：第70页，此课件共81页哦第71页，此课件共81页哦图图5.3 5.3 原方程的残差图原方程的残差图 从从残残差差图图