时间序列分析实验报告1.docx

1、时间序列分析实验报告1时间序列分析课 程 实 验 报 告 项目名称：数据集建立&时间序列的预处理组员姓名：李菲指导教师：牛宪华完成日期：201 年 3 月 16 日 练习题：1. Input语句数据输入格式有：列表方式或自由格式、列方式、格式化方式、命名方式，分别给出程序例子和运行结果，并分析运行结果。（1） 自由格式data lianxi1_1;input name$ sex$ age;cards;naw 22 ejniw 23 dhf 23 husi 21 huh 24;proc print data=lianxi1_1;run;（2） 列方式data lianxi1_2;input na

2、me$2-5 sex$6-8 age 9-12;cards;naw 22 ejniw 23 dhf 23 husi 21 huh 24;proc print data=lianxi1_2;run;（3） 格式化方式data lianxi1_3;input name$5. sex$ age 2.;cards;naw 22 ejniw 23 dhf 23 husi 21 huh 24;proc print data=lianxi1_3;run;（4） 命名方式data lianxi1_4;Input id name=$ age=;Cards;200012 name=marry age=152000

3、15 age=16 name=join200011 age=16 name=smith;proc print data=lianxi1_4;run;分析：由以上结果及程序可知，自由方式的input后必须加“”表示连续输入，而下面的数据可以随意排列；列方式中每个列名后必须加这个列是下面数据的第几个字节到第几个字节，并且输入的数据要对应整齐；格式化方式中在列名后加“数字.”表示这个列取到第几位（数字表示位数），输入的数据也要一行一行对齐不能多加几列上去；而命名方式在列名后要加“=”，并且在输入数据时要把对应的数据都写成“列名=数据”的形式，输入时可以与自由方式一样进行输入，不一定要整齐。2. 下面

4、已给出相应的数据步，请分别对数据集sample进行如下操作：(1) 使用sort过程，按照x变量进行降序排列，并输出查看结果；(2) 使用print过程中的by、pageby、sumby语句，查看输出结果，并分析by、pageby、sumby语句的含义。data sample2;input x y z ;cards;1 1 1 1 2 3 1 7 8 9 5 4 3 4 2 45 5 6 6 7 1 4 3 5 6 2;（1） 排序data sample2;input x y z ;cards;1 1 1 1 2 3 1 7 8 9 5 4 3 4 2 45 5 6 6 7 1 4 3 5 6

5、 2;proc sort data=sample2;by descending x;proc print data=sample2;run;（2）by语句data sample2;input x y z ;cards;1 1 1 1 2 3 1 7 8 9 5 4 3 4 2 45 5 6 6 7 1 4 3 5 6 2;proc sort data=sample2;by x;proc print data=sample2;by x;run;分析：by语句就是针对不同的x对y、z进行分组。例如x=1时，y、z有4组。pageby语句data sample2;input x y z ;cards

6、;1 1 1 1 2 3 1 7 8 9 5 4 3 4 2 45 5 6 6 7 1 4 3 5 6 2;proc sort data=sample2;by x;proc print data=sample2;by x;pageby x;run;分析：pageby语句需要和by语句一起使用，他就是将by进行的分组再分页表示出来。sumby语句data sample2;input x y z ;cards;1 1 1 1 2 3 1 7 8 9 5 4 3 4 2 45 5 6 6 7 1 4 3 5 6 2;proc sort data=sample2;by x;proc print dat

7、a=sample2;by x;sumby x;run;分析：sumby语句和by语句一起使用就是在by语句的基础上把每组的y、z的和求出来，如x=1时的结果，而下面的因为每组只有一组数据所以就省略的求和的结果，在最下名时也求出了所有y、z的和。实训题：一、 在D盘根目录下使用自己学号后三位定义一个SAS库标记，并在此库下用DATA步把下列数据创建成一个student数据集，再按身高降序输出。结果要求：1.完整程序；2.给出新定义的SAS库的截图；3.给出排序后student数据集的输出结果。NameAgeSexHeight(cm)Weight(kg)Mike23Male17062Alice21

8、Female16148James21Male18080Katie22Female15750Susan20Female16855Tim22Male17575二、已知某公司2003-2005年computer和MP3销量如下表所示，对数据做如下操作：1. 为上述数据建立一个临时SAS数据集company；2. 分别创建产品com和P3的SAS数据集，计算3年间computer和MP3各自的总销量，并输出该结果；结果要求：1.完整程序；2.给出第二问的输出结果。yearproductsale2003computer20002003mp314002004computer14002004mp313002

9、005computer16002005mp322001、程序：libname _102 D:;data _102.student;input name$ age sex$ height weight;cards;Mike 23 male 170 62Alice 21 female 161 48James 21 male 180 80Katie 22 female 157 50Susan 20 female 168 55Tim 22 male 175 75;proc sort data=_102.student;by descending height;proc print data=_102.

10、student;run;输出结果： 2、程序：data company;input year pruduct$ sale;cards;2003 computer 20002003 mp3 14002004 computer 14002004 mp3 13002005 computer 16002005 mp3 2200;data com;set company;where pruduct=computer;data p3;set company;where pruduct=mp3;proc print data=com;sum sale;proc print data=p3;sum sale;

11、run;输出结果：3、下图分别为间隔函数intnx(month,01jan2005d,_n_k)中k取-2、-1、0、1、2时的结果图：分析：由上图可知k的取值代表开始的观测指针，-1代表当前参展时间开始，-2时是当前时间向前，0、1、2表示从-1的时间依次向后一个时间作为开始值。P34页：习题3、程序：data xiti3;input rain;month=intnx(month,1jan1945d,_n_-1);format month date.;cards;69.3 80.0 40.9 74.9 84.6 101.1 225.0 95.3 100.6 48.3 144.5 128.33

12、8.4 52.3 68.6 37.1 148.6 218.7 131.6 112.8 81.8 31.0 47.5 70.196.8 61.5 55.6 171.7 220.5 119.4 63.2 181.6 73.9 64.8 166.9 48.0137.7 80.5 105.2 89.9 174.8 124.0 86.4 136.9 31.5 35.3 112.3 143.0160.8 97.0 80.5 62.5 158.2 7.6 165.9 106.7 92.2 63.2 26.2 77.052.3 105.4 144.3 49.5 116.1 54.1 148.6 159.3 8

13、5.3 67.3 112.8 59.4;proc arima data=xiti3;identify var=rain;run;运行结果：分析：由上图分析后可知该序列的自相关系数一直都比较小，始终控制在2倍的标准差范围内，基本可以认为该序列自始至终都在零轴附近波动，故该序列为平稳时间序列。分析：由上图可知当延迟阶数为18时p值为0.1352，大于显著性水平0.05，接受原假设，即该序列为纯随机性序列。习题5、程序：data xiti5;input sale;month=intnx(month,1jan2000d,_n_-1);format month date.;cards;153 187 2

14、34 212 300 221 201 175 123 104 85 78134 175 243 227 298 256 237 165 124 106 87 74145 203 189 214 295 220 231 174 119 85 67 75117 178 149 178 248 202 162 135 120 96 90 63;proc gplot data=xiti5;plot sale*month=1;symbol c=red v=star i=spline;proc arima data=xiti5;identify var=sale;run;输出结果：分析：由该序列的时序图可知该序列具有明显的周期性，其一年为周期，且每年基本无递增或递减趋势，故该序列不是平稳序列。分析：由上图分析可知，该序列的自相关系数开始从1减到0再到负轴直至接近于-1再慢慢向正轴的1靠近，具有明显的周期性，故该序列为非平稳序列。和时序图的变化基本是吻合的。分析：由上图可知p值均小于0.0001，故也小于显著性水平0.05，即拒绝原假设，故该序列为非白噪声序列，即为非纯随机性序列。