data newlib.a11;
input a1 @@;
day=intnx( 'day','01mar2007'd,_n_-1);
format day yymmdd10.;
cards;
330.45 330.97 331.64 332.87 333.61 333.55
331.90 330.05 328.58 328.31 329.41 330.63
331.63 332.46 333.36 334.45 334.82 334.32
333.05 330.87 329.24 328.87 330.18 331.50
332.81 333.23 334.55 335.82 336.44 335.99
334.65 332.41 331.32 330.73 332.05 333.53
334.66 335.07 336.33 337.39 337.65 337.57
336.25 334.39 332.44 332.25 333.59 334.76
335.89 336.44 337.63 338.54 339.06 338.95
337.41 335.71 333.68 333.69 335.05 336.53
337.81 338.16 339.88 340.57 341.19 340.87
;
run;
proc gplot data=newlib.a11;
symbol i=spline v=star h=2 c=red;
plot a1*day;
run;
proc arima data=newlib.a11;
identify var=a1(1) nlag=24;
run;
estimate p=1 ;
run;
forecast lead=6 interval=day id=day out=out;
run;
proc print data=out;
run;