SAS中文论坛

 找回密码
 立即注册

扫一扫,访问微社区

查看: 713|回复: 2
打印 上一主题 下一主题

各位英雄,帮帮忙吧~~~

[复制链接]

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
楼主
 楼主| 发表于 2009-3-13 20:28:06 | 只看该作者

各位英雄,帮帮忙吧~~~

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]
请问该如何解决?
谢谢
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
沙发
 楼主| 发表于 2009-3-14 09:15:41 | 只看该作者

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.
回复 支持 反对

使用道具 举报

49

主题

76

帖子

1462

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
1462
板凳
 楼主| 发表于 2009-3-15 13:23:50 | 只看该作者

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;
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|SAS中文论坛  

GMT+8, 2026-2-5 10:29 , Processed in 0.076692 second(s), 19 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表