试验ARIMA模型分析

1、第七章ARIMA模型分析一、实验目的了解AR,MA以及ARIMA模型的特点，了解三者之间的区别联系，以及AR与MA的转换，掌握如何利用自相关系数和偏自相关系数对ARIMA模型进行识别，利用最小二乘法等方法对ARIMA模型进行估计，利用信息准则对估计的ARIMA模型进行诊断，以及如何利用ARIMA模型进行预测。掌握在实证研究如何运用Eviews软件进行ARIMA模型的识别、诊断、估计和预测。二、基本概念所谓ARIMA模型，是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA模型根据原序列是否平稳以及回归中所含部分的不同，包

2、括移动平均过程(MA)、自回归过程(AR)、自回归移动平均过程(ARMA)以及ARIMA过程。在ARIMA模型的识别过程中，我们主要用到两个工具：自相关函数(简称ACF),偏自相关函数(简称PACF)以及它们各自的相关图(即ACF、PACF相对于滞后长度描图)。对于一个序列Yt来说，它的第j阶自相关系数(记作j)定义为它的j阶自协方差除以它的方差，即j=j/0，它是关于j的函数，因此我们也称之为自相关函数，通常记ACF(j)o偏自相关函数PACF(j)度量了消除中间滞后项影响后两滞后变量之间的相关关系。三、实验内容及要求1、实验内容：根据1991年1月2005年1月我国货币供应量(广义货币M2

3、)的月度时间数据来说明在Eviews软件中如何利用B-J方法论建立合适的ARIMA(p,d,q)模型，并利用此模型进行数据的预测。2、实验要求：(1)深刻理解上述基本概念；(2)思考：如何通过观察自相关，偏自相关系数及其图形，利用最小二乘法，以及信息准则建立合适的ARIMA模型；如何利用ARIMA模型进行预测；(3)熟练掌握相关Eviews操作。四、实验指导1、ARIMA模型的识别(1)导入数据打开Eviews软件，选择File菜单中的New-WorMle选项，出现WorkfileRange对话框，在Workfilefrequency框中选择Monthly,在Startdate和Enddate

4、椎中分别输入1991:01”和2005:01，然后单击“OK；选择“FileW中的“Import-ReadText-Lotus-Excel选项，找到要导入的名为EX6.2.xls的Excel文档，单击打开出现ExcelSpreadsheetImport对话框并在其中输入相关数据名称(M2),再单击“OK完成数据导入。(2)模型的识别首先利用ADF检验，确定d值，判断M2序列为2阶非平稳过程(由于具体操作方法我们在第五章中予以说明，此处略)，即d的值为2,将两次差分后得到的平稳序列命名为W2;下面我们来看W2的自相关、偏自相关函数图。打开W2序列，点击View”Correlogram菜单，会弹出

5、如图71所示的窗口，图71自相关形式设定我们选择滞后项数为36,然后点击“OK,就得到了W2的自相关函数图和偏自相关函数图,如图72所示。AutocorrelationPartialCorrelationAC=PACQ-StatProti1111:1c11.1II111II11111111II11N11111il1i1I IIl1llIIli1lI:aL1_II匚11n1I11iEII111IIIII111111llii:11IIII1II111111IIn11IIII1111II11IIIIII1IIII1-0.701-070183.J370.00020.141-0.608BS.B140.0

6、0030.277-0.143100.040,0004-0.453-0295135.53000050.371-0.224159.500,00060.135-0.156162.690,0007-0.0S5-0.095163970.00080.143-0.304167.590,0009-O.aiO0.052167.610,00010-0.1420029171.240,000110.143-0.121174.960.00012-0.011-0.015174.980,00013-0.1150.122177420.000140.136-0.050180.960,00015-001

7、0,00016-0.Q140054132.050000170.066-0.130182.630,00013-0.023-0.0971S272O.OOC19-0.0300.057182.090,000200.028-0.162183.0J0,000210.042-0.092183.380,00022-0.140-0.1041S7.190,000230.2030.042195.660,00024-0.1540.070200.350,00025-0.0190.0462口口420.00026口.15-0,108205.200,00027-0.175-0.103211.410,000280.117-0.

8、019214.170,00029-0.002OOG22U.170,00030-0.069-0.032215.200000310.131O,OB5219770.00032-0.0810.067221160.00033-0.0Q8-0015221.320,000340.075-0.118222.500,00035-0B120.02B222.53O.OOC3B00550DOB223.190,000图72W2自相关函数图和偏自相关函数图从W2的自相关函数图和偏自相关函数图中我们可以看到，他们都是拖尾的，因此可设定为ARMA过程。W2的自相关函数1-5阶都是显著的，并且从第6阶开始下降很大，数值也不太显

9、著，因此我们先设定q值为5。W2的偏自相关函数1-2阶都很显著，并且从第3阶开始下降很大，因此我们先设定p的值为2,于是对于序列W2,我们初步建立了ARMA(2,5)模型。2、模型的估计点击QuickEstimateEquation,会弹出如图73所示的窗口，在aEquationSpecification空白栏中键入“W2CMA(1)MA(2)MA(3)MA(4)MA(5)AR(1)AR(2)”，在“EstimationSettings中选择”“L-LeastSquares(NLSandARMA)”，然后“OK；得到如图74所示的估计结果。图73回归方程设定FView?-Eguatxon:UK

10、TTTLEDTorkfile;WTITLEDIFileFrintNun*FrezeForacastStatsResideEdit.OfcjeelsViewTjfocsQuickOtiansIPiTtdowHelpDependentVariatle:W2Method:LeastSquaresDate:03/26/05Time:23:26Sample(adjustecT):1991,052005:01Includedobservations:165afteradjustingendpaintsConvergenceachieved白R曰i4Citerations=iackcast:1990:121

11、991:D4variabeCoefficientStd.ErrorStatisticProb.YI.i,，1.H(2aG(3(4(5crraaaaaAAMMMMM14.4260410784D6133772001829-10071580.05496918.322280.0000-0.8399380.046282-1B.14S35c.oooc-0.412644C.037434-4719462C.00000.1427600.0914731.560S840.1206-D.360579O.D5O160-14,30273u.oooc303361C口365西3.503370000060.21S9640039

12、64724425200.0157R-squared0.805271Meandependentvai1634363Adjusted艮-总qunr白日0.796590S.D.dependentvar2309544S.E.ofnegrassiori1041.632Akaikeinfocriterion1B7E223Sumsquaredresid17口EMBSchvwgrzcriterion16.932B3Loglikelihood-137EM4F-statistic92,74953Durbin-Watsonstat2053S93ProbfF-staiistic)o.oooooo图74ARMA(2,5

13、)回归结果可以看到，除常数项外，其它解释变量的系数估计值在15%的显著性水平下都是显著的。3、模型的诊断点击View”ResidualtestGorfelogram-Q-statistics,在弹出的窗口中选择滞后阶数为36,点击“0匕就可以得到Q统计量，此时为30.96,p值为0.367,因此不能拒绝原假设，可以认为模型较好的拟合了数据。我们再来看是否存在一个更好的模型。我们的做法是增加模型的滞后长度，然后根据信息值来判断。表5-1是我们试验的几个p,q值白AIC信息值。表7-1不同p,q值的AIC信息值p234222333444q555678678678AI16.716.716.716.7

14、16.716.716.716.716.716.716.716.7C857667789598可以看到，根据AIC信息值，我彳门应选择p=3、q=5或p=4、q=6,但是按照后者建立的模型中有的解释变量的系数估计值是不显著的，而按照前者建立的模型其解释变量的系数值都是显著的(如图75所示)，因此我们最终建立的模型是ARMA(3,5)。ri.nEViews-EqualzLDn=UNTITLEDTorlfile：tTHTITLEDlFileEd.1+01?jaCt.EV1-fewFrQua0口卫士1Lw1Helpin|Fro%|Qbj*ctTrint|Eftim&t?|口口工tllUfi也工epend

15、entVariable:Method:LeastSquaresDat&:03/27/05Tune;0D:2GSample(adjusted)：1991:0G2006:01Includedobsenraiicns:164afteradjustingendpointsCorvergenceachievedafter24ileralionsBackcasi：1991:011991.05VanbleCoefticientStd.Errort-StatisticProb.U.490798.6215171600770O.OS48-1.7061990074288-22.967350OCOO-1.520646

16、0.099913-15.29S80O.OCOO-0.6067520.072167过407574O.OCOO0.29510100409927.2006720.0000-2196470065464-33502090.0010-0.6406090.05700311,23621D.OCOO-0.3934270.03B2791027778D.OCOO0.4792220.0E7O9471425D?oncnoR-squared614916Meandependentvar16,18107AdjustedR-siuared0.805363S.D.dependentvar2316.617S.E.ofregress

17、ion1022037Akaikeinfocriterion16,75030Sumsquaredresid162E-HJ6Schwarzcriterioin1G,92041Loglikelihood-1364.524F-startistic05.30712Durbin-Watsonstat2.012600Prob(F*tatistic)o.ooocoo图75ARMA(3,5)回归结果4、模型的预测点击Forecast,会弹出如图76所示的窗口。在Eviews中有两种预测方式：Dynamic和“Static，前者是根据所选择的一定的估计区间，进行多步向前预测；后者是只滚动的进行向前一步预测，即每预

18、测一次，用真实值代替预测值，加入到估计区间，再进行向前一步预测。我们首先用前者来估计2003年1月至U2005年1月的W2，在Samplerangeforforecast”空白栏中键入“2003:012005:01（如图76所示）,选择Dynamic”，其他的一些选项诸如预测序列的名称、以及输出结果的形式等，我们可以根据目的自行选择，不再介绍，点击“OK；得到如图77所示的预测结果。图76ARMA（3,5）模型预测设定ForecastW2FActual:W2Forecastsample:2003:012005;Includedobservations:25RootVeanSquaredError1509,882MeanAbsoluteError1141.053MeanAbs.PercentError122,7764TheilIneqvalilyUoeffkient0.320159BiasProportion。000010VarianceProportion0.436438CovarianceProportion0561552W2F-SE图77Dynamic预测方式结果图中实线代表的是W2的预测值，两条虚线则提供了2倍标准差的置