经济时间序列分各种模型分析

1、 目 录实验一 分析太阳黑子数序列3实验二 模拟AR模型4实验三 模拟MA模型和ARMA模型6实验四 分析化工生产量数据8实验五 模拟ARIMA模型和季节ARIMA模型10实验六 分析美国国民生产总值的季度数据13实验七 分析国际航线月度旅客总数数据16实验八 干预模型的建模19实验九 传递函数模型的建模22实验十 回归与时序相结合的建模25太阳黑子年度数据28美国国民收入数据29化工生产过程的产量数据30国际航线月度旅客数据30洛杉矶臭氧每小时读数的月平均值数据31煤气炉数据35芝加哥某食品公司大众食品周销售数据37牙膏市场占有率周数据39某公司汽车生产数据44加拿大山猫数据44 实验一 分

2、析太阳黑子数序列一、实验目的：了解时刻序列分析的差不多步骤，熟悉SAS/ETS软件使用方法。二、实验内容：分析太阳黑子数序列。三、实验要求：了解时刻序列分析的差不多步骤，注意各种语句的输出结果。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。创建名为exp1的SAS数据集，即在窗中输入下列语句：保存此步骤中的程序，供以后分析使用（只需按工具条上的保存按钮然后填写完提问后就能够把这段程序保存下来即可）。绘数据与时刻的关系图，初步识不序列，输入下列程序：ods html;ods listing close;run;提交程序，在graph窗口中观看序列，能够看出此

3、序列是均值平稳序列。识不模型，输入如下程序。提交程序，观看输出结果。初步识不序列为AR(2)模型。可能和诊断。输入如下程序：提交程序，观看输出结果。假设通过了白噪声检验，且模型合理，则进行预测。进行预测，输入如下程序：提交程序，观看输出结果。退出SAS系统，关闭计算机。总程序：data exp1;infile D:exp1.txt;input a1 ;year=intnx(year,1jan1742d,_n_-1);format year year4.;proc print;run;ods html;ods listing close;proc gplot data=exp1 ; symbol

4、 i=spline v=dot h=1 cv=red ci=green w=1; plot a1*year/autovref lvref=2 cframe=yellow cvref=black ; title 太阳黑子数序列;run; proc arima data=exp1; identify var=a1 nlag=24 minic p=(0:5) q=(0:5); estimate p=3; forecast lead=6 interval=year id=year out=out;run;proc print data=out;run;选取拟合模型的规则:1.模型显著有效(残差检验为白

5、噪声)2.模型参数尽可能少3.结合自相关图和偏自相关图以及minic条件(BIC信息量最小原则),选取显著有效的参数 实验二 模拟AR模型实验目的：熟悉各种AR模型的样本自相关系数和偏相关系数的特点，为理 论学习提供直观的印象。实验内容：随机模拟各种AR模型。实验要求：记录各AR模型的样本自相关系数和偏相关系数，观看各种序列 图形，总结AR模型的样本自相关系数和偏相关系数的特点实验时刻：2小时。实验软件：SAS系统。实验步骤1、开机进入SAS系统。模拟实根情况，模拟过程。在edit窗中输入如下程序： 4、观看输出的数据，输入如下程序，并提交程序。观看样本自相关系数和偏相关系数，输入输入如下程序

6、，并提交程序。作为作业把样本自相关系数和偏相关系数记录下来。可能模型参数，并与实际模型的系数进行对比，即输入如下程序，并提交。模拟虚根情况，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟AR(3)模型，模拟过程。重复步骤3-7即可(但部分程序需要修改，请读者自己完成).10、回到graph窗口观看各种序列图形的异同11、退出SAS系统,关闭计算机.总程序:title;data a; x1=0.5; x2=0.5; do i=-50 to 250; a=rannor(32565); x=a-0.6*x1+0.4*x2; x2=x1; x1=x; output; end;

7、 run;proc print data=a; var x; proc gplot data=a;symbol i=spline c=red; plot x*i/haxis=-50 to 255 by 20; run;quit;proc arima data=a; identify var=x nlag=10 minic p=(0:3) q=(0:3) outcov=exp1; estimate p=2 noint;run; proc gplot data=exp1;symbol i=needle width=6;plot corr*lag;run;proc gplot data=exp1;s

8、ymbol i=needle width=6;plot partcorr*lag;run; 实验三 模拟MA模型和ARMA模型实验目的：熟悉各种MA模型和ARMA模型的样本自相关系数和偏相关系数 的特点，为理论学习提供直观的印象。实验内容：随机模拟各种MA模型和ARMA模型。实验要求：记录各MA模型和ARMA模型的样本自相关系数和偏相关系数， 观看各序列的异同，总结MA模型和ARMA模型的样本自相关系 数和偏相关系数的特点实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。2、模拟情况，模拟过程。在edit窗中输入如下程序：4、观看输出的数据序列，输入如下程序，并提交程序。观看

9、样本自相关系数和偏相关系数，输入输入如下程序，并提交程序。可能模型参数，并与实际模型的系数进行对比，即输入如下程序，并提交。模拟情况，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟情况，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟情况，模拟过程。重复步骤3-7即可（但部分程序需要修改，请读者自己完成）。模拟ARMA模型，模拟过程。重复步骤3-7即可(但部分程序需要修改，请读者自己完成).回到graph窗口观看各种序列图形的异同。退出SAS系统,关闭计算机.总程序:data a; a1=0; a2=0; do n=1to 250; a=ran

10、nor(32565); x=a+0.65*a1+0.24*a2; a2=a1; a1=a; output; end; run;proc gplot data=a;symbol i=spline h=1 w=1;plot x*n /haxis=-10 to 260 by 10;run;proc arima data=a; identify var=x nlag=10 minic p=(0:3) q=(0:3) outcov=exp1; estimate q=2 noint;run; proc gplot data=exp1;symbol1 i=needle c=red;plot corr*lag

11、=1;run;proc gplot data=exp1;symbol2 i=needle c=green;plot partcorr*lag=2;run;quit; 实验四 分析化工生产量数据实验目的：进一步熟悉时刻序列建模的差不多步骤，掌握用SACF及SPACF定 模型的阶的方法。实验内容：分析化工生产过程的产量序列。实验要求：掌握ARMA模型建模的差不多步骤，初步掌握数据分析技巧。写出 实验报告。实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。创建名为exp2的SAS数据集，即在窗中输入下列语句：保存此步骤中的程序，供以后分析使用（只需按工具条上的保存按钮然后填写完提问

12、后就能够把这段程序保存下来即可）。绘数据与时刻的关系图，初步识不序列，输入下列程序：提交程序，在graph窗口中观看序列，能够看出此序列是均值平稳序列。识不模型，输入如下程序。提交程序，观看输出结果，发觉二阶样本自相关系数和一阶的样本偏相关系数都在2倍的标准差之外，那么我们首先作为一阶AR模型可能，输入如下程序：提交程序，观看输出结果，发觉残差能通过白噪声检验，但它的二阶的样本偏相关系数比较大，那么我们考虑二阶AR模型。输入如下程序：提交程序，观看输出结果，发觉残差样本自相关系数和样本偏相关系数都 在2倍的标准差之内。且能通过白噪声检验。比较两个模型的AIC和SBC, 发觉第二个模型的AIC和

13、SBC都比第一个的小，故我们选择第二个模型为 我们的结果。记录参数可能值，写出模型方程式。进行预测，输入如下程序：提交程序，观看输出结果。退出SAS系统，关闭计算机。data exp2; infile D:exp1.txt; input x ; n=_n_;proc print;run;proc gplot data=exp2; symbol i=join v=star h=2 ci=green cv=red; plot x*n/vref=30 50 70 cvref=red lvref=2 ;run;proc arima data=exp2;identify var=x nlag=12 mi

14、nic p=(0:3) q=(0:3);estimate plot p=1;forecast lead=2 out=out;run;quit; 实验五 模拟ARIMA模型和季节ARIMA模型实验目的：熟悉各种ARIMA模型的样本自相关系数和偏相关系数的特点， 区不各种ARIMA模型的图形，为理论学习提供直观的印象。实验内容：随机模拟各种ARIMA模型。实验要求：记录各ARIMA模型的样本自相关系数和偏相关系数观看各序列 图形的异同，总结ARIMA模型的样本自相关系数和偏相关系数 的特点实验时刻：2小时。实验软件：SAS系统。实验步骤开机进入SAS系统。2、模拟ARIMA(0,1,1)过程，模拟

15、过程。创建数据集，在edit窗中输入如下程序： 4、观看输出的数据序列，输入如下程序：。 5、提交程序，在Graph窗口中观看图形。6、观看样本自相关系数和偏相关系数，输入输入如下程序： 提交程序，发觉自相关系数成缓慢下降的趋势，讲明要做差分运算，做一阶差分运算，输入如下程序：提交程序，观看样本自相关系数与样本偏相关系数，发觉自相关系数1阶截尾，故推断差分后序列为MA(1)模型。进行模型参数可能，输入如下程序：提交程序，并观看残差图，发觉模型拟合完全。10、写出模型的方程，并与真实模型对比。title;data a; x1=0.9; a1=0; do n=0 to 250; a=rannor(

16、32565); x=x1+a-0.8*a1; x1=x; a1=a; output; end; run;proc gplot data=a; symbol i=join v=dot h=1 ci=green cv=red; plot x*n/vref=-2 1 4 cvref=red lvref=2 haxis=-10 to 260 by 10;run;proc arima data=a; identify var=x nlag=10 minic p=(0:3) q=(0:3) outcov=exp1; run; proc gplot data=exp1;symbol1 i=needle c=

17、red;plot corr*lag=1;run;proc gplot data=exp1;symbol2 i=needle c=green;plot partcorr*lag=2;run;proc arima data=a; identify var=x(1) nlag=24 minic p=(0:3) q=(0:3); /*一阶差分x(1)*/run;estimate q=1 plot noint;run;quit;11、模拟ARIMA(1,1,0)模型，模拟过程。重复步骤 3-10即可（但部分程序需要修改，请读者自己完成）。模拟模型， 模拟模型, 即模型。13、创建数据集，在edit窗中输

18、入如下程序：绘序列图，输入如下程序：提交程序，到graph窗口中观看序列图形。初步识不模型，输入如下程序：提交程序，观看样本自相关系数和样本偏相关系数。做季节差分和一阶差分除掉季节因子和趋势因子，输入如下程序：提交程序，观看样本自相关系数和样本偏相关系数，确定模型阶数。可能模型参数，输入如下程序：提交程序，观看残差的样本自相关系数和样本偏相关系数，看是否通过 了白噪声检验。写出模型方程式，并与真实模型对比。回到graph窗口观看各种序列图形的异同。退出SAS系统,关闭计算机. data c; x1=0.9;x2=0;x3=0;x4=0;x5=0;x6=0;x7=0; x8=0;x9=0;x10

19、=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=0 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=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

20、;a6=a5;a5=a4;a4=a3;a3=a2;a2=a1;a1=a; output; end; run;proc gplot data=c; symbol i=join v=dot h=1 ci=green cv=red; plot x*n/vref=-20 1 10 cvref=red lvref=2 haxis=-10 to 260 by 10;run;proc arima data=c; identify var=x nlag=20 minic p=(0:3) q=(0:3); run; identify var=x(1,12) nlag=36 minic p=(0:3) q=(0:

21、3) outcov=exp1;run;estimate q=(1)(12) method=cls noint;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;quit; 实验六 分析美国国民生产总值的季度数据一、实验目的：进一步学习数据分析技巧，进一步了解ARIMA模型。二、实验内容：47年1季度到96年3季度美国国民生产总值的季度数据。三、实验要求：写出分析报告。四、

22、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp3的SAS数据集，输入如下程序：保存上述程序，供以后分析使用（只需按工具条上的保存按钮，然后填写 完提问后就能够把这段程序保存下来）。绘序列图，输入如下程序：观看图形，发觉图形成指数函数上升形式，故做对数变换，输入如下程序：绘变换后序列图，输入如下程序：提交程序，到graph窗口中观看变换后的序列图，能够看出它成直线上升趋势。对序列做初步识不，输入如下程序：提交程序，观看样本自相关系数，可看出有缓慢下降趋势，结合我们观看的图形，我们明白要对序列做差分运算，作一阶差分，输入如下程序：提交程序，观看样

23、本自相关系数，可看出样本自相关系数5步后是截尾的,那么确定为MA(5)模型,进行参数可能,输入如下程序:提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟合充分。且MA1,3 , MA1,4的T值较小，讲明参数显著为0，除掉这两项重新进行可能，输入如下程序：提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟合充分，且残差标准误与前一可能相差专门小，故以此结果为我们所要的结果，依此结果写出方程式。进行预测，预测美国以后2年的每季国民生产总值。输入如下程序：提交程序，并把预测值记录下来。退出SAS系统，关闭计算机。data exp3; infile C:Documents

24、and SettingsAdministrator桌面exp1.txt; input gnp; date=intnx(qtr,1jan47d,_n_-1); format date yyqc.;run;proc gplot data=exp3; symbol i=join w=2 ci=green; /*w：线的大小 h：点的大小 l：线的类型*/ plot gnp*date=1;run;data lexp; set exp3; lgnp=log(gnp);run;proc gplot data=lexp; symbol2 i=spline c=red w=2; plot lgnp*date=

25、2 ;run;proc arima data=lexp; identify var=lgnp(1) nlag=12 minic p=(0:3) q=(0:3);run;estimate q=5;run;estimate q=(1,2,5) plot;run;forecast lead=6 interval=qtr id=date out=results;run;data results2;set results;gnp=exp(lgnp);l95=exp(l95);u95=exp(u95);forecast=exp(forecast);keep date gnp l95 forecast u9

26、5;run;proc print data=results2;var date gnp l95 forecast u95;where date=1jan96d;run;proc gplot data=results2;plot gnp*date=1 forecast*date=2 l95*date=3 u95*date=3/overlay legend;symbol1 v=dot cv=red i=none h=1 w=1;symbol2 i=join ci=green;symbol3 i=spline ci=black l=2;run;quit; 实验七 分析国际航线月度旅客总数数据实验目的

27、：熟悉运用SAS建立模型的方法，进一步 了解模型的特征。二、实验内容：1949年1月至1960年12月国际航线月度旅客总数数据。三、实验要求：写出分析报告。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp4的SAS数据集，输入如下程序：绘序列图，输入如下程序：提交程序，观看图形，发觉图形有专门强的季节性，且成指数函数上升形式，故做对数变换，输入如下程序：绘变换后序列图，输入如下程序：提交程序，到graph窗口中观看变换后的序列图，能够看出它总的趋势成直线上升，且有专门强的季节性。对序列做初步识不，输入如下程序：提交程序，观看样本自相关系数和

28、偏相关系数，可看出样本自相关系数有缓慢下降趋势，偏相关系数在1步，13步，25步较大，我们作一步一阶差分，输入如下程序：提交程序，观看样本自相关系数和偏相关系数，发觉样本自相关系数在12步，24步，36步特不大,而偏相关系数在12步特不大，那么我们再做12步的一阶差分，输入如下程序：10、提交程序，观看样本自相关系数和偏相关系数，发觉样本自相关系数在1步，12步特不大,而偏相关系数看不出有特不的规律，我们可确定模型的MA因子为。11、进行参数可能，输入如下程序：、提交程序，观看输出结果，可看出模型通过了白噪声检验，讲明模型拟 合充分，故以此结果为我们所要的结果，依此结果写出方程式。13、进行预

29、测，输入如下程序： 14、提交程序，认真观看预测的结果有什么规律，考虑什么缘故有如此的规律？15、变换预测值，以猎取原度量下的预测值，输入如下程序：绘预测和置信限的散点图，输入如下程序：提交程序，观看图形。退出SAS系统，关闭计算机。data exp4; infile C:Documents and SettingsAdministrator桌面exp1.txt; input air; date=intnx(month,1jan49d,_n_-1); format date monyy.;run;proc gplot data=exp4;symbol1 i=join v=dot c=red;p

30、lot air*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;proc arima data=lair; identify var=lair nlag=36;run;identify var=lair(1) nlag=24;run;identify var=lair(1,12) nlag=24 minic;run;estimate q=(1)(12) noconstant method=uls plot;

31、run;forecast lead=3 interval=month id=date out=b;run;proc print data=b;run;data c; set b; air=exp(lair); forecast=exp(forecast+std*std/2); l95=exp(l95); u95=exp(u95);run;proc print data=c;run;symbol1 I=none v=star r=1 c=red;symbol2 I=join v=plus r=1 c=green;symbol3 I=join v=none l=3 r=1 c=blue;proc

32、gplot data=c;where data=1jan59d;plot air*date=1 forecast*date=2 l95*date=3 u95*date=3/ overlay haxis=1jan59d to 1mar61d by month;run; 实验八 干预模型的建模实验目的：掌握干预模型的分析方法，进一步熟悉ARIMA过程的使用方法。二、实验内容：1955年1月至1972年12月洛杉矶月平均臭氧数据。三、实验要求：写出实验报告，掌握干预模型的建模方法。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp5的SAS数据集，

33、输入如下程序：或者输入如下程序： data exp5; input ozone ; date=intnx(month,1jan55d,_n_-1);format date monyy.;month=month(date);year=year(date);x1=year=1960;summer=(5month1965);winter=(year1965)-summer; cards; 只输入 ozone 一栏的数据 ； run;保存上述程序，供以后分析使用（只需按工具条上的保存按钮，然后填写 完提问后就能够把这段程序保存下来）。4、绘序列图，输入如下程序：5、提交程序，观看图形，发觉图形有专门强

34、的季节性和缓慢下降的趋势。初步识不模型，输入如下程序：提交程序，观看样本自相关系数和偏相关系数，可看出样本自相关系数在1步，12步，24步，36步都较大，且具有周期性，偏相关系数在1步最大，我们作季节差分，输入如下程序：提交程序，观看样本自相关系数和偏相关系数，发觉样本自相关系数在 1步，12步较大,而偏相关系数在1步，12步，24步都较大，且呈现拖尾 现象，我们可确定模型的MA因子为。9、进行参数可能，输入如下程序：10、提交程序，观看输出结果，可看出模型不是专门洁净，且不能通过白噪声检验。我们能够做残差序列图，观看残差的特性，输入如下程序：11、进行预测，输入如下程序：12、提交程序，观看

35、图形，可看出前面一段时期的残差比后面的要大。 13、我们考察修建高速公路后，是否对臭氧有显著性阻碍，输入如下程序：14、提交程序，观看输出结果，发觉模型的标准差,AIC,SBC都变小了专门多， 且x1的阻碍显著。考虑什么缘故要对x1进行季节差分？15、我们再来考察汽车装上尾气过滤器，是否对臭氧有显著性阻碍，输入如 下程序：提交程序，观看输出结果，发觉模型的标准差,AIC,SBC都变小了，且 模型差不多上通过了白噪声检验，同时x1,summer的阻碍显著,而winter 的阻碍不显著。考虑什么缘故不对summer和 winter进行差分？ 17、进行预测值，输入如下程序： 注：如此的预测是x1,

36、summer,winter已知的预测。 18、提交程序，观看预测值。 19、退出SAS系统，关闭计算机。data exp5; infile C:Documents and SettingsAdministrator桌面exp1.txt; input n ozone x1 summer winter; date=intnx(month,1jan55d,_n_-1); format date monyy.;run;proc gplot data=exp5; symbol1 i=join v=dot c=red; plot ozone*date=1;run;proc arima data=exp5;

37、identify var=ozone nlag=36;run;identify var=ozone(12) nlag=36 minic;run;estimate q=(1)(6)(12) method=cls plot;run;forecast lead=36 interval=month id=date out=b noprint;run;proc gplot data=b;symbol I=spline v=dot c=red;plot residual*date;run;proc arima data=exp5; identify var=ozone(12) crosscorr=(x1(

38、12) noprint; estimate q=(1)(6)(12) input=(x1) noconstant method=ml itprint plot;run;proc arima data=exp5; identify var=ozone(12) crosscorr=(x1(12) summer winter) noprint; estimate q=(1)(12) input=(x1 summer winter) noconstant method=ml itprint plot;run;forecast lead=12 id=date interval=month out=c;r

39、un;proc gplot data=c;symbol I=spline v=dot c=red;plot residual*date;run;quit; 实验九 传递函数模型的建模实验目的：熟悉传递函数模型的建模方法。二、实验内容：煤气炉数据。三、实验要求：写出实验报告，总结传递函数模型的建模的一般步骤。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp6的SAS数据集，输入如下程序：保存上述程序，供以后分析使用（只需按工具条上的保存按钮，然后填写 完提问后就能够把这段程序保存下来）。4、绘序列图，输入如下程序：5、提交程序，认真观看两序列

40、图形，看两者有何联系。6、先观看和的相关情况，看是否要做差分，输入如下程序：7、提交程序，观看的自相关和互相关系数，发觉都专门快的衰减，表明不 要做差分运算。识不输入序列, 输入如下程序：提交程序，观看的自相关和偏相关系数，能够看到偏相关系数是3步 截尾的。10、对拟合AR(3)模型，看是否充分，输入如下程序：11、提交程序，观看输出结果，可看到模型通过了白噪声检验，讲明拟合效果不错，把拟合的方程式写出来。12、观看预白噪声化后的两序列的互相关系数，输入如下程序： 13、提交程序，观看样本自相关系数和偏相关系数和互相关系数，我们能够 初步识不传递函数模型为（2，2，3）(考虑:什么缘故?)，即

41、： 14、进行参数可能，并查看残差的相关情况，输入如下程序：15、提交程序，观看输出结果，能够看到残差的偏相关系数是2步 截尾的。那么模型可识不为： 16、进行参数可能，输入如下程序：提交程序，观看输出结果，可看到专门小，且模型通过了白噪声检验，那么我们除掉这一项,再进行可能，输入如下程序：提交程序，观看输出结果，可看到模型通过了白噪声检验提交程序，讲明模型拟合充分,请写出方程式。19、进行预测，输入如下程序： 20、提交程序，观看预测结果。 21、退出SAS系统，关闭计算机。data exp6; infile C:Documents and SettingsAdministrator桌面ex

42、p1.txt; input x y; t=_n_;run;proc gplot data=exp6; symbol1 i=spline c=red; symbol2 i=spline c=green; plot x*t=1 y*t=2;run;proc arima data=exp6;identify var=y crosscorr=(x) nlag=12;run;proc arima data=exp6; identify var=x nlag=12 minic;run;estimate p=(1,2,3,5) noint plot ;run;estimate input=(3$(1,2)/

43、(1,2)x) plot;run;estimate p=2 input=(3$(1,2)/(1,2)x) plot;run;estimate p=2 input=(3$(1,2)/(1)x) plot;run;forecast lead=6 ;run; 实验十 回归与时序相结合的建模实验目的：熟悉回归与时序相结合的建模方法。二、实验内容：芝加哥某食品公司大众食品周销售数据。三、实验要求：写出实验报告，总结回归与时序相结合的建模的一般步骤。四、实验时刻：2小时。五、实验软件：SAS系统。六、实验步骤1、开机进入SAS系统。2、建立名为exp7的SAS数据集，输入如下程序：保存上述程序，供以后分析

44、使用（只需按工具条上的保存按钮，然后填写 完提问后就能够把这段程序保存下来）。4、首先只分析销售额的数据，不加回归项。绘序列图，输入如下程序：5、提交程序，认真观看序列图形。6、初步识不模型，输入如下程序：7、提交程序，观看y1的相关系数，发觉偏相关系数是4阶截尾的，那么 我们初步识不为AR(4)模型，进行参数可能，并观看残差相关系数。输入 如下程序：8、提交程序，观看输出结果，可看到模型拟合得依旧比较好。然后能够实验其他一些模型，最后依照AIC和BIC准则，我们最后选定 模型为: 下面我们开始加入回归项，首先我们绘四个序列的图形。输入如下程序：提交程序，观看这四个序列有什么特点。绘y1对y2

45、、y3、y4的散点图，输入如下程序：提交程序，观看他们的相关性，可看出y1和y2负相关，y1和y3正相 关，而y1和y4仿佛不相关。做纯回归分析，输入如下程序：提交程序，观看输出结果，可看到y4的系数接近于零，我们除掉这一 项再做回归，并观看残差的相关系数，输入如下程序：提交程序，观看输出结果，可看到残差不是白躁声。我们把残差用ARMA 模型拟合，输入如下程序：提交程序，观看输出结果,能够看出模型拟合比较充分，且y4、MA1,1、 AR1,3的系数接近如零，除掉这几项，再观看，输入如下程序：提交程序，观看输出结果,能够看出模型拟合比较充分，且残差的标准 误和前一模型没有多大变化，且AIC 和B

46、IC也比前一模型小，故我们 就选择这一模型，把这一结果记录下来。下面我们来看看残差对预测值,y2,y3的关系图。输入如下程序：提交程序，观看图形，可看出残差对y2,y3还不是十分充分，我们加入y2,y3的滞后一阶，看结果有什么变化，输入如下程序： 23、提交程序，观看输出结果，并与原来结果比较，看是否有进步。(没进步) 24、进行预测，输入如下程序： 25、提交程序，观看预测结果。 26、退出SAS系统，关闭计算机。data exp7; infile C:Documents and SettingsAdministrator桌面exp1.txt; input y1 y2 y3 y4 ; dat

47、e=intnx(week,14sep91d,_n_-1); format date date9.;run;proc gplot data=exp7; symbol1 i=spline c=red; plot y1*date=1;run;proc arima data=exp7; identify var=y1 nlag=15 minic;run;estimate p=(1,4) plot;run;proc gplot data=exp7;symbol3 i=spline c=green;plot y1*date=3 y2*date=3 y3*date=3 y4*date=3;run;proc

48、plot data=exp7; plot y1*y2=* y1*y3=* y1*y4=*;run;proc arima data=exp7; identify var=y1 crosscorr=(y2 y3 y4) noprint; estimate input=(y2 y3 y4);run;identify var=y1 crosscorr=(y2 y3);estimate input=(y2 y3) plot;run;identify var=y1 crosscorr=(y2 y3 y4);estimate p=4 q=3 input=(y2 y3 y4) plot;run;identif

49、y var=y1 crosscorr=(y2 y3);estimate p=(1,2,4) q=(2,3) input=(y2 y3 ) plot;run;forecast lead=6 id=date interval=week out=a noprint;run;data b;merge exp7 a;run;proc gplot data=b;plot residual*forecast=* residual*y2=* residual*y3=*;symbol c=red;run;data d; set exp7; y21=lag(y2); y31=lag(y3); run;proc a

50、rima data=d; identify var=y1 crosscorr=(y2 y21 y3 y31) noprint; estimate p=(1,2,4) q=(2,3) input=(y2 y21 y3 y31) plot;run;forecast lead=6 ;run;附数据： 太阳黑子年度数据(1742-1957) 美国国民收入数据(1947第一季度到1996第三季度) （顺序是横向排列） 化工生产过程的产量数据（顺序是横向排列） 国际航线月度旅客数据 (1949.01-1960.12,单位:千人)（顺序是横向排列） 112 118 132 129 121 135 148 1

51、48 136 119 104 118115 126 141 135 125 149 170 170 158 133 114 140145 150 178 163 172 178 199 199 184 162 146 166171 180 193 181 183 218 230 242 209 191 172 194196 196 236 235 229 243 264 272 237 211 180 201204 188 235 227 234 264 302 293 259 229 203 229242 233 267 269 270 315 364 347 312 274 237 278

52、284 277 317 313 318 374 413 405 355 306 271 306315 301 356 348 355 422 465 467 404 347 305 336340 318 362 348 363 435 491 505 404 359 310 337360 342 406 396 420 472 548 559 463 407 362 405417 391 419 461 472 535 622 606 508 461 390 432 洛杉矶臭氧每小时读数的月平均值数据 (1955-1972) 1 2.7 0.0 0.0 0.0 2 2.0 0.0 0.0 0.

53、0 3 3.6 0.0 0.0 0.0 4 5.0 0.0 0.0 0.0 5 6.5 0.0 0.0 0.0 6 6.1 0.0 0.0 0.0 7 5.9 0.0 0.0 0.0 8 5.0 0.0 0.0 0.0 9 6.4 0.0 0.0 0.0 10 7.4 0.0 0.0 0.0 11 8.2 0.0 0.0 0.0 12 3.9 0.0 0.0 0.0 13 4.1 0.0 0.0 0.0 14 4.5 0.0 0.0 0.0 15 5.5 0.0 0.0 0.0 16 3.8 0.0 0.0 0.0 17 4.8 0.0 0.0 0.0 18 5.6 0.0 0.0 0.0 1

54、9 6.3 0.0 0.0 0.0 20 5.9 0.0 0.0 0.0 21 8.7 0.0 0.0 0.0 22 5.3 0.0 0.0 0.0 23 5.7 0.0 0.0 0.0 24 5.7 0.0 0.0 0.0 25 3.0 0.0 0.0 0.0 26 3.4 0.0 0.0 0.0 27 4.9 0.0 0.0 0.0 28 4.5 0.0 0.0 0.0 29 4.0 0.0 0.0 0.0 30 5.7 0.0 0.0 0.0 31 6.3 0.0 0.0 0.0 32 7.1 0.0 0.0 0.0 33 8.0 0.0 0.0 0.0 34 5.2 0.0 0.0 0

55、.0 35 5.0 0.0 0.0 0.0 36 4.7 0.0 0.0 0.0 37 3.7 0.0 0.0 0.0 38 3.1 0.0 0.0 0.0 39 2.5 0.0 0.0 0.0 40 4.0 0.0 0.0 0.0 41 4.1 0.0 0.0 0.0 42 4.6 0.0 0.0 0.0 43 4.4 0.0 0.0 0.0 44 4.2 0.0 0.0 0.0 45 5.1 0.0 0.0 0.0 46 4.6 0.0 0.0 0.0 47 4.4 0.0 0.0 0.0 48 4.0 0.0 0.0 0.0 49 2.9 0.0 0.0 0.0 50 2.4 0.0 0

56、.0 0.0 51 4.7 0.0 0.0 0.0 52 5.1 0.0 0.0 0.0 53 4.0 0.0 0.0 0.0 54 7.5 0.0 0.0 0.0 55 7.7 0.0 0.0 0.0 56 6.3 0.0 0.0 0.0 57 5.3 0.0 0.0 0.0 58 5.7 0.0 0.0 0.0 59 4.8 0.0 0.0 0.0 60 2.7 0.0 0.0 0.0 61 1.7 1.0 0.0 0.0 62 2.0 1.0 0.0 0.0 63 3.4 1.0 0.0 0.0 64 4.0 1.0 0.0 0.0 65 4.3 1.0 0.0 0.0 66 5.0 1

57、.0 0.0 0.0 67 5.5 1.0 0.0 0.0 68 5.0 1.0 0.0 0.0 69 5.4 1.0 0.0 0.0 70 3.8 1.0 0.0 0.0 71 2.4 1.0 0.0 0.0 72 2.0 1.0 0.0 0.0 73 2.2 1.0 0.0 0.0 74 2.5 1.0 0.0 0.0 75 2.6 1.0 0.0 0.0 76 3.3 1.0 0.0 0.0 77 2.9 1.0 0.0 0.0 78 4.3 1.0 0.0 0.0 79 4.2 1.0 0.0 0.0 80 4.2 1.0 0.0 0.0 81 3.9 1.0 0.0 0.0 82 3

58、.9 1.0 0.0 0.0 83 2.5 1.0 0.0 0.0 84 2.2 1.0 0.0 0.0 85 2.4 1.0 0.0 0.0 86 1.9 1.0 0.0 0.0 87 2.1 1.0 0.0 0.0 88 4.5 1.0 0.0 0.0 89 3.3 1.0 0.0 0.0 90 3.4 1.0 0.0 0.0 91 4.1 1.0 0.0 0.0 92 5.7 1.0 0.0 0.0 93 4.8 1.0 0.0 0.0 94 5.0 1.0 0.0 0.0 95 2.8 1.0 0.0 0.0 96 2.9 1.0 0.0 0.0 97 1.7 1.0 0.0 0.0

59、98 3.2 1.0 0.0 0.0 99 2.7 1.0 0.0 0.0 100 3.0 1.0 0.0 0.0 101 3.4 1.0 0.0 0.0 102 3.8 1.0 0.0 0.0 103 5.0 1.0 0.0 0.0 104 4.8 1.0 0.0 0.0 105 4.9 1.0 0.0 0.0 106 3.5 1.0 0.0 0.0 107 2.5 1.0 0.0 0.0 108 2.4 1.0 0.0 0.0 109 1.6 1.0 0.0 0.0 110 2.3 1.0 0.0 0.0 111 2.5 1.0 0.0 0.0 112 3.1 1.0 0.0 0.0 11

60、3 3.5 1.0 0.0 0.0 114 4.5 1.0 0.0 0.0 115 5.7 1.0 0.0 0.0 116 5.0 1.0 0.0 0.0 117 4.6 1.0 0.0 0.0 118 4.8 1.0 0.0 0.0 119 2.1 1.0 0.0 0.0 120 1.4 1.0 0.0 0.0 121 2.1 1.0 0.0 0.0 122 2.9 1.0 0.0 0.0 123 2.7 1.0 0.0 0.0 124 4.2 1.0 0.0 0.0 125 3.9 1.0 0.0 0.0 126 4.1 1.0 0.0 0.0 127 4.6 1.0 0.0 0.0 12