时间序列分析第三章平稳时间序列分析.docx
《时间序列分析第三章平稳时间序列分析.docx》由会员分享，可在线阅读，更多相关《时间序列分析第三章平稳时间序列分析.docx（24页珍藏版）》请在冰点文库上搜索。
时间序列分析第三章平稳时间序列分析
应用时间序列分析实验报告
实验名称第三章平稳时间序列分析
根据所得图像，对序列进行平稳性检验。
时序图就是一个平面二维坐标图，通常横轴表示时间，纵
轴表示序列取值。
时序图可以直观地帮助我们掌握时间序列的一些基本分布特征。
根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界的特点。
如果观察序列的时序图，显示出该序列有明显的趋势性或周期性，那它通常不是平稳序列。
从图上可以看出，数值围绕在0附近随机波动，没有明显或周期,
其本可以视为平稳序列，时序图显示该序列波动平稳。
procarimadata=example3_1;
dentifyvar=xnlag=8;
run;
图一
图二样本自相关图
图四样本偏自相关图
图五纯随机检验图
实验结果分析：
（1）由图一我们可以知道序列样本的序列均值为，标准差为，观察值个数为84个。
（2）根据图二序列样本的自相关图我们可以知道该图横轴表示自相关系数，综轴表示延迟时期数，
用水平方向的垂线表示自相关系数的大小。
我们发现样本自相关图延迟3阶之后，自相关系数都落入2
倍标准差范围以内，而且自相关系数向衰减的速度非常快，延迟5阶之后自相关系数即在值附近波动。
这是一个短期相关的样本自相关图。
所以根据样本自相关图的相关性质，可以认为该序列平稳。
（3）根据图五的检验结果我们知道，在各阶延迟下LB检验统计量的P值都非常小（＜）,所以我们
可以以很大的把握（置信水平>%断定该序列样本属于非白噪声序列。
procarimadata=example3_1;
identifyvar=xnlag=8minicp=（0:
5）q=（0:
5）;
run;
IDENTIFY命令输出的最小信息量结果
某个观察值序列通过序列预处理，可以判定为平稳非白噪声序列，就可以利用ARMA莫型对该序列
建模。
建模的基本步骤如下：
A:
求出该观察值序列的样本自相关系数（ACF和样本偏自相关系数（PACF的值。
B：
根据样本自相关系数和偏自相关系数的性质，选择适当地ARMA（p，q）模型进行拟合。
C:
估计模型中未知参数的值。
D:
检验模型有效性。
如果拟合模型不通过检验，转向步骤B,重新选择模型再拟合。
E:
模型优化。
如果拟合模型通过检验，仍然转向步骤B,充分考虑各种可能，建立多个拟合模型，
从所有通过检验中选择最优模型。
F:
利用拟合模型，预测序列的将来走势。
为了尽量避免因个人经验不足导致的模型识别问题，SAS系统还提供了相对最优模型识别。
最后一
条信息显示，在自相关延迟阶数小于等于5，移动平均延迟阶数也小于等于5的所有ARMR（p,q）模型
中，BIC信息量相对最小的是ARMFR0,4）模型，即MA（4）模型。
需要注意的是，MINIC只给出一定范围内SBC最小的模型定阶结果，但该模型的参数未必都能通过参数检验，即经常会出现MINIC给出的模型阶数依然偏高的情况。
estimateq=4;run;
本例参数估计输出结果显示均值MU不显著（t的检验统计量的P值为），其他参数均显著（t检验
统计量的P值均小于），所以选择NOINT选项，除去常数项，再次估计未知参数的结果，即可输入第二
条ESTIMATE命令：
estimateq=4noint;
run;
参数估计部分输出结果如图六所示：
图六ESTIMATE命令消除常数项之后的输出结果显然四个未知参数均显著。
拟合统计量的值
SBC信息
这部分输出五个统计量的值，由上到下分别是方差估计值、标准差估计值、AIC信息量、
量及残差个数，如图七所示:
图七ESTIMATE命令输出的拟合统计量的值
系数相关阵
这部分输出各参数估计值的相关阵，如图八所示:
图八ESTIMATE命令输出的系数相关阵
残差自相关检验结果
这部分的输出格式（图九）和序列自相关系数白噪声检验部分的输出结果一样。
本例中由于延迟各阶的LB统计量的P值均显著大于a（a=）,所以该拟合模型显著成立。
图九ESTIMATE命令输出的残差自相关检验结果
拟合模型的具体形式
ESTIMATE命令输出的拟合模型的形式
序列预测
forecastlead=5id=timeout=results;
run;
其中，lead是指定预测期数；id是指定时间变量标识；out是指定预测后的结果存入某个数据集。
该命令运行后输出结果如下：
FORECAS命令输出的预测结果
95%勺置信
该输出结果从左到右分别为序列值的序号、预测值、预测值的标准差、95%勺置信下限、
上限。
利用存储在临时数据集RESULT里的数据，我们还可以绘制漂亮的拟合预测图，相关命令如下:
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;
run;
输出图像如下：
拟合效果图
注：
图中，S号代表序列的观察值；连续曲线代表拟合序列曲线；虚线代表拟合序列的95%上下置信限。
所谓预测就是要利用序列以观察到的样本值对序列在未来某个时刻的取值进行估计。
目前对平稳序列最常用的预测方法是线性最小方差预测。
线性是指预测值为观察值序列的线性函数，最小方差是指预测方差达到最小。
在预测图上可以看到，数据围绕一个范围内波动，即说明未来的数值变化时平稳的。
二、课后习题
第十七题：
根据某城市过去63年中每年降雪量数据（单位：
mm得：
（书本P94）
程序：
dataexample17_1;
inputx@@;
time=_n_;
cards;
97110
;
procgplotdata=example171;
plotx*time=1;
symbolc=redi=joinv=star;
run;
proc
arimadata=example171;
identifyvar=xnlag=15minicp=（0:
5）q=（0:
5）;
run;
estimate
p=1;
run;
estimaterun;
p=1noin;
forecast
lead=5id=timeout=results;
run;
procgplotdata=results;
plotx*time=1forecast*time=2l95*time=3u95*time=3/overlaysymbol1c=blacki=nonev=start;
symbol2c=redi=joinv=none;
symbol3c=greeni=joinv=nonel=32;
run;
（1）判断该序列的平稳性与纯随机性
该序列的时序图如下（图a）
由时序图显示过去63年中每年降雪量数据围绕早70mm附近随机波动，没有明显趋势或周期,
基本可以看成平稳序列，为了稳妥起见，做了如下自相关图（图b）
AutocorreIIati
0
553.6闕
1.00000
0
1
16^.516
0・$闸5
■
常瘴出HC弗舟
1259阴
2
163.778
0,29679
■
0.137290
3
2D.G77781
0.03724
0.1470EB
70.149&06
D.12669
附常舟a
IL147208
5
^.374730
D.00BSB
J.14B923
e
4.甜日
a.14前腮
7
-39.057&37
-.07054
屮
0.I4S94S
站
-SB.4B6733
-.ID6E80
J.14J474
9
.044Q0
屮
0.143935
1C
9.1G723S
0.01475
V.W018G
11
31.453617
0.05661
*:
0.1502IS
12
42.071531
0.07680
a.1DODSO
13
ED.2EG9D4
O.D^O7B
M■
a.1E11B7
14
63.271771
0.11J27
a.ieeoao
11
65.555370
D.1IS3B
J.IG33&7
LagCov«ltlance
Uorreliticr»
StdError
narksIwostiririarderrors
时序图就是一个平面二维坐标图，通常横轴表示时间，纵轴表示序列取值。
时序图可以直观地帮助我们掌握时间序列的一些基本分布特征。
根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界的特点。
如果观察序列的时序图，显示出该序列有明显的趋势性或周期性，那它通常不是平稳序列。
样本的自相关图我们可以知道该图横轴表示自相关系数，综轴表示延迟时期数，用水平方向的
垂线表示自相关系数的大小。
我们发现样本自相关图延迟2阶之后，自相关系数都落入2倍标准差范围
以内，
自相关图显示该序列自相关系数一直都比较小，1阶开始控制在2倍的标准差范围以内，可以认为
该序列自始自终都在零轴附近波动，这是随即性非常强的平稳时间序列。
纯随机性检验见下图：
（图c）
4ut«：
QiirrelationCheckforWhlieNoIse
To
Chk
Pr>
Sqjftrs
DF
ChiSq
■■hutocorre1ations
6
13.28
6
u』朋卩
0^06
0.286
0.037
（J.127
0.008
0.008
12
U.88
IE
C.2474
-0.071
-C.066
-B.060
Q-01&
0.057
0.C76
图c
根据图c的检验结果我们知道，在6阶延迟下LB检验统计量的P值显著小于，所以我们可以以很
大的把握（置信水平>95%断定这个拟合模型的残差序列属于非白噪声序列。
（2）如果序列平稳且非白躁声，选择适当模型拟合该序列的发展。
模型识别如下图（图d）
MnimumImformationCriterion
Error3t?
riesmodeIARC10）
MiniihuhTable★玄丨“士EICC1,0）=E.91EG7C
图d
假如某个观察值序列通过序列预处理，可以判定为平稳非白噪声序列，就可以利用ARMA莫型对该
序列建模。
建模的基本步骤如下：
1:
求出该观察值序列的样本自相关系数（ACF和样本偏自相关系数（PACF的值。
2:
根据样本自相关系数和偏自相关系数的性质，选择适当地ARMA（p,q）模型进行拟合。
3:
估计模型中未知参数的值。
4:
检验模型有效性。
如果拟合模型不通过检验，转向步骤B，重新选择模型再拟合。
5:
模型优化。
如果拟合模型通过检验，仍然转向步骤B,充分考虑各种可能，建立多个拟合模型，
从所有通过检验中选择最优模型。
6:
利用拟合模型，预测序列的将来走势。
最后一条信息显示，在自相数迟阶数小于等于5，移动平均延迟阶数也小于等于5的所有ARMAp,q）
模型中，BIC信息量相对最小的是ARMA（1，0）模型，既AR
（1）模型。
它们的自相关系数都呈现出拖尾性和呈指数衰减到零值附近的性质。
自相关系数是按负指数单调收敛到零；
利用拟合模型，预测该城市未来5年的降雪量•
由
（2）可以知道该模型是AR
（1）模型；
预测结果如下图（图e）
riheARIMIAProcedure
Forecast®farriabIx
由图得未来5（64-68年）的降雪量分别为、、
18.某地区连续74年的谷物产量（单位：
千吨）
procgplotdata=example18_1;
plotx*time=1;
symbolc=redi=joinv=star;
run;
procarimadata=example18_1;
identifyvar=xnlag=18minicp=（0:
5）q=（0:
5）;run;
estimateq=1;
run;
forecastlead=5id=timeout=results;
run;
procgplotdata=results;
plotx*time=1forecast*time=2I95*time=3u95*time=3/overlay
symbollc=blacki=nonev=start;symbol2c=redi=joinv=none;
symbol3c=greeni=joinv=nonel=32;run;
（1）判断该序列的平稳性与纯随机性
该序列的时序图如下（图f）
时序图就是一个平面二维坐标图，通常横轴表示时间，纵轴表示序列取值。
时序图可以直观地帮助我们掌握时间序列的一些基本分布特征。
根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界的特点。
如果观察序列的时序图，显示出该序列有明显的趋势性或周期性，那它通常不是平稳序列。
由时序图显示过去74年中每年谷物产量数据围绕早千吨附近随机波动，没有明显趋势或周期,基本可以看成平稳序列，为了稳妥起见，做了如下自相关图（图g）
图g
样本的自相关图我们可以知道该图横轴表示自相关系数，综轴表示延迟时期数，用水平方向的垂
线表示自相关系数的大小。
我们发现样本自相关图延迟2阶之后，自相关系数都落入2倍标准差范围以
内，自相关图显示该序列自相关系数一直都比较小，1阶开始控制在2倍的标准差范围以内，可以认为
该序列自始自终都在零轴附近波动，这是随即性非常强的平稳时间序列。
纯随机性检验见下图：
（图h）
i+.ccarrelaiionChe&k「nrWhiteNciss
Chi-Pr>
SquareDPLhiSq
图h
根据图h的检验结果我们知道，在各阶延迟下LB检验统计量的P值显著小于，所以我们可以以很
大的把握（置信水平>95%断定这个拟合模型的残差序列属于非白噪声序列。
选择适当模型拟合该序列的发展。
如果序列平稳且非白躁声，选折适当模型拟合序列的发展
模型识别如下图（图
網0
也1
MAt
NA3
MA4
MA5
-2.72095
-2.75017
-2.7^0C2
-2-B77B9
-2,89896
-2.
61S06
2.02ODQ
E.G1621
2.70125
2.74130
2.892^7
G3O72
-2.78813
-2.7052
-E.68383
-2,63746
-z
*5844
75662
-2,73983
-2.66167
-2.A4734
-2,81778
-2
，站孙
-2.79089
-2.G4U1
-2,C143G
-3.68906
-2.
&^777
-2+63881
-2.64Z7f
-2.53317
-£.06433
-E.S400B
-2.
50554
The：
APIMAProceduro
Minimile1nform^tinnCri+srinn
Errcrserieshodel:
（6）
MinimunTablcLluo：
DIC（1-0
2,02003
u®alZ345s
aRRRRRRLAAAAAA
假如某个观察值序列通过序列预处理，可以判定为平稳非白噪声序列，就可以利用ARMA莫型对该
序列建模。
建模的基本步骤如下：
A:
求出该观察值序列的样本自相关系数（ACF和样本偏自相关系数（PACF的值。
B：
根据样本自相关系数和偏自相关系数的性质，选择适当地ARMA（p,q）模型进行拟合。
C:
估计模型中未知参数的值。
D:
检验模型有效性。
如果拟合模型不通过检验，转向步骤B，重新选择模型再拟合。
E:
模型优化。
如果拟合模型通过检验，仍然转向步骤B,充分考虑各种可能，建立多个拟合模型，
从所有通过检验中选择最优模型。
F:
利用拟合模型，预测序列的将来走势。
最后一条信息显示，在自相数迟阶数小于等于5，移动平均延迟阶数也小于等于5的所有ARMAp,q）
模型中，BIC信息量相对最小的是ARMA（1，0）模型，既AR
（1）模型。
它们的自相关系数都呈现出拖尾性和呈指数衰减到零值附近的性质。
自相关系数是按负指数单调收敛到零；
利用拟合模型，预测该地区未来5年的谷物产量，预测结果如下图（图j）
由
（2）可知，该模型为AR
（1）模型；
SAS羸歸20的年（18月16E
TheARIMAPro亡色臼Li广色
Forecast孝forvariabIex
Obc
Fcrocwot
SidError
師篝CbniidoncoLimiic
7&
0.?
fc43
0.290（11.禹CM
75
D・2937
0,27611-
77
0.6513
D.2S3?
0.2?
611.4275
?
0
o.nr.iit
nr?
ti07
n.p7&11.4S75
7S
n.?
SFI7
0.P7011.4775
图j未来5年的谷物产量一次为，，，。
19.现有201个连续的生产记录
procgplotdata=example19_1;plotx*time=1;
symbolc=redi=joinv=star;
run;
procarimadata=example19_1;
identify
var=xnlag=24minicp=（0:
5）q=（0:
5）;
run;
estimate
q=1;]
run;
forecast
lead=5id=timeout=results;
run;
procgplot
data=results;
plotx*time=1forecast*time=2I95*time=3u95*time=3/overlay
symbol1c=blacki=nonev=start;
symbol2
c=redi=joinv=none;
symbol3c=greeni=joinv=nonel=32;
run;
（1）判断该序列的平稳性与纯随机性
该序列的时序图如下（图k）
由时序图显示过去201个连续的生产记录数据围绕早84附近随机波动，没有明显趋势或周期,
基本可以看成平稳序列，为了稳妥起见，做了如下自相关图（图I）
LcleCavariblnee
Autocorrelaiions
S.387843
i.Donno
»iXrJ*^jLp■le-rrlt""■ib■»■1■«irsp^■■■!
■
0
-2.504863
-1.010039
-.28528
0.070535
1202?
s**
0.076568
-0^98701
-.04743
■弗
0.0774SB
0.305540
ft.10705
H叽
0.077002
-0.601629
-・0SB44
0.078374
1,324528
0.15772
船啊
0.078950
-0.486367
-.05799
O.0SO502
0.11«74S
U.191414
U.090/10
-0.B18014
-.06160
:
诽
0.090722
0.390769
0.D4GG3
Q・OQOBE?
0.00IE909
Q.00019
0.0010^0
-0.E15242
-.06135
0.081090
O.OOSSEflE
0.00111
0.081321
Q*76JJ624
ojaaaa
粧：
0.081921
-10.487332
-.05523
BSB
0.061B2S
0.682647
ft.1）6541
Hi.
0.002030
-0.46W9
-.0E58E
n.002282
-0.476164
-.0E668
■昭
0.092467
1.103007
ft.13134
0.022骷0
-Q//IQZV3
-.D85S3
■
（LUS3692
0.B80072
0.03684
Hi*,
0.00^125
-IQ.084-0IE
00418
■
0.094701
-1.（11G329
121C2
・3fe弗
D.0S4702
0.44304?
Q.Q5286
犒・
0.005556
ntiarkstwostandarderrors
CorreIation
-190765432101234567091
SldlError
时序图就是一个平面二维坐标图，通常横轴表示时间，纵轴表示序列取值。
时序图可以直观地帮助我们掌握时间序列的一些基本分布特征。
根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界的特点。
如果观察序列的时序图，显示出该序列有明显的趋势性或周期性，那它通常不是平稳序列。
样本的自相关图我们可以知道该图横轴表示自相关系数，综轴表示延迟时期数，用水平方向的
垂线表示自相关系数的大小。
我们发现样本自相关图延迟1阶之后，自相关系数都落入2
升级会员 