多元时间序列建模分析.doc

应用时间序列分析实验报告实验名称多元时间序列建模分析姓名学号班级实验地点实验日期指导教师实验目的：1、熟悉单位根检验；2、掌握ARIMAX模型建模。涉及实验的相关情况介绍（包含使用软件或实验设备等情况）：SAS、excel表格、word。实验内容：1、 我国1950-2008年进出口总额数据（单位：亿元）如表6-15所示。表6-15年份出口总额进口总额19502021.3195124.235.3195227.137.5195334.846.119544044.7195548.761.1195655.753195754.55019586761.7195978.171.2196063.365.1196147.743196247.133.819635035.7196455.442.1196563.155.319666661.1196758.853.4196857.650.9196959.847.2197056.856.1197168.552.4197282.9641973116.9103.61974139.4152.81975143147.41976134.8129.31977139.7132.81978167.6187.41979211.7242.91980271.2298.81981367.6367.71982413.8357.51983438.3421.81984580.5620.51985808.91257.819861082.11498.3198714701614.219881766.72055.1198919562199.919902985.82574.319913827.13398.719924676.34443.319935284.85986.2199410421.89960.1199512451.811048.1199612576.411557.4199715160.711806.5199815223.611626.1199916159.813736.5200020634.418638.8200122024.420159.2200226947.924430.3200336287.934195.6200449103.346435.8200562648.154273.7200677594.663376.9200793455.673284.62008100394.979526.5（1）使用单位根检验，分别考察进口总额和出口总额序列的平稳。（2）分别对进口总额序列和出口总额数据拟合模型。（3）考察这两个序列是否具有协整关系。（4）如果这两个序列具有协整关系，请建立适当模型拟合它们之间的相关关系。（5）构造该协整模型的误差修正模型。实验过程记录（含程序、数据记录及分析和实验结果等）：时序图如下：单位根检验输出结果如下：序列x的单位根检验结果：序列y的单位根检验结果：序列y和序列x之间的相关图如下：残差序列自相关图： 自相关图显示。延迟6阶之后自相关系数都在2倍标准差范围之内，可以认为残差序列平稳。对残差序列进行2阶自相关单位根检验，检验结果显示残差序列显著平稳，如下图：残差序列单位根检验结果：残差序列平稳，说明序列Y与序列X之间具有协整关系，我可以大胆的在这两个序列之间建立回归模型而不必担心虚假回归问题。考察残差序列白噪声检验结果，如下图：残差序列白噪声检验结果：输出结果显示，延迟各阶LB统计量的P值都大于显著水平0.05，可以认为残差序列为白噪声检验结果，结束分析。出口序列拟合的模型为：lnxtARIMA(1,1,0),具体口径为：进口序列拟合的模型为 lnytARIMA(1,1,0) ，具体口径为：lnyt和lnxt具有协整关系。协整模型为：误差修正模型为：SAS程序如下：data example6_4;input x y;t=_n_;cards;195020.021.3195124.235.3195227.137.5195334.846.1195440.044.7195548.761.1195655.753.0195754.550.0195867.061.7195978.171.2196063.365.1196147.743.0196247.133.8196350.035.7196455.442.1196563.155.3196666.061.1196758.853.4196857.650.9196959.847.2197056.856.1197168.552.4197282.964.01973116.9103.61974139.4152.81975143.0147.41976134.8129.31977139.7132.81978167.6187.41979211.7242.91980271.2298.81981367.6367.71982413.8357.51983438.3421.81984580.5620.51985808.91257.819861082.11498.319871470.01614.219881766.72055.119891956.02199.919902985.82574.319913827.13398.719924676.34443.319935284.85986.2199410421.89960.1199512451.811048.1199612576.411557.4199715160.711806.5199815223.611626.1199916159.813736.5200020634.418638.8200122024.420159.2200226947.924430.3200336287.934195.6200449103.346435.8200562648.154273.7200677594.663376.9200793455.673284.62008100394.979526.5run;proc gplot;plot x*t=1 y*t=2/overlay;symbol1 c=black i=join v=none;symbol2 c=red i=join v=none w=2 l=2;run;proc arima data=example6_4;identify var=x stationarity=(adf=1);identify var=y stationarity=(adf=1);run;proc arima;identify var=y crrosscorr=x;estimate methed=ml input=x plot;forecast lead=0 id=t out=out;proc aima data=out;identify varresidual st