基于SAS软件的时间序列实验的代码.docx

1、基于SAS软件的时间序列实验的代码基于SAS软件的时间序列实验的代码 实 验 指 南 实验一 分析太阳黑子数序列3实验二 模拟AR模型4实验三 模拟MA模型和ARMA模型6实验四 分析化工生产量数据8实验五 模拟ARIMA模型和季节ARIMA模型10实验六 分析美国国民生产总值的季度数据13实验七 分析国际航线月度旅客总数数据16实验八 干预模型的建模19实验九 传递函数模型的建模22实验十 回来与时序相结合的建模25太阳黑子年度数据28美国国民收入数据29化工生产过程的产量数据30国际航线月度旅客数据30洛杉矶臭氧每小时读数的月平均值数据31煤气炉数据35芝加哥某食品公司大众食品周销售数据3

2、7牙膏市场占有率周数据39某公司汽车生产数据44加拿大山猫数据44 实验一 分析太阳黑子数序列实验目的：了解时刻序列分析的差不多步骤，熟悉SAS/ETS软件使用方法。二、实验内容：分析太阳黑子数序列。三、实验要求：了解时刻序列分析的差不多步骤，注意各种语句的输出结果。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。创建名为exp1的SAS数据集，即在窗中输入下列语句：data exp1;input a1 ;year=intnx(year,1jan1742d,_n_-1);format year year4.;cards; 输入太阳黑子数序列(见附表)run;

3、储存此步骤中的程序，供以后分析使用（只需按工具条上的储存按钮然后填写完提咨询后就能够把这段程序储存下来即可）。绘数据与时刻的关系图，初步识别序列，输入下列程序：proc gplot data=exp1; symbol i=spline v=star h=2 c=green; plot a1*year;run;提交程序，在graph窗口中观看序列，能够看出此序列是均值平稳序列。识别模型，输入如下程序。 proc arima data=exp1; identify var=a1 nlag=24; run;提交程序，观看输出结果。初步识别序列为AR(3)模型。估量和诊断。输入如下程序： estima

4、te p=3; run;提交程序，观看输出结果。假设通过了白噪声检验，且模型合理，则进行推测。进行推测，输入如下程序：forecast lead=6 interval=year id=year out=out;run;proc print data=out;run;提交程序，观看输出结果。退出SAS系统，关闭运算机。 实验二 模拟AR模型实验目的：熟悉各种AR模型的样本自有关系数和偏有关系数的特点，为理 论学习提供直观的印象。实验内容：随机模拟各种AR模型。实验要求：记录各AR模型的样本自有关系数和偏有关系数，观看各种序列 图形，总结AR模型的样本自有关系数和偏有关系数的特点实验时刻：2小时。

5、实验软件：SAS系统。实验步骤1、开机进入SAS系统。模拟实根情形，模拟过程。在edit窗中输入如下程序： data a; x1=0.5; x2=0.5; n=-50; do i=-50 to 250; a=rannor(32565); x=a-0.6*x1+0.4*x2; x2=x1; x1=x; n=n+1; if i0 then output; end; run; 4、观看输出的数据，输入如下程序，并提交程序。 proc print data=a; var x; proc gplot data=a; symbol i=spline c=red; plot x*n; run;观看样本自有关

6、系数和偏有关系数，输入输入如下程 序，并提交程序。 proc arima data=a; identify var=x nlag=10 outcov=exp1; run; proc gplot data=exp1; symbol i=needle width=6; plot corr*lag; run; proc gplot data=exp1; symbol i=needle width=6; plot partcorr*lag; run;作为作业把样本自有关系数和偏有关系数记录下来。估量模型参数，并与实际模型的系数进行对比，即输入如下程序，并提交。 proc arima data=a; i

7、dentify var=x nlag=10 ; run; estimate p=2; run;模拟虚根情形，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟AR(3)模型，模拟过程。重复步骤3-7即可(但部分程序需要修改，请读者自己完成).10、回到graph窗口观看各种序列图形的异同11、退出SAS系统,关闭运算机. 实验三 模拟MA模型和ARMA模型实验目的：熟悉各种MA模型和ARMA模型的样本自有关系数和偏有关系数 的特点，为理论学习提供直观的印象。实验内容：随机模拟各种MA模型和ARMA模型。实验要求：记录各MA模型和ARMA模型的样本自有关系数和偏有关系数，

8、观看各序列的异同，总结MA模型和ARMA模型的样本自有关系 数和偏有关系数的特点实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。2、模拟情形，模拟过程。在edit窗中输入如下程序： data a; a1=0; a2=0; do n=-50 to 250; a=rannor(32565); x=a+0.65*a1+0.24*a2; a2=a1; a1=a; if n0 then output; end; run; 4、观看输出的数据序列，输入如下程序，并提交程序。 proc gplot data=a;symbol i=spline;plot x*n; run;5、观看样本自有

9、关系数和偏有关系数，输入输入如下程 序，并提交程序。 proc arima data=a; identify var=x nlag=10 outcov=exp1; run; proc gplot data=exp1; symbol1 i=needle c=red; plot corr*lag=1; run; proc gplot data=exp1; symbol2 i=needle c=green; plot partcorr*lag=2; run;作为作业把样本自有关系数和偏有关系数记录下来。估量模型参数，并与实际模型的系数进行对比，即输入如下程序，并提交。 proc arima data

10、=a; identify var=x nlag=10 ; run; estimate q=2; run;模拟情形，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟情形，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟情形，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟ARMA模型，模拟过程。重复步骤3-7即可(但部分程序需要修改，请读者自己完成).回到graph窗口观看各种序列图形的异同。退出SAS系统,关闭运算机. 实验四 分析化工生产量数据实验目的：进一步熟悉时刻序列建模的差不多步骤，把握用SACF及SPACF定 模

11、型的阶的方法。实验内容：分析化工生产过程的产量序列。实验要求：把握ARMA模型建模的差不多步骤，初步把握数据分析技巧。写出 实验报告。实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。创建名为exp2的SAS数据集，即在窗中输入下列语句：data exp2; input x ; n=_n_; cards; 输入化工生产产量数据序列（见附表) ；run;储存此步骤中的程序，供以后分析使用（只需按工具条上的储存按钮然后填写完提咨询后就能够把这段程序储存下来即可）。绘数据与时刻的关系图，初步识别序列，输入下列程序：proc gplot data=exp2; symbol i=spl

12、ine v=star h=2 c=green; plot x*n;run;提交程序，在graph窗口中观看序列，能够看出此序列是均值平稳序列。识别模型，输入如下程序。 proc arima data=exp2; identity var=x nlag=12; run;提交程序，观看输出结果，发觉二阶样本自有关系数和一阶的样本偏有关系数都在2倍的标准差之外，那么我们第一作为一阶AR模型估量，输入如下程序： estimate plot p=1; run;提交程序，观看输出结果，发觉残差能通过白噪声检验，但它的二阶的样本偏有关系数比较大，那么我们考虑二阶AR模型。输入如下程序： estimate p

13、lot p=2; run;提交程序，观看输出结果，发觉残差样本自有关系数和样本偏有关系数都 在2倍的标准差之内。且能通过白噪声检验。比较两个模型的AIC和SBC, 发觉第二个模型的AIC和SBC都比第一个的小，故我们选择第二个模型为 我们的结果。记录参数估量值，写出模型方程式。进行推测，输入如下程序： forecast lead=12 out=out; run; proc print data=out; run;提交程序，观看输出结果。退出SAS系统，关闭运算机。 实验五 模拟ARIMA模型和季节ARIMA模型实验目的：熟悉各种ARIMA模型的样本自有关系数和偏有关系数的特点， 区别各种ARI

14、MA模型的图形，为理论学习提供直观的印象。实验内容：随机模拟各种ARIMA模型。实验要求：记录各ARIMA模型的样本自有关系数和偏有关系数观看各序列 图形的异同，总结ARIMA模型的样本自有关系数和偏有关系数 的特点实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。2、模拟ARIMA(0,1,1)过程，模拟过程。创建数据集，在edit窗中输入如下程序： data a; x1=0.9; a1=0; do n=-50 to 250; a=rannor(32565); x=x1+a-0.8*a1; x1=x; a1=a; if n0 then output; end; run; 4

15、、观看输出的数据序列，输入如下程序：。 proc gplot data=a;symbol i=spline; plot x*n; run; 5、提交程序，在Graph窗口中观看图形。6、观看样本自有关系数和偏有关系数，输入输入如下程序： proc arima data=a; identify var=x nlag=10 outcov=exp1; run; proc gplot data=exp1; symbol1 i=needle c=red; plot corr*lag=1; run; proc plot data=exp1; symbol2 i=needle c=green; plot p

16、artcorr*lag=2; run;提交程序，发觉自有关系数成缓慢下降的趋势，讲明要做差分运算，做一阶差分运算，输入如下程序：proc arima data=a; identity var=x(1) nlag=24; run;提交程序，观看样本自有关系数与样本偏有关系数，发觉自有关系数1阶截尾，故判定差分后序列为MA(1)模型。进行模型参数估量，输入如下程序：estimate q=1 plot;run;提交程序，并观看残差图，发觉模型拟合完全。10、写出模型的方程，并与真实模型对比。11、模拟ARIMA(1,1,0)模型，模拟过程。重复步骤 3-10即可（但部分程序需要修改，请读者自己完成）

17、。模拟模型， 模拟模型, 即模型。13、创建数据集，在edit窗中输入如下程序： data c; x1=0.9;x2=0;x3=0;x4=0;x5=0;x6=0;x7=0; x8=0;x9=0;x10=0;x11=0;x12=0;x13=0; a1=0;a2=0;a3=0;a4=0;a5=0;a6=0;a7=0; a8=0;a9=0;a10=0;a11=0;a12=0;a13=0; do n=-50 to 250; a=rannor(12345); x=x1+x12-x13+a-0.4*a1-0.6*a12+0.24*a13; x13=x12;x12=x11;x11=x10;x10=x9;x9

18、=x8;x8=x7; x7=x6;x6=x5;x5=x4;x4=x3;x3=x2;x2=x1;x1=x; a13=a12;a12=a11;a11=a10;a10=a9;a9=a8;a8=a7; a7=a6;a6=a5;a5=a4;a4=a3;a3=a2;a2=a1;a1=a; if n0 then output; end; run;绘序列图，输入如下程序： proc gplot data=c; symbol i=spline c=red; plot x*n; run;提交程序，到graph窗口中观看序列图形。初步识别模型，输入如下程序： proc arima data=c; identify

19、var=x nlag=36; run;提交程序，观看样本自有关系数和样本偏有关系数。做季节差分和一阶差分除掉季节因子和趋势因子，输入如下程序： identify var=x(1,12) nlag=36; run;提交程序，观看样本自有关系数和样本偏有关系数，确定模型阶数。估量模型参数，输入如下程序： estimate q=(1)(12) method=uls plot; run;提交程序，观看残差的样本自有关系数和样本偏有关系数，看是否通过 了白噪声检验。写出模型方程式，并与真实模型对比。回到graph窗口观看各种序列图形的异同。退出SAS系统,关闭运算机. 实验六 分析美国国民生产总值的季度

20、数据一、实验目的：进一步学习数据分析技巧，进一步了解ARIMA模型。二、实验内容：47年1季度到96年3季度美国国民生产总值的季度数据。三、实验要求：写出分析报告。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp3的SAS数据集，输入如下程序：data exp3; input gnp; date=intnx(qtr,1jan47d,_n_-1); format date yyqc.; cards; 输入美国国民生产总值的数据 ；run; 注：Intnx函数按间隔递增日期，Intnx函数运算某个区间通过若干区间间 隔之后的间隔的开始日期或日期

21、时刻值，其中开始间隔内的一个日期或 日期时刻值给出。 Intnx函数的格式如下： Intnx(interval,from,n)储存上述程序，供以后分析使用（只需按工具条上的储存按钮，然后填写 完提咨询后就能够把这段程序储存下来）。绘序列图，输入如下程序： proc gplot data=exp3; symbol1 i=spline; plot gnp*date=1; run;观看图形，发觉图形成指数函数上升形式，故做对数变换，输入如下程序：data lexp; set exp3; lgnp=log(gnp);run;绘变换后序列图，输入如下程序：proc gplot data=lexp; sy

22、mbol2 i=spline c=red; plot lgnp*date=2;run;提交程序，到graph窗口中观看变换后的序列图，能够看出它成直线上升趋势。对序列做初步识别，输入如下程序：proc arima data=lexp; identify var=lgnp nlag=12; run;提交程序，观看样本自有关系数，可看出有缓慢下降趋势，结合我们观看的图形，我们明白要对序列做差分运算，作一阶差分，输入如下程序：identify var=lgnp(1) nlag=12;run;提交程序，观看样本自有关系数，可看出样本自有关系数5步后是截尾的,那么确定为MA(5)模型,进行参数估量,输入

23、如下程序:estimate q=5 plot;run;提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟合充分。且MA1,3 , MA1,4的T值较小，讲明参数明显为0，除掉这两项重新进行估量，输入如下程序：estimate q=(1,2,5) plot;run;提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟合充分，且残差标准误与前一估量相差专门小，故以此结果为我们所要的结果，依此结果写出方程式。进行推测，推测美国以后2年的每季国民生产总值。输入如下程序： forcast lead=6 interval=qtr id=date out=results; run; d

24、ata results; set results; gnp=exp(lgnp); l95=exp(l95); u95=exp(u95); forecast=exp(forecast+std*std/2); run; proc print data=results; var date forcast; where date=1jan96d; run;提交程序，并把推测值记录下来。退出SAS系统，关闭运算机。 实验七 分析国际航线月度旅客总数数据实验目的：熟悉运用SAS建立模型的方法，进一步 了解模型的特点。二、实验内容：19497年1月至1960年12月国际航线月度旅客总数数据。三、实验要求：写

25、出分析报告。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp4的SAS数据集，输入如下程序：data exp4; input air; date=intnx(month,1jan49d,_n_-1); format date monyy.; cards; 输入国际航线月度旅客总数数据 ；run;储存上述程序，供以后分析使用（只需按工具条上的储存按钮，然后填写 完提咨询后就能够把这段程序储存下来）。绘序列图，输入如下程序： proc gplot data=exp4; symbol1 i=spline v=dot c=red; plot air

26、*date=1; run;提交程序，观看图形，发觉图形有专门强的季节性，且成指数函数上升形式，故做对数变换，输入如下程序：data lair; set exp4; lair=log(air);run;绘变换后序列图，输入如下程序：proc gplot data=lair; symbol2 i=spline c=green; plot lair*date=2;run;提交程序，到graph窗口中观看变换后的序列图，能够看出它总的趋势成直线上升，且有专门强的季节性。对序列做初步识别，输入如下程序：proc arima data=lair; identify var=lair nlag=36; ru

27、n;提交程序，观看样本自有关系数和偏有关系数，可看出样本自有关系数有缓慢下降趋势，偏有关系数在1步，13步，25步较大，我们作一步一阶差分，输入如下程序：identity var=lair(1) nlag=36;run;提交程序，观看样本自有关系数和偏有关系数，发觉样本自有关系数在12步，24步，36步专门大,而偏有关系数在12步专门大，那么我们再做12步的一阶差分，输入如下程序：identify var=lair(1,12) nlag=36;run;10、提交程序，观看样本自有关系数和偏有关系数，发觉样本自有关系数在1步，12步专门大,而偏有关系数看不出有专门的规律，我们可确定模型的MA因子为。11、进行参数估量，输入如下程序： estimate q=(1)（12）noconstant method=uls plot; run;、提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟