《时间序列统计讲义》 非平稳序列的随机分析

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第五章,非平稳序列的随机分析,本章结构,差分运算,ARIMA,模型,Auto-Regressive,模型,异方差的性质,方差齐性变化,条件异方差模型,5.1 差分运算,差分运算的实质,差分方式的选择,过差分,差分运算的实质,差分方法是一种非常简便、有效的确定性信息提取方法,Cramer,分解定理在理论上保证了适当阶数的差分一定可以充分提取确定性信息,差分运算的实质是使用自回归的方式提取确定性信息,差分方式的选择,序列蕴含着显著的线性趋势，一阶差分就可以实现趋势平稳,序列蕴含着曲线趋势，通常低阶二阶或三阶差分就可以提取出曲线趋势的影响,对于蕴含着固定周期的序列进行步长为周期长度的差分运算，通常可以较好地提取周期信息,【例,1.1】1964,年,1999,年中国纱年产量序列蕴含着一个近似线性的递增趋势。对该序列进行一阶差分运算,考察差分运算对该序列线性趋势信息的提取作用,差分前后时序图,原序列时序图,差分后序列时序图,尝试提取,1950,年,1999,年北京市民用车辆拥有量序列,的确定性信息,差分后序列时序图,一阶差分,二阶差分,差分运算提取,1962,年,1,月,1975,年,12,月平均每头奶牛的月产奶量序列中的确定性信息,差分后序列时序图,一阶差分,1阶12步差分,过差分,足够多次的差分运算可以充分地提取原序列中的非平稳确定性信息,但过度的差分会造成有用信息的浪费,假设序列如下,考察一阶差分后序列和二阶差分序列 的平稳性与方差,比较,一阶差分,平稳,方差小,二阶差分过差分,平稳,方差大,5.2 ARIMA模型,ARIMA模型结构,ARIMA模型性质,ARIMA模型建模,ARIMA模型预测,疏系数模型,季节模型,ARIMA模型结构,使用场合,差分平稳序列拟合,模型结构,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,随机游走模型( random walk),模型结构,模型产生典故,Karl Pearson(1905)在?自然?杂志上提问：假如有个醉汉醉得非常严重，完全丧失方向感，把他放在荒郊野外，一段时间之后再去找他，在什么地方找到他的概率最大呢？,ARIMA模型的平稳性,ARIMA(p,d,q),模型共有,p+d,个特征根，其中,p,个在单位圆内，,d,个在单位圆上。所以当 时,ARIMA(p,d,q),模型非平稳。,ARIMA(0,1,0),时序图,ARIMA模型的方差齐性,时，原序列方差非齐性,d,阶差分后，差分后序列方差齐性,ARIMA模型建模步骤,获,得,观,察,值,序,列,平稳性,检验,差分,运算,Y,N,白噪声,检验,Y,分,析,结,束,N,拟合,ARMA,模型,对,1952,年,1988,年中国农业实际国民收入指数序列建模,一阶差分序列时序图,一阶差分序列自相关图,一阶差分后序列白噪声检验,延迟阶数,统计量,P,值,6,15.33,0.0178,12,18.33,0.1060,18,24.66,0.1344,拟合ARMA模型,偏自相关图,建模,定阶,ARIMA(0,1,1),参数估计,模型检验,模型显著,参数显著,ARIMA模型预测,原那么,最小均方误差预测原理,Green函数递推公式,预测值,ARIMA(1,1,1)模型为,且,求 的95的置信区间,预测值,等价形式,计算预测值,计算置信区间,Green函数值,方差,95置信区间,例5.6续：,对中国农业实际国民收入指数序列做为期,10,年的预测,疏系数模型,ARIMA(p,d,q),模型是指,d,阶差分后自相关最高阶数为,p，,移动平均最高阶数为,q,的模型，通常它包含,p+q,个独立的未知系数：,如果该模型中有部分自相关系数 或部分移动平滑系数 为零，即原模型中有部分系数省缺了，那么该模型称为疏系数模型。,疏系数模型类型,如果只是自相关部分有省缺系数，那么该疏系数模型可以简记为,为非零自相关系数的阶数,如果只是移动平滑部分有省缺系数，那么该疏系数模型可以简记为,为非零移动平均系数的阶数,如果自相关和移动平滑部分都有省缺，可以简记为,对,1917,年,1975,年美国,23,岁妇女每万人生育率序列建模,一阶差分,自相关图,偏自相关图,建模,定阶,ARIMA(1,4),1,0),参数估计,模型检验,模型显著,参数显著,季节模型,简单季节模型,乘积季节模型,简单季节模型,简单季节模型是指序列中的季节效应和其它效应之间是加法关系,简单季节模型通过简单的趋势差分、季节差分之后序列即可转化为平稳，它的模型结构通常如下,拟合,1962,1991,年德国工人季度失业率序列,差分平稳,对原序列作一阶差分消除趋势，再作,4,步差分消除季节效应的影响，差分后序列的时序图如下,白噪声检验,延迟阶数,统计量,P,值,6,43.84,0.0001,12,51.71,0.0001,18,54.48,0.0001,差分后序列自相关图,差分后序列偏自相关图,模型拟合,定阶,ARIMA(1,4),(1,4),0),参数估计,模型检验,残差白噪声检验,参数显著性检验,延迟,阶数,统计量,P,值,待估,参数,统计量,P,值,6,2.09,0.7191,5.48,0.0001,12,10.99,0.3584,-3.41,0.0001,拟合效果图,乘积季节模型,使用场合,序列的季节效应、长期趋势效应和随机波动之间有着复杂地相互关联性，简单的季节模型不能充分地提取其中的相关关系,构造原理,短期相关性用低阶,ARMA(p,q),模型提取,季节相关性用以周期步长,S,为单位的,ARMA(P,Q),模型提取,假设短期相关和季节效应之间具有乘积关系，模型结构如下,例5.10 :拟合1948,1981,年美国女性月度失业率序列,差分平稳,一阶、12步差分,差分后序列自相关图,差分后序列偏自相关图,简单季节模型拟合结果,延迟阶数,拟合模型残差白噪声检验,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,结果,拟合模型均不显著,乘积季节模型拟合,模型定阶,ARIMA(1,1,1)(0,1,1),12,参数估计,模型检验,残差白噪声检验,参数显著性检验,延迟,阶数,统计量,P,值,待估,参数,统计量,P,值,6,4.50,0.2120,-4.66,0.0001,12,9.42,0.4002,23.03,0.0001,18,20.58,0.1507,-6.81,0.0001,结果,模型显著,参数均显著,乘积季节模型拟合效果图,5.3 Auto-Regressive模型,构造思想,首先通过确定性因素分解方法提取序列中主要的确定性信息,然后对残差序列拟合自回归模型，以便充分提取相关信息,Auto-Regressive模型结构,对趋势效应的常用拟合方法,自变量为时间t的幂函数,自变量为历史观察值,对季节效应的常用拟合方法,给定季节指数,建立季节自回归模型,使用,Auto-Regressive,模型分析,1952,年,1988,年中国农业实际国民收入指数序列。,时序图显示该序列有显著的线性递增趋势，但没有季节效应，所以考虑建立如下结构的,Auto-Regressive,模型,趋势拟合,方法一：,变量为时间,t,的幂函数,方法二：变量为一阶延迟序列值,趋势拟合效果图,残差自相关检验,检验原理,回归模型拟合充分，残差的性质,回归模型拟合得不充分，残差的性质,Durbin-Waston检验DW检验,假设条件,原假设：残差序列不存在一阶自相关性,备择假设：残差序列存在一阶自相关性,DW统计量,构造统计量,DW统计量和自相关系数的关系,DW统计量的判定结果,正,相,关,相,关,性,待,定,不相关,相,关,性,待,定,负,相,关,0,4,2,例续,检验第一个确定性趋势模型,残差序列的自相关性。,DW检验结果,检验结果,检验结论,检验结果显示残差序列高度正自相关。,DW,统计量的值,P,值,0.1378,1.42,1.53,0.0001,Durbin h,检验,DW,统计量的缺陷,当回归因子包含延迟因变量时，,残差序列的,DW,统计量是一个有偏统计量。在这种场合下使用,DW,统计量容易产生残差序列正自相关性不显著的误判,Durbin h,检验,检验第二个确定性趋势模型,残差序列的自相关性。,Dh检验结果,检验结果,检验结论,检验结果显示残差序列高度正自相关。,Dh,统计量的值,P,值,2.8038,0.0025,残差序列拟合,确定自回归模型的阶数,参数估计,模型检验,对第一个确定性趋势模型,的残差序列,进行拟合,残差序列自相关图,残差序列偏自相关图,模型拟合,定阶,AR(2),参数估计方法,极大似然估计,最终拟合模型口径,第二个AutoRegressive模型的拟合结果,三个拟合模型的比较,模型,AIC,SBC,ARIMA(0,1,1),模型：,249.3305,252.4976,AutoRegressive,模型一：,260.8454,267.2891,AutoRegressive,模型二：,250.6317,253.7987,5.4 异方差的性质,异方差的定义,如果随机误差序列的方差会随着时间的变化而变化，这种情况被称作为异方差,异方差的影响,忽视异方差的存在会导致残差的方差会被严重低估，继而参数显著性检验容易犯纳伪错误，这使得参数的显著性检验失去意义，最终导致模型的拟合精度受影响。,异方差直观诊断,残差图,残差平方图,残差图,方差齐性残差图,递增型异方差残差图,残差平方图,原理,残差序列的方差实际上就是它平方的期望。,所以考察残差序列是否方差齐性，主要是考察残差平方序列是否平稳,直观考察美国,1963,年,4,月,1971,年,7,月短期国库券的月度收益率序列的方差齐性。,一阶差分后残差图,一阶差分后残差平方图,异方差处理方法,假如异方差函数具体形式，进行方差齐性变化,假如不知异方差函数的具体形式，拟合条件异方差模型,5.5 方差齐性变换,使用场合,序列显示出显著的异方差性，且方差与均值之间具有某种函数关系,其中： 是某个函数,处理思路,尝试寻找一个转换函数 ，使得经转换后的变量满足方差齐性,转换函数的确定原理,转换函数 在 附近作一阶泰勒展开,求转换函数的方差,转换函数的确定,常用转换函数的确定,假定,转换函数的确定,例,对美国,1963,年,4,月,1971,年,7,月短期国库券的月度收益率序列使用方差齐性变换方法进行分析,假定,函数变换,对数序列时序图,一阶差分后序列图,白噪声检验,延迟阶数,LB,统计量,P,值,6,3.58,0.7337,12,10.82,0.5441,18,21.71,0.2452,拟合模型口径及拟合效果图,5.6 条件异方差模型,ARCH,模型,GARCH,模型,GARCH,模型的变体,EGARCH,模型,IGARCH,模型,GARCH-M,模型,AR-GARCH,模型,ARCH模型,假定,原理,通过构造残差平方序列的自回归模型来拟合异方差函数,ARCH(q)模型结构,GARCH 模型结构,使用场合,ARCH,模型实际上适用于异方差函数短期自相关过程,GARCH,模型实际上适用于异方差函数长期自相关过程,模型结构,GARCH模型的约束条件,参数非负,参数有界,EGARCH模型,IGARCH模型,GARCH-M模型,AR-GARCH模型,GARCH模型拟合步骤,回归拟合,残差自相关性检验,异方差自相关性检验,ARCH模型定阶,参数估计,正态性检验,例,使用条件异方差模型拟合某金融时间序列。,回归拟合,拟合模型,参数估计,参数显著性检验,P,值0.0001,参数高度显著,残差自相关性检验,残差序列,DW,检验结果,拟合残差自回归模型,方法：逐步回归,模型口径,异方差自相关检验,Portmantea Q检验,拉格朗日乘子LM检验,Portmantea Q,检验,假设条件,检验统计量,检验结果,拒绝原假设,接受原假设,LM,检验,假设条件,检验统计量,检验结果,拒绝原假设,接受原假设,例,ARCH模型拟合,定阶：GARCH(1,1),参数估计：极大似然估计,拟合模型口径：AR(2)-GARCH(1,1),模型检验,检验方法：正态性检验,假设条件：,检验统计量,检验结果,拒绝原假设,接受原假设,例,P,值,AR(2)-GARCH(1,1),模型显著成立,拟合效果图,