data score1;/*建立数据集*/
input car@@;
date=intnx('month','1jan07'd,_n_-1);
format date monyy.;
cards;
457085 331735 460839 461752 418692 432254 400209 419869 481016 415031 491995 544639
552812 397899 578484 503868 469345 491227 406125 387672 465885 455306 441098 514179
488746 474207 607472 658934 651909 708827 677754 717420 837848 777801 859530 948548
;
run;
proc gplot data=score1;/*绘制散点图*/
symbol1 i=join v=star;
plot car*date=1;
run;
data score2;/*差分*/
set score1;
dif1=dif(car);
run;
proc gplot data=score2;/*差分图*/
plot car*date dif1*date;
symbol v=star c=red i=join;
run;
proc arima data=score2;/*序列的初步识别*/
identify var=dif1 nlag=12;
run;
proc arima data=score2;/*参数估计、假设检验*/
identify var=dif1;
estimate p=1 q=1 method=ml;
forecast id=date lead=5 out=out;
run;
proc gplot data=score2;/*图像*/
plot dif1*date=2 forecast*date=3 l95*date=4 u95*date=4/overlay;
symbol2 v=star i=none c=black ;
symbol3 v=none i=join c=red w=2;
symbol4 v=none i=join c=green l=2 ;
run;
各位大侠,我第一次学计量,参考资料很不方便,这是我根据根据网上的一些程序改编的ARIMA预测的过程,希望大家帮我看看运行的结果,指出我应该改进的地方;还有预测程序我不会写,哪位大哥帮我弄下,万分感谢!