ata a;
input x@@;
difx=dif(x);
difxx=dif4(dif(x));
t=intnx ('quarter','1jan1998'd,_n_-1 );
format t yyq4.;
cards;
180.1        433         305.7        226.9        240.3        509.2        281.2        264.2        284         623.9        476.3        392.2
463.9        903.5        541.4        358.4        615.7        1296.2        689.6        500.3        985.3        1760        928.1        1095.4
1450.3        2463.5        1696.2        1395.5        1741.3        2923.5        1257.7        1511.7        1929.8        2871.9        1562.9        1522.8
2060.1        3360.5        1962.5        1433.6        2285.6        3235.5        1843.6        1152.5        1740.4        2864.7        2026.1        1439.1
2264.7        3711.9        2626.9        2303        3594.6        4770.5        3002.5        2420.4        3312.4        4473.8        2795.2        2653.6
3500.7        4443.9        3456.8        2655.5
;
proc gplot data=a;
plot x*t difx*t difxx*t;
symbol c=black v=star i=join;
run;
proc arima data=a;
identify var=x(1,4) nlag=8 minic p=(0:5) q=(0:5);
estimate q=1;
forecast lead=8 id=t out=results;
proc gplot data=results;
plot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay;
symbol1 c=black v=star i=none;
symbol2 c=red v=none i=join;
symbol3 c=green v=none i=join l=2;
run;