《时间序列分析第五章上机指导》由会员分享,可在线阅读,更多相关《时间序列分析第五章上机指导(15页珍藏版)》请在金锄头文库上搜索。
1、上机指导第五章5.8.1 拟合ARIMA模型由于ARMA模型是ARIMA模型的一种特例,所以在SAS系统中这两种模型的拟合都放在了ARIMA过程中。我们已经在第3章进行了ARMA模型拟合时介绍了ARIMA过程的基本命令格式。再次以临时数据集example5_1的数据为例介绍ARIMA模型拟合与ARMA模型拟合的不同之处。 data example5_1; input x; difx=dif(x); t=_n_; cards;1.05 -0.84 -1.42 0.20 2.81 6.72 5.40 4.38 5.52 4.46 2.89 -0.43 -4.86 -8.54 -11.54 -16.
2、22-19.41 -21.61 -22.51 -23.51 -24.49 -25.54 -24.06 -23.44-23.41 -24.17 -21.58 -19.00 -14.14 -12.69 -9.48 -10.29-9.88 -8.33 -4.67 -2.97 -2.91 -1.86 -1.91 -0.80 proc gplot; plot x*t difx*t; symbol v=star c=black i=join;run;输出时序图显示这是一个典型的非平稳序列。如图5-49所示图5-49 序列x时序图考虑对该序列进行1阶差分运算,同时考察查分后序列的平稳性,在原程序基础上添加相
3、关命令,程序修改如下: data example5_1; input x; difx=dif(x); t=_n_; cards;1.05 -0.84 -1.42 0.20 2.81 6.72 5.40 4.38 5.52 4.46 2.89 -0.43 -4.86 -8.54 -11.54 -16.22-19.41 -21.61 -22.51 -23.51 -24.49 -25.54 -24.06 -23.44-23.41 -24.17 -21.58 -19.00 -14.14 -12.69 -9.48 -10.29-9.88 -8.33 -4.67 -2.97 -2.91 -1.86 -1.
4、91 -0.80 proc gplot; plot x*t difx*t; symbol v=star c=black i=join; proc arima; identify var=x(1); estimate p=1; forecast lead=5 id=t ; run;语句说明:(1)DATA步中的命令“difx=dif(x);”,这是指令系统对变量x进行1阶差分,差分后的序列值赋值给变量difx。其中dif()是差分函数,假如要差分的变量名为x,常见的几种差分表示为:1阶差分:dif(x)2阶差分:dif(dif(x)k步差分:difk(x)(2)我们在GPLOT过程中添加绘制了一
5、个时序图“difx*t”,这是为了直观考察1阶差分后序列的平稳性。所得时序图如图5-50所示。图5-50 序列difx时序图时序图显示差分后序列difx没有明显的非平稳特征。(3)“identify var=x(1);”,使用该命令可以识别查分后序列的平稳性、纯随机性和适当的拟合模型阶数。其中x(1)表示识别变量x的1阶差分后序列。SAS支持多种形式的差分序列识别:var=x(1),表示识别变量x的1阶查分后序列xt;var=x(1,1),表示识别变量x的2阶查分后序列2xt;var=x(k),表示识别变量x的k步差分后序列kxt;var=x(k,s),表示识别变量x的k步差分后,再进行s步查
6、分后序列skxt。识别部分的输出结果显示1阶查分后序列difx为平稳非白噪声序列,且具有显著的自相关系数不截尾、偏自相关系数1截尾的性质。(4)“estimate p=1;”对1阶差分后序列xt拟合AR(1)模型。输出拟合结果显示常数项不显著,添加或修改估计命令如下:estimate p=1 nonit;这是命令系统不要常数项拟合AR(1)模型,拟合结果显示模型显著且参数显著。如图5-51所示。图5-51 序列difx模型拟合结果输出结果显示,序列xt的拟合模型为ARIMA(1,1,0)模型,模型口径为:xt=艾普龙t/1-0.66933B或等阶记为:xt=1.66933xt-1-0.6693
7、3xt-2+艾普龙t(5)“forecast lead=5 id=t;”,利用拟合模型对序列xt作5期预测。一、 建立数据集,绘制时序图data example5_2; input x; lagx=lag(x); t=_n_; cards; 3.03 8.46 10.22 9.80 11.96 2.83 8.43 13.77 16.18 16.84 19.57 13.26 14.78 24.48 28.16 28.27 32.62 18.44 25.25 38.36 43.70 44.46 50.66 33.01 39.97 60.17 68.12 68.84 78.15 49.84 62.2
8、3 91.49 103.20 104.53 118.18 77.8894.75 138.36 155.68 157.46 177.69 117.15 ; proc gplot data=example5_2;plot x*t=1;symbol1 c=black i=join v=star; run; 输出时序图如5-52所示。 图5-52 序列x时序图时序图显示,序列X有一个明显的随时间线性递增的趋势,同时又有一定规律性的波动,所以不妨考虑使用误差自回归模型拟合该序列的发展。2、 因变量关于时间的回归模型proc autoreg data=example5_2;model x=t/ dwpro
9、b;run;语句说明:(1) “proc autoreg data=example5_2;”指令SAS系统对临时数据集example5_2进行回归程序分析。(2) “model x=t/ dwprob;”指令SAS系统以变量t作为自变量,变量x作为因变量,建立线性模型:并给出残差序列 DW检验统计量的分为点。本例中,序列x关于变量t的线性回归模型最小二乘估计输出结果如图5-53所示。图5-53 序列x关于变量t的线性回归模型最小二乘估计结果本例输出结果显示,DW统计量的值等于0.7628,输出概率显示残差序列显著正相关。所以应该考虑对残差序列拟合自相关模型,修改AUTOREG程序如下:proc
10、 autoreg data=example5_2; model x=t/nlag=5 backstep method=ml; run;Model语句是指令系统对线性回归模型 的残差序列显示延迟5阶的自相关图,并拟合延迟5阶自相关模型,特别注意,SAS输出的自回归模型结构为:即输出的自相关回归参数值与我们习惯定义的自回归参数值相差一个负号。 由于自相关延迟阶数的确定是由我们尝试选择的,所以nlag得阶数通常会指得大一些。这就导致残差自回归模型中可能有部分参数不显著,因而添加逐步回归选项backstep,指令系统使用逐步回归的方法筛选出显著自相关因子,并使用极大似然的方法进行参数估计。输出如下四方
11、面的结果:1因变量说明如图5-54所示:图5-54 因变量说明2普通最小二乘估计结果该部分输出信息包括差平方和(SSE)、自由度(DFE)、均方误差(MSE)、根号均方误差(Root MSE)、SBC信息量。回归部分相关系数平方(regress rsquare)、总的相关系数平方(totel rsquare),DW统计量(durbin watson)及所有待参数的自由度、估计值、标准差、t值和统计量的P值。如图5-55所示。图5-55 普通最小二乘估计结果3回归误差分析该部分共输出四方面的信息:残差序列自相关图、逐步回归消除的不显著项报告、初步均方误差(MSE)、自回归参数估计值。本例该部分输
12、出结果如图5-56所示。图5-56 自回归误差分析输出结果本例输出的残差序列自相关图显示残差序列有非常显著的1阶正相关性。逐步回归消除报告显示除了1阶的序列值显著自相关外,延迟其他阶数的序列值均不具有显著的自相关性,因此延迟2阶-5阶的自相关项被剔除。初步均方误差为4、最终拟合模型该部分包括三方面的汇总信息:收敛状况、极大似然估计结果和回归系数估计。本例该部分输出结果如图5-57所示。图5-57 最终拟合模型输出结果本例得到最终拟合模型为:为了得到直观的拟合效果,我们可以利用OUTPUT命令将拟合结果存入SAS数据集中,并对输出结果作图,相关命令如下:proc autoreg data=example5_2; model x=t/nlag=5 backstep method=ml; output out=out p=xp pm=trend;proc autoreg data=example5_2; model x=t/nlag=5 backstep method=ml noint; output out=out p=xp pm=trend; proc gplot data=out;