时间序列分析 第三章平稳时间序列分析Word格式.docx
《时间序列分析 第三章平稳时间序列分析Word格式.docx》由会员分享,可在线阅读,更多相关《时间序列分析 第三章平稳时间序列分析Word格式.docx(18页珍藏版)》请在冰豆网上搜索。
procgplotdata=example3_1;
plotx*time=1;
symbolc=redi=joinv=star;
run;
建立该数据集,绘制该序列时序图得:
根据所得图像,对序列进行平稳性检验.时序图就是一个平面二维坐标图,通常横轴暗示时间,纵轴暗示序列取值.时序图可以直观地帮手我们掌握时间序列的一些基天职布特征.
根据平稳时间序列均值、方差为常数的性质,平稳序列的时序图应该显示出该序列始终在一个常数值附近随机摆荡,而且摆荡的范围有界的特点.如果观察序列的时序图,显示出该序列有明显的趋势性或周期性,那它通常不是平稳序列.从图上可以看出,数值围绕在0附近随机摆荡,没有明显或周期,其本可以视为平稳序列,时序图显示该序列摆荡平稳.
procarimadata=example3_1;
identifyvar=xnlag=8;
图一
图二样本自相关图
图三样本逆自相关图
图四样本偏自相关图
图五纯随机检验图
实验结果分析:
(1)由图一我们可以知道序列样本的序列均值为-0.06595,标准差为1.561613,观察值个数为84个.
(2)根据图二序列样本的自相关图我们可以知道该图横轴暗示自相关系数,综轴暗示延迟时期数,用水平方向的垂线暗示自相关系数的年夜小.我们发现样本自相关图延迟3阶之后,自相关系数都落入2倍标准差范围以内,而且自相关系数向0.03衰减的速度非常快,延迟5阶之后自相关系数即在0.03值附近摆荡.这是一个短时间相关的样本自相关图.所以根据样本自相关图的相关性质,可以认为该序列平稳.
(3)根据图五的检验结果我们知道,在各阶延迟下LB检验统计量的P值都非常小(<
0.0001),所以我们可以以很年夜的掌控(置信水平>
99.999%)判定该序列样本属于非白噪声序列.
identifyvar=xnlag=8minicp=(0:
5)q=(0:
5);
IDENTIFY命令输出的最小信息量结果
某个观察值序列通过序列预处置,可以判定为平稳非白噪声序列,就可以利用ARMA模型对该序列建模.建模的基本步伐如下:
A:
求出该观察值序列的样本自相关系数(ACF)和样本偏自相关系数(PACF)的值.
B:
根据样本自相关系数和偏自相关系数的性质,选择适本地ARMA(p,q)模型进行拟合.
C:
估计模型中未知参数的值.
D:
检验模型有效性.如果拟合模型欠亨过检验,转向步伐B,重新选择模型再拟合.
E:
模型优化.如果拟合模型通过检验,仍然转向步伐B,充沛考虑各种可能,建立多个拟合模型,从所有通过检验中选择最优模型.
F:
利用拟合模型,预测序列的将来走势.
为了尽量防止因个人经验缺乏招致的模型识别问题,SAS系统还提供了相对最优模型识别.最后一条信息显示,在自相关延迟阶数小于即是5,移动平均延迟阶数也小于即是5的所有ARMR(p,q)模型中,BIC信息量相对最小的是ARMR(0,4)模型,即MA(4)模型.
需要注意的是,MINIC只给出一定范围内SBC最小的模型定阶结果,但该模型的参数未必都能通过参数检验,即经常会呈现MINIC给出的模型阶数依然偏高的情况.
estimateq=4;
本例参数估计输出结果显示均值MU不显著(t的检验统计量的P值为0.9968),其他参数均显著(t检验统计量的P值均小于0.00001),所以选择NOINT选项,除去常数项,再次估计未知参数的结果,即可输入第二条ESTIMATE命令:
estimateq=4noint;
参数估计部份输出结果如图六所示:
图六ESTIMATE命令消除常数项之后的输出结果
显然四个未知参数均显著.
拟合统计量的值
这部份输出五个统计量的值,由上到下分别是方差估计值、标准差估计值、AIC信息量、SBC信息量及残差个数,如图七所示:
图七ESTIMATE命令输出的拟合统计量的值
系数相关阵
这部份输出各参数估计值的相关阵,如图八所示:
图八ESTIMATE命令输出的系数相关阵
残差自相关检验结果
这部份的输出格式(图九)和序列自相关系数白噪声检验部份的输出结果一样.本例中由于延迟各阶的LB统计量的P值均显著年夜于a(a=0.05),所以该拟合模型显著成立.
图九ESTIMATE命令输出的残差自相关检验结果
拟合模型的具体形式
ESTIMATE命令输出的拟合模型的形式
序列预测
forecastlead=5id=timeout=results;
其中,lead是指定预测期数;
id是指按时间变量标识;
out是指定预测后的结果存入某个数据集.
该命令运行后输出结果如下:
FORECAST命令输出的预测结果
该输出结果从左到右分别为序列值的序号、预测值、预测值的标准差、95%的置信下限、95%的置信上限.
利用存储在临时数据集RESULTS里的数据,我们还可以绘制漂亮的拟合预测图,相关命令如下:
procgplotdata=results;
plotx*time=1forecast*time=2l95*time=3u95*time=3/overlay;
symbol1c=blacki=nonev=start;
symbol2c=redi=joinv=none;
symbol3c=greeni=joinv=nonel=32;
输出图像如下:
拟合效果图
注:
图中,S号代表序列的观察值;
连续曲线代表拟合序列曲线;
虚线代表拟合序列的95%上下置信限.
所谓预测就是要利用序列以观察到的样本值对序列在未来某个时刻的取值进行估计.目前对平稳序列最经常使用的预测方法是线性最小方差预测.线性是指预测值为观察值序列的线性函数,最小方差是指预测方差到达最小.在预测图上可以看到,数据围绕一个范围内摆荡,即说明未来的数值变动时平稳的.
二、课后习题
第十七题:
根据某城市过去63年中每年降雪量数据(单元:
mm)得:
(书本P94)
法式:
dataexample17_1;
98.355.566.178.4120.597110
procgplotdata=example17_1;
procarimadata=example17_1;
identifyvar=xnlag=15minicp=(0:
estimatep=1;
estimatep=1noin;
(1)判断该序列的平稳性与纯随机性
该序列的时序图如下(图a)
图a
由时序图显示过去63年中每年降雪量数据围绕早70mm附近随机摆荡,没有明显趋势或周期,基本可以看成平稳序列,为了稳妥起见,做了如下自相关图(图b)
图b
时序图就是一个平面二维坐标图,通常横轴暗示时间,纵轴暗示序列取值.时序图可以直观地帮手我们掌握时间序列的一些基天职布特征.根据平稳时间序列均值、方差为常数的性质,平稳序列的时序图应该显示出该序列始终在一个常数值附近随机摆荡,而且摆荡的范围有界的特点.如果观察序列的时序图,显示出该序列有明显的趋势性或周期性,那它通常不是平稳序列.
样本的自相关图我们可以知道该图横轴暗示自相关系数,综轴暗示延迟时期数,用水平方向的垂线暗示自相关系数的年夜小.我们发现样本自相关图延迟2阶之后,自相关系数都落入2倍标准差范围以内,
自相关图显示该序列自相关系数一直都比力小,1阶开始控制在2倍的标准差范围以内,可以认为该序列自始自终都在零轴附近摆荡,这是随即性非常强的平稳时间序列.
纯随机性检验见下图:
(图c)
图c
根据图c的检验结果我们知道,在6阶延迟下LB检验统计量的P值显著小于0.05,所以我们可以以很年夜的掌控(置信水平>
95%)判定这个拟合模型的残差序列属于非白噪声序列.
(2)如果序列平稳且非白躁声,选择适当模型拟合该序列的发展.
模型识别如下图(图d)
图d
假如某个观察值序列通过序列预处置,可以判定为平稳非白噪声序列,就可以利用ARMA模型对该序列建模.建模的基本步伐如下:
1:
2:
3:
4:
5:
6:
最后一条信息显示,在自相数迟阶数小于即是5,移动平均延迟阶数也小于即是5的所有ARMA(p,q)模型中,BIC信息量相对最小的是ARMA(1,0)模型,既AR
(1)模型.它们的自相关系数都呈现出拖尾性和呈指数衰减到零值附近的性质.自相关系数是按负指数单调收敛到零;
利用拟合模型,预测该城市未来5年的降雪量.
由
(2)可以知道该模型是AR
(1)模型;
预测结果如下图(图e)
由图得未来5(64-68年)的降雪量分别为103.6820mm、97.7270mm、92.1139mm、86.8232mm、81.8365mm.
18.某地域连续74年的谷物产量(单元:
千吨)
dataexample18_1;
procgplotdata=example18_1;
procarimadata=example18_1;
identifyvar=xnlag=18minicp=(0:
estimateq=1;
symbol2c=red