SAS讲义十二课非平稳序列随机分析研究

个人收集整理 仅供参考学习

OBS DATE XLOG FORECAST STD L95 U95 RESIDUAL 140 213 6.40688 6.43874 0.044559 6.35141 6.52607 -0.03186 141 244 6.01127 6.24129 0.044559 6.15396 6.32863 -0.23003 142 274 6.13340 5.99662 0.044559 5.90928 6.08395 0.13678 143 305 5.96615 5.93671 0.044559 5.84938 6.02405 0.02943 144 335 6.06843 6.05720 0.044559 5.96987 6.14453 0.01123 145 366 . 6.09873 0.044559 6.01139 6.18606 . 146 397 . 6.04148 0.049764 5.94395 6.13902 . 147 425 . 6.15726 0.054474 6.05049 6.26403 . 148 456 . 6.18807 0.058808 6.07281 6.30333 . 149 486 . 6.22160 0.062843 6.09843 6.34477 . 150 517 . 6.35699 0.066635 6.22638 6.48759 . 151 547 . 6.49668 0.070222 6.35904 6.63431 . 152 578 . 6.49172 0.073635 6.34740 6.63604 . 153 609 . 6.20952 0.076896 6.05880 6.36023 . 154 639 . 6.19788 0.080025 6.04104 6.35473 . 155 670 . 6.05147 0.083035 5.88873 6.21422 . 156 700 . 6.15565 0.085941 5.98720 6.32409 . l95=exp(l95); u95=exp(u95); proc print data=arimad04; run;

proc gplot data=arimad04;

where date>='1jan57'd;

plot x*date forecast*date l95*date u95*date

/ove

29 / 33

个人收集整理 仅供参考学习

rlay vaxis=axis1 haxis

30 / 33

个人收集整理 仅供参考学习

=axis2 href='31dec60'd

31 / 33

个人收集整理 仅供参考学习

;

symbol1 i=join v=C h=2.5 l=1 font=swissb

c=red;

symbol2 i=join v=F h=3 l=1 font=swissb

c=blue;

symbol3 i=join l=1 font=swissb c=green; symbol4 i=join l=1 font=swissb c=green; axis1 label=('Passenger') order=(250 to 800

by 50);

axis2 label=('Month') order=('1jan57'd

to '1jan62'd by year);

format date monyy. ; title1 'Forecast Chart'; title2 'C--x'; title3 'F--forecast'; title4 'None--u95 and l95'; run;

32 / 33