sumgawing
前面的没问题,就是最后画图那不行
data sasuser.example6_1;
input x t;
cards;
100 1952
115.1 1953
116 1954
120.3 1955
128.7 1956
131.8 1957
126.2 1958
127.2 1959
113.8 1960
99.3 1961
102 1962
118.7 1963
132 1964
144.3 1965
146.5 1966
150 1967
148.6 1968
152.9 1969
158.3 1970
161.8 1971
177.1 1972
182.8 1973
186.8 1974
194 1975
202.9 1976
210.4 1977
217.3 1978
228.6 1979
243.1 1980
268.8 1981
269.8 1982
276.7 1983
296.2 1984
320.4 1985
343.7 1986
373.8 1987
407 1988
400.6 1989
430.5 1990
470.6 1991
543.1 1992
591.3 1993
601.2691 1994
629.0835 1995
644.7291 1996
667.111 1997
718.1765 1998
;
proc gplot data=sasuser.example6_1;
plot x*t;
symbol i=join;
run;
data sasuser.example6_1;
set sasuser.example6_1;
dif=dif(x);
proc gplot data= sasuser.example6_1;
plot dif*t;
symbol c=red v=star i=join;
run;
proc arima data= sasuser.example6_1;
identify var=x(1);
run;
proc arima data= sasuser.example6_1;
identify var=x(1);
estimate p=(1 5) noint;
forecast lead=5 id=t out=out;
run;
proc gplot data=sasuser.example6_1;
plot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay;
symbol1 c=black i=join v=diamond;
symbol2 c=red i=join v=dot;
symbol3 c=green i=join v=star l=3 h=3;
run;