《平稳时间序列模型的建立教材P》由会员分享，可在线阅读，更多相关《平稳时间序列模型的建立教材P（98页珍藏版）》请在金锄头文库上搜索。
1、1、建模流程 （有限长度）时序样本模型识别与定阶模型参数估计模型适用性检验模型优化 2、基本前提 平稳序列Xt 零均值序列EXt=0,建模步骤,流程图,平 稳 非 白 噪 声 序 列,计 算 样 本 相 关 系 数,模型 识别,参数 估计,模型 检验,模 型 优 化,序 列 预 测,Y,N,一、平稳性检验 二、纯随机性检验 三、计算样本自相关函数 四、关于非零均值的平稳序列,第一节 时间序列的预处理,本章所介绍的是对零均值平稳序列建立ARMA模型，因此，在对实际的序列进行模型识别之前，应首先检验序列是否平稳，若序列非平稳，应先通过适当变换将其化为平稳序列，然后再进行模型识别.,序列的非平稳包括
2、均值非平稳和方差非平稳. 均值非平稳序列平稳化的方法：差分变换. 方差非平稳序列平稳化的方法：对数变换、平方根变换等. 序列平稳性的检验方法和手段主要有：序列趋势图、自相关图、单位根检验、非参数检验方法等等.,一、平稳性检验图检验方法,（一）时序图检验 根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界、无明显趋势及周期特征. （二）自相关图检验 平稳序列通常具有短期相关性.该性质用自相关函数来描述就是随着延迟期数的增加，平稳序列的自相关函数会很快地衰减向零.,例题,例1 检验1964年1999年中国纱年产量序列的平稳性 例
3、2 检验1962年1月1975年12月平均每头奶牛月产奶量序列的平稳性 例3 检验1949年1998年北京市每年最高气温序列的平稳性,例1 时序图,例1 自相关图,例2 时序图,例2 自相关图,例3 时序图,例3 自相关图,（一）纯随机序列的定义 纯随机序列也称为白噪声序列，它满足如下两条性质,二、纯随机性检验,（二）纯随机性检验,检验原理 假设条件 检验统计量 判别原则 应用举例,Barlett定理 如果一个时间序列是纯随机的，得到一个观察期数为 的观察序列，那么该序列的延迟非零期的样本自相关系数将近似服从均值为零，方差为序列观察期数倒数的正态分布,1、检验原理,2、假设条件,原假设：延迟期
4、数小于或等于 期的序列值之间相互独立 备择假设：延迟期数小于或等于 期的序列值之间有相关性,3、检验统计量,Q统计量 LB统计量,4、判别原则,拒绝原假设 当检验统计量大于 分位点，或该统计量的P值小于 时，则可以以 的置信水平拒绝原假设，认为该序列为非白噪声序列 接受原假设 当检验统计量小于 分位点，或该统计量的P值大于 时，则认为在 的置信水平下无法拒绝原假设，即不能显著拒绝序列为纯随机序列的假定,例4、标准正态白噪声序列纯随机性检验,样本自相关图,5、应用举例,检验结果,由于P值显著大于显著性水平 ，所以该序列不能拒绝纯随机的原假设.,例5、对1950年1998年北京市城乡居民定期储蓄所
5、占比例序列的平稳性与纯随机性进行检验,自相关图,白噪声检验结果,三、计算样本相关函数,样本自相关函数,样本偏自相关函数,四、关于非零均值的平稳序列,非零均值的平稳序列有两种处理方法： 设xt为一非零均值的平稳序列，且有E(xt)= 方法一: 用样本均值 作为序列均值的估计，建模前先对序列作如下处理： 令 然后对零均值平稳序列wt建模.,方法二 在模型识别阶段对序列均值是否为零不予考虑，而在参数估计阶段，将序列均值作为一个参数加以估计. 以一般的ARMA(p,q)为例说明如下：,将上式展开得：,此时，所要估计的未知参数有p+q+1个.,第二节 模型识别与定阶,一、模型识别 二、模型定阶,一、模型
6、识别,基本原则,零均值平稳序列模型识别的主要根据是序列的自相关函数(ACF)和偏自相关函数(PACF)的特征. 若序列xt的偏自相关函数 在kp以后截尾，即kp 时， ，而且它的自相关函数 拖尾，则可判断此序列是AR(p)序列.,若序列xt的自相关函数 在kq以后截尾，即kq 时， ，而且它的偏自相关函数 拖尾，则可判断此序列是MA(q)序列. 若序列xt的自相关函数、偏相关函数都呈拖尾形态，则可断言此序列是ARMA序列. 若序列的自相关函数和偏自相关函数不但都不截尾，而且至少有一个下降趋势势缓慢或呈周期性衰减，则可认为它也不是拖尾的，此时序列是非平稳序列，应先将其转化为平稳序列后再进行模型识
7、别.,模型定阶的困难,因为由于样本的随机性，样本的相关系数不会呈现出理论截尾的完美情况，本应截尾的 或 仍会呈现出小值振荡的情况 由于平稳时间序列通常都具有短期相关性，随着延迟阶数 ， 与 都会衰减至零值附近作小值波动 ？当 或 在延迟若干阶之后衰减为小值波动时，什么情况下该看作为相关系数截尾，什么情况下该看作为相关系数在延迟若干阶之后正常衰减到零值附近作拖尾波动呢？,二、模型定阶,样本相关系数的近似分布,Barlett Quenouille,1、经验定阶方法,95的置信区间,模型定阶的经验方法 如果样本(偏)自相关系数在最初的p阶明显大于两倍标准差范围，而后几乎95的自相关系数都落在2倍标准
8、差的范围以内，而且通常由非零自相关系数衰减为小值波动的过程非常突然.这时，通常视为(偏)自相关系数截尾.截尾阶数为p.,例1,上海延中实业股票数据识别（一阶差分后）,平均每日生产汽车废品数据的识别(n=45),美国女性失业月数据识别（差分后）,上海延中实业股份有限公司是上海首家向社会公开发行股票的企业. 1985年1月底发行股票500万元，其中由上海延中复印工业公司出资30万元.上海延中实业股票收盘价基本反映了沪市股票的大致走向.总观测期n619，先作出原序列的样本自相关函数和样本偏相关函数，其结果见表1和图1.,上海延中实业股票数据识别（一阶差分后）,表1 延中股票的样本自相关和样本偏自相关
9、函数值,美国1961年1月至1985年12月间女性失业月人数时间序列,美国女性失业月数据识别（差分后）,2、残差方差图定阶法,（1）基本思想 如果拟合的模型阶数与真正阶数不符合，则模型的残差平方和SSE必然偏大，残差方差将比真正模型的残差方差大。 如果是不足拟合，那么逐渐增加模型阶数，模型的残差方差会渐减少，直到残差方差达到最小。 如果是过度拟合，此时逐渐少模型阶数，模型残差方差分逐渐下降，直到残差方差达到最小。,（2）残差方差的估计公式,注：式中 “实际观察值个数”是指拟合模型时实际使用 的观察值项数，即经过平稳化后的有效样本容量。 设原序列有n个样本，若建立的模型中有含有自回归AR部分,
10、且阶数为p，则实际观察值个数为n-p个。若没有AR部分，则 实际观察值个数即为n个。 模型的参数个数指模型中所含的参数个数，如：若是不 带常数项的ARMA(p,q)模型，参数个数为p+q个，若带有 常数项，则参数个数为p+q+1个。,用Eviews建立ARMA模型后，可直接得到剩余平方和SSE(Sum squared resid) 输出结果中也可直接得到残差标准差： S.E.of regression ，此项的平方即为残差方差。因此，对不同的模型残差方差进行比较，直接比较此项既可。,例：以磨轮剖面数据为例，分别建立适应性模 型，输出结果见图示，从中选择最佳模型。,三个模型残差方差比较,3、F检
11、验定阶法,基本思想（以一般情形和ARMA(p,q)模型为例） 先对数据拟合ARMA(p,q)模型(假设不含常数项)，设其残差平方和为Q0，再对数据拟合 较低阶的模型ARMA(p-m,q-s)，设其残差平方和为Q1。 建立原假：,在原假设成立的条件下有：,于是计算统计量F，在给定的显著性水平下。 若FF ，则拒绝原假设，说明两模型差异是显著的， 此时模型阶数存在升高的可能性。 若FF ，此不能拒绝原假设，说明两模型差异不显著， 此时模型阶数存在降低的可能性。,注：F检验定阶法的应用条件：两模型中有一个为合适模型。,4、最佳准则函数定阶法,最佳准则函数法，即确定出一个准则函数，该函数既要考虑某一模
12、型拟合时对原始数据的接近程度，同时又要考虑模型中所含待定参数的个数。 建模时，使准则函数达到极小的是最佳模型。,4.1 赤池的AIC准则和BIC准则 4.1.1 AIC 准则(Akaike iformationcriterion) AIC准则是1973年由赤池（Akaike）提出，此准则是对FPE准则(用来判别AR模型的阶数是否合适)的推广，用来识别ARMA模型的阶数。,AIC准则函数为：,式中，M为模型中参数的个数。,AIC的简化式为：,式中： 是残差方差 的极大似然估计值。,Eviews输出的Akaike info criterion与上述形式略有差别(参见Eviews help)，其定义
13、为：,其中：n是实际观察值的个数。,4.1.2 BIC准则 柴田（Shibata）1976年证明AIC有过分估计自回归参数的倾向，于是Akaike又提出了AIC方法的贝叶斯扩展，即BIC。 BIC准则函数为：,式中：C为常数。余同前。,4.2 施瓦茨(Schwarz)的SC准则 此准则1978年由Schwarz提出，被称为SBC(Schwartzs Bayesian criterion)。 准则函数：,简化式为：,同样Eviews输出的结果与上形式略有差别，其定义为：,准则函数使用注意,1、当样本量趋于无穷时，用AIC准则挑选的最佳模型的阶数往往比真实模型阶数高，而用SBC准则确定的最佳模型的
14、阶数往往与真实模型的阶数相一致。 2、样本量不是很大时，SBC准则的定阶效果不及AIC。,一、矩估计 二、极大似然估计 三、最小二乘估计,第三节 模型参数估计,一、矩估计,原理 样本自相关系数估计总体自相关系数 样本一阶均值估计总体均值，样本方差估计总体方差,例1 求AR(2)模型系数的矩估计,AR(2)模型 Yule-Walker方程 矩估计（Yule-Walker方程的解）,例2 求MA(1)模型系数的矩估计,MA(1)模型 方程 矩估计,例3 求ARMA(1,1)模型系数的矩估计,ARMA(1,1)模型 方程 矩估计,对矩估计的评价,优点 估计思想简单直观 不需要假设总体分布 计算量小（
15、低阶模型场合） 缺点 信息浪费严重 只用到了p+q个样本自相关系数信息，其他信息都被忽略 估计精度差 通常矩估计方法被用作极大似然估计和最小二乘估计迭代计算的初始值,二、极大似然估计,原理 在极大似然准则下，认为样本来自使该样本出现概率最大的总体。因此未知参数的极大似然估计就是使得似然函数（即联合密度函数）达到最大的参数值,似然方程,由于 和 都不是 的显式表达式。因而似然方程组实际上是由p+q+1个超越方程构成，通常需要经过复杂的迭代算法才能求出未知参数的极大似然估计值,对极大似然估计的评价,优点 极大似然估计充分应用了每一个观察值所提供的信息，因而它的估计精度高 同时还具有估计的一致性、渐
16、近正态性和渐近有效性等许多优良的统计性质 缺点 需要假定总体分布,三、最小二乘估计,原理 使残差平方和达到最小的那组参数值即为最小二乘估计值,条件最小二乘估计,实际中最常用的参数估计方法 假设条件 残差平方和方程 解法 迭代法,对最小二乘估计的评价,优点 最小二乘估计充分应用了每一个观察值所提供的信息，因而它的估计精度高 条件最小二乘估计方法使用率最高 缺点 需要假定总体分布,例4,确定1950年1998年北京市城乡居民定期储蓄比例序列拟合模型的口径 拟合模型：AR(1) 估计方法：极大似然估计 模型口径,例5,确定美国科罗拉多州某一加油站连续57天的OVERSHORTS序列拟合模型的口径 拟合模型：MA(1) 估计方法：条件最小二乘估计 模型口径,例6,确定1880-1985全球气表平均温度改变值差分序列拟合模型的口径 拟合模型：ARMA(1,1) 估计方法：条件
