应 用 时 间 序 列 分 析 实 验 报 告
实验名称 第二章 时间序列的预处理
一、上机练习
2.4.1绘制时序图
data example2_1;
input price1 price2;
time=intnx('month','01jul2004'd,_n_-1); format time date.; cards; 12.85 15.21 13.29 14.23 12.41 14.69 15.21 13.27 14.23 16.75 13.56 15.33 ;
proc gplot data=example2_1;
plot price1*time=1 price2*time=2/overlay; symbol1 c=black v=star i=join; symbol2 c=red v=circle i=spline; run;
语句说明：
（1）“proc gplot data=example2_1；”是告诉系统，下面准备对临时数据集example2_1中的数据绘图。
（2）“plot price1*time=1 price2*time=2/overlay;”是要求系统要绘制两条时序曲线。 （3）“symbol1 c=black v=star i=join;”，symbol语句是专门指令绘制的格式。
输出的时序图见下图：
1
两时间序列重叠显示时序图 2.4.2 平稳性与纯随机性检验 1、平稳性检验 为了判断序列是否平稳，除了需要考虑时序图的性质，还需要对自相关图进行检验。SAS系统ARIMA过程中的IDENTIFY语句可以提供非常醒目的自相关图。 data example2_2; input freq@@; year=intnx ('year','1jan1970'd,_n_-1); format year year4.; cards; 97 154 137.7 149 164 157 188 204 179 210 202 218 209 204 211 206 214 217 210 217 219 211 233 316 221 239 215 228 219 239 224 234 227 298 332 245 357 301 389 ; proc arima data=example2_2; identify var=freq; run; 语句说明： （1）“proc arima data=example2_2;”是告诉系统，下面要对临时数据集example2_2中的数据进行ARIMA程序分析。 （2）“identify var=freq;”是对指令变量freq的某些重要性质进行识别。 2
执行本例程序，IDENTIFY语句输出的描述性信息如下： 这部分给出了分析变量的名称、序列均值、标准差和观察值个数。 IDENTIFY语句输出结果的第二部分分为自相关图，本例获得的样本自相关见下图。 序列FREQ样本自相关图 其中： Lag——延迟阶数。 Covariance——延迟阶数给定后的自协方差函数。 Correlation——自相关系数的标准差。 “.”——2倍标准差范围。 2、纯随机性检验 为了判断序列是否有分析价值，我们必须对序列进行纯随机性检验，即白噪声检验。在IDENTIFY输出结果的最后一部分信息就是白噪声检验结果。本例中白噪声检验输出结果如下： 其中： To Lag——延迟阶数。 检验结果显示，在6阶延迟下LB检验统计量的P值非常小（<0.0001），所以我们可以以很大的把握（置信水平>99.999%）断定该序列属于非白噪声序列。 3
二、课后习题
2.1975-1980年夏威夷岛莫那罗亚火山（Mauna loa）每月释放的co2数据如下（单位：ppm），见表2-7.
330.45 331.9 331.63 333.05 332.81 334.65 334.66 336.25 335.89 337.41 337.81 339.25
330.97 330.05 332.46 330.87 333.23 332.41 335.07 334.39 336.44 335.71 338.16 337.19
331.64 328.58 333.36 329.24 334.55 331.32 336.33 332.44 337.63 333.68 339.88 335.49
332.87 328.31 334.45 328.87 335.82 330.73 337.39 332.25 338.54 333.69 340.57 336.63
333.61 329.41 334.82 330.18 336.44 332.05 337.65 333.59 339.06 335.05 341.19 337.74
333.55 330.63 334.32 331.5 335.99 333.53 337.57 334.76 338.95 336.53 340.87 338.36
（1）绘制序列时序图，并判断该系列是否平稳。
实验程序：
data example2_1; input ppm@@;
time=intnx('month','01jan1975'd,_n_-1); format time date.; cards; 330.45 330.97 331.64 332.87 333.61 333.55 331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 ;
proc gplot data=example2_1; plot ppm*time=1;
symbol1 c=black v=star i=join;
run;
4
实验结果：
实验分析体会：
时序图给我们的提供的信息非常明确，夏威夷岛莫那罗亚火山（Mauna loa）每月释放的co2时间序列图有明显的递增趋势，所以它不是平稳序列。
?k(k?1,2,（2）计算该序列的样本自相关系数?,24)。
实验程序：
data example2_1; input ppm@@;
time=intnx('month','01jan1975'd,_n_-1); format time date.; cards; 330.45 330.97 331.64 332.87 333.61 333.55 331.90 330.05 328.58 328.31 329.41 330.63 331.63 332.46 333.36 334.45 334.82 334.32 333.05 330.87 329.24 328.87 330.18 331.50 332.81 333.23 334.55 335.82 336.44 335.99 334.65 332.41 331.32 330.73 332.05 333.53 334.66 335.07 336.33 337.39 337.65 337.57 336.25 334.39 332.44 332.25 333.59 334.76 335.89 336.44 337.63 338.54 339.06 338.95 337.41 335.71 333.68 333.69 335.05 336.53 337.81 338.16 339.88 340.57 341.19 340.87 339.25 337.19 335.49 336.63 337.74 338.36 ;
5