时间序列分析第二章-时间序列的预处理3936

应用时间序列分析实验报告实验名称第二章时间序列的预处理一、上机练习2.4.1绘制时序图dataexample2_1;inputprice1price2;time=intnx('month','01jul2004'd,_n_-1);formattimedate.;cards;12.8515.2113.2914.2312.4114.6915.2113.2714.2316.7513.5615.33;procgplotdata=example2_1;plotprice1*time=1price2*time=2/overlay;symbol1c=blackv=stari=join;symbol2c=redv=circlei=spline;run;语句说明：（1）“procgplotdata=example2_1；”是告诉系统，下面准备对临时数据集数据绘图。example2_1中的（2）“plotprice1*time=1price2*time=2/overlay;”是要求系统要绘制两条时序曲线。（3）“symbol1c=blackv=stari=join;”，symbol语句是专门指令绘制的格式。输出的时序图见下图：1两时间序列重叠显示时序图2.4.2平稳性与纯随机性检验1、平稳性检验为了判断序列是否平稳，除了需要考虑时序图的性质，还需要对自相关图进行检验。SAS系统ARIMA过程中的IDENTIFY语句可以提供非常醒目的自相关图。dataexample2_2;inputfreq@@;year=intnx('year','1jan1970'd,_n_-1);formatyearyear4.;cards;97154137.7149164157188204179210202218209204211206214217210217219211233316221239215228219239224234227298332245357301389;procarimadata=example2_2;identifyvar=freq;run;语句说明：（1）“procarimadata=example2_2;”是告诉系统，下面要对临时数据集example2_2中的数据进行ARIMA程序（2）“identifyvar=freq;”是对指令变量freq的某些重要性质进行分析。识别。2执行本例程序，IDENTIFY语句输出的描述性信息如下：这部分给出了分析变量的名称、序列均值、标准差和观察值个数。IDENTIFY语句输出结果的第二部分分为自相关图，本例获得的样本自相关见下图。序列FREQ样本自相关图其中：Lag——延迟阶数。Covariance——延迟阶数给Correlation——自相关“.”——2倍标准差2、纯随机性检验定后的自协方差函数。系数的标准差。范围。为了判断序列是否有分析价值，我们必须对序列进行纯随机性检验，即白噪声检验。在IDENTIFY输出结果的最后一部分信息就是白噪声检验结果。本例中白噪声检验输出结果如下：其中：ToLag——延迟阶数。在6阶延迟下LB检验统计量的P值非常小（<0.0001），所以我们可以以很大的把握（置信水平>99.999%）断定该序列检验结果显示，属于非白噪声序列。3二、课后习题co2.1975-1980年夏威夷岛莫那罗亚火山（Maunaloa）每月释放的数据如下（单位：ppm），见表2-7.2330.45331.9330.97330.05332.46330.87333.23332.41335.07334.39336.44335.71338.16337.19331.64328.58333.36329.24334.55331.32336.33332.44337.63333.68339.88335.49332.87333.61329.41334.82330.18336.44332.05337.65333.59339.06335.05341.19337.74333.55330.63334.32331.5328.31334.45328.87335.82330.73337.39332.25338.54333.69340.57336.63331.63333.05332.81334.65334.66336.25335.89337.41337.81339.25335.99333.53337.57334.76338.95336.53340.87338.36（1）绘制序列时序图，并判断该系列是否平稳。实验程序：dataexample2_1;inputppm@@;time=intnx('month','01jan1975'd,_n_-1);formattimedate.;cards;330.45330.97331.64332.87333.61333.55331.90330.05328.58328.31329.41330.63331.63332.46333.36334.45334.82334.32333.05330.87329.24328.87330.18331.50332.81333.23334.55335.82336.44335.99334.65332.41331.32330.73332.05333.53334.66335.07336.33337.39337.65337.57336.25334.39332.44332.25333.59334.76335.89336.44337.63338.54339.06338.95337.41335.71333.68333.69335.05336.53337.81338.16339.88340.57341.19340.87339.25337.19335.49336.63337.74338.36;procgplotdata=example2_1;plotppm*time=1;symbol1c=blackv=stari=join;run;4实验结果：实验分析体会：co时间序2时序图给我们的提供的信息非常明确，夏威夷岛莫那罗亚火山（列图有明显的递增趋势，所以它不是平稳序列。Maunaloa）每月释放的ˆ(k1,2,,24)。（2）计算该序列的样本自相关系数k实验程序：dataexample2_1;inputppm@@;time=intnx('month','01jan1975'd,_n_-1);formattimedate.;cards;330.45330.97331.64332.87333.61333.55331.90330.05328.58328.31329.41330.63331.63332.46333.36334.45334.82334.32333.05330.87329.24328.87330.18331.50332.81333.23334.55335.82336.44335.99334.65332.41331.32330.73332.05333.53334.66335.07336.33337.39337.65337.57336.25334.39332.44332.25333.59334.76335.89336.44337.63338.54339.06338.95337.41335.71333.68333.69335.05336.53337.81338.16339.88340.57341.19340.87339.25337.19335.49336.63337.74338.36;5procarimadata=example2_1;identifyvar=ppm;run;实验结果:Correlation0.907510.721710.512520.349820.24690.203090.210210.264290.364330.484720.584560.601980.518410.368560.206710.081380.00135-0.03248（3）绘制该样本自相关图，并解释该图形。6自相关图显示序列子相关系数长期位于零轴的一边，这是具有单调趋势序列的典型特征，同时自相关图呈现出明显的正弦波动规律，这是具有周期变化规律的非平稳序列的典型特征。自相关图显示出来的这两个性质和该序列时序图显示的带长期递增趋势的周期性质是非常吻合的。3.1945-1950年费城月度降雨量数据如下（单位：mm）实验程序：dataexample2_3;inputfreq@@;time=intnx('month','1jan1945'd,_n_-1);formattimedate.;cards;69.380.040.974.984.6101.1225.095.3100.648.3144.5128.338.452.368.637.1148.6218.7131.6112.881.831.047.570.196.861.555.6171.7220.5119.463.2181.673.964.8166.948.0137.780.5105.289.9174.8124.086.4136.931.535.3112.3143.0160.897.080.562.5158.27.6165.9106.792.263.226.277.052.3105.4114.349.5116.154.1148.6159.385.367.3112.859.4;procarimadata=example2_3;identifyvar=freq;run;自相关图：7ˆ(k1,2,,24)。（1）计算该序列的样本自相关系数k从上面的自相关图可以看出样本的自相关系数为Correlation0.06005-0.04326-0.09752-0.21647-0.13151-0.057300.03337-0.09036-0.002320.025200.169910.02973-.16785-.15233-.183190.080968（2）判断该序列的平稳性。如下图是该序列的时序图：根据序列图可以知道，图上可以看出该序列在一个常值附近上下波动，且不具有周期性，判断该序列为平稳序列。（3）判断该序列的纯随机性。本序列的检验结果如下：由于P值显著大于显著性水平0.05，所以该序列不能拒绝纯随机的原假设。因而可以认为费城月度降雨量的变动属于纯随机波动。95.表2-9数据是某公司在2000-2003年期间每月的销售量。（1）绘制该序列时序图及样本自相关图。实验程序：dataexample2_3;inputnumber@@;time=intnx('month','1jan2000'd,_n_-1);formattimeyymmdd10.;cards;15318723421230022120117512310485781341752432272982562371651241068774145203189214295220231174119856775117178149178248202162135120969063;procgplotdata=example2_3;plotnumber*time=1;symbol1c=blackv=stari=join;procarimadata=e