2022年时间序列和ARIMA模型 .pdf
实验五 ARIMA 模型的概念和构造一、实验目的了解 AR，MA以及 ARIMA模型的特点，了解三者之间的区别联系，以及AR与 MA的转换，掌握如何利用自相关系数和偏自相关系数对ARIMA模型进行识别，利用最小二乘法等方法对ARIMA模型进行估计，利用信息准则对估计的ARIMA模型进行诊断，以及如何利用ARIMA模型进行预测。掌握在实证研究如何运用Eviews 软件进行ARIMA模型的识别、诊断、估计和预测。二、基本概念所谓 ARIMA模型，是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA 模型根据原序列是否平稳以及回归中所含部分的不同，包括移动平均过程（MA）、自回归过程（AR）、自回归移动平均过程（ARMA）以及 ARIMA过程。在 ARIMA模型的识别过程中，我们主要用到两个工具：自相关函数(简称 ACF），偏自相关函数(简称 PACF)以及它们各自的相关图（即ACF、PACF相对于滞后长度描图）。对于一个序列来说，它的第 j 阶自相关系数(记作 )定义为它的j 阶自协方差除以它的方差，即，它是关于j 的函数，因此我们也称之为自相关函数，通常记ACF(j)。偏自相关函数 PACF(j)度量了消除中间滞后项影响后两滞后变量之间的相关关系。三、实验内容及要求1、实验内容：根据 1991 年 1 月 2005 年 1 月我国货币供应量（广义货币M2）的月度时间数据来说明在 Eviews3.1 软件中如何利用B-J 方法论建立合适的ARIMA（p,d,q）模型，并利用此模型进行数据的预测。2、实验要求：（1）深刻理解上述基本概念；（2）思考：如何通过观察自相关，偏自相关系数及其图形，利用最小二乘法，以及信息准则建立合适的ARIMA模型；如何利用ARIMA模型进行预测；（3）熟练掌握相关Eviews 操作。四、实验指导1、ARIMA模型的识别（1）导入数据打开 Eviews 软件，选择“File”菜单中的“New-Workfile”选项，出现“Workfile Range”对话框，在“Workfile frequency”框中选择“Monthly”，在“Start date”和“End date”框中分别输入“1991:01”和“2005:01”，然后单击“OK”，选择“File”菜单中的“Import-Read Text-Lotus-Excel”选项，找到要导入的名为EX6.2.xls的 Excel 文档，单击“打开”出现“Excel Spreadsheet Import”对话框并在其中输入相关数据名称(M2)，再单击“OK”完成数据导入。（2）模型的识别首先利用ADF检验，确定 d 值，判断 M2序列为 2 阶非平稳过程（由于具体操作方法我们在第五章中予以说明，此处略），即 d 的值为 2，将两次差分后得到的平稳序列命名为W2；下面我们来看W2的自相关、偏自相关函数图。打开 W2序列，点击“View”“Correlogram”菜单，会弹出如图51 所示的窗口，名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 4 页 -图 51 自相关形式设定我们选择滞后项数为36，然后点击“OK”，就得到了W2的自相关函数图和偏自相关函数图，如图 52 所示。图 52 W2 自相关函数图和偏自相关函数图从 W2的自相关函数图和偏自相关函数图中我们可以看到，他们都是拖尾的，因此可设定为 ARMA 过程。W2的自相关函数1-5 阶都是显著的，并且从第6 阶开始下降很大，数值也不太显著，因此我们先设定q 值为 5。W2的偏自相关函数1-2 阶都很显著，并且从第3 阶开始下降很大，因此我们先设定 p 的值为 2，于是对于序列W2，我们初步建立了ARMA(2,5)模型。2、模型的估计点击“Quick”“Estimate Equation”，会弹出如图53 所示的窗口，在“Equation Specification”空白栏中键入“W2 C MA(1)MA(2)MA(3)MA(4)MA(5)AR(1)AR(2)”，在“Estimation Settings”中选择“LS-Least Squares(NLS and ARMA)”，然后“OK”，得到如图5 4所示的估计结果。图 53 回归方程设定图 54 ARMA(2,5)回归结果可以看到，除常数项外，其它解释变量的系数估计值在15%的显著性水平下都是显著的。3、模型的诊断点击“View”“Residual test”“Correlogram-Q-statistics”，在弹出的窗口中选择滞后阶数为36，点击“Ok”，就可以得到Q统计量，此时为30.96，p 值为 0.367，因此不能拒绝原假设，可以认为模型较好的拟合了数据。我们再来看是否存在一个更好的模型。我们的做法是增加模型的滞后长度，然后根据信息值来判断。表5-1 是我们试验的几个p,q值的 AIC 信息值。表 5-1 不同 p,q 值的 AIC 信息值p 2 3 4 2 2 2 3 3 3 4 4 4 q 5 5 5 名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 4 页 -6 7 8 6 7 8 6 7 8 AIC 16.78 16.75 16.77 16.76 16.76 16.77 16.77 16.78 16.79 16.75 16.79 16.78 可以看到，根据AIC 信息值，我们应选择p=3、q=5 或 p=4、q=6，但是按照后者建立的模型中有的解释变量的系数估计值是不显著的，而按照前者建立的模型其解释变量的系数值都是显著的（如图55 所示），因此我们最终建立的模型是ARMA(3,5)。图 55 ARMA(3,5)回归结果4、模型的预测点击“Forecast”，会弹出如图56 所示的窗口。在 Eviews 中有两种预测方式：“Dynamic”和“Static”，前者是根据所选择的一定的估计区间，进行多步向前预测；后者是只滚动的进行向前一步预测，即每预测一次，用真实值代替预测值，加入到估计区间，再进行向前一步预测。我们首先用前者来估计2003年 1月到 2005 年 1月的 W2，在“Sample range for forecast”空白栏中键入“2003:01 2005:01”（如图 56 所示），选择“Dynamic”，其他的一些选项诸如预测序列的名称、以及输出结果的形式等，我们可以根据目的自行选择，不再介绍，点击“OK”，得到如图57 所示的预测结果。图 56 ARMA(3,5)模型预测设定图 57 Dynamic预测方式结果图中实线代表的是W2的预测值，两条虚线则提供了2 倍标准差的置信区间。可以看到，正如我们在前面所讲的，随着预测时间的增长，预测值很快趋向于序列的均值（接近0）。图名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 4 页 -的右边列出的是评价预测的一些标准，如平均预测误差平方和的平方根（RMSE），Theil不相等系数及其分解。可以看到，Theil不相等系数为0.82，表明模型的预测能力不太好，而对它的分解表明偏误比例很小，方差比例较大，说明实际序列的波动较大，而模拟序列的波动较小，这可能是由于预测时间过长。下面我们再利用“Static”方法来估计2004 年 1 月 2005 年 1 月的 W2，（操作过程略），我们可以得到如图58 所示的结果。从图中可以看到，“Static”方法得到的预测值波动性要大；同时，方差比例的下降也表明较好的模拟了实际序列的波动，Theil不相等系数为0.62，其中协方差比例为0.70，表明模型的预测结果较理想。图 58 Static预测方式结果名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 4 页 -
