第二章：上机指导

一、课本上的上机指导示例程序：

1. 识别序列的平稳性和白噪声属性

data example2_1;
input  freq@@;
year=intnx('year','1jan1970'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 245 357 301 389
;
proc gplot;
plot freq*year;
symbol v=square c=red i=join;
proc arima data= example2_1;
identify var=freq nlag=22;
run;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

二、课本中的例题：复现程序及数据

1. 课本43页，例2.7

goptions vsize=6.8cm hsize=10cm;
data a;
do time=-50 to 1000 by 1;
noise=rannor(12345);
if time>0 then output;
end;
proc gplot;
plot noise*time;
symbol v=none i=join c=red;
proc arima data=a;
identify var=noise;
run;

 2. 课本29页，例题2_3:

goptions vsize=6.8cm hsize=10cm;
data a;
input percent;
year=intnx('year','01jan1978'd,_n_-1);
format year year4.;
cards;
23.9
21.6
21.6
22.0
21.8
22.4
24.8
28.7
29.1
29.6
30.5
32.1
31.5
33.7
34.8
33.7
33.6
32.9
32.8
34.2
36.2
37.8
39.0
40.5
41.5
41.2
40.4
40.5
40.9
41.9
41.8
43.4
43.2
43.4
44.6
;
proc gplot;
plot percent*year=1;
symbol1 c=red v=star i=join;
run;