经济时间序列分_各种模型分析报告.docx
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《经济时间序列分_各种模型分析报告.docx》由会员分享,可在线阅读,更多相关《经济时间序列分_各种模型分析报告.docx(49页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
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系统。2、 创建名为exp1的SAS数据集,即在窗中输入下列语句:3、 保存此步骤中的程序,供以后分析使用(只需按工具条上的保存按钮然后填写完提问后就可以把这段程序保存下来即可)。4、 绘数据与时间的关系图,初步识别序列,输入下列程序:ods html;ods listing close;5、 run;提交程序,在graph窗口中观察序
3、列,可以看出此序列是均值平稳序列。6、 识别模型,输入如下程序。7、 提交程序,观察输出结果。初步识别序列为AR(2)模型。8、 估计和诊断。输入如下程序:9、 提交程序,观察输出结果。假设通过了白噪声检验,且模型合理,则进行预测。10、 进行预测,输入如下程序:11、 提交程序,观察输出结果。12、 退出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 clos
4、e;proc gplot data=exp1 ; symbol 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=ou
5、t;run;选取拟合模型的规则:1.模型显著有效(残差检验为白噪声)2.模型参数尽可能少3.结合自相关图和偏自相关图以与minic条件(BIC信息量最小原则),选取显著有效的参数实验二模拟AR模型一、 实验目的:熟悉各种AR模型的样本自相关系数和偏相关系数的特点,为理论学习提供直观的印象。二、 实验容:随机模拟各种AR模型。三、 实验要求:记录各AR模型的样本自相关系数和偏相关系数,观察各种序列图形,总结AR模型的样本自相关系数和偏相关系数的特点四、 实验时间:2小时。五、 实验软件:SAS系统。六、 实验步骤1、开机进入SAS系统。2、 模拟实根情况,模拟过程。3、 在edit窗中输入如下程
6、序:4、观察输出的数据,输入如下程序,并提交程序。观察样本自相关系数和偏相关系数,输入输入如下程序,并提交程序。作为作业把样本自相关系数和偏相关系数记录下来。5、 估计模型参数,并与实际模型的系数进行对比,即输入如下程序,并提交。6、 模拟虚根情况,模拟过程。重复步骤3-7即可(但部分程序需要修改,请读者自己完成)。7、 模拟AR(3)模型,模拟过程。重复步骤3-7即可(但部分程序需要修改,请读者自己完成).10、回到graph窗口观察各种序列图形的异同11、退出SAS系统,关闭计算机.总程序:title;data a; x1=0.5; x2=0.5; do i=-50 to 250; a=r
7、annor(32565); x=a-0.6*x1+0.4*x2; x2=x1; x1=x; output; end; 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
8、 i=needle width=6;plot corr*lag;run;proc gplot data=exp1;symbol i=needle width=6;plot partcorr*lag;run;实验三 模拟MA模型和ARMA模型一、 实验目的:熟悉各种MA模型和ARMA模型的样本自相关系数和偏相关系数的特点,为理论学习提供直观的印象。二、 实验容:随机模拟各种MA模型和ARMA模型。三、 实验要求:记录各MA模型和ARMA模型的样本自相关系数和偏相关系数,观察各序列的异同,总结MA模型和ARMA模型的样本自相关系 数和偏相关系数的特点四、 实验时间:2小时。五、 实验软件:SAS系
9、统。六、 实验步骤1、 开机进入SAS系统。2、模拟情况,模拟过程。3 在edit窗中输入如下程序:4、观察输出的数据序列,输入如下程序,并提交程序。5、 观察样本自相关系数和偏相关系数,输入输入如下程序,并提交程序。6、 估计模型参数,并与实际模型的系数进行对比,即输入如下程序,并提交。7、 模拟情况,模拟过程。重复步骤3-7即可(但部分程序需要修改,请读者自己完成)。8、 模拟情况,模拟过程。重复步骤3-7即可(但部分程序需要修改,请读者自己完成)。9、 模拟情况,模拟过程。重复步骤3-7即可(但部分程序需要修改,请读者自己完成)。10、 模拟ARMA模型,模拟过程。重复步骤3-7即可(但
10、部分程序需要修改,请读者自己完成).11、 回到graph窗口观察各种序列图形的异同。12、 退出SAS系统,关闭计算机.总程序:data a; a1=0; a2=0; do n=1to 250; a=rannor(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
11、) q=(0:3) outcov=exp1; estimate q=2 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;实验四 分析化工生产量数据一、 实验目的:进一步熟悉时间序列建模的基本步骤,掌握用SACF与SPACF定模型的阶的方法。二、 实验容:分析化工生产过程的产量序列。三、 实验要求:掌握ARMA模型建模的基本步骤,初步掌握数据分析技
12、巧。写出实验报告。四、 实验时间:2小时。五、 实验软件:SAS系统。六、 实验步骤1、 开机进入SAS系统。2、 创建名为exp2的SAS数据集,即在窗中输入下列语句:3、 保存此步骤中的程序,供以后分析使用(只需按工具条上的保存按钮然后填写完提问后就可以把这段程序保存下来即可)。4、 绘数据与时间的关系图,初步识别序列,输入下列程序:5、 提交程序,在graph窗口中观察序列,可以看出此序列是均值平稳序列。6、 识别模型,输入如下程序。7、 提交程序,观察输出结果,发现二阶样本自相关系数和一阶的样本偏相关系数都在2倍的标准差之外,那么我们首先作为一阶AR模型估计,输入如下程序:8、 提交程
13、序,观察输出结果,发现残差能通过白噪声检验,但它的二阶的样本偏相关系数比较大,那么我们考虑二阶AR模型。输入如下程序:9、 提交程序,观察输出结果,发现残差样本自相关系数和样本偏相关系数都 在2倍的标准差之。且能通过白噪声检验。比较两个模型的AIC和SBC, 发现第二个模型的AIC和SBC都比第一个的小,故我们选择第二个模型为 我们的结果。10、 记录参数估计值,写出模型方程式。11、 进行预测,输入如下程序:12、 提交程序,观察输出结果。13、 退出SAS系统,关闭计算机。data exp2; infile D:exp1.txt; input x ; n=_n_;proc print;ru
14、n;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 minic p=(0:3) q=(0:3);estimate plot p=1;forecast lead=2 out=out;run;quit;实验五 模拟ARIMA模型和季节ARIMA模型一、 实验目的:熟悉各种ARIMA模型的样本自相关系数和偏相关系数的特点,区别各种ARIMA模型
15、的图形,为理论学习提供直观的印象。一、 实验容:随机模拟各种ARIMA模型。二、 实验要求:记录各ARIMA模型的样本自相关系数和偏相关系数观察各序列图形的异同,总结ARIMA模型的样本自相关系数和偏相关系数的特点三、 实验时间:2小时。四、 实验软件:SAS系统。五、 实验步骤1. 开机进入SAS系统。2. 2、模拟ARIMA(0,1,1)过程,模拟过程。3. 创建数据集,在edit窗中输入如下程序:4、观察输出的数据序列,输入如下程序:。 5、提交程序,在Graph窗口中观察图形。6、观察样本自相关系数和偏相关系数,输入输入如下程序: 提交程序,发现自相关系数成缓慢下降的趋势,说明要做差分
16、运算,做一阶差分运算,输入如下程序:7、 提交程序,观察样本自相关系数与样本偏相关系数,发现自相关系数1阶截尾,故判断差分后序列为MA(1)模型。进行模型参数估计,输入如下程序:8、 提交程序,并观察残差图,发现模型拟合完全。10、写出模型的方程,并与真实模型对比。title;data a; x1=0.9; a1=0; do n=0 to 250; a=rannor(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; plo
17、t 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=red;plot corr*lag=1;run;proc gplot data=exp1;symbol2 i=needle c=green;plot partcorr*lag=2;run;proc arima data=a; id
18、entify 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即可(但部分程序需要修改,请读者自己完成)。12 模拟模型, 模拟模型, 即模型。13、创建数据集,在edit窗中输入如下程序:14、 绘序列图,输入如下程序:15、 提交程序,到graph窗口中观察序列图形。16、 初步识别模型,输入如下程序:17、 提交程序,观察样本自相关系数和样本偏相关系数。18、 做季节差分和一阶差分除掉季节因子
19、和趋势因子,输入如下程序:19、 提交程序,观察样本自相关系数和样本偏相关系数,确定模型阶数。20、 估计模型参数,输入如下程序:21、 提交程序,观察残差的样本自相关系数和样本偏相关系数,看是否通过 了白噪声检验。写出模型方程式,并与真实模型对比。22、 回到graph窗口观察各种序列图形的异同。23、 退出SAS系统,关闭计算机. 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=
20、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;a6=a5;a5=a4;a4=a3;a3=a2;a2=a1;a1=a; output; end; run;proc gplot data=c;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 经济 时间 序列 各种 模型 分析 报告
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内