SAS經(jīng)濟(jì)時(shí)間序列分-各種模型分析

1、PAGE 1PAGE 61 目 錄實(shí)驗(yàn)一 分析太陽黑子數(shù)序列3實(shí)驗(yàn)二 模擬AR模型4實(shí)驗(yàn)三 模擬MAA模型和和ARMMA模型型6實(shí)驗(yàn)四 分析化化工生產(chǎn)產(chǎn)量數(shù)據(jù)據(jù)8實(shí)驗(yàn)五 模擬ARRIMAA模型和和季節(jié)AARIMMA模型型10實(shí)驗(yàn)六 分析美美國國民民生產(chǎn)總總值的季季度數(shù)據(jù)據(jù)113實(shí)驗(yàn)七 分析國國際航線線月度旅旅客總數(shù)數(shù)數(shù)據(jù)16實(shí)驗(yàn)八 干預(yù)模模型的建建模199實(shí)驗(yàn)九 傳遞函函數(shù)模型型的建模模222實(shí)驗(yàn)十 回歸與與時(shí)序相相結(jié)合的的建模25太陽黑子子年度數(shù)數(shù)據(jù)28美國國民民收入數(shù)數(shù)據(jù)29化工生產(chǎn)產(chǎn)過程的的產(chǎn)量數(shù)數(shù)據(jù)30國際航線線月度旅旅客數(shù)據(jù)據(jù)300洛杉磯臭臭氧每小小時(shí)讀數(shù)數(shù)的月平平均值數(shù)數(shù)據(jù)31煤氣

2、爐數(shù)數(shù)據(jù)35芝加哥某某食品公公司大眾眾食品周周銷售數(shù)數(shù)據(jù)37牙膏市場場占有率率周數(shù)據(jù)據(jù)399某公司汽汽車生產(chǎn)產(chǎn)數(shù)據(jù)444加拿大山山貓數(shù)據(jù)據(jù)444 實(shí)驗(yàn)一一 分析析太陽黑黑子數(shù)序序列一、實(shí)驗(yàn)驗(yàn)?zāi)康模毫私鈺r(shí)時(shí)間序列列分析的的基本步步驟，熟熟悉SAAS/EETS軟軟件使用用方法。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：分析太太陽黑子子數(shù)序列列。三、實(shí)驗(yàn)驗(yàn)要求：了解時(shí)時(shí)間序列列分析的的基本步步驟，注注意各種種語句的的輸出結(jié)結(jié)果。四、實(shí)驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。創(chuàng)建名為為expp1的SASS數(shù)據(jù)集集，即在在窗中輸輸入下列列語句：保存此步步驟中的的程序，供以后后

3、分析使使用（只只需按工工具條上上的保存存按鈕然然后填寫寫完提問問后就可可以把這這段程序序保存下下來即可可）。繪數(shù)據(jù)與與時(shí)間的的關(guān)系圖圖，初步步識別序序列，輸輸入下列列程序：ods htmml;ods lisstinng cclosse;run;提交程程序，在在graaph窗窗口中觀觀察序列列，可以以看出此此序列是是均值平平穩(wěn)序列列。識別模型型，輸入入如下程程序。提交程序序，觀察察輸出結(jié)結(jié)果。初初步識別別序列為為AR(2)模型。估計(jì)和診診斷。輸輸入如下下程序：提交程序序，觀察察輸出結(jié)結(jié)果。假假設(shè)通過過了白噪噪聲檢驗(yàn)驗(yàn)，且模模型合理理，則進(jìn)進(jìn)行預(yù)測測。進(jìn)行預(yù)測測，輸入入如下程程序：提交程序序，觀察

4、察輸出結(jié)結(jié)果。退出SAAS系統(tǒng)統(tǒng)，關(guān)閉閉計(jì)算機(jī)機(jī)?？偝绦颍篸ataa exxp1;infiile D:exxp1.txtt;inpuut aa1 ;yearr=inntnxx(yyearr,1jaan17742d,_n_-1);formmat yeaar yyearr4.;procc prrintt;ruun;ods htmml;ods lisstinng cclosse;procc gpplott daata=expp1 ; syymbool ii=spplinne vv=doot hh=1 cv=redd cii=grreenn w=1; pllot a1*yeaar/aautoovree

5、f llvreef=22 cfframme=yyelllow cvrref=blaack ; tiitlee 太太陽黑子子數(shù)序列列;run; prooc aarimma ddataa=exxp1; iddenttifyy vaar=aa1 nnlagg=244 miinicc p=(0:5) q=(0:55); esstimmatee p=3; fooreccastt leead=6 iinteervaal=yyearr idd=yeear outt=ouut;run;procc prrintt daata=outt;run;選取擬合合模型的的規(guī)則:1.模型型顯著有有效(殘殘差檢驗(yàn)驗(yàn)為白噪噪聲

6、)2.模型型參數(shù)盡盡可能少少3.結(jié)合合自相關(guān)關(guān)圖和偏偏自相關(guān)關(guān)圖以及及minnic條條件(BBIC信信息量最最小原則則),選選取顯著著有效的的參數(shù) 實(shí)驗(yàn)二二 模擬ARR模型實(shí)驗(yàn)?zāi)康牡模菏煜はじ鞣NAAR模型型的樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)的的特點(diǎn)，為理 論學(xué)學(xué)習(xí)提供供直觀的的印象。實(shí)驗(yàn)內(nèi)容容：隨機(jī)機(jī)模擬各各種ARR模型。實(shí)驗(yàn)要求求：記錄錄各ARR模型的的樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，觀觀察各種種序列 圖形形，總結(jié)結(jié)AR模型型的樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)的的特點(diǎn)實(shí)驗(yàn)時(shí)間間：2小小時(shí)。實(shí)驗(yàn)軟件件：SAAS系統(tǒng)統(tǒng)。實(shí)驗(yàn)步驟驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。模擬實(shí)根根情況，模擬過過程。

7、在ediit窗中中輸入如如下程序序： 4、觀察輸輸出的數(shù)數(shù)據(jù)，輸輸入如下下程序，并提交交程序。觀察樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)，輸入輸輸入如下下程序，并提交交程序。作為作作業(yè)把樣樣本自相相關(guān)系數(shù)數(shù)和偏相相關(guān)系數(shù)數(shù)記錄下下來。估計(jì)模型型參數(shù)，并與實(shí)實(shí)際模型型的系數(shù)數(shù)進(jìn)行對對比，即即輸入如如下程序序，并提提交。模擬虛根根情況，模擬過過程。重重復(fù)步驟驟3-7即可可（但部部分程序序需要修修改，請請讀者自自己完成成）。模擬ARR(3)模型，模擬過過程。重重復(fù)步驟驟3-7即可可(但部分分程序需需要修改改，請讀讀者自己己完成).10、回回到grraphh窗口觀觀察各種種序列圖圖形的異異同11、退退出S

8、AAS系統(tǒng)統(tǒng),關(guān)閉計(jì)計(jì)算機(jī).總程序:titlle;dataa a; x1=00.5; x2=00.5; do ii=-550 tto 2250; a=raannoor(3325665); x=a-0.66*x11+0.4*xx2; x2=xx1; x1=xx; outpput; end; run;procc prrintt daata=a; var x; procc gpplott daata=a;symbbol i=sspliine c=rred; plott x*i/hhaxiis=-50 to 2555 byy 200; run;quitt;procc arrimaa daata=a;

9、idenntiffy vvar=x nnlagg=100 miinicc p=(0:3) q=(0:33) ooutccov=expp1; estiimatte pp=2 noiint;run; procc gpplott daata=expp1;symbbol i=nneeddle widdth=6;plott coorr*lagg;run;procc gpplott daata=expp1;symbbol i=nneeddle widdth=6;plott paartccorrr*laag;run; 實(shí)驗(yàn)三三 模擬擬MA模型型和ARRMA模模型實(shí)驗(yàn)?zāi)康牡模菏煜はじ鞣NMMA模型型和ARRMA

10、模模型的樣樣本自相相關(guān)系數(shù)數(shù)和偏相相關(guān)系數(shù)數(shù) 的特特點(diǎn)，為為理論學(xué)學(xué)習(xí)提供供直觀的的印象。實(shí)驗(yàn)內(nèi)容容：隨機(jī)機(jī)模擬各各種MAA模型和和ARMMA模型型。實(shí)驗(yàn)要求求：記錄錄各MAA模型和和ARMMA模型型的樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)， 觀察察各序列列的異同同，總結(jié)結(jié)MA模型型和ARRMA模模型的樣樣本自相相關(guān)系 數(shù)和和偏相關(guān)關(guān)系數(shù)的的特點(diǎn)實(shí)驗(yàn)時(shí)間間：2小小時(shí)。實(shí)驗(yàn)軟件件：SAAS系統(tǒng)統(tǒng)。實(shí)驗(yàn)步驟驟開機(jī)進(jìn)入入SASS系統(tǒng)。2、模擬擬情況，模擬過過程。在ediit窗中中輸入如如下程序序：4、觀察察輸出的的數(shù)據(jù)序序列，輸輸入如下下程序，并提交交程序。觀察樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)，輸入輸

11、輸入如下下程序，并提交交程序。估計(jì)模型型參數(shù)，并與實(shí)實(shí)際模型型的系數(shù)數(shù)進(jìn)行對對比，即即輸入如如下程序序，并提提交。模擬情況況，模擬擬過程。重復(fù)步步驟3-7即可可（但部部分程序序需要修修改，請請讀者自自己完成成）。模擬情況況，模擬擬過程。重復(fù)步步驟3-7即可可（但部部分程序序需要修修改，請請讀者自自己完成成）。模擬情況況，模擬擬過程。重復(fù)步步驟3-7即可可（但部部分程序序需要修修改，請請讀者自自己完成成）。模擬ARRMA模模型，模模擬過程程。重復(fù)復(fù)步驟33-7即可可(但部分分程序需需要修改改，請讀讀者自己己完成).回到grraphh窗口觀觀察各種種序列圖圖形的異異同。退出SAAS系統(tǒng)統(tǒng),關(guān)閉計(jì)計(jì)

12、算機(jī).總程序:dataa a; a1=00; a2=00; do nn=1tto 2250; a=raannoor(3325665); x=a+0.665*aa1+00.244*a22; a2=aa1; a1=aa; outpput; end; run;procc gpplott daata=a;symbbol i=sspliine h=11 w=1;plott x*n /haxxis=-100 too 2660 bby 110;run;procc arrimaa daata=a; idenntiffy vvar=x nnlagg=100 miinicc p=(0:3) q=(0:33) oo

13、utccov=expp1; estiimatte qq=2 noiint;run; procc gpplott daata=expp1;symbbol11 i=neeedlee c=redd;plott coorr*lagg=1;run;procc gpplott daata=expp1;symbbol22 i=neeedlee c=greeen;plott paartccorrr*laag=22;run;quitt; 實(shí)驗(yàn)四四 分析析化工生生產(chǎn)量數(shù)數(shù)據(jù)實(shí)驗(yàn)?zāi)康牡模哼M(jìn)一一步熟悉悉時(shí)間序序列建模模的基本本步驟，掌握用用SACCF及SPAACF定定 模型型的階的的方法。實(shí)驗(yàn)內(nèi)容容：分析析化工生生產(chǎn)

14、過程程的產(chǎn)量量序列。實(shí)驗(yàn)要求求：掌握握ARMMA模型型建模的的基本步步驟，初初步掌握握數(shù)據(jù)分分析技巧巧。寫出出 實(shí)驗(yàn)驗(yàn)報(bào)告。實(shí)驗(yàn)時(shí)間間：2小小時(shí)。實(shí)驗(yàn)軟件件：SAAS系統(tǒng)統(tǒng)。實(shí)驗(yàn)步驟驟開機(jī)進(jìn)入入SASS系統(tǒng)。創(chuàng)建名為為expp2的SAAS數(shù)據(jù)據(jù)集，即即在窗中中輸入下下列語句句：保存此步步驟中的的程序，供以后后分析使使用（只只需按工工具條上上的保存存按鈕然然后填寫寫完提問問后就可可以把這這段程序序保存下下來即可可）。繪數(shù)據(jù)與與時(shí)間的的關(guān)系圖圖，初步步識別序序列，輸輸入下列列程序：提交程序序，在ggrapph窗口口中觀察察序列，可以看看出此序序列是均均值平穩(wěn)穩(wěn)序列。識別模型型，輸入入如下程程序。提

15、交程序序，觀察察輸出結(jié)結(jié)果，發(fā)發(fā)現(xiàn)二階階樣本自自相關(guān)系系數(shù)和一一階的樣樣本偏相相關(guān)系數(shù)數(shù)都在22倍的標(biāo)標(biāo)準(zhǔn)差之之外，那那么我們們首先作作為一階階AR模型型估計(jì)，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，發(fā)發(fā)現(xiàn)殘差差能通過過白噪聲聲檢驗(yàn)，但它的的二階的的樣本偏偏相關(guān)系系數(shù)比較較大，那那么我們們考慮二二階ARR模型。輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，發(fā)發(fā)現(xiàn)殘差差樣本自自相關(guān)系系數(shù)和樣樣本偏相相關(guān)系數(shù)數(shù)都 在在2倍的的標(biāo)準(zhǔn)差差之內(nèi)。且能通通過白噪噪聲檢驗(yàn)驗(yàn)。比較較兩個(gè)模模型的AAIC和和SBCC, 發(fā)發(fā)現(xiàn)第二二個(gè)模型型的AIIC和SBCC都比第第一個(gè)的的小，故故我們選選擇第二二個(gè)

16、模型型為 我我們的結(jié)結(jié)果。記錄參數(shù)數(shù)估計(jì)值值，寫出出模型方方程式。進(jìn)行預(yù)測測，輸入入如下程程序：提交程序序，觀察察輸出結(jié)結(jié)果。退出SAAS系統(tǒng)統(tǒng)，關(guān)閉閉計(jì)算機(jī)機(jī)。dataa exxp2; innfille D:expp1.ttxt; innputt x ; n=_n_;procc prrintt;ruun;procc gpplott daata=expp2; syymbool ii=jooin v=sstarr h=2 cci=ggreeen ccv=rred; pllot x*nn/vrref=30 50 70 cvrref=redd lvvreff=2 ;run;procc arrimaa

17、 daata=expp2;idenntiffy vvar=x nnlagg=122 miinicc p=(0:3) q=(0:33);estiimatte pplott p=1;foreecasst lleadd=2 outt=ouut;run;quitt; 實(shí)驗(yàn)五五 模擬擬ARIIMA模模型和季季節(jié)ARRIMAA模型實(shí)驗(yàn)?zāi)康牡模菏煜はじ鞣NAARIMMA模型型的樣本本自相關(guān)關(guān)系數(shù)和和偏相關(guān)關(guān)系數(shù)的的特點(diǎn)， 區(qū)別別各種AARIMMA模型型的圖形形，為理理論學(xué)習(xí)習(xí)提供直直觀的印印象。實(shí)驗(yàn)內(nèi)容容：隨機(jī)機(jī)模擬各各種ARRIMAA模型。實(shí)驗(yàn)要求求：記錄錄各ARRIMAA模型的的樣本自自相關(guān)系系數(shù)和偏偏相

18、關(guān)系系數(shù)觀察察各序列列 圖形形的異同同，總結(jié)結(jié)ARIIMA模模型的樣樣本自相相關(guān)系數(shù)數(shù)和偏相相關(guān)系數(shù)數(shù) 的特特點(diǎn)實(shí)驗(yàn)時(shí)間間：2小小時(shí)。實(shí)驗(yàn)軟件件：SAAS系統(tǒng)統(tǒng)。實(shí)驗(yàn)步驟驟開機(jī)進(jìn)入入SASS系統(tǒng)。2、模擬擬ARIIMA(0,11,1)過程，模擬過過程。創(chuàng)建數(shù)據(jù)據(jù)集，在在ediit窗中中輸入如如下程序序： 4、觀察輸輸出的數(shù)數(shù)據(jù)序列列，輸入入如下程程序：。 5、提交程程序，在在Graaph窗窗口中觀觀察圖形形。6、觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，輸輸入輸入入如下程程序： 提交程程序，發(fā)發(fā)現(xiàn)自相相關(guān)系數(shù)數(shù)成緩慢慢下降的的趨勢，說明要要做差分分運(yùn)算，做一階階差分運(yùn)運(yùn)算，輸輸入如下下程序：提

19、交程序序，觀察察樣本自自相關(guān)系系數(shù)與樣樣本偏相相關(guān)系數(shù)數(shù)，發(fā)現(xiàn)現(xiàn)自相關(guān)關(guān)系數(shù)11階截尾尾，故判判斷差分分后序列列為MAA(1)模型。進(jìn)行模模型參數(shù)數(shù)估計(jì)，輸入如如下程序序：提交程序序，并觀觀察殘差差圖，發(fā)發(fā)現(xiàn)模型型擬合完完全。10、寫寫出模型型的方程程，并與與真實(shí)模模型對比比。titlle;dataa a; x1=00.9; a1=00; do nn=0 to 2500; a=raannoor(3325665); x=x11+a-0.88*a11; x1=xx; a1=aa; outpput; end; run;procc gpplott daata=a; syymbool ii=jooin

20、v=ddot h=11 cii=grreenn cvv=reed; pllot x*nn/vrref=-2 1 44 cvvreff=reed llvreef=22 haaxiss=-110 tto 2260 by 10;run;procc arrimaa daata=a; idenntiffy vvar=x nnlagg=100 miinicc p=(0:3) q=(0:33) ooutccov=expp1; run; procc gpplott daata=expp1;symbbol11 i=neeedlee c=redd;plott coorr*lagg=1;run;procc gppl

21、ott daata=expp1;symbbol22 i=neeedlee c=greeen;plott paartccorrr*laag=22;run;procc arrimaa daata=a; iddenttifyy vaar=xx(1) nllag=24 minnic p=(0:33) qq=(00:3); /*一一階差分分x(11)*/run;estiimatte qq=1 ploot nnoinnt;run;quitt;11、模模擬ARRIMAA(1,1,00)模型型，模擬擬過程。重復(fù)步步驟 3-10即可可（但部部分程序序需要修修改，請請讀者自自己完成成）。模擬模型型， 模擬擬模型,

22、 即模型。13、創(chuàng)創(chuàng)建數(shù)據(jù)據(jù)集，在在ediit窗中中輸入如如下程序序：繪序列圖圖，輸入入如下程程序：提交程序序，到ggrapph窗口口中觀察察序列圖圖形。初步識別別模型，輸入如如下程序序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)和樣樣本偏相相關(guān)系數(shù)數(shù)。做季節(jié)差差分和一一階差分分除掉季季節(jié)因子子和趨勢勢因子，輸入如如下程序序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)和樣樣本偏相相關(guān)系數(shù)數(shù)，確定定模型階階數(shù)。估計(jì)模型型參數(shù)，輸入如如下程序序：提交程序序，觀察察殘差的的樣本自自相關(guān)系系數(shù)和樣樣本偏相相關(guān)系數(shù)數(shù)，看是是否通過過 了白白噪聲檢檢驗(yàn)。寫寫出模型型方程式式，并與與真實(shí)模模型對比比?；氐絞rraphh

23、窗口觀觀察各種種序列圖圖形的異異同。退出SAAS系統(tǒng)統(tǒng),關(guān)閉計(jì)計(jì)算機(jī). daata c; x11=0.9;xx2=00;x33=0;x4=0;xx5=00;x66=0;x7=0; x88=0;x9=0;xx10=0;xx11=0;xx12=0;xx13=0; a11=0;a2=0;aa3=00;a44=0;a5=0;aa6=00;a77=0; a88=0;a9=0;aa10=0;aa11=0;aa12=0;aa13=0; doo n=0 tto 2250; a=rrannnor(123345); x=xx1+xx12-x133+a-0.44*a11-0.6*aa12+0.224*aa13; x

24、133=x112;xx12=x111;x111=xx10;x100=x99;x99=x88;x88=x77; x7=x6;x6=x5;x5=x4;x4=x3;x3=x2;x2=x1;x1=x; a133=a112;aa12=a111;a111=aa10;a100=a99;a99=a88;a88=a77; a7=a6;a6=a5;a5=a4;a4=a3;a3=a2;a2=a1;a1=a; ouutpuut; ennd; ruun;procc gpplott daata=c; syymbool ii=jooin v=ddot h=11 cii=grreenn cvv=reed; pllot x*n

25、n/vrref=-200 1 10 cvrref=redd lvvreff=2 haxxis=-100 too 2660 bby 110;run;procc arrimaa daata=c; idenntiffy vvar=x nnlagg=200 miinicc p=(0:3) q=(0:33); run; idenntiffy vvar=x(11,122) nnlagg=366 miinicc p=(0:3) q=(0:33) outtcovv=exxp1;run;estiimatte qq=(11)(112) metthodd=clls nnoinnt;run;procc gpplott

26、 daata=expp1;symbbol11 i=neeedlee c=redd;plott coorr*lagg=1;run;procc gpplott daata=expp1;symbbol22 i=neeedlee c=greeen;plott paartccorrr*laag=22;run;quitt; 實(shí)驗(yàn)六六 分析析美國國國民生產(chǎn)產(chǎn)總值的的季度數(shù)數(shù)據(jù)一、實(shí)驗(yàn)驗(yàn)?zāi)康模哼M(jìn)一步步學(xué)習(xí)數(shù)數(shù)據(jù)分析析技巧，進(jìn)一步步了解AARIMMA模型型。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：47年年1季度度到966年3季季度美國國國民生生產(chǎn)總值值的季度度數(shù)據(jù)。三、實(shí)驗(yàn)驗(yàn)要求：寫出分分析報(bào)告告。四、實(shí)驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)

27、軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。2、建立立名為eexp33的SASS數(shù)據(jù)集集，輸入入如下程程序：保存上述述程序，供以后后分析使使用（只只需按工工具條上上的保存存按鈕，然后填填寫 完完提問后后就可以以把這段段程序保保存下來來）。繪序列圖圖，輸入入如下程程序：觀察圖形形，發(fā)現(xiàn)現(xiàn)圖形成成指數(shù)函函數(shù)上升升形式，故做對對數(shù)變換換，輸入入如下程程序：繪變換后后序列圖圖，輸入入如下程程序：提交程序序，到ggrapph窗口口中觀察察變換后后的序列列圖，可可以看出出它成直直線上升升趨勢。對序列列做初步步識別，輸入如如下程序序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)，可可看出有有緩慢

28、下下降趨勢勢，結(jié)合合我們觀觀察的圖圖形，我我們知道道要對序序列做差差分運(yùn)算算，作一一階差分分，輸入入如下程程序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)，可可看出樣樣本自相相關(guān)系數(shù)數(shù)5步后后是截尾尾的,那那么確定定為MAA(5)模型,進(jìn)行參參數(shù)估計(jì)計(jì),輸入入如下程程序:提交程序序，觀察察輸出結(jié)結(jié)果，可可看出模模型通過過了白噪噪聲檢驗(yàn)驗(yàn)，說明明模型擬擬合充分分。且MMA1,3 , MAA1,44的T值較小小，說明明參數(shù)顯顯著為00，除掉掉這兩項(xiàng)項(xiàng)重新進(jìn)進(jìn)行估計(jì)計(jì)，輸入入如下程程序：提交程序序，觀察察輸出結(jié)結(jié)果，可可看出模模型通過過了白噪噪聲檢驗(yàn)驗(yàn)，說明明模型擬擬合充分分，且殘殘差標(biāo)準(zhǔn)準(zhǔn)誤與前前一估計(jì)

29、計(jì)相差很很小，故故以此結(jié)結(jié)果為我我們所要要的結(jié)果果，依此此結(jié)果寫寫出方程程式。進(jìn)行預(yù)測測，預(yù)測測美國未未來2年年的每季季國民生生產(chǎn)總值值。輸入入如下程程序：提交程序序，并把把預(yù)測值值記錄下下來。退出SAAS系統(tǒng)統(tǒng)，關(guān)閉閉計(jì)算機(jī)機(jī)。dataa exxp3; innfille C:Doccumeentss annd SSetttinggsAAdmiinisstraatorr桌面面exxp1.txtt; innputt gnnp; daate=inttnx(qttr,1jjan447dd,_nn_-11); foormaat ddatee yyyqc.;run;procc gpplott daata

30、=expp3; syymbool ii=jooin w=2 cci=ggreeen; /*ww：線的的大小 h：點(diǎn)點(diǎn)的大小小 l：線的類類型*/ pllot gnpp*daate=1;run;dataa leexp; sett exxp3; lgnnp=llog(gnpp);run;procc gpplott daata=lexxp; syymbool2 i=sspliine c=rred w=22; pllot lgnnp*ddatee=2 ;run;procc arrimaa daata=lexxp; iddenttifyy vaar=llgnpp(1) nllag=12 minnic p

31、=(0:33) qq=(00:3);run;estiimatte qq=5;run;estiimatte qq=(11,2,5) ploot;run;foreecasst lleadd=6 inttervval=qtrr idd=daate outt=reesullts;run;dataa reesullts22;set ressultts;gnp=expp(lggnp);l95=expp(l995);u95=expp(u995);foreecasst=eexp(forrecaast);keepp daate gnpp l995 fforeecasst uu95;run;procc prrin

32、tt daata=ressultts2;var datte ggnp l955 fooreccastt u995;wherre ddatee=1jaan966d;run;procc gpplott daata=ressultts2;plott gnnp*ddatee=1 forrecaast*datte=22 l995*ddatee=3 u955*daate=3/ooverrlayy leegennd;symbbol11 v=dott cvv=reed ii=noone h=11 w=1;symbbol22 i=joiin cci=ggreeen;symbbol33 i=spllinee ci

33、i=bllackk l=2;run;quitt; 實(shí)驗(yàn)七七 分析析國際航航線月度度旅客總總數(shù)數(shù)據(jù)據(jù)實(shí)驗(yàn)?zāi)康牡模菏煜はみ\(yùn)用SSAS建建立模型型的方法法，進(jìn)一一步 了解解模型的的特征。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：19449年11月至119600年122月國際際航線月月度旅客客總數(shù)數(shù)數(shù)據(jù)。三、實(shí)驗(yàn)驗(yàn)要求：寫出分分析報(bào)告告。四、實(shí)驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。2、建立立名為eexp44的SAAS數(shù)據(jù)據(jù)集，輸輸入如下下程序：繪序列圖圖，輸入入如下程程序：提交程序序，觀察察圖形，發(fā)現(xiàn)圖圖形有很很強(qiáng)的季季節(jié)性，且成指指數(shù)函數(shù)數(shù)上升形形式，故故做對數(shù)數(shù)變換，

34、輸入如如下程序序：繪變換后后序列圖圖，輸入入如下程程序：提交程序序，到ggrapph窗口口中觀察察變換后后的序列列圖，可可以看出出它總的的趨勢成成直線上上升，且且有很強(qiáng)強(qiáng)的季節(jié)節(jié)性。對對序列做做初步識識別，輸輸入如下下程序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，可可看出樣樣本自相相關(guān)系數(shù)數(shù)有緩慢慢下降趨趨勢，偏偏相關(guān)系系數(shù)在11步，113步，25步步較大，我們作作一步一一階差分分，輸入入如下程程序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，發(fā)發(fā)現(xiàn)樣本本自相關(guān)關(guān)系數(shù)在在12步步，244步，336步特特別大,而偏相相關(guān)系數(shù)數(shù)在122步特別別大，那那么我們們再做112步的

35、的一階差差分，輸輸入如下下程序：10、提提交程序序，觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，發(fā)發(fā)現(xiàn)樣本本自相關(guān)關(guān)系數(shù)在在1步，112步特特別大,而偏相相關(guān)系數(shù)數(shù)看不出出有特別別的規(guī)律律，我們們可確定定模型的MA因因子為。11、進(jìn)進(jìn)行參數(shù)數(shù)估計(jì)，輸入如如下程序序：、提交程程序，觀觀察輸出出結(jié)果，可看出出模型通通過了白白噪聲檢檢驗(yàn)，說說明模型型擬 合充分分，故以以此結(jié)果果為我們們所要的的結(jié)果，依此結(jié)結(jié)果寫出出方程式式。13、進(jìn)進(jìn)行預(yù)測測，輸入入如下程程序： 144、提交交程序，仔細(xì)觀觀察預(yù)測測的結(jié)果果有什么么規(guī)律，思考為為什么有有這樣的的規(guī)律？15、變變換預(yù)測測值，以以獲取原原度量下下的預(yù)測測值

36、，輸輸入如下下程序：繪預(yù)測和和置信限限的散點(diǎn)點(diǎn)圖，輸輸入如下下程序：提交程序序，觀察察圖形。退出SAAS系統(tǒng)統(tǒng)，關(guān)閉閉計(jì)算機(jī)機(jī)。dataa exxp4; innfille C:Doccumeentss annd SSetttinggsAAdmiinisstraatorr桌面面exxp1.txtt; innputt aiir; daate=inttnx(moonthh,1jaan499d,_n_-1); foormaat ddatee moonyyy.;run;procc gpplott daata=expp4;symbbol11 i=joiin vv=doot cc=reed;plott ai

37、ir*ddatee=1;run;dataa laair; sett exxp4; laiir=llog(airr);run;procc gpplott daata=laiir; syymbool2 i=sspliine c=ggreeen; pllot laiir*ddatee=2;run;procc arrimaa daata=laiir; iddenttifyy vaar=llairr nllag=36;run;idenntiffy vvar=laiir(11) nnlagg=244;run;idenntiffy vvar=laiir(11,122) nnlagg=244 miinicc;

38、run;estiimatte qq=(11)(112) nocconsstannt mmethhod=ulss pllot;run;foreecasst lleadd=3 inttervval=monnth id=datte oout=b;run;procc prrintt daata=b;run;dataa c; sett b; airr=exxp(llairr); forrecaast=expp(fooreccastt+sttd*sstd/2); l955=exxp(ll95); u955=exxp(uu95);run;procc prrintt daata=c;run;symbbol11

39、 I=nonne vv=sttar r=11 c=redd;symbbol22 I=joiin vv=pllus r=11 c=greeen;symbbol33 I=joiin vv=noone l=33 r=1 cc=bllue;procc gpplott daata=c;wherre ddataa=1jaan599d;plott aiir*ddatee=1 forrecaast*datte=22 l995*ddatee=3 u955*daate=3/ ooverrlayy haaxiss=11jann59d tto 1maar611d by monnth;run; 實(shí)驗(yàn)八八 干預(yù)預(yù)模型的的

40、建模實(shí)驗(yàn)?zāi)康牡模赫莆瘴崭深A(yù)模模型的分分析方法法，進(jìn)一一步熟悉悉ARIIMA過過程的使使用方法法。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：19555年11月至119722年122月洛杉杉磯月平平均臭氧氧數(shù)據(jù)。三、實(shí)驗(yàn)驗(yàn)要求：寫出實(shí)實(shí)驗(yàn)報(bào)告告，掌握握干預(yù)模模型的建建模方法法。四、實(shí)驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。2、建立立名為eexp55的SAAS數(shù)據(jù)據(jù)集，輸輸入如下下程序：或者輸入入如下程程序： datta eexp55; inpput ozzonee ; datte=iintnnx(monnth,1jaan555d,_n_-1);formmat datte

41、 mmonyyy.;montth=mmontth(ddatee);yearr=yeear(datte);x1=yyearr=119600;summmer=(5monnth19965);wintter=(yeear19665)-summmerr; carrds; 只輸入入 ozzonee 一欄欄的數(shù)據(jù)據(jù) ； runn;保存上述述程序，供以后后分析使使用（只只需按工工具條上上的保存存按鈕，然后填填寫 完完提問后后就可以以把這段段程序保保存下來來）。4、繪序序列圖，輸入如如下程序序：5、提交交程序，觀察圖圖形，發(fā)發(fā)現(xiàn)圖形形有很強(qiáng)強(qiáng)的季節(jié)節(jié)性和緩緩慢下降降的趨勢勢。初步識別別模型，輸入如如下程序序：提

42、交程序序，觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，可可看出樣樣本自相相關(guān)系數(shù)數(shù)在1步步，122步，224步，36步步都較大大，且具具有周期期性，偏偏相關(guān)系系數(shù)在11步最大大，我們們作季節(jié)節(jié)差分，輸入如如下程序序：提交程序序，觀察察樣本自自相關(guān)系系數(shù)和偏偏相關(guān)系系數(shù)，發(fā)發(fā)現(xiàn)樣本本自相關(guān)關(guān)系數(shù)在在 11步，112步較較大,而而偏相關(guān)關(guān)系數(shù)在在1步，12步步，244步都較較大，且且呈現(xiàn)拖拖尾 現(xiàn)現(xiàn)象，我我們可確確定模型型的MAA因子為為。9、進(jìn)行行參數(shù)估估計(jì)，輸輸入如下下程序：10、提提交程序序，觀察察輸出結(jié)結(jié)果，可可看出模模型不是是很干凈凈，且不不能通過過白噪聲聲檢驗(yàn)。我我們可以以做殘差差序列圖

43、圖，觀看看殘差的的特性，輸入如如下程序序：11、進(jìn)進(jìn)行預(yù)測測，輸入入如下程程序：12、提提交程序序，觀察察圖形，可看出出前面一一段時(shí)期期的殘差差比后面面的要大大。 13、我們考考察修建建高速公公路后，是否對對臭氧有有顯著性性影響，輸入如如下程序序：14、提提交程序序，觀察察輸出結(jié)結(jié)果，發(fā)發(fā)現(xiàn)模型型的標(biāo)準(zhǔn)準(zhǔn)差,AAIC,SBCC都變小小了很多多， 且x11的影響響顯著。思考為為什么要要對x11進(jìn)行季季節(jié)差分分？15、我我們再來來考察汽汽車裝上上尾氣過過濾器，是否對對臭氧有有顯著性性影響，輸入如如 下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，發(fā)發(fā)現(xiàn)模型型的標(biāo)準(zhǔn)準(zhǔn)差,AAIC,SBCC都變小小了，且且

44、模型型基本上上通過了了白噪聲聲檢驗(yàn)，并且xx1,ssummmer的的影響顯顯著,而而winnterr 的影影響不顯顯著。思思考為什什么不對對summmerr和 wiinteer進(jìn)行行差分？ 177、進(jìn)行行預(yù)測值值，輸入入如下程程序： 注：這樣的的預(yù)測是是x1,summmerr,wiinteer已知知的預(yù)測測。 18、提提交程序序，觀察察預(yù)測值值。 199、退出出SASS系統(tǒng)，關(guān)閉計(jì)計(jì)算機(jī)。dataa exxp5; innfille C:Doccumeentss annd SSetttinggsAAdmiinisstraatorr桌面面exxp1.txtt; innputt n ozoone x

45、1 summmerr wiinteer; daate=inttnx(moonthh,1jaan555d,_n_-1); foormaat ddatee moonyyy.;run;procc gpplott daata=expp5; symmboll1 ii=jooin v=ddot c=rred; ploot oozonne*ddatee=1;run;procc arrimaa daata=expp5;idenntiffy vvar=ozoone nlaag=336;run;idenntiffy vvar=ozoone(12) nllag=36 minnic;run;estiimatte qq

46、=(11)(66)(112) metthodd=clls pplott;run;foreecasst lleadd=366 innterrvall=moonthh idd=daate outt=b nopprinnt;run;procc gpplott daata=b;symbbol I=sspliine v=ddot c=rred;plott reesidduall*daate;run;procc arrimaa daata=expp5; iddenttifyy vaar=oozonne(112) croossccorrr=(xx1(112) noopriint; esstimmatee q

47、=(1)(6)(122) iinpuut=(x1) nooconnstaant metthodd=mll ittpriint ploot;run;procc arrimaa daata=expp5; iddenttifyy vaar=oozonne(112) croossccorrr=(xx1(112) summmerr wiinteer) noopriint; esstimmatee q=(1)(122) iinpuut=(x1 summmerr wiinteer) nocconsstannt meethood=mml iitprrintt pllot;run;foreecasst llea

48、dd=122 idd=daate inttervval=monnth outt=c;run;procc gpplott daata=c;symbbol I=sspliine v=ddot c=rred;plott reesidduall*daate;run;quitt; 實(shí)驗(yàn)九九 傳遞遞函數(shù)模模型的建建模實(shí)驗(yàn)?zāi)康牡模菏煜は鬟f函函數(shù)模型型的建模模方法。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：煤氣爐爐數(shù)據(jù)。三、實(shí)驗(yàn)驗(yàn)要求：寫出實(shí)實(shí)驗(yàn)報(bào)告告，總結(jié)結(jié)傳遞函函數(shù)模型型的建模模的一般般步驟。四、實(shí)驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。2、建立立名為eexp66的SAAS數(shù)

49、據(jù)據(jù)集，輸輸入如下下程序：保存上述述程序，供以后后分析使使用（只只需按工工具條上上的保存存按鈕，然后填填寫 完完提問后后就可以以把這段段程序保保存下來來）。4、繪序序列圖，輸入如如下程序序：5、提交交程序，仔細(xì)觀觀察兩序序列圖形形，看兩兩者有何何聯(lián)系。6、先觀觀察和的相關(guān)關(guān)情況，看是否否要做差差分，輸輸入如下下程序：7、提交交程序，觀察的的自相關(guān)關(guān)和互相相關(guān)系數(shù)數(shù)，發(fā)現(xiàn)現(xiàn)都很快快的衰減減，表明明不 要要做差分分運(yùn)算。識別輸入入序列, 輸入入如下程程序：提交程序序，觀察察的自相相關(guān)和偏偏相關(guān)系系數(shù)，可可以看到到偏相關(guān)關(guān)系數(shù)是是3步 截截尾的。10、對對擬合ARR(3)模型，看是否否充分，輸入如如

50、下程序序：11、提提交程序序，觀察察輸出結(jié)結(jié)果，可可看到模模型通過過了白噪噪聲檢驗(yàn)驗(yàn)，說明明擬合效效果不錯(cuò)，把擬合合的方程程式寫出出來。12、觀觀察預(yù)白白噪聲化化后的兩兩序列的的互相關(guān)關(guān)系數(shù)，輸入如如下程序序： 133、提交交程序，觀察樣樣本自相相關(guān)系數(shù)數(shù)和偏相相關(guān)系數(shù)數(shù)和互相相關(guān)系數(shù)數(shù)，我們們可以 初初步識別別傳遞函函數(shù)模型型為（22，2，3）(思考:為什么么?)，即： 144、進(jìn)行行參數(shù)估估計(jì)，并并查看殘殘差的相相關(guān)情況況，輸入入如下程程序：15、提提交程序序，觀察察輸出結(jié)結(jié)果，可可以看到到殘差的的偏相關(guān)關(guān)系數(shù)是是2步 截截尾的。那么模模型可識識別為： 16、進(jìn)進(jìn)行參數(shù)數(shù)估計(jì)，輸入如如下程

51、序序：提交程序序，觀察察輸出結(jié)結(jié)果，可可看到很很小，且且模型通通過了白白噪聲檢檢驗(yàn)，那那么我們們除掉這這一項(xiàng),再進(jìn)行行估計(jì)，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，可可看到模模型通過過了白噪噪聲檢驗(yàn)驗(yàn)提交程程序，說說明模型型擬合充充分,請請寫出方方程式。19、進(jìn)進(jìn)行預(yù)測測，輸入入如下程程序： 200、提交交程序，觀察預(yù)預(yù)測結(jié)果果。 211、退出出SASS系統(tǒng)，關(guān)閉計(jì)計(jì)算機(jī)。dataa exxp6; innfille C:Doccumeentss annd SSetttinggsAAdmiinisstraatorr桌面面exxp1.txtt; innputt x y; t=_n_;run

52、;procc gpplott daata=expp6; symmboll1 ii=spplinne cc=reed; symmboll2 ii=spplinne cc=grreenn; ploot xx*t=1 yy*t=2;run;procc arrimaa daata=expp6;idenntiffy vvar=y ccrossscoorr=(x) nllag=12;run;procc arrimaa daata=expp6; iddenttifyy vaar=xx nllag=12 minnic;run;estiimatte pp=(11,2,3,55) nnoinnt pplott ;

53、run;estiimatte iinpuut=(3$(1,22)/(1,22)x) pllot;run;estiimatte pp=2 inpput=(3$(1,2)/(1,2)xx) pplott;run;estiimatte pp=2 inpput=(3$(1,2)/(1)x) ploot;run;foreecasst lleadd=6 ;run; 實(shí)實(shí)驗(yàn)十 回歸與與時(shí)序相相結(jié)合的的建模實(shí)驗(yàn)?zāi)康牡模菏煜はせ貧w與與時(shí)序相相結(jié)合的的建模方方法。二、實(shí)驗(yàn)驗(yàn)內(nèi)容：芝加哥哥某食品品公司大大眾食品品周銷售售數(shù)據(jù)。三、實(shí)驗(yàn)驗(yàn)要求：寫出實(shí)實(shí)驗(yàn)報(bào)告告，總結(jié)結(jié)回歸與與時(shí)序相相結(jié)合的的建模的的一般步步驟。四、實(shí)

54、驗(yàn)驗(yàn)時(shí)間：2小時(shí)時(shí)。五、實(shí)驗(yàn)驗(yàn)軟件：SASS系統(tǒng)。六、實(shí)驗(yàn)驗(yàn)步驟1、開機(jī)機(jī)進(jìn)入SSAS系系統(tǒng)。2、建立立名為eexp77的SAAS數(shù)據(jù)據(jù)集，輸輸入如下下程序：保存上述述程序，供以后后分析使使用（只只需按工工具條上上的保存存按鈕，然后填填寫 完完提問后后就可以以把這段段程序保保存下來來）。4、首先先只分析析銷售額額的數(shù)據(jù)據(jù)，不加加回歸項(xiàng)項(xiàng)。繪序序列圖，輸入如如下程序序：5、提交交程序，仔細(xì)觀觀察序列列圖形。6、初步步識別模模型，輸輸入如下下程序：7、提交交程序，觀察yy1的相相關(guān)系數(shù)數(shù)，發(fā)現(xiàn)現(xiàn)偏相關(guān)關(guān)系數(shù)是是4階截截尾的，那么 我我們初步步識別為為AR(4)模模型，進(jìn)進(jìn)行參數(shù)數(shù)估計(jì)，并觀察察殘差

55、相相關(guān)系數(shù)數(shù)。輸入入 如如下程序序：8、提交交程序，觀察輸輸出結(jié)果果，可看看到模型型擬合得得還是比比較好。然后可以以實(shí)驗(yàn)其其他一些些模型，最后根根據(jù)AIIC和BICC準(zhǔn)則，我們最最后選定定 模型型為: 下面我們們開始加加入回歸歸項(xiàng)，首首先我們們繪四個(gè)個(gè)序列的的圖形。輸入如如下程序序：提交程序序，觀察察這四個(gè)個(gè)序列有有什么特特點(diǎn)。繪y1對對y2、y3、y4的散散點(diǎn)圖，輸入如如下程序序：提交程序序，觀察察他們的的相關(guān)性性，可看看出y11和y2負(fù)相相關(guān)，yy1和y3正相相 關(guān)，而y11和y4好象象不相關(guān)關(guān)。做純回歸歸分析，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，可可看到y(tǒng)y4的系系數(shù)接近近于

56、零，我們除除掉這一一 項(xiàng)再再做回歸歸，并觀觀察殘差差的相關(guān)關(guān)系數(shù)，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果，可可看到殘殘差不是是白躁聲聲。我們們把殘差差用ARRMA 模型型擬合，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果,可以看看出模型型擬合比比較充分分，且yy4、MA11,1、 ARR1,33的系數(shù)數(shù)接近如如零，除除掉這幾幾項(xiàng)，再再觀察，輸入如如下程序序：提交程序序，觀察察輸出結(jié)結(jié)果,可以看看出模型型擬合比比較充分分，且殘殘差的標(biāo)標(biāo)準(zhǔn) 誤和和前一模模型沒有有多大變變化，且且AICC 和BICC也比前前一模型型小，故故我們 就選選擇這一一模型，把這一一結(jié)果記記錄下來來。下面我們們來看

57、看看殘差對對預(yù)測值值,y22,y33的關(guān)系系圖。輸輸入如下下程序：提交程序序，觀察察圖形，可看出出殘差對對y2,y3還還不是十十分充分分，我們們加入y2,yy3的滯滯后一階階，看結(jié)結(jié)果有什什么變化化，輸入入如下程程序： 23、提提交程序序，觀察察輸出結(jié)結(jié)果，并并與原來來結(jié)果比比較，看看是否有有進(jìn)步。(沒進(jìn)進(jìn)步) 244、進(jìn)行行預(yù)測，輸入如如下程序序： 255、提交交程序，觀察預(yù)預(yù)測結(jié)果果。 266、退出出SASS系統(tǒng)，關(guān)閉計(jì)計(jì)算機(jī)。dataa exxp7; innfille C:Doccumeentss annd SSetttinggsAAdmiinisstraatorr桌面面exxp1.tx

58、tt; innputt y11 y22 y33 y44 ; daate=inttnx(weeek,114seep911d,_n_-1); foormaat ddatee daate99.;run;procc gpplott daata=expp7; symmboll1 ii=spplinne cc=reed; ploot yy1*ddatee=1;run;procc arrimaa daata=expp7; iddenttifyy vaar=yy1 nnlagg=155 miinicc;run;estiimatte pp=(11,4) pllot;run;procc gpplott daata

59、=expp7;symbbol33 i=spllinee c=greeen;plott y11*daate=3 yy2*ddatee=3 y3*datte=33 y44*daate=3;run;procc pllot datta=eexp77; ploot yy1*yy2=* y1*y3=* y11*y44=*;run;procc arrimaa daata=expp7; ideentiify varr=y11 crrossscorrr=(y2 y3 y4) noopriint; esttimaate inpput=(y22 y33 y44);run;idenntiffy vvar=y1 cro

60、ossccorrr=(yy2 yy3);estiimatte iinpuut=(y2 y3) pllot;run;idenntiffy vvar=y1 croossccorrr=(yy2 yy3 yy4);estiimatte pp=4 q=33 innputt=(yy2 yy3 yy4) ploot;run;idenntiffy vvar=y1 croossccorrr=(yy2 yy3);estiimatte pp=(11,2,4) q=(2,33) iinpuut=(y2 y3 ) pplott;run;foreecasst lleadd=6 id=datte iinteervaal=w