SAS中文论坛

标题: 江湖告急！！！ [打印本页]

作者: shiyiming 时间: 2009-3-15 17:39
标题: 江湖告急！！！
WARNING: The intervals on the axis labeled date are not evenly spaced.
WARNING: No minor tick marks will be drawn because major tick increments have been specified in uneven or unordered intervals.
NOTE: 12 observation(s) contained a MISSING value for the cpi * date request.
NOTE: 1 observation(s) contained a MISSING value for the FORECAST * date request.
各位，我做出来个时间序列的预测，可是预测值在和真实值反映的一张图的时候总是要比真实值延后一期，并且出现上面的提示
请问该如何解决？
thanks for tangyongsas advice,and my code is here

data zlh;
input cpi@@;
date=intnx('month','1jan98'd,_n_-1);
format date monyy.;
cards;
100.3 99.9 100.7 99.7 99 98.7 98.6 98.6 98.5 98.9 98.8 99
98.8 98.7 98.2 97.8 97.8 97.9 98.6 98.7 99.2 99.4 99.1 99
99.8 100.7 99.8 99.7 100.1 100.5 100.5 100.3 100 100 101.3 101.5
101.2 100 100.8 101.6 101.7 101.4 101.5 101 99.9 100.2 99.7 99.7
99 100 99.2 98.7 98.9 99.2 99.1 99.3 99.3 99.2 99.3 99.6
100.4 100.2 100.9 101 100.7 100.3 100.5 100.9 101.1 101.8 103 103.2
103.2 102.1 103 103.8 104.4 105 105.3 105.3 105.2 104.3 102.8 102.4
101.9 103.9 102.7 101.8 101.8 101.6 101.8 101.3 100.9 101.2 101.3 101.6
101.9 100.9 100.8 101.2 101.4 101.5 101 101.3 101.5 101.4 101.9 102.8
102.2 102.7 103.3 103 103.4 104.4 105.6 106.5 106.2 106.5 106.9 106.5
107.1 108.7 108.3 108.5 107.7 107.1 106.3 104.9 104.6 104 102.4
;
run;
proc gplot data=Zlh;
symbol1 i=join v=dot c=red;
plot cpi*date=1;
run;
data cpi;
set zlh;
difcpi=dif(cpi);
proc gplot;
plot cpi*date difcpi*date;
symbol v=star c=black i=join;
run;
proc arima;
identify var=cpi(1) nlag=24;
estimate p=(1)q=(1)(12) noint ;
forecast lead=12 interval=month id=date out=b;
run;
data cpi;
set b;
proc gplot data=cpi;
symbol1 I=join v=star r=1 c=red;
symbol2 I=join v=plus r=1 c=green;
proc gplot data=cpi;
where date>='1jan98'd;
plot cpi*date=1 forecast*date=2/
overlay haxis='1jan98'd to '1jan10'd by year;
run;

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)