上机试验-时间序列分析.doc
文本预览下载声明
上机实验
第一章 数据导入
一、SAS操作界面
1．菜单栏、工具栏、窗口
2．菜单栏
3．工具栏
4．窗口
编辑窗口、运行窗口、结果输出窗口
二、创建时间序列数据集
1．数据集及创建
SAS数据集使用二级命名制——数据库名.数据集名。数据库包括存放临时数据集的数据库——WORK数据库，当退出本次操作后，该库中的数据集将被全部删除；永久数据库——SASUSER数据库，数据集不会因为退出SAS系统而丢失。
创建WORK数据集
Data exa1;
Input time monyy7. price;
Format time monyy5. ;
Cards;
Jan2005 12.3
Feb2005 12.5
Mar2005 13.1
Apr2005 12.7
May2005 13.2
;
Run;
创建SASUSER数据集
Data sasuser.exam1;
查看数据集
（1）使用PRINT查看数据集
Proc print data=exa1;
Run;
（2）使用资源管理器查看数据集
2．直接导入外部数据文件转换成SAS数据集
三、时间序列数据集的处理
1．间隔函数的使用intnx
Data exam2;
Input price;
Time=intnx(‘month’,’01jan2005’d, _n_-1);
Format time monyy5.;
Cards;
12.3
12.5
12.7
13.1
13.4
12.3
;
Run;
Proc print data=exam2;
Run;
2.序列变换
Data exam3;
Input price;
Lprice=log(price);
Time=intnx(‘month’,’01jan2005’d, _n_-1);
Format time monyy5.;
Cards;
12.3
12.5
12.7
13.1
13.4
12.3
;
Run;
Proc print data=exam3;
Run;
3.调用子集
Data exam4;
Set exam3;
Keep time lprice;
Where time=’01mar2005’d;（或者Where ‘01may2005’=time=’01mar2005’d;）
Proc print data=exam4;
Run;
4.缺失值插值
Data exam5;
Input price;
Time=intnx (‘month’,’01jan2005’d,_n_-1);
Format time monyy5.;
Cards;
12.3
12.5
12.7
13.1
.
13.4
12.3
;
Run;
Proc print data=exam5;
Run;
生成具有缺失值的exam5，用插补充生成exam6
Proc expand data=exam5 out=exam6;
Id time;
Proc print data=exam5;
Proc print data=exam6;
Run;
第二章 时间序列的预处理
一、绘制时序图
data exam2_1;
input price1 price2;
time=intnx(month,01jul2004d,_n_-1);
format time monyy5.;
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
;
run;
proc print data=exam2_1;
run;
proc gplot data=exam2_1;
plot price1*time=1 price2*time=2/overlay;
symbol1 c=red i=spline v=star;
symbol2 c=green i=join v=cycle;
run;
二、平稳性检验和纯随机性检验
利用IDENTIFY语句进行检验。IDENTIFY语句执行后，SAS会输出五个方面的信息：
（1）分析变量的描述性统计；
（2）样本自相关图；
（3）样本偏自相关图；
（4）样本逆自相关图；
（5）纯随机检验结果。
Data example2_2;
Input freq@@;
Year=intnx(‘year’,’01jan1970’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 24
显示全部