1、應(yīng)用時間序列分析第4 章答案河南大學(xué) :姓名 : 汪寶班級 : 七班學(xué)號 :1122314451班級序號 :68我國 1949 年 2008 年年末人口總數(shù) ( 單位 : 萬人 ) 序列如表 48 所示 ( 行數(shù)據(jù) ) 選擇適當(dāng)?shù)哪Ｐ蛿M合該序列的長期數(shù)據(jù)，并作 5 期預(yù)測。解：具體解題過程如下： （本題代碼我是做一問寫一問的）觀察時序圖 :datawangbao4_5;inputx;time=1949 +_n_-1;cards;54167 55196 56300 57482 58796 60266 61465 6282864653 65994 67207 66207 65859 67295 6

2、9172 7049972538 74542 76368 78534 80671 82992 85229 8717789211 90859 92420 93717 94974 96259 97542 98705100072 101654 103008 104357 105851 107507 109300 111026112704 114333 115823 117171 118517 119850 121121 122389123626 124761 125786 126743 127627 128453 129227 129988130756 131448 132129 132802;pro

3、cgplotdata=wangbao4_5;plotx*time=1 ;symbol1c =blackv =stari =join;run ;分析 : 通過時序圖 , 我可以發(fā)現(xiàn)我國 1949 年 2008 年年末人口總數(shù) ( 隨時間的變化呈現(xiàn)出線性變化 . 故此時我可以用線性模型擬合序列的發(fā)展X t =a+bt +Itt=1,2,3,60tt2其中， I t 為隨機波動；Xt =a+b就是消除隨機波動的影響之后該序列的長期趨勢。進行線性模型擬合：procautoregdata=wangbao4_5;modelx=time;outputout =outp =wangbao4_5_cup;run

4、 ;procgplotdata=out;plotx*time=1 wangbao4_5_cup*time=symbol2c =redv =nonei =joinw= 2run ;2/ overlay l = 3;分析：由上面輸出結(jié)果可知：兩個參數(shù)的p 值明顯小于0.05，即這兩個參數(shù)都是具有顯著非零，：模型檢驗又因為Regress R-square=total R-square=0.9931,即擬合度達到99.31% 所以用這個模型擬合的非常好。：結(jié)論所以本題擬合的模型為：X t =-2770828+1449t+I t t=1,2,3,60E(It )=0,var(It )= 26：作 5 期

5、預(yù)測procforecastdata=wangbao4_5 method=stepar trend=2 lead= 5out=out outfull outest=est;id t;varx;procgplotdata=out;plot x*time=_type_/href= 2008;symbol1i =nonev =starc =black;symbol2i =joinv =nonec =red;symbol3i =joinv =nonec =blackl =2 ;symbol4i =joinv =nonec =blackl =2 ;run;6：愛荷華州1948-1979 年非農(nóng)產(chǎn)品季度收

6、入數(shù)據(jù)如表49 所示（行數(shù)據(jù)），選擇適當(dāng)?shù)哪Ｐ蛿M合該序列的長期趨勢。解：具體做題過程如下： （本題代碼我是做一問寫一問的）1、繪制時序圖datawangbao4_6;inputx;time=_n_;cards;601 604 620 626 641 642 645 655 682 678 692 707736 753 763 775 775 783 794 813 823 826 829 831830 838 854 872 882 903 919 937 927 962 975 9951001 1013 1021 1028 1027 1048 1070 1095 1113 1143 1154

7、 11731178 1183 1205 1208 1209 1223 1238 1245 1258 1278 1294 13141323 1336 1355 1377 1416 1430 1455 1480 1514 1545 1589 16341669 1715 1760 1812 1809 1828 1871 1892 1946 1983 2013 20452048 2097 2140 2171 2208 2272 2311 2349 2362 2442 2479 25282571 2634 2684 2790 2890 2964 3085 3159 3237 3358 3489 3588

8、3624 3719 3821 3934 4028 4129 4205 4349 4463 4598 4725 48274939 5067 5231 5408 5492 5653 5828 5965;procgplotdataplotx*time;symbolc =blackrun ;=wangbao4_6;v =stari =join;分析 ; 可知時序圖顯示該序列有明顯的曲線遞增趨勢。嘗試使用修正指數(shù)型模型進行迭代擬合：x tabc t +?t ，t=1，2， 1282、擬合模型procnlinmethod=gauss;modelx=a+b*c*time;parametersa= 0.1b=

9、 0.1c= 1.1 ;outputpredicted=xhatout =out;run ;NLIN 過程輸出以下六方面信息：（1）迭代過程（2）收斂狀況（本次迭代收斂）（3）估計信息摘要4）主要統(tǒng)計量5）參數(shù)信息摘要得到的擬合模型為：xt604.8112.21.0307 ttt=1，2， 128（6）近似相關(guān)矩陣3、擬合效果為了直觀看出擬合效果， 我們可以將原序列值和擬合值聯(lián)合作圖：procgplotdata=out;plotx*t=1 xhat*t=2 / overlay;symbol1c =blackv =stari =join;symbol2c =redv =nonei =join;分

10、析：由上圖圖我們可以看出，原序列值和擬合值很接近，擬合效果較好。綜合以上的分析，我們可以選擇模型：xt 604.8 112.2 1.0307 tt 來擬合該序列的長期趨勢。擬合效果很不錯。某城市 1980 年 1 月至 1995 年 8 月每月屠宰生豬的數(shù)量（單位：頭）如表 4 11 所示（行數(shù)據(jù)），選擇適當(dāng)?shù)哪Ｐ蛿M合該序列的發(fā)展， 并預(yù)測 1995 年9月至 1997 年9月該城市的生豬屠宰量。解： 具體解題過程如下：（本題代碼我是做一問寫一問的）atawangbao4_8;inputx;time=_n_;cards;76378 71947 33873 96428 105084 95741

11、110647 100311 94133 10305590595 101457 76889 81291 91643 96228 102736 100264 103491 9702795240 91680 101259 109564 76892 85773 95210 93771 98202 97906100306 94089 102680 77919 93561 117062 81225 88357 106175 91922104114 109959 97880 105386 96479 97580 109490 110191 90974 98981107188 94177 115097 113

12、696 114532 120110 93607 110925 103312 120184103069 103351 111331 106161 111590 99447 101987 85333 86970 10056189543 89265 82719 79498 74846 73819 77029 78446 86978 7587869571 75722 64182 77357 63292 59380 78332 72381 55971 6975085472 70133 79125 85805 81778 86852 69069 79556 88174 6669872258 73445 7

13、6131 86082 75443 73969 78139 78646 66269 7377680034 70694 81823 75640 75540 82229 75345 77034 78589 7976975982 78074 77588 84100 97966 89051 93503 84747 74531 9190081635 89797 81022 78265 77271 85043 95418 79568 103283 9577091297 101244 114525 101139 93866 95171 100183 103926 102643 10838797077 9090

14、1 90336 88732 83759 99267 73292 78943 94399 9293790130 91055 106062 103560 104075 101783 93791 102313 82413 83534109011 96499 102430 103002 91815 99067 110067 101599 97646 104930 88905 89936 106723 84307 114896 106749 87892 100506 ;procgplotdata=wangbao4_8;plotx*time=1 ;symbol1c =redi =joinv =star;r

15、un ;procarimadata=wangbao4_8;identifyvar =x;run ;1：時序圖與平穩(wěn)性判別分析：上圖是數(shù)據(jù)對應(yīng)的時序圖，從圖上曲線分析來看，數(shù)據(jù)并沒有周期性或者趨向性規(guī)律， 并且每月的生豬的屠宰量大約在 80000 上下波動。所以由該序列圖我可以認(rèn)為它是個平穩(wěn)的數(shù)列。即可以用第三章的 AR模型或 MA模型或 ARMA模型進行擬合。但是為了穩(wěn)妥起見，我還需要利用自相關(guān)圖進一步輔助識別。具體如下：自相關(guān)圖：分析：由上圖可知：樣本自相關(guān)圖中的自相關(guān)系數(shù)在延遲 4 階之后幾乎全部落入 2 個標(biāo)準(zhǔn)差范圍之內(nèi)，并且向零衰減的速度還是比較快的。所以我認(rèn)為該序列是平穩(wěn)的序列。由時

16、序圖與自相關(guān)圖可知其是平穩(wěn)序列。 故可以用第三章的AR模型或 MA模型或 ARMA模型進行擬合分析：由上面數(shù)據(jù)可知： 由于 p 值顯著小于 0.05 ，故可以否定原假設(shè)（ H0），接受備選假設(shè) (H1) ，即我可以認(rèn)為該序列是平穩(wěn)的非純隨機性序列。 這樣就說明了我們可以根據(jù)歷史信息預(yù)測未來月的生豬屠宰量。分析：觀察自相關(guān)圖和偏自相關(guān)圖，從這兩圖來看，偏自相關(guān)圖是拖尾，而自相關(guān)系數(shù)是拖尾的。 因而我們可以用 ARMA 模型進行擬合。但是為了穩(wěn)妥起見，我還需要利用計算機進行相對最優(yōu)定階。：相對最優(yōu)定階：identifyvar =xnlag=18minicp =( 0 : 5)q=( 0: 5 );

17、run ;分析：從上圖可以看出，在眾多模型中，ARMA模型的 BIC信息量是最小的是 ARMA（4 ，5），因而我們接下來會采用ARMA模型來進行擬合分析，這和我們?nèi)斯ゎA(yù)測的相吻合。：參數(shù)估計：estimatep=4 q=5;run ;具體輸出結(jié)果如下圖 :該輸出形式等價于xt =(1-1.21457B+0.70228B 2-0.04985B 3-0.41243B 4) ?t 故該模型為： xt =+( （ B）/ (B) ?t =(3) 序列預(yù)測(1995年9月至1997年9月)forecastlead=24id =timeout =wangbaoyc;run ;分析：以上是我們對數(shù)據(jù)進行了 2年24 期的