SAS建立時(shí)間序列模型

1、SAS軟件簡(jiǎn)介 SAS系統(tǒng)系統(tǒng)是由美國(guó)是由美國(guó)SAS軟件研究所開(kāi)發(fā)的用于軟件研究所開(kāi)發(fā)的用于決策支持的大型集成信息系統(tǒng)，是數(shù)決策支持的大型集成信息系統(tǒng)，是數(shù)據(jù)處理和統(tǒng)計(jì)領(lǐng)域的國(guó)際標(biāo)準(zhǔn)軟件之據(jù)處理和統(tǒng)計(jì)領(lǐng)域的國(guó)際標(biāo)準(zhǔn)軟件之一，廣泛應(yīng)用于金融、醫(yī)藥衛(wèi)生、生一，廣泛應(yīng)用于金融、醫(yī)藥衛(wèi)生、生產(chǎn)、運(yùn)輸、通訊、政府、教育和科研產(chǎn)、運(yùn)輸、通訊、政府、教育和科研等領(lǐng)域。等領(lǐng)域。應(yīng)用應(yīng)用SAS軟件建立時(shí)間序列模型軟件建立時(shí)間序列模型 準(zhǔn)備工作：建立一個(gè)時(shí)間序列數(shù)據(jù)集準(zhǔn)備工作：建立一個(gè)時(shí)間序列數(shù)據(jù)集SAS語(yǔ)句：語(yǔ)句：Data 數(shù)據(jù)集名；數(shù)據(jù)集名； Input 序號(hào)（序號(hào)（year or month）變量名）變量名

2、 ; Cards；/（輸入數(shù)據(jù)，按（輸入數(shù)據(jù)，按input格式逐個(gè)輸入數(shù)格式逐個(gè)輸入數(shù)據(jù)，以分號(hào)結(jié)束）；據(jù)，以分號(hào)結(jié)束）； Proc print data數(shù)據(jù)集名；數(shù)據(jù)集名；/輸出數(shù)據(jù)表輸出數(shù)據(jù)表Run；例例1 磨輪剖面資料磨輪剖面資料data li; input x ; cards;13.5 4.0 4.0 4.5 3.0 3.0 10.0 10.2 9.0 10.0 8.5 7.0 10.5 7.5 7.0 10.5 9.5 7.0 12.0 13.5 12.5 15.0 13.0 11.0 9.0 10.5 10.5 11.5 10.5 9.0 8.2 8.5 9.2 8.5 10.0

3、14.5 13.0 2.0 6.0 6.0 11.0 9.5 12.5 13.8 12.0 12.0 12.0 13.0 12.0 14.0 14.5 13.5 12.3 7.0 7.0 7.0 6.5 12.5 15.0 12.5 11.6 11.0 10.0 8.5 3.0 11.5 11.5 11.5 11.0 9.0 2.5 7.0 6.0 6.6 14.0 11.0 9.0 6.5 4.0 6.0 12.0 11.0 12.0 12.5 12.5 13.6 13.0 8.0 6.5 6.8 6.0 7.2 10.2 8.0 7.5 11.0 11.8 11.8 6.5 8.0 9.0

4、 8.0 8.0 9.0 9.5 10.0 9.0 12.0 13.5 13.8 15.0 12.5 11.0 11.5 14.5 11.5 11.8 13.0 15.0 14.5 13.0 9.0 11.0 9.0 10.0 14.0 13.5 3.0 2.2 6.0 8.0 9.0 9.0 9.0 7.0 6.0 6.5 7.0 7.5 8.5 9.0 9.5 10.0 11.5 11.2 12.5 11.6 8.0 7.0 6.0 6.0 6.0 9.0 12.0 13.5 13.0 3.5 1.8 1.6 7.5 8.0 7.9 11.6 12.5 10.5 8.0 9.0 11.6

5、11.8 12.6 10.2 10.0 5.0 7.0 -1.0 0.0 0.0 3.0 11.0 12.0 12.2 11.0 8.0 7.0 5.5 10.0 11.5 7.0 4.0 7.0 7.0 10.0 9.0 8.0 10.0 13.0 10.0 6.5 11.0 13.0 13.0 14.0 13.0 12.5 12.0 9.0 8.5 7.0 8.5 10.0 8.0 4.0 3.0 10.0 13.0 13.0 13.0 12.0 11.0 11.0 11.0 14.5 14.0 14.0 13.5 10.0 9.5 10.0 12.5 10.0 9.0 9.0 4.0 3

6、.0 6.0 5.0 7.0 6.0 5.0 8.510.5 11.1 11.0 10.0 11.2 8.0 2.5 5.0 13.2 14.0;Proc print data=li;run;data seriesg; input x ; xlog=log(x); date=intnx(month,31dec1948d,_n_); format date monyy.; cards; 112 118 132 129 121 135 148 148 136 119 104 118 115 126 141 126 126 149 170 170 158 133 114 140 145 150 17

7、8 163 172 178 199 199 184 162 146 166 171 180 193 181 183 218 230 242 209 191 172 194 196 196 236 235 229 243 264 272 237 211 180 201 204 188 235 227 234 264 302 233 259 229 203 229 242 233 267 269 270 316 364 347 312 274 237 278 284 277 317 313 318 374 413 405 355 306 271 306 315 301 356 348 356 42

8、2 465 467 404 347 305 336 340 318 362 348 363 435 491 505 404 359 310 337 360 342 406 396 420 472 548 559 463 407 362 405 417 391 419 461 472 535 622 606 508 461 390 432；Proc print data= seriesg;run;例例2 國(guó)航客票數(shù)國(guó)航客票數(shù)(Airline)數(shù)據(jù)數(shù)據(jù)SAS的建模步驟的建模步驟 SAS建模根據(jù)建模根據(jù)BoxJenkins建模方法建模方法,主要包括三個(gè)階段：主要包括三個(gè)階段：模型識(shí)別階段（包括模型定

9、階）模型識(shí)別階段（包括模型定階）模型參數(shù)估計(jì)階段（包括模型檢驗(yàn)）模型參數(shù)估計(jì)階段（包括模型檢驗(yàn)）模型的預(yù)測(cè)階段模型的預(yù)測(cè)階段第一階段第一階段: 模型的識(shí)別模型的識(shí)別 平穩(wěn)性模型識(shí)別平穩(wěn)性模型識(shí)別首先判定時(shí)間序列數(shù)據(jù)是否為平穩(wěn)隨機(jī)數(shù)首先判定時(shí)間序列數(shù)據(jù)是否為平穩(wěn)隨機(jī)數(shù)據(jù)，據(jù)，（一）通過(guò)（一）通過(guò)時(shí)間序列數(shù)據(jù)趨勢(shì)圖時(shí)間序列數(shù)據(jù)趨勢(shì)圖判別。判別。Sas語(yǔ)句：語(yǔ)句：symbol1 i=join v=star;proc gplot data=seriesg;plot x*date=1/haxis=1jan49d to 1jan61d by year;run;（二）通過(guò)自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的截（二）通

10、過(guò)自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的截尾性識(shí)別模型尾性識(shí)別模型“IDENTIFY”語(yǔ)句語(yǔ)句通過(guò)通過(guò)SAS軟件，運(yùn)行程序如下：軟件，運(yùn)行程序如下：proc arima data=數(shù)據(jù)集數(shù)據(jù)集 identify var=變量名變量名 nlag=時(shí)間間隔個(gè)數(shù)時(shí)間間隔個(gè)數(shù)run; 計(jì)算出自相關(guān)系數(shù)計(jì)算出自相關(guān)系數(shù)ACF, 逆自相關(guān)系數(shù)逆自相關(guān)系數(shù)SIACF, 偏自相關(guān)系數(shù)偏自相關(guān)系數(shù)PACF和互相關(guān)系數(shù)。和互相關(guān)系數(shù)。根據(jù)樣本自相關(guān)系數(shù)根據(jù)樣本自相關(guān)系數(shù)ACF和偏相關(guān)系數(shù)和偏相關(guān)系數(shù)PACF的形態(tài)來(lái)識(shí)別模型類(lèi)別。的形態(tài)來(lái)識(shí)別模型類(lèi)別。 如果序列的樣本自相關(guān)系數(shù)在如果序列的樣本自相關(guān)系數(shù)在q步后截尾，則是步后截

11、尾，則是MA序列，如果偏相關(guān)系數(shù)在序列，如果偏相關(guān)系數(shù)在p步后截尾，則是步后截尾，則是AR序列。序列。如果都不截尾，只是按負(fù)指數(shù)衰減或以阻尼正弦波形式如果都不截尾，只是按負(fù)指數(shù)衰減或以阻尼正弦波形式趨于零（即是拖尾的），則應(yīng)判斷為趨于零（即是拖尾的），則應(yīng)判斷為ARMA序列，但是序列，但是不能確定階次。不能確定階次。 若序列的樣本自相關(guān)和偏相關(guān)系數(shù)都不截尾，而且若序列的樣本自相關(guān)和偏相關(guān)系數(shù)都不截尾，而且至少有一個(gè)不是拖尾，即下降趨勢(shì)很慢，不能被負(fù)指數(shù)至少有一個(gè)不是拖尾，即下降趨勢(shì)很慢，不能被負(fù)指數(shù)函數(shù)所控制，或是不具有下降的趨勢(shì)而是周期變化，那函數(shù)所控制，或是不具有下降的趨勢(shì)而是周期變化，那

12、么我們便認(rèn)為序列具有增長(zhǎng)趨勢(shì)或季節(jié)性變化，是非平么我們便認(rèn)為序列具有增長(zhǎng)趨勢(shì)或季節(jié)性變化，是非平穩(wěn)序列?？蓱?yīng)用提取趨勢(shì)性和季節(jié)性的方法，對(duì)數(shù)據(jù)進(jìn)穩(wěn)序列?？蓱?yīng)用提取趨勢(shì)性和季節(jié)性的方法，對(duì)數(shù)據(jù)進(jìn)行處理，就是主要通過(guò)差分等變換將非平穩(wěn)序列變成一行處理，就是主要通過(guò)差分等變換將非平穩(wěn)序列變成一個(gè)平穩(wěn)序列。個(gè)平穩(wěn)序列。非平穩(wěn)序列的平穩(wěn)化非平穩(wěn)序列的平穩(wěn)化 若序列是非平穩(wěn)的，下面是通過(guò)差分變換變成一個(gè)平穩(wěn)若序列是非平穩(wěn)的，下面是通過(guò)差分變換變成一個(gè)平穩(wěn)序列。序列。SAS的程序?yàn)橐浑A差分的程序?yàn)橐浑A差分 變量名（變量名（1） identify var=變量（變量（1） nlog=N ； run; 若一階差

13、分是平穩(wěn)的，對(duì)差分序列建模，觀測(cè)若一階差分是平穩(wěn)的，對(duì)差分序列建模，觀測(cè)ACF、PACF的變化趨勢(shì)，初步給出的階數(shù)。因?yàn)檩斎霐?shù)據(jù)是序的變化趨勢(shì)，初步給出的階數(shù)。因?yàn)檩斎霐?shù)據(jù)是序列的有限樣本，所以由輸入序列計(jì)算出樣本自相關(guān)系數(shù)是列的有限樣本，所以由輸入序列計(jì)算出樣本自相關(guān)系數(shù)是逼近產(chǎn)生序列的理論自相關(guān)系數(shù)。這意味著樣本自相關(guān)系逼近產(chǎn)生序列的理論自相關(guān)系數(shù)。這意味著樣本自相關(guān)系數(shù)不能夠恰好等于任何模型的理論自相關(guān)系數(shù)，并且可能數(shù)不能夠恰好等于任何模型的理論自相關(guān)系數(shù)，并且可能會(huì)具有一種或多種不同的模型的理論自相關(guān)系數(shù)相似的類(lèi)會(huì)具有一種或多種不同的模型的理論自相關(guān)系數(shù)相似的類(lèi)型。若一階差分序列仍不平

14、穩(wěn)，重復(fù)以上過(guò)程，（二階差型。若一階差分序列仍不平穩(wěn)，重復(fù)以上過(guò)程，（二階差分，三階差分等等）直到差分序列平穩(wěn)。分，三階差分等等）直到差分序列平穩(wěn)。 由于大多數(shù)時(shí)間序列是非平穩(wěn)的，所以必須在進(jìn)行建模之由于大多數(shù)時(shí)間序列是非平穩(wěn)的，所以必須在進(jìn)行建模之前，就應(yīng)該把它們變換為平穩(wěn)序列。前，就應(yīng)該把它們變換為平穩(wěn)序列。 對(duì)數(shù)變換對(duì)數(shù)變換 ylog=log(y) 差分差分 identify var=變量名（變量名（1） /1為一階差分為一階差分季節(jié)性也會(huì)導(dǎo)致非平穩(wěn)，若數(shù)據(jù)是月記錄序列季節(jié)性也會(huì)導(dǎo)致非平穩(wěn)，若數(shù)據(jù)是月記錄序列 identify var=數(shù)據(jù)名（數(shù)據(jù)名（12） identify var=數(shù)

15、據(jù)名（數(shù)據(jù)名（1，1）/二階差分二階差分 identify var=數(shù)據(jù)名（數(shù)據(jù)名（1，12） var=數(shù)據(jù)名（數(shù)據(jù)名（1，12）是當(dāng)前相鄰時(shí)刻的差分以及相隔）是當(dāng)前相鄰時(shí)刻的差分以及相隔12個(gè)時(shí)刻的差分，對(duì)于差分階數(shù)的指定和差分間隔長(zhǎng)度的指?jìng)€(gè)時(shí)刻的差分，對(duì)于差分階數(shù)的指定和差分間隔長(zhǎng)度的指定是沒(méi)有限制的。定是沒(méi)有限制的。例例2 國(guó)航客票數(shù)國(guó)航客票數(shù)(Airline)數(shù)據(jù)數(shù)據(jù)proc arima data=seriesg; identify var=xlog(1,12) nlag=15; run; 差分不僅影響用于差分不僅影響用于IDENTIFY語(yǔ)句輸出的語(yǔ)句輸出的序列，而且應(yīng)用于任何隨后的序

16、列，而且應(yīng)用于任何隨后的ESTIMATE和和FORECAST語(yǔ)句。語(yǔ)句。ESTIMATE語(yǔ)句對(duì)差分序語(yǔ)句對(duì)差分序列擬模型，列擬模型，F(xiàn)ORECAST語(yǔ)句預(yù)測(cè)差分并自動(dòng)語(yǔ)句預(yù)測(cè)差分并自動(dòng)把差分加起來(lái)以取消有把差分加起來(lái)以取消有IDENTIFY語(yǔ)句指定的語(yǔ)句指定的差分操作。差分操作。第一階段第一階段(IDENTIFY)的輸出的輸出描述統(tǒng)計(jì)量：描述統(tǒng)計(jì)量：N E(X) D(X) (X) 自相關(guān)系數(shù)散點(diǎn)圖自相關(guān)系數(shù)散點(diǎn)圖 ACF 自相關(guān)系數(shù)圖表及序列的當(dāng)前值和過(guò)去值的相關(guān)自相關(guān)系數(shù)圖表及序列的當(dāng)前值和過(guò)去值的相關(guān)程度，圖中以圖像的形式顯示相關(guān)系數(shù)的值。程度，圖中以圖像的形式顯示相關(guān)系數(shù)的值。偏相關(guān)系數(shù)

17、偏相關(guān)系數(shù) PACF 與自相關(guān)系數(shù)圖格式相同與自相關(guān)系數(shù)圖格式相同逆自相關(guān)系數(shù)逆自相關(guān)系數(shù) SIACF 在在ARIMA建模中，樣本遞自相關(guān)系數(shù)和建模中，樣本遞自相關(guān)系數(shù)和PACF起起大體一樣的作用，但大體一樣的作用，但SIACF在指出子集和季節(jié)自在指出子集和季節(jié)自回歸模型時(shí)效果優(yōu)于回歸模型時(shí)效果優(yōu)于PACF。對(duì)偶模型的自相關(guān)系數(shù)稱(chēng)作原模型的逆自相關(guān)系對(duì)偶模型的自相關(guān)系數(shù)稱(chēng)作原模型的逆自相關(guān)系數(shù)。數(shù)。BtBtWaBtBtZa白噪聲檢驗(yàn)白噪聲檢驗(yàn)卡方檢驗(yàn)卡方檢驗(yàn) H0 ：直到某一給定時(shí)間間隔的樣本自相關(guān)系數(shù)沒(méi)：直到某一給定時(shí)間間隔的樣本自相關(guān)系數(shù)沒(méi)有顯著不為零的有顯著不為零的.(Xt為白噪聲，獨(dú)立

18、的隨機(jī)擾動(dòng)）為白噪聲，獨(dú)立的隨機(jī)擾動(dòng)） 如果對(duì)所有時(shí)間間隔，該零假設(shè)成立，則沒(méi)有需要如果對(duì)所有時(shí)間間隔，該零假設(shè)成立，則沒(méi)有需要建模的信息，也不需要建立建模的信息，也不需要建立ARIMA模型模型. 被檢查的時(shí)間間隔個(gè)數(shù)依賴(lài)于被檢查的時(shí)間間隔個(gè)數(shù)依賴(lài)于=選項(xiàng)選項(xiàng) 對(duì)前對(duì)前N-2個(gè)自相關(guān)系數(shù)的檢驗(yàn)個(gè)自相關(guān)系數(shù)的檢驗(yàn)P值。值。 P=0.005 拒絕拒絕 H0 （拒絕為白噪聲，（拒絕為白噪聲，P=0時(shí)，時(shí)， Xt高度自相關(guān)）高度自相關(guān)） P=0.005 接受接受 H0 (即對(duì)所有時(shí)間間隔，自相關(guān)即對(duì)所有時(shí)間間隔，自相關(guān)系數(shù)為零，說(shuō)明沒(méi)有建模信息，不必要做下去了）系數(shù)為零，說(shuō)明沒(méi)有建模信息，不必要做下去

19、了）第二階段：估計(jì)和診斷檢驗(yàn)階段第二階段：估計(jì)和診斷檢驗(yàn)階段 時(shí)間序列時(shí)間序列Wt由由IDENTIFY語(yǔ)句識(shí)別并且由語(yǔ)句識(shí)別并且由ESTIMATE語(yǔ)句處理語(yǔ)句處理, 即在完成可能的模型即在完成可能的模型識(shí)別后，開(kāi)始估計(jì)和診斷檢查階段識(shí)別后，開(kāi)始估計(jì)和診斷檢查階段. 估計(jì)模型估計(jì)模型AR(P)estimate p=?; (根據(jù)偏相關(guān)系數(shù)的截尾點(diǎn)根據(jù)偏相關(guān)系數(shù)的截尾點(diǎn))run;ESTIMATE語(yǔ)句的功能語(yǔ)句的功能是用模型擬合數(shù)據(jù)，并打印是用模型擬合數(shù)據(jù)，并打印出參數(shù)估計(jì)值和診斷統(tǒng)計(jì)量，指出模型對(duì)數(shù)據(jù)的出參數(shù)估計(jì)值和診斷統(tǒng)計(jì)量，指出模型對(duì)數(shù)據(jù)的擬合優(yōu)度。擬合優(yōu)度。ESTIMATE的輸出的輸出 參數(shù)估

20、計(jì)表參數(shù)估計(jì)表：估計(jì)方法：估計(jì)方法：METHOD=選擇不同的估計(jì)方法（選擇不同的估計(jì)方法（條件最條件最小二乘法估計(jì)小二乘法估計(jì)，極大似然估計(jì)，條件、無(wú)條件估計(jì)，極大似然估計(jì)，條件、無(wú)條件估計(jì)，線性或非線性估計(jì)）線性或非線性估計(jì)）參數(shù)估計(jì)值參數(shù)估計(jì)值(提供：估計(jì)值，標(biāo)準(zhǔn)差，提供：估計(jì)值，標(biāo)準(zhǔn)差，t比值比值)t比值比值：關(guān)于參數(shù)估計(jì)值的顯著性檢驗(yàn)（近似值）。：關(guān)于參數(shù)估計(jì)值的顯著性檢驗(yàn)（近似值）。當(dāng)觀測(cè)序列的長(zhǎng)度很短，并且被估計(jì)參數(shù)的個(gè)數(shù)相對(duì)當(dāng)觀測(cè)序列的長(zhǎng)度很短，并且被估計(jì)參數(shù)的個(gè)數(shù)相對(duì)于序列昌都很長(zhǎng)時(shí)，于序列昌都很長(zhǎng)時(shí)，t統(tǒng)計(jì)量的近似效果統(tǒng)計(jì)量的近似效果很差。很差。均值項(xiàng)均值項(xiàng)MU常數(shù)項(xiàng)常數(shù)項(xiàng)：

21、Constant Estimate 模型的常數(shù)項(xiàng)可以表模型的常數(shù)項(xiàng)可以表示為均值項(xiàng)示為均值項(xiàng)MU和自回歸參數(shù)的函數(shù)。和自回歸參數(shù)的函數(shù)。擬合優(yōu)度統(tǒng)計(jì)量表擬合優(yōu)度統(tǒng)計(jì)量表Variance Estimate 殘差序列的方差殘差序列的方差Std Error Estimate 方差估計(jì)值的平方根方差估計(jì)值的平方根AIC和和BIC兩個(gè)信息準(zhǔn)則兩個(gè)信息準(zhǔn)則 （貝葉斯準(zhǔn)則）（貝葉斯準(zhǔn)則） 其中其中L是一個(gè)近似值，以是一個(gè)近似值，以AIC的絕對(duì)值越小擬合優(yōu)度越好的絕對(duì)值越小擬合優(yōu)度越好.參數(shù)估計(jì)值的相關(guān)系數(shù)表參數(shù)估計(jì)值的相關(guān)系數(shù)表 判斷其變量之間互線性可能影響結(jié)果的程度判斷其變量之間互線性可能影響結(jié)果的程度.

22、 如果兩個(gè)參數(shù)估計(jì)值高度相關(guān)如果兩個(gè)參數(shù)估計(jì)值高度相關(guān)(相關(guān)系數(shù)很高相關(guān)系數(shù)很高), 可以可以考慮模型從模型中去掉一個(gè)參數(shù)對(duì)應(yīng)的變量?？紤]模型從模型中去掉一個(gè)參數(shù)對(duì)應(yīng)的變量。AIC=-2ln(L)+2RBIC(SBC)=-2ln(L)+ln(n)R殘差自相關(guān)系數(shù)的檢驗(yàn)殘差自相關(guān)系數(shù)的檢驗(yàn) 對(duì)于殘差序列的卡方檢驗(yàn)統(tǒng)計(jì)量將指對(duì)于殘差序列的卡方檢驗(yàn)統(tǒng)計(jì)量將指明殘差是否不相關(guān)，或是否包含可以被更明殘差是否不相關(guān)，或是否包含可以被更復(fù)雜模型利用的附加信息。復(fù)雜模型利用的附加信息。H0：無(wú)自相關(guān)：無(wú)自相關(guān) P值值=1jan58d;plot x*date=1 forecast*date=2 l95*date

23、=3 u95*date=3 / overlay haxis=1jan58d to 1jan62d by year; run;例例2 國(guó)航客票數(shù)國(guó)航客票數(shù)(Airline)時(shí)間序列分析時(shí)間序列分析 Data 數(shù)據(jù)集名數(shù)據(jù)集名; input x ; date=intnx(month, datamonthyear(開(kāi)始時(shí)刻開(kāi)始時(shí)刻),_n_-1) /注：生成周期日期值，從第一個(gè)觀測(cè)日期起所需的增量數(shù)注：生成周期日期值，從第一個(gè)觀測(cè)日期起所需的增量數(shù) cards; /輸入數(shù)據(jù)，輸入數(shù)據(jù)，建立數(shù)據(jù)集建立數(shù)據(jù)集(至少五年觀測(cè)值）。至少五年觀測(cè)值）。 Run； Proc x11 data=數(shù)據(jù)集名數(shù)據(jù)集名；

24、 monthly date=date; var=變量名變量名； tables d11(最終的季節(jié)調(diào)整過(guò)的序列）最終的季節(jié)調(diào)整過(guò)的序列） /注：注：D部表格給出了各成分的終估計(jì)，表部表格給出了各成分的終估計(jì)，表d10為季節(jié)因子終估為季節(jié)因子終估計(jì)，表計(jì)，表d12為趨勢(shì)起伏的終估計(jì)，表為趨勢(shì)起伏的終估計(jì)，表d13為不規(guī)則項(xiàng)終估計(jì)，為不規(guī)則項(xiàng)終估計(jì)， run; / SAS軟件軟件X11過(guò)程的實(shí)現(xiàn)過(guò)程的實(shí)現(xiàn)先驗(yàn)星期權(quán)重與交易日回歸先驗(yàn)星期權(quán)重與交易日回歸 若序列分離出長(zhǎng)期趨勢(shì)以后若序列分離出長(zhǎng)期趨勢(shì)以后, 若存在星期構(gòu)若存在星期構(gòu)成因素的影響成因素的影響,需要對(duì)序列中進(jìn)行需要對(duì)序列中進(jìn)行先驗(yàn)星期權(quán)重先驗(yàn)星期權(quán)重與交易日回歸分析與交易日回歸分析。歷法效應(yīng)：歷法效應(yīng)：假定歷法效應(yīng)既有先驗(yàn)周期權(quán)重也有假定歷法效應(yīng)既有先驗(yàn)周期權(quán)重也有其他影響，需要同時(shí)使用以下語(yǔ)句：其他影響，需要同時(shí)使用以下語(yǔ)句：Proc x11 data=數(shù)據(jù)集名數(shù)據(jù)集名； monthly date=date tdregr=adjust; var=變量名變量名； tables a1 a4 b15 b16 c14 c15 c18 d11; pdweights mon=1 tue=1 wed=1 thu=2 fri=2 sat=2;(sun=0) ou