标题: 各位英雄，帮帮忙吧~~~ [打印本页] 作者: shiyiming 时间: 2009-3-13 20:28 标题: 各位英雄，帮帮忙吧~~~ 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.
[b:1nrezqfa]各位，我做出来个时间序列的预测，可是预测值在和真实值反映的一张图的时候总是要比真实值延后一期，并且出现上面的提示[/b:1nrezqfa]
请问该如何解决？
谢谢作者: shiyiming 时间: 2009-3-14 09:15 标题: Re: 各位英雄，帮帮忙吧~~~ Give your code. then maybe we can try to figure it out.
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.
these warning is not a big deal, if the actual data are not evenly spaced.作者: shiyiming 时间: 2009-3-15 13:23 标题: Re: 各位英雄，帮帮忙吧~~~ thanks for your advice,and my code is here,and my personal opinion of the model is qutie well,because the result of the forecast compared with the real figure seems suitable:
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;