应用时间序列实验报告(6)

实验二程序： data example4_2;

input huoyunliang@@;

difhuoyunliang=dif(huoyunliang);

time=intnx ('year','01JAN1949'd,_n_-1); format time monyy7.; cards;

5589 9983 11083 13217 16131 19288 19376 24605 27421 38109 54410 67219 44988 35261 36418 41786 49100 54951 43089 42095 53120 68132 76471 80873 83111 78772 88955 84066 95309 110119 111893 111279 107673 113495 118784 124074 130709 135635 140653 144948 151489 150681 152893 157627 162794 163216 165982 171024 172149 164309 167554 178581 193189 204956 224248 249017 269296 288224 314237 330354 ;

/*proc print data=example4_2;*/ /*proc gplot;*/

/*plot difhuoyunliang*time=2;*/ /*symbol1 v=star c=black i=join;*/ /*symbol2 c=black i=join v=star;*/ proc arima data=example4_2;

identify var=difhuoyunliang(1) nlag=8 minic p=(0:5) q=(0:5); estimate q=2;

forecast lead=5 id=time

23

实验三程序： data example4_3; input death@@;

time=intnx ('month','01JAN1973'd,_n_-1); format time monyy7.; cards;

9007 8106 8928 9137 10017 10826 11317 10744 9713 9938 9161 8927 7750 6981 8038 8422 8714 9512 10120 9823 8743 9129 8710 8680 8162 7306 8124 7870 9387 9556 10093 9620 8285 8433 8160 8034 7717 7461 7776 7925 8634 8945 10078 9179 8037 8488 7874 8647 7792 6957 7726 8106 8890 9299 10625 9302 8314 8850 8265 8796 7836 6892 7791 8129 9115 9434 10484 9827 9110 9070 8633 9240 ;

/*proc gplot data=example4_3;*/ /*plot death*time=1;*/

/*symbol1 c=black i=join v=star;*/ proc autoreg data=example4_3; /*model death=time/ dwprob;*/

model death=time/ nlag=5 backstep method=ml noint ; output out =out p=xp pm=trend; proc gplot data=out;

plot death*time=2 xp*time=3 trend*time=4 / overlay; symbol2 v=star i=none c=blak;

symbol3 v=none i=join c=red w=2 l=3; symbol4 v=none i=join c=green w=2; run;

24

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库应用时间序列实验报告(6)在线全文阅读。