时间序建模分析

目录l 1、ARIMA模型 1.1 模型的适用条件与构建过程 1.2 EVIEWS操作简单说明 1.3 模型构建实例l 2、季节时间序列模型 2.1 确定性季节时间序列模型 2.2 随机性季节时间序列模型时间序列的预处理： 拿到一个时间序列后，首先要对它的平稳性和纯随机性进行检验，这两个重要的检验称为序列的预处理。根据检验的结果可以将序列分为不同的类型，对不同类型的序列采取不同的分析方法。时间序列的基本类型：时间序列平稳时间序列非平稳时间序列平稳白噪声序列平稳非白噪声序列确定性时序分析随机性时序分析没有分析价值模型拟合（常用ARMA模型）长期趋势循环波动季节性变化平稳性检验纯随机性检验随机波动ARIMA模型残差自回归模型条件异方差模型平稳性检验方法：l 图检验方法 l 构造检验统计量时序图检验自相关图检验 主观色彩较强单位根检验平稳非平稳有明显趋势或周期性，则为非平稳随着延迟期数增加，自相关系数会很快衰减向零反之，自相关系数衰减向零的速度较慢纯随机性检验方法：构造检验统计量大样本场合大，小样本场合Q统计量LB统计量否则，认为该序列为纯随机序列对Q统计量修正若P值非常小（0.05）则认为该序列属于非白噪声序列检验结果（有分析价值）（无分析价值）平稳非白噪声序列建模步骤：平稳非白噪声序列预测序列将来的走势计算ACF,PACFARMA模型识别估计模型中未知参数的值模型优化模型检验NYARIMA模型建模流程：获得观察值序列拟合ARMA模型差分运算分析结束平稳性检验白噪声检验NYNYEVIEWS 操作数据录入画图自相关和偏自相关图建立方程Q检验预测例：某国1980年至1993年GNP平减指数的季节时间序列，共56个观测值，见下表19801981198219831984198519861987198819891990199119921993表5.1 某国GNP平减指数季度资料189.8994.498.72102.95110.72122.88130.12136.8145.12158.6171.94190.01203.98212.87291.0795.799.42104.75113.48124.44131.3139.01148.89161.85176.46193.03206.77214.25391.7996.52100.25106.53116.42126.68132.89141.03152.02165.12180.24197.7208.53215.89493.0397.39101.54108.74119.79128.99134.99143.24155.38168.05185.13201.69210.27218.21年/季该序列时序图（1.1）和自相关图（1.2）如下：图（1.1）图（1.2）该图显示有明显的长期趋势 自相关系数随延迟期数的增加，衰减向零的速度相当缓慢，且后期有反向递增趋势序列非平稳序列GNP的单位根检验结果：检验t统计量的值是0.325604，大于各个显著性水平下的临界值，所以不能拒绝原假设。也就是说，序列GNP存在单位根，因此，是非平稳的。一阶差分后的时序图与自相关图：图（1.3）图（1.4）时序图仍显示有长期趋势自相关系数向零衰减的速度依然较慢一阶差分序列仍不平稳一阶差分序列D(GNP)的单位根检验结果：检验t统计量的值是-1.929760，大于各个显著性水平下的临界值，所以不能拒绝原假设。也就是说，一阶差分序列D(GNP)存在单位根，因此，一阶差分序列也是非平稳的。2阶差分时序图与自相关图：图（1.5）图（1.6） 差分序列在零附近波动，无明显趋势或周期自相关系数在零值附近波动认为2阶差分序列平稳二阶差分序列的单位根检验：检验t统计量的值是-3.709559，小于各个显著性水平下的临界值，所以拒绝原假设。也就是说，二阶差分序列不存在单位根。二阶差分序列平稳。对平稳的2阶差分序列进行白噪声检验： 在显著性水平为0.05的条件下，延迟期数为6和12时，Q统计量的P值均小于0.052阶差分序列为非白噪声序列 结合前面分析，认为该序列为2阶差分平稳非白噪声序列，可考虑建立ARIMA模型根据2阶差分序列的自相关图ACF和偏自相关图PACF的特点，判断阶数进行建模:可以尝试用ARMA(2,2) ARMA(3,2) ARMA(3,3);也就是说，对原序列GNP尝试用ARIMA(2,2,2) ARIMA(3,2,2) ARIMA(3,2,3)进行拟合，首先建立ARIMA(2,2,2)如下：C与MA(1)系数的T检验显示：由于P值均大于0.05，故接受原假设，即二者系数显著为零，所以剔除模型ARiMA(2,2,2)：d(gnp,2) ar(1) ar(2) c ma(1) ma(2)剔除C与MA（1）:ARIMA(2,2,(2) : d(gnp,2) ar(1) ar(2) ma(2)可供选用模型一模型参数均通过检验建立ARIMA(3,2,2)如下：ARIMA(3,2,2):d(gnp,2) ar(1) ar(2) ar(3) ma(1) ma(2)AR(3)系数未通过检验，予以剔除结果和前述模型相同建立ARIMA(3,2,3):命令为：d(gnp,2) ar(1) ar(2) ar(3) ma(1) ma(2) ma(3)可供选用模型二模型适用性检验：模型ARIMA(2,2,(2)模型ARIMA(3,2,3)通过对模型的适用性检验，左侧拟合模型中的残差白噪声检验显示延迟6阶，12阶，18阶的残差序列属于白噪声序列，模型ARIMA(2,2,(2)显著有效，对序列适应性更强。因此，选用该模型作为最终拟合模型。模型预测结果：2220.868001ttt（1-B） （1+0.328913B+0.806248B）X2220.868001ttt（1-B） （1+0.328913B+0.806248B ）X拟合曲线对比： 拟合曲线与原序列曲线十分接近，直观来看，拟合效果较好！预测值的比较原始值 ARIMA(2,2,(2) ARIMA(3,2,3)212.87212.01211.69214.87215.51216.01215.89216.08214.91218.21217.32219.0693Q193Q293Q393Q4l对我国1990年1月至1997年12月工业总产值的月度资料（1990年为不变价格）共有96个观测值进行时间序列拟合，并对1998年工业总产值进行预测。1990年年1月至月至1997年年12月我国工业总产值月我国工业总产值单位：亿元单位：亿元一阶差分一阶差分二阶差分二阶差分 具体操作具体操作利用指数平滑法对以上图形进行拟合实际值实际值预测值预测值预测误差预测误差最终预测值最终预测值季节指数季节指数指数平滑预测值指数平滑预测值 123,1,11,1,1 120,1,10,1,1