非平稳性序列的随机分析

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第五章,非平稳序列的随机分析,1,本章结构,差分运算,ARIMA,模型,Auto-Regressive,模型,异方差的性质,方差齐性变化,条件异方差模型,2,5.1 差分运算,差分运算的实质,差分方式的选择,过差分,3,差分运算的实质,差分方法是一种非常简便、有效的确定性信息提取方法,Cramer,分解定理在理论上保证了适当阶数的差分一定可以充分提取确定性信息,差分运算的实质是使用自回归的方式提取确定性信息,4,差分方式的选择,序列蕴含着显著的线性趋势，一阶差分就可以实现趋势平稳,序列蕴含着曲线趋势，通常低阶（二阶或三阶）差分就可以提取出曲线趋势的影响,对于蕴含着固定周期的序列进行步长为周期长度的差分运算，通常可以较好地提取周期信息,5,例5.1,【例,1.1】1964,年,1999,年中国纱年产量序列蕴含着一个近似线性的递增趋势。对该序列进行一阶差分运算,考察差分运算对该序列线性趋势信息的提取作用,6,差分前后时序图,原序列时序图,差分后序列时序图,7,例5.2,尝试提取,1950,年,1999,年北京市民用车辆拥有量序列,的确定性信息,8,差分后序列时序图,一阶差分,二阶差分,9,例5.3,差分运算提取,1962,年,1,月,1975,年,12,月平均每头奶牛的月产奶量序列中的确定性信息,10,差分后序列时序图,一阶差分,1阶12步差分,11,过差分,足够多次的差分运算可以充分地提取原序列中的非平稳确定性信息,但过度的差分会造成有用信息的浪费,12,例5.4,假设序列如下,考察一阶差分后序列和二阶差分序列 的平稳性与方差,13,比较,一阶差分,平稳,方差小,二阶差分（过差分）,平稳,方差大,14,5.2 ARIMA模型,ARIMA模型结构,ARIMA模型性质,ARIMA模型建模,ARIMA模型预测,疏系数模型,季节模型,15,ARIMA模型结构,使用场合,差分平稳序列拟合,模型结构,16,ARIMA 模型族,d=0,ARIMA(p,d,q)=ARMA(p,q),P=0,ARIMA(P,d,q)=IMA(d,q),q=0,ARIMA(P,d,q)=ARI(p,d),d=1,P=q=0,ARIMA(P,d,q)=random walk model,17,随机游走模型( random walk),模型结构,模型产生典故,Karl Pearson,(,1905,),在自然杂志上提问,：,假如有个醉汉醉得非常严重，完全丧失方向感，把他放在荒郊野外，一段时间之后再去找他，在什么地方找到他的概率最大呢？,18,ARIMA模型的平稳性,ARIMA(p,d,q),模型共有,p+d,个特征根，其中,p,个在单位圆内，,d,个在单位圆上。 所以当 时,ARIMA(p,d,q),模型非平稳。,例5.5,ARIMA(0,1,0),时序图,19,ARIMA模型的方差齐性,时，原序列方差非齐性,d,阶差分后，差分后序列方差齐性,20,ARIMA模型建模步骤,获,得,观,察,值,序,列,平稳性,检验,差分,运算,Y,N,白噪声,检验,Y,分,析,结,束,N,拟合,ARMA,模型,21,例5.6,对,1952,年,1988,年中国农业实际国民收入指数序列建模,22,一阶差分序列时序图,23,一阶差分序列自相关图,24,一阶差分后序列白噪声检验,延迟阶数,统计量,P,值,6,15.33,0.0178,12,18.33,0.1060,18,24.66,0.1344,25,拟合ARMA模型,偏自相关图,26,建模,定阶,ARIMA(0,1,1),参数估计,模型检验,模型显著,参数显著,27,ARIMA模型预测,原则,最小均方误差预测原理,Green,函数递推公式,28,预测值,29,例5.7,已知,ARIMA(1,1,1),模型为,且,求 的,95,的置信区间,30,预测值,等价形式,计算预测值,31,计算置信区间,Green函数值,方差,95置信区间,32,例5.6续：,对中国农业实际国民收入指数序列做为期,10,年的预测,33,疏系数模型,ARIMA(p,d,q),模型是指,d,阶差分后自相关最高阶数为,p，,移动平均最高阶数为,q,的模型，通常它包含,p+q,个独立的未知系数：,如果该模型中有部分自相关系数 或部分移动平滑系数 为零，即原模型中有部分系数省缺了，那么该模型称为疏系数模型。,34,疏系数模型类型,如果只是自相关部分有省缺系数，那么该疏系数模型可以简记为,为非零自相关系数的阶数,如果只是移动平滑部分有省缺系数，那么该疏系数模型可以简记为,为非零移动平均系数的阶数,如果自相关和移动平滑部分都有省缺，可以简记为,35,例5.8,对,1917,年,1975,年美国,23,岁妇女每万人生育率序列建模,36,一阶差分,37,自相关图,38,偏自相关图,39,建模,定阶,ARIMA(1,4),1,0),参数估计,模型检验,模型显著,参数显著,40,季节模型,简单季节模型,乘积季节模型,41,简单季节模型,简单季节模型是指序列中的季节效应和其它效应之间是加法关系,简单季节模型通过简单的趋势差分、季节差分之后序列即可转化为平稳，它的模型结构通常如下,42,例5.9,拟合,1962,1991,年德国工人季度失业率序列,43,差分平稳,对原序列作一阶差分消除趋势，再作,4,步差分消除季节效应的影响，差分后序列的时序图如下,44,白噪声检验,延迟阶数,统计量,P,值,6,43.84,0.0001,12,51.71,0.0001,18,54.48,0.0001,45,差分后序列自相关图,46,差分后序列偏自相关图,47,模型拟合,定阶,ARIMA(1,4),(1,4),0),参数估计,48,模型检验,49,拟合效果图,50,乘积季节模型,使用场合,序列的季节效应、长期趋势效应和随机波动之间有着复杂地相互关联性，简单的季节模型不能充分地提取其中的相关关系,构造原理,短期相关性用低阶,ARMA(p,q),模型提取,季节相关性用以周期步长,S,为单位的,ARMA(P,Q),模型提取,假设短期相关和季节效应之间具有乘积关系，模型结构如下,51,例5.10 :拟合1948,1981,年美国女性月度失业率序列,52,差分平稳,一阶、12步差分,53,差分后序列自相关图,54,差分后序列偏自相关图,55,简单季节模型拟合结果,延迟阶数,拟合模型残差白噪声检验,AR(1,12),MA(1,2,12),ARMA(1,12),(1,12),值,P,值,值,P,值,值,P,值,6,14.58,0.0057,9.5,0.0233,15.77,0.0004,12,16.42,0.0883,14.19,0.1158,17.99,0.0213,结果,拟合模型均不显著,56,乘积季节模型拟合,模型定阶,ARIMA(1,1,1)(0,1,1),12,参数估计,57,模型检验,58,乘积季节模型拟合效果图,59,5.3 Auto-Regressive模型,构造思想,首先通过确定性因素分解方法提取序列中主要的确定性信息,然后对残差序列拟合自回归模型，以便充分提取相关信息,60,Auto-Regressive模型结构,61,对趋势效应的常用拟合方法,自变量为时间t的幂函数,自变量为历史观察值,62,对季节效应的常用拟合方法,给定季节指数,建立季节自回归模型,63,例5.6续,使用,Auto-Regressive,模型分析,1952,年,1988,年中国农业实际国民收入指数序列。,时序图显示该序列有显著的线性递增趋势，但没有季节效应，所以考虑建立如下结构的,Auto-Regressive,模型,64,趋势拟合,方法一：,变量为时间,t,的幂函数,方法二：变量为一阶延迟序列值,65,趋势拟合效果图,66,残差自相关检验,检验原理,回归模型拟合充分，残差的性质,回归模型拟合得不充分，残差的性质,67,Durbin-Waston,检验（,DW,检验）,假设条件,原假设：残差序列不存在一阶自相关性,备择假设：残差序列存在一阶自相关性,68,DW统计量,构造统计量,DW统计量和自相关系数的关系,69,DW统计量的判定结果,正,相,关,相,关,性,待,定,不相关,相,关,性,待,定,负,相,关,0,4,2,70,例,5.6,续,检验第一个确定性趋势模型,残差序列的自相关性。,71,DW检验结果,检验结果,检验结论,检验结果显示残差序列高度正自相关。,72,Durbin h,检验,DW,统计量的缺陷,当回归因子包含延迟因变量时，,残差序列的,DW,统计量是一个有偏统计量。在这种场合下使用,DW,统计量容易产生残差序列正自相关性不显著的误判,Durbin h,检验,73,例5.6续,检验第二个确定性趋势模型,残差序列的自相关性。,74,Dh检验结果,检验结果,检验结论,检验结果显示残差序列高度正自相关。,Dh,统计量的值,P,值,2.8038,0.0025,75,残差序列拟合,确定自回归模型的阶数,参数估计,模型检验,76,例5.6续,对第一个确定性趋势模型,的残差序列,进行拟合,77,残差序列自相关图,78,残差序列偏自相关图,79,模型拟合,定阶,AR(2),参数估计方法,极大似然估计,最终拟合模型口径,80,例5.6,第二个AutoRegressive模型的拟合结果,81,三个拟合模型的比较,模型,AIC,SBC,ARIMA(0,1,1),模型：,249.3305,252.4976,AutoRegressive,模型一：,260.8454,267.2891,AutoRegressive,模型二：,250.6317,253.7987,82,