1、對1950-2009年的新疆社會消費品零售總額的時間序列分析與預測利用1950-2009年的新疆社會消費品零售總額（記為：save，單位：萬元）的時間序列數(shù)據(jù)進行分析，建立時間序列ARIMA模型，并預測未來10年的社會消費品零售總額。表1 1950-2009年的新疆社會消費品零售總額年份社會消費品品零售總總額年份社會消費品品零售總總額年份社會消費品品零售總總額年份社會消費品品零售總總額1950219200196511297701980293599019952536447519512902331966121344919813407339199629535597195236646619671295

2、330198236413331997310411971953431988196812297711983413322419983275221019545221661969131311819844614339199934739958195561379919701323006198557384422000374499991956714644197113795581986638988120014063448719578557881972143411619877239113200244288871195892490019731546776198888698862003421166801959110522

3、61974158033519899814997200456365520196011905591975168488619901043004120056402000019611067880197618137771991121511802006733200001962105455419771934557199213824452200785750000196310083371978218866519931683773720081041550000196410540061979247799619941971008620091177553000數(shù)據(jù)來源：新疆疆統(tǒng)計年年鑒20010，新新疆五十十年模型應用d

4、ata a;inputt daate cosst; cardss; 1950 2199201951 2900231952 3666461953 4311981954 5222161955 6133791956 7144641957 8555781958 9244901959 110052661960 119905991961 106678001962 105545441963 100083771964 105540661965 112297001966 121134991967 129953001968 122297111969 131131881970 132230661971 137795

5、881972 143341661973 154467661974 158803551975 168848661976 181137771977 193345771978 218886551979 247779661980 293359001981 340073991982 364413331983 413332441984 461143991985 573384221986 638898111987 723391331988 886698661989 981149771990 1044304411991 1211518801992 1388245521993 1688373371994 197

6、7108861995 2533647751996 2955359971997 3100419971998 3277521101999 3477395582000 3744499992001 4066348872002 4422887712003 4211168802004 5633652202005 6400200002006 7333200002007 8577500002008 10441500002009 1177753300;procpprinntdaata=a; Run;程序說明：這段程程序是錄錄入19950年年到20009年年的新疆疆社會消消費品零零售總額額的數(shù)據(jù)據(jù)。data a;

7、set wworkk.a;procggplootdaata=a;plotccostt*daate;symboolv=doot ii=jooin c=bllackk l=1w=2; run;程序說明：這段程程序首先先建立數(shù)數(shù)據(jù)集aa，然后后調用ggploot過程程作圖來來顯示新新疆社會會消費品品零售總總額趨勢勢圖，結結果如下下圖。新疆社會消消費品零零售總額額的時序序圖從圖中發(fā)現(xiàn)現(xiàn)，隨著著時間的的增加波波動越來來越大，因因而此序序列是非非平穩(wěn)的的，從而而有必要要對原數(shù)數(shù)據(jù)進行行函數(shù)變變換，我我們這里里對原數(shù)數(shù)據(jù)取對對數(shù)，程程序如下下：data a;set wworkk.a;lc=loog(ccos

8、tt);procggplootdaata=a;plot lc*datte;symboolv=doot ii=jooin c=bllackk l=1w=2; run;新疆社會消消費品零零售總額額取對數(shù)數(shù)后的時時序圖上圖取對數(shù)數(shù)后的序序列還是是不平穩(wěn)穩(wěn)，但序序列的原原始數(shù)據(jù)據(jù)趨勢圖圖只能大大概的看看出此序序列的平平穩(wěn)性，為為了進一一步了解解新疆社社會消費費品零售售總額取取對數(shù)后后序列（lnin）的性質，下面給出lc序列的ACF，程序如下：procaarimmadaata=a;identtifyyvarr=lcc nllag=30essacffp=(0:6) q=(0:6) miiniccp=(0:

9、6) q=(0:6);run;新疆社會消消費品零零售總額額取對數(shù)數(shù)后的AACF從這個圖上上可以看看出，自自相關函函數(shù)緩慢慢衰減，從從而表明明序列l(wèi)lc是非非平穩(wěn)的的，這和和我們從從lc的的序列圖圖上觀察察的結果果是相同同的。因因此需要要對變量量lniin進行行一階差差分操作作并對差差分后的的序列進進行平穩(wěn)穩(wěn)性識別別，程序序如下：identtifyyvarr=lcc(1) nllag=30essacffp=(0:8) q=(0:8) miiniccp=(0:6) q=(0:6);run;識別過程結結果會給給出三個個可能不不同的模模型，分分別對這這三個模模型進行行估計，已已得到擬擬合最好好模型。模

10、型一差分后的對對數(shù)序列列的自相相關函數(shù)數(shù)圖差分后的對對數(shù)序列列的偏自相關關函數(shù)圖圖由圖可知，建建議對一一階差分分后的llc序列列擬合AARMAA（2,2）過過程，接接著進行行模型的的估計檢檢驗過程程，程序序如下：estimmateep=2q=2metthodd=mll;run;極大似然估估計結果果（顯著著性檢驗驗）結果發(fā)現(xiàn)，系系數(shù)除了了常數(shù)項項外其他他都顯著著為0，此此模型擬擬合效果果差，排排除這種種模型的的可能。模型二ESACFF的建議議模型通過ESAACF，建建議對一一階差分分后的llc序列列擬合AARMAA（5,1）過過程，接接著進行行模型的的估計檢檢驗過程程，程序序如下：estimma

11、teep=5q=1metthodd=mll;run;lc序列估估計模型型ariima（11,1,1）結結果此模型擬合合效果挺挺好，可可以考慮慮選擇這這個模型型。模型三最小信息準準則建議議模型由最小信息息準則，擬擬合模型型ARMMA（22，0），程程序如下下：estimmateep=2metthodd=mll;run;lc序列估估計模型型ariima（2,1,0）結果結果不僅各各系數(shù)顯顯著不為為0，而而且殘差差序列經(jīng)經(jīng)驗證也也是白噪噪聲序列列，符合合要求。AIC與SSBC所所計算出出來的值值越小,則代表表模型的的配適度度越佳，故故兩個模模型比較較，選擇擇ARIIMA(2,11,0)模型。所以最終

12、估估計結果果對應的的模型為為：下面利用我我們擬合合的模型型進行預預測，我我們預測測未來110年的的新疆社社會消費費品零售售總額，預預測程序序如下：foreccasttleaad=10idd=daate ouut=ppredd1;run;程序說明：調用aarimma過程程中的fforeecasst語句句來進行行預測，lleadd是指向向后預測測的步數(shù)數(shù)，idd指定變變量用來來識別觀觀測的周周期，iinteervaal指定定觀測間間隔的時時間，不不寫明的的時候默默認值為為1，iid和iinteervaal一般般一起使使用，將將預測結結果輸出出到數(shù)據(jù)據(jù)集prred11中，預預測結果果如下：對lc未來來五年的的預測結結果變換預測值值，以獲獲取原度度量下的的預測值值，程序序如下：data w;set ppredd1;y=expp(lcc);foreccastt1=eexp(forrecaast);run;procpprinntdaata=w;run;程序說明：對取對對數(shù)后的的序列l(wèi)lninn再進行行expp函數(shù)變變換，以以得到原原序列，同同時預測測十年的的值。YearCost6120101354770188.5