《自回归移动平均模型.ppt》由会员分享，可在线阅读，更多相关《自回归移动平均模型.ppt（115页珍藏版）》请在三一办公上搜索。
1、第四章 时间序列计量经济学模型的理论与方法,第一节 随机时间序列的特征第二节 随机时间序列分析模型第三节 协整分析与误差修正模型第四节 向量自回归模型,4.1 随机时间序列的特征,一、随机时间序列模型简介二、趋势平稳与差分平稳三、时间序列平稳性的检验,一、随机时间序列模型简介,一个标有时间脚标的随机变量序列被称为时间序列(time series)。前提假设：时间序列是由某个随机过程(Stochastic process)生成的。即，假定序列X1,X2,XT 的每一个数值都是从一个概率分布中随机得到。当收集到一个时间序列数据集时，就得到该随机过程的一个可能结果或实现(realization)。,
2、假定某个时间序列是由某一随机过程生成，即假定时间序列Xt的每一个数值都是从一个概率分布中随机得到，如果时间序列Xt 满足：1）均值E(Xt)=是与时间t 无关的常数；2）方差Var(Xt)=2是与时间t 无关的常数；3）协方差Cov(Xt,Xt+k)=k 是只与时期间隔k 有关，与时间t 无关的常数；则称该随机时间序列是平稳的（stationary)，而该随机过程是一平稳随机过程（stationary stochastic process）。,1.时间序列的平稳性,经典计量模型的数学基础是极限法则，以独立随机抽样为样本，如果模型设定正确，模型随机误差项满足极限法则和由极限法则导出的基本假设，继
3、而进行的参数估计和统计推断是可靠的。以时间序列数据为样本，破坏了随机抽样的假定，则经典计量模型的数学基础能否被满足成为一个重要问题。对照极限法则和时间序列的平稳性条件研究发现，如果模型设定正确，并且所有时间序列是平稳的，时间序列的平稳性可以替代随机抽样假定，模型随机误差项仍然满足极限法则。,2.平稳性与经典回归,3.白噪声和随机游走,由定义知：白噪声序列是平稳的。,一个最简单的随机时间序列是一具有零均值同方差的独立同分布序列：Xt=t，t N(0,2),该序列常被称为是一个白噪声(white noise)。,另一个简单的随机时间列序被称为随机游走(random walk)，该序列由如下随机过程
4、生成：Xt=Xt-1+t 这里，t 是一个白噪声,t N(0,2)。,该序列 同均值，但方差不同：E(Xt)=E(Xt-1)X1=X0+1 X2=X1+2=X0+1+2 Xt=X0+1+2+t var(Xt)=t2，Xt的方差与时间 t 有关，而非常数，因此随机游走是非平稳序列。,4.齐次非平稳过程,如果一个时间序列是非平稳的，经过一次或多次差分后成为平稳序列，产生这样的非平稳序列的随机过程称为齐次随机过程。原序列转化为平稳序列所需的差分次数称为齐次的阶数。,对随机游走序列Xt取一阶差分(first difference)：,由于 t 是一个白噪声，则序列Xt 是平稳的。这提示我们如果一个时间
5、序列是非平稳的，常常可以通过取差分的方法形成平稳序列。,如果Yt 是一阶齐次非平稳过程，则序列：Wt=Yt Yt-1=Yt 就是平稳的。如果Yt 是二阶齐次非平稳过程，则序列：Wt=Yt Yt-1=2Yt 就是平稳的。,5.单整与非单整,如果一个时间序列经过一次差分变成平稳序列，也称原序列是1阶单整(integrated of 1)序列,记为I(1)过程。如果经过d 次差分后变成平稳序列,则称原序列是d 阶单整(integrated of d),记为I(d)。I(0)代表平稳时间序列。多次差分无法变为平稳的时间序列称为非单整的(non-integrated)。,随机时间序列Yt 的自相关函数(
6、autocorrelation function,ACF)：k=k/0 自相关函数是关于滞后期k的递减函数。对一个随机过程只有一个实现(样本),因此，只能计算样本自相关函数(Sample autocorrelation function)：,6.自相关函数、Q统计量,为了检验自相关函数的某个数值 k 是否为0，可以用Bartlett的研究结果：如果时间序列由白噪声生成，则对所有k 0，k N(0,1/T),为了检验所有k 0的自相关函数 k 都为0的联合假设，可以采用Box-Pierce的Q 统计量：,Q 统计量近似地服从自由度为k 的 分布。如果计算出Q 值大于显著性水平 下的临界值，就有1
7、-的把握拒绝所有k(k 0)同时为0的原假设。,1.确定性时间趋势 描述非平稳经济时间序列一般有两种方法，一种方法是包含一个确定性时间趋势：(*)其中 ut 是平稳序列；a+t 是线性趋势函数。这种过程也称为趋势平稳的，因为如果从式(*)中减去 a+t，结果是一个平稳过程。,二、趋势平稳与差分平稳随机过程,一般时间序列可能存在一个非线性函数形式的确定性时间趋势，例如可能存在多项式趋势：(*)t=1,2,T 同样可以除去这种确定性趋势，然后分析和预测去势后的时间序列。对于中长期预测而言，能准确地给出确定性时间趋势的形式很重要。如果 Yt 能够通过去势方法排除确定性趋势，转化为平稳序列，称为退势平
8、稳过程。,2.差分平稳过程 非平稳序列中有一类序列可以通过差分运算，得到具有平稳性的序列，考虑下式(*)也可写成：(*),其中 a 是常数,ut 是一个白噪声序列。式(*)的差分序列是含漂移 a 的随机游走，说明 yt 的差分序列yt是平稳序列。（*）式中L表示滞后算子。,实际上，以往讨论的回归方程的序列自相关问题暗含着残差序列是一个平稳序列。因为如果残差序列是一个非平稳序列，则说明因变量除了能被解释变量解释的部分以外，其余的部分变化仍然不规则，随着时间的变化有越来越大的偏离因变量均值的趋势，这样的模型是不能够用来预测未来信息的。,残差序列是一个非平稳序列的回归被称为伪回归，这样的一种回归有可
9、能拟合优度、显著性水平等指标都很好，但是由于残差序列是一个非平稳序列，说明了这种回归关系不能够真实的反映因变量和解释变量之间存在的均衡关系，而仅仅是一种数字上的巧合而已。伪回归的出现说明模型的设定出现了问题，有可能需要增加解释变量或者减少解释变量，抑或是把原方程进行差分，以使残差序列达到平稳。一个可行的办法是先把一个非平稳时间序列通过某种变换化成一个平稳序列。,一个平稳的时间序列在图形上往往表现出一种围绕其均值不断波动的过程；而非平稳序列则往往表现出在不同的时间段具有不同的均值（如持续上升或持续下降）。,1.平稳性检验的图示判断,三、时间序列的平稳性检验,平稳时间序列与非平稳时间序列图,单位根
10、检验（unit root test）是普遍应用的一类检验时间序列平稳性的方法，以ADF检验最为常用。(1)DF检验 我们已知道，随机游走序列Yt=Yt-1+t 是非平稳的，其中t 是白噪声。序列可看成是随机模型Yt=Yt-1+t 中参数=1时的情形。,2.平稳性的单位根检验,也就是说，对式 Yt=Yt-1+t（*）回归，如果确实发现=1，就说随机变量Yt有一个单位根。,（*）式可变成差分形式：Yt=(-1)Yt-1+t=Yt-1+t(*),检验（*）式是否存在单位根=1，也可通过（*）式判断是否有=0。,一般地:,检验一个时间序列Yt的平稳性，可通过检验带有截距项的一阶自回归模型 Yt=+Yt
11、-1+t（*）中的参数是否小于1。,或者：检验其等价变形式 Yt=+Yt-1+t（*）中的参数是否小于0。,（*）式中的参数 1或=1时，时间序列是非平稳的；对应于（*）式，则是 0或=0。,针对（*）式 Yt=+Yt-1+t 零假设 H0：=0，即原序列存在单位根；备择假设 H1：0，即原序列是平稳的；,上述检验可通过OLS法下的t 检验完成。Dicky和Fuller于1976年提出了这一情形下t 统计量服从的分布（这时的t 统计量称为 统计量），即DF分布（见下表）。,DF分布临界值表,通过OLS法估计 Yt=+Yt-1+t 计算t 统计量的值，与DF分布表中给定显著性水平下的临界值比较：
12、如果：t 临界值（左尾单侧检验），则拒绝原假设H0：=0，认为时间序列不存在单位根，是平稳的。,DF检验的问题：在上述使用 Yt=+Yt-1+t 对时间序列进行平稳性检验中，实际上假定时间序列是由一阶自回归过程AR(1)生成的，并且随机误差项是白噪声。为了保证DF检验中随机误差项的白噪声特性，Dicky和Fuller对DF检验进行了扩充，形成了ADF（Augment Dickey-Fuller）检验。,(2)ADF检验,ADF检验是通过以下3个模型完成的：,3个模型检验的原假设都是：H0：=0，即存在一单位根，备择假设：H1:0。,同时估计出上述3个模型的适当形式，然后通过ADF临界值表检验零
13、假设H0：=0。1）只要其中有一个模型的检验结果拒绝了零假设，就可以认为时间序列是平稳的；2）当3个模型的检验结果都不能拒绝零假设时，则认为时间序列是非平稳的。检验原理与DF检验相同。Dicky和Fuller推导了3个模型所使用的ADF分布临界值表。ADF检验也可判断时间序列的单整阶数。,ADF检验过程：,例1：检验19782000年间中国支出法GDP时间序列的平稳性及单整性。,1）经过偿试，模型3取了2阶滞后：,通过拉格朗日乘数检验对随机误差项的自相关性进行检验：LM(1)=0.92，LM(2)=4.16，小于5%显著性水平下自由度分别为1与2的2分布的临界值，可见不存在自相关性。,从 看，
14、t临界值(查ADF分布表)，不能拒绝存在单位根的零假设。,2）经试验，模型2中滞后项取2阶：,LM检验表明模型残差不存在自相关性。从GDPt-1的参数值看，其 t 统计量为正值，大于临界值(查ADF分布表)，不能拒绝存在单位根的零假设。,3）经试验，模型1中滞后项取2阶：,LM检验表明模型残差项不存在自相关性，因此模型的设定是正确的。从GDPt-1的参数值看，其t统计量为正值，大于临界值(查ADF分布表)，不能拒绝存在单位根的零假设。结论：根据ADF检验结果，可断定中国支出法核算的GDP时间序列是非平稳的。,支出法GDP时间序列的平稳性ADF检验模型3结果：,Eviews中，GDP平稳性ADF
15、检验结果：,Eviews中，GDP平稳性ADF检验结果（续）：,4）中国支出法GDP的单整性。,经过试算，发现中国支出法GDP是1阶单整的，适当的检验模型为：,支出法GDP时序一阶差分后的平稳性ADF检验模型3结果：,结论：根据ADF检验结果，可断定中国支出法核算的GDP的一阶差分序列是平稳的，即I(1)。,Eviews中，GDP序列ADF检验模型3的检验结果：,例2：检验关于人均居民消费与人均国内生产总值这两时间序列的平稳性及单整性。,1)对中国人均国内生产总值GDPP来说，经过偿试，三个模型的适当形式分别为：,模型3：,ADF检验过程：,模型2：,模型1：,3个模型中参数的估计值的t统计量
16、均大于各自的临界值，因此不能拒绝存在单位根的零假设。结论:人均国内生产总值(GDPP)是非平稳的。经过进一步检验发现，人均国内生产总值(GDPP)和人均居民消费(CONSPP)都是二阶单整序列，I(2),Eviews中GDPP序列ADF检验给出的模型3的检验结果：,2）对于人均居民消费CONSP时间序列来说，3个模型的适当形式为：,模型3：,模型2：,3个模型中参数CONSPt-1的t统计量的值均比ADF临界值表中各自的临界值大，不能拒绝该时间序列存在单位根的假设。结论：可判断人均居民消费序列CONSP是非平稳的。,模型1：,Eviews中CONSP序列ADF检验给出的模型3的检验结果：,Ev
17、iews中CONSP序列ADF检验给出的模型2的检验结果：,Eviews中CONSP序列ADF检验给出的模型1的检验结果：,中国人均居民消费与人均国内生产总值的单整性：,经过试算，发现中国人均国内生产总值GDPP是2阶单整的，适当的检验模型为：,CONSP也是2阶单整的，适当的检验模型为：,4.2 随机时间序列分析模型,一、模型的一般形式及其适用性二、模型的平稳性条件三、模型的识别四、模型的参数估计五、模型的检验,随机时间序列模型（Time Series Modeling）一般形式为：Xt=F(Xt-1,Xt-2,t)建立具体的时间序列模型的三个问题：(1)模型的具体形式(2)时序变量的滞后期
18、(3)随机扰动项的结构,一、随机时间序列模型的一般形式及适用性,例如，取线性方程、一期滞后以及白噪声随机扰动项（t=t），模型将是一个1阶自回归过程AR(1)：Xt=Xt-1+t(t 特指白噪声),一般的，p阶自回归过程AR(p)为：Xt=1Xt-1+2Xt-2+pXt-p+t(*),(1)如果随机扰动项是一个白噪声(t=t)，则称(*)式为一纯AR(p)过程(pure AR(p)process)。,(2)如果t不是一个白噪声，通常认为它是一个q阶的移动平均(moving average)过程MA(q)：t=t 1t-1 2t-2 qt-q 该式给出了一个纯MA(q)过程(pure MA(q)
19、process)。,一般的p阶自回归过程AR(p)是：Xt=1Xt-1+2Xt-2+pXt-p+t(*),将纯AR(p)与纯MA(q)结合，得到一个一般的自回归移动平均(autoregressive moving average)过程ARMA(p,q)：,Xt=1Xt-1+2Xt-2+pXt-p+t-1t-1-2t-2-qt-q,ARMA(p,q)：,该式表明：（1）一个随机时间序列可以通过一个自回归移动平均过程生成，即该序列可以由其自身的过去或滞后值以及随机扰动项来解释。（2）如果该序列是平稳的，即它的行为并不会随着时间的推移而变化，那么我们就可以通过该序列过去的行为来预测未来。,Xt=1X
20、t-1+2Xt-2+pXt-p+t-1t-1-2t-2-qt-q,经典回归模型的问题：（1）经典的计量经济学模型是以因果关系为基础，且具有一定的模型结构，因此也常称为结构式模型（structural model）。（2）然而，如果Xt波动的主要原因可能是我们无法解释的因素，则利用结构式模型来解释Xt的变动就比较困难或不可能。,时间序列分析模型的适用性,在这些情况下，采用另一条预测途径：通过时间序列的历史数据,得出关于其过去行为的有关结论，进而对时间序列未来行为进行推断。随机时间序列分析模型，就是要通过序列过去的变化特征来预测未来的变化趋势。,1.AR(p)模型的平稳性条件,如果一个p阶自回归模
21、型AR(p)生成的时间序列是平稳的，就说该AR(p)模型是平稳的，否则，就说该AR(p)模型是非平稳的。,二、随机时间序列模型的平稳性条件,考虑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 equation)。可以证明，如果该特征方程的所有根在单位圆外（根的模大于1），
22、则AR(p)模型是平稳的。,例，AR(1)模型的平稳性条件,对1阶自回归模型AR(1):,由于Xt 仅与t 相关,因此,E(Xt-1t)=0。如果该模型平稳,则有E(Xt2)=E(Xt-12),从而上式可变换为：,在平稳条件下,该方差是一非负的常数,从而有|1,AR(1)的特征方程：,的根为 z=1/AR(1)稳定,即|1,意味着特征根大于1,根的模大于1。,对高阶自回归模型AR(p)：,(1)AR(p)模型稳定的必要条件是：1+2+p1(2)由于i(i=1,2,p)可正可负,AR(p)模型稳定的充分条件是：|1|+|2|+|p|1,对于移动平均模型MA(q)：Xt=t-1t-1-2t-2-q
23、t-q 其中，t 是一个白噪声，于是,2.MA(q)模型的平稳性,当滞后期大于q时，Xt 的自协方差系数为0。因此：有限阶移动平均模型总是平稳的。,q+1=cov(Xt,Xt-q-1)=E(XtXt-q-1)=E(t-1t-1-2t-2-qt-q)*(t-q-1-1t-q-2-2t-q-3-qt-2q-1)=0,由于ARMA(p,q)模型是AR(p)模型与MA(q)模型的组合：Xt=1Xt-1+2Xt-2+pXt-p+t-1t-1-2t-2-qt-q,3.ARMA(p,q)模型的平稳性,而MA(q)模型总是平稳的，因此ARMA(p,q)模型的平稳性取决于AR(p)部分的平稳性。当AR(p)部分
24、平稳时，则该ARMA(p,q)模型是平稳的，否则，不是平稳的。,（1）一个平稳的时间序列总可以找到生成它的平稳的随机过程或模型；（2）一个非平稳的随机时间序列通常可以通过差分的方法将它变换为平稳的，对差分后平稳的时间序列也可找出对应的平稳随机过程或模型。如果将一个非平稳时间序列通过d次差分,将它变为平稳的，然后用一个平稳的ARMA(p,q)模型作为它的生成模型，则我们就说该原始时间序列是一个自回归单整移动平均（autoregressive integrated moving average）时间序列，记为ARIMA(p,d,q)。,4.ARIMA(p,d,q)模型,随机时间序列模型的识别，就是
25、对于一个平稳的随机时间序列，找出生成它的合适的随机过程或模型，即判断该时间序列是遵循一纯AR过程、还是遵循一纯MA过程或ARMA过程。所使用的工具主要是时间序列的自相关函数（autocorrelation function，ACF）及偏自相关函数（partial autocorrelation function，PACF）。,三、随机时间序列模型的识别,1.AR(p)过程,(1)自相关函数ACF 1阶自回归模型AR(1)：Xt=Xt-1+t 的k阶滞后自协方差为：,k=1,2,AR(1)模型的自相关函数为：,k=1,2,由AR(1)的稳定性知|1,因此,k时,ACF呈指数形衰减,直到零。这种现
26、象称为拖尾或称AR(1)有无穷记忆（infinite memory）。注意，0时，呈振荡衰减状。,一般地，p阶自回归模型AR(p):Xt=1Xt-1+2Xt-2+pXt-p+t,k期滞后协方差为:,从而有自相关函数:,无论k有多大，k 的计算均与其1到p阶滞后的自相关函数有关，因此呈拖尾状。如果AR(p)是稳定的,则|k|递减且趋于零。,其中：zi=1/i是AR(p)特征方程(z)=0的特征根，由AR(p)平稳的条件知，|i|1;因此，当zi均为实数根时，k呈几何型衰减（单调或振荡）；当存在虚数根时，则一对共扼复根构成通解中的一个阻尼正弦波项，k呈正弦波衰减。,事实上，自相关函数,是一p阶差分
27、方程，其通解为,(2)偏自相关函数,Xt与Xt-k间的偏自相关函数(partial autocorrelation,简记为PACF)是消除了中间变量Xt-1，Xt-k+1 带来的间接相关后的直接相关性，它是在已知序列值Xt-1，Xt-k+1的条件下，Xt 与Xt-k 间关系的度量。,Xt 与Xt-k 之间的条件相关性。其相关程度用偏自相关系数k,k 度量。在 k 阶滞后下估计偏自相关系数的计算公式如下（*）其中：rk 是在 k 阶滞后时的自相关系数估计值。（*）这是偏自相关系数的一致估计。,66,要得到k,k的更确切的估计，需要进行回归 t=1,2,T（*）因此滞后 k 阶的偏自相关系数是当
28、Xt 对 Xt-1，Xt-k 作回归时 Xt-k 的系数。称之为偏相关是因为它度量了k 期间距的相关而不考虑 k-1 期的相关。,从Xt中去掉Xt-1的影响，则只剩下随机扰动项t，显然它与Xt-2无关，因此Xt与Xt-2的偏自相关系数为零，记为,在AR(1)中，,同样，在AR(p)过程中，对所有的kp，Xt与Xt-k间的偏自相关系数为零。AR(p)的一个主要特征是：kp时，k*=Corr(Xt,Xt-k)=0 即k*在p以后截尾。,AR(p)随机时间序列的识别原则：若Xt的偏自相关函数在p以后截尾,即kp时，k*=0，而它的自相关函数k是拖尾的，则此序列是自回归AR(p)序列。,AR(1)过程
29、，时序图=0.2,AR(1)过程，时序图,=-0.2,AR(1)过程，时序图,=1.02,对MA(1)过程,2.MA(q)过程,它的自协方差系数：,于是，MA(1)过程的自相关函数为：,可见，当k1时，k0，即Xt与Xt-k不相关，MA(1)自相关函数是截尾的。,MA(1)过程可以等价地写成t 关于无穷序列Xt，Xt-1，的线性组合的形式：,(*)是一个AR()过程，它的偏自相关函数非截尾但却趋于零，因此MA(1)的偏自相关函数是非截尾但却趋于零的。,其自协方差系数为:,一般地，q阶移动平均过程MA(q),相应的自相关函数为:,当kq时，Xt与Xt-k不相关,k=0,即存在截尾现象,这是MA(
30、q)的一个特征。可以根据自相关系数是否从某一点开始一直为0来判断MA(q)模型的阶数。MA(q)过程的偏自相关函数是非截尾但趋于零的。,在实际识别时，由于样本自相关函数rk 是总体自相关函数k的一个估计，由于样本的随机性，当k q 时，rk不会全为0，而是在0的上下波动。但可以证明,当k q时,rk 服从如下渐近正态分布:rkN(0,1/n)式中n表示样本容量。,则有95.5%的把握判断原时间序列在q之后截尾。,因此，如果计算的rk满足,MA(q)模型的识别规则：若随机序列的自相关函数截尾，即自q以后,k=0(kq)：而它的偏自相关函数拖尾,则此序列是移动平均MA(q)序列。,ARMA(p,q
31、)的自相关函数，可以看作MA(q)的自相关函数和AR(p)的自相关函数的混合。当p=0时，它具有截尾性质；当q=0时，它具有拖尾性质；当p、q都不为0时，它具有拖尾性质;从识别上看，通常：ARMA(p,q)过程的偏自相关函数(PACF)可能在p阶滞后前有几项明显的尖柱(spikes)，但从p阶滞后项开始逐渐趋向于零；而它的自相关函数(ACF)则是在q阶滞后前有几项明显的尖柱,从q阶滞后项开始逐渐趋向于零。,3.ARMA(p,q)过程,ARMA(p,q)模型的ACF与PACF理论模式,ACF PACF,模型2:Xt=-0.7Xt-1+t,模型3:Xt=t 0.7t,例1中GDP是一阶单整的，GD
32、P是平稳序列：,AR(p)、MA(q)、ARMA(p,q)模型的估计方法较多，大体上分为3类：利用自相关函数的直接估计 矩估计 最小二乘估计,结构阶数,模型识别,确定,估计,参数,四、随机时间序列模型的估计,1.AR(p)模型的Yule Walker方程估计,在AR(p)模型的识别中，曾得到,利用k=-k，得到如下方程组：,此方程组称为Yule Walker方程组。该方程组建立了AR(p)模型的模型参数1,2,p与自相关函数1,2,p的关系。,一般地，p阶自回归模型AR(p):Xt=1Xt-1+2Xt-2+pXt-p+t,k期滞后协方差为:,从而有自相关函数:,利用实际时间序列提供的信息，首先
33、求得自相关函数的估计值,然后利用Yule Walker方程组，求解模型参数的估计值,由于,于是,从而可得2的估计值,2.MA(q)模型的矩估计,将MA(q)模型的自协方差函数中的各个量用估计量代替，得到：,(*),首先求得自协方差函数的估计值，(*)是一个包含(q+1)个待估参数 的非线性方程组,可以用直接法或迭代法求解。,常用的迭代方法有线性迭代法和Newton-Raphsan迭代法。,（1）MA(1)模型的直接算法,对于MA(1)模型，（*）式相应地写成,于是有解,由于参数估计有两组解，可根据可逆性条件|1|1来判断选取一组。,于是,（2）MA(q)模型的迭代算法,对于q1的MA(q)模型
34、，一般用迭代算法估计参数。由（*）式得,代入（*）式，计算出第一次迭代值,（*）,第二步，将第一次迭代值代入（*）式，计算出第二次迭代值,按此反复迭代下去，直到第m步的迭代值与第m-1步的迭代值相差不大时（满足一定的精度），便停止迭代，并用第m步的迭代结果作为（*）的近似解。,3.ARMA(p,q)模型的矩估计,在ARMA(p,q)中共有(p+q+1)个待估参数1,2,p与1,2,q以及2，其估计量计算步骤及公式如下：第一步，估计1,2,p,第二步，改写模型,求1,2,q以及2的估计值.,将模型:,改写为：,令,于是(*)可以写成：,构成一个MA模型。按照估计MA模型参数的方法,可以得到1,2
35、,q以及2的估计值。,4.AR(p)的最小二乘估计,假设模型AR(p)的参数估计值已经得到，即有,残差的平方和为：,(*),根据最小二乘原理，所要求的参数估计值是下列方程组的解：,即,j=1,2,p(*),解该方程组，就可得到待估参数的估计值。,注意，在上述模型的平稳性、识别与估计的讨论中，ARMA(p,q)模型中均未包含常数项。,如果包含常数项，该常数项并不影响模型的原有性质，因为通过适当的变形，可将包含常数项的模型转换为不含常数项的模型。以一般的ARMA(p,q)模型为例，对含有常数项的模型：,方程两边同减/(1-1-p)，则可得到,其中,(*),如果估计的ARMA(p,q)模型正确，残差
36、应代表一白噪声序列，否则说明模型的识别与估计有误，需重新识别与估计。,1.残差项的白噪声检验,五、模型的检验,实际检验时，主要检验残差序列是否存在自相关。可用Q 统计量进行2检验来检验是否拒绝残差序列为白噪声的零假设。,2.AIC与SBC模型选择标准,有可能存在不止一组(p,q)值都能通过识别检验。对可能的适当的模型，存在着模型的“简洁”与模型的拟合优度的权衡选择问题。,常用的模型选择的判别标准有：AIC与SBC,注意：在不同模型间进行比较时，必须选取相同的时间段。,中国支出法GDP是非平稳的，但它的一阶差分是平稳的，即支出法GDP是I(1)时间序列。可以对经过一阶差分后的GDP建立适当的AR
37、MA(p,q)模型。记GDP经一阶差分后的新序列为GDPD1，该新序列的样本自相关函数图与偏自相关函数图如下：Eviews,例，中国支出法GDP的ARMA(p,q)模型估计。,图形：样本自相关函数图形呈正弦线型衰减波，而偏自相关函数图形则在滞后两期后迅速趋于0。因此可初步判断该序列满足2阶自回归过程AR(2)。,可以认为：偏自相关函数是截尾的。再次验证了一阶差分后的GDP满足AR(2)随机过程。,自相关函数与偏自相关函数的函数值：自相关函数具有明显的拖尾性；偏自相关函数值在k2以后，,设序列GDPD1的模型形式为：,有如下Yule Walker 方程：,解为：,用OLS法回归的结果为：,(1)
38、,GDPD1序列AR(2)模型Eviews的OLS法参数估计结果：,有时，在用回归法时，也可加入常数项。本例中加入常数项的回归OLS估计结果为：,模型检验,Eviews可以给出三模型的残差项的自相关系数及Q 检验值。,模型1与模型3的残差项接近于一白噪声，但模型2存在4阶滞后相关问题，Q统计量的检验也得出模型2拒绝所有自相关系数为零的假设。因此：模型1与3可作为描述中国支出法GDP一阶差分序列的随机生成过程。,GDPD1序列模型1残差序列的ACF与Q统计量：,GDPD1序列模型2残差序列的ACF与Q统计量：,GDPD1序列模型3残差序列的ACF与Q统计量：,用建立的AR(2)模型对中国支出法GDP进行外推预测。,模型1可作如下展开：,如果已知t-1、t-2、t-3期的GDP时，就可对第t期的GDP作出外推预测。模型3的预测式与此相类似,只多出一项常数项。,对2001年中国支出法GDP的预测结果（亿元）预测值 实际值 误差 模型1 95469 95933-0.48%模型3 97160 95933 1.28%,
