西安市社会固定资产总投资拟合数据

1、统计预测与决策课程论文题 目基于ARMA模型的西安市社会固定资产总投资时间序列分析与预测 学生姓名谭嘉宝学生学号13610801150962专 业经济统计学班 级统本统计1302班提交日期二一六年五月基于ARMA模型的我国TIIFA时间序列分析与预测摘要:本文分析了西安市1989-2009年的社会固定资产总投时间序列，在将该时间序列平稳化的基础上，建立自回归移动平均模（ARMA）,从中得出西安市全社会固定资产投资序列的变化规律，并且预测2010、2011年社会固定资产总投资的数值。关键词：时间序列预测；TIIFA ；ARMA模型 1. 前 言全社会固定资产投资包括基本建设投资、更新改造投资、国

2、有单位其他固定资产投资、房地产开发投资、城镇集体固定资产投资、联营经济、股份制经济、外商投资经济、港澳台投资经济及其他经济类型的固定资产投资，农村集体5万元以上固定资产投资，城镇工矿区私人建房投资和国防、人防基本建设投资。全社会固定资产投资按登记注册类型可分为：国有、集体、个体、联营、股份制、外商、港澳台商、其他等。按照管理渠道可分为：基本建设、更新改造、房地产开发和其他固定资产投资四个部分。按经济统计可划分第一产业是指农、林、牧、渔业。第二产业是指采矿业、制造业、电力、燃气及水的生产和供应业、建筑业。第三产业是指除第一、二产业以外的其他行业。第三产业包括：交通运输、仓储和邮政业，信息传输、计

3、算机服务和软件业，批发和零售业，住宿和餐饮业，金融业，房地产业，租赁和商务服务业，科学研究、技术服务和地质勘查业，水利、环境和公共设施管理业，居民服务和其他服务业，教育，卫生、社会保障和社会福利业，文化、体育和娱乐业，公共管理和社会组织，国际组织。2. ARMA模型2.1 ARMA模型概述ARMA模全称为自回归移动平均模型(Auto-regressive Moving Average Model，简称 ARMA)是研究时间序列的重要方法。其在经济预测过程中既考虑了经济现象在时间序列上的依存性, 又考虑了随机波动的干扰性, 对经济运行短期趋势的预测准确率较高, 是近年应用比较广泛的方法之一。AR

4、MA模型是由美国统计学家G.E.P.Box 和 G.M.Jenkins在20世纪70年代提出的著名时序分析模型,即自回归移动平均模型。ARMA模型有自回归模型AR(q)、移动平均模型MR(q)、自回归移动平均模型ARMA(p,q) 3种基本类型。其中ARMA(p,q）自回归移动平均模型，模型可表示为：其中，为自回归模型的阶数，为移动平均模型的介数；表示时间序列在时刻的值；为自回归系数；表示移动平均系数；表示时间序列在时期的误差或偏差。2.2 ARMA模型建模流程首先用ARMA模型预测要求序列必须是平稳的，也就是说，在研究的时间范围内研究对象受到的影响因素必须基本相同。若所给的序列并非稳定序列，

5、则必须对所给的序列做预处理，使其平稳化，然后用ARMA模型建模。建模的基本步骤如下：（1）求出该观察值序列的样本自相关系数（ACF)和样本偏相关（PACF)的值。（2）根据样本自相关系数和偏自相关系数的性质选择适当的模型进行拟合。 （3）估计模型中未知参数的值。 （4）检验模型的有效性。如果拟合模型通不过检验，转向步骤（2），重新选择模型再拟合。 （5）模型优化。如果拟合模型通过检验，仍然转向步骤（2)，充分考虑各种可能，建立多个拟合模型，从所有通过检验的拟合模型中选择最优模型。（6）利用拟合模型，预测序列的将来走势。 3. 西安市社会固定资产总投资时间序列模型的建立3.1 数据的预处理 本文

6、选取了西安市1989-2009年的社会固定资产总投资数据作为时间序列观察值。对此时间序列做时序图如图1所示：图1 西安市社会固定资产总投资时序图由时间序列的时序图可以发现西安社会固定资产总投资随时间的增长是呈指数趋势。因此，对原始序列作对数变换并作出其时序图如图1所示：图1 取对数后的西安市社会固定资产总投资时序图通过观察取对数后的社会固定资产总投资时序图，发现经过处理后的序列具有趋势性。由于固定资产总投资带有很强的趋势成分, 而我们的目的主要是利用ARMA 模型对其周期成分进行分析, 因此需要对此类的数据先进行消除趋势性的处理, 然后建立ARMA模型。拿到观察值序列之后，无论是采用确定性时序

7、分析方法还是随机时序分析方法，分析的第一步都是要通过有效的手段提取信息中所蕴含的确定性信息。在Box和Jenkins在Time Series Analysis Forecasting and Control一书中特别强调差分方法的使用，他们使用大量的案例分析证明差分方法是一种非常简便有效的确定性信息提取方法。实践中，我们会根据序列的不同特点选择合适的差分方式，常见情况有以下三种；序列蕴含着显著的线性趋势，一阶差分就可以实现趋势平稳。序列蕴含着曲线趋势，通常低阶（2阶或3阶）差分就可以提取出曲线趋势的影响。蕴含固定周期的序列，一般进行步长为周期长度的差分运算就可以较好地提取周期信息。从理论上而言

8、，足够多次的差分运算可以充分地提取原序列中的非平稳确定性信息。但应当注意的是，差分运算的阶数并不是越多越好。差分运算是一种对信息的提取加工过程，每次差分都会有信息的损失，在实际中差分运算的阶数要适当，应当避免过差分。观察时序图2，可使用一阶差分就可以提取序列的足够信息。做一阶差分后，做其序列图3如下：图3 一阶差分后对数西安市社会固定资产总投资时序图从图（3）可以观察得出，序列大致趋于平稳。为了进一步检验序列是否真正平稳，在此使用Eviews统计软件对已转换进行平稳性检验。对时间序列的平稳性有两种检验方法，一种是根据时序图和自相关图显示的特征作出判断的图检验方法；一种是构造检验统计量进行假设检

9、验的方法。目前最常用的平稳性统计检验方法是单位根检验（unit root test）。使用单位根检验法对变换数据进行检验得出检验结果如表1所示：表1 DLOGTIIFA一阶差分单位根检验结合图3与表1，结果表明序列logTIIFA经过一阶差分之后序列平稳。3.2 模型的识别与选择计算出样本自相关系数和偏相关系数的值之后，我们主要是根据它们表现出来的性质，选择适当的ARMA模型拟合观察值序列。这个过程实际上就是要根据样本自相关系数和偏相关系数的性质估计自相关阶数 和移动平均阶数，因此模型的识别过程也成为定阶过程。一般ARMA模型定阶的基本原则如表2示：表2 ARMA(p,q)模型选择原则ACFP

10、ACF模型定阶拖尾p阶截尾AR(p)模型q阶截尾拖尾MA(q)模型拖尾拖尾ARMA(p,q)模型利用Eviews统计软件对差分数据进行操作，可得样本自相关系数和偏相关系数图如图4所示： 图4 差分序列自相关系数与偏相关系数图通过对一阶差分的对数序列的自相关系数和偏相关系数图的分析观察，可以知道模型大致可选取两种模型。第一种，自相关系数为拖尾，而偏相关系数为一阶截尾。此时选取模型可以为AR（1）模型。第二种，自相关二阶截尾，而偏相关系数为一阶截尾。此时选取模型可以为AR(1）模型。3.3 参数估计 选择拟合好后的模型之后，下一步就是要利用序列的观察值确定该模型的口径，即估计模型中未知参数的值。对

11、于一个非中心化ARMA（p,q）模型，有 式中，该模型共含个未知参数：。对于未知参数的估计方法有三种：矩估计极大似然估计和最小二乘估计。其中本文使用最小二乘估计法对序列进行参数估计。在ARMA(p,q)模型场合，记 残差项为： 残差平方和为： 是残差平方和达到最小的那组参数值即为的最小估计值。使用Eviews统计软件操作可得序列两种可能的参数估计图如图5、6所示：Dependent Variable: LNXMethod: Least SquaresDate: 05/30/16 Time: 19:34Sample (adjusted): 1990 2009Included observatio

12、ns: 20 after adjustmentsConvergence achieved after 4 iterationsVariableCoefficientStd. Errort-StatisticProb.  C-1.1486634.756488-0.2414940.8119AR(1)1.0281380.01759458.435250.0000R-squared0.994756    Mean dependent var6.582884Adjusted R-squared0.994465   &

13、#160;S.D. dependent var1.212894S.E. of regression0.090237    Akaike info criterion-1.878122Sum squared resid0.146568    Schwarz criterion-1.778549Log likelihood20.78122    F-statistic3414.678Durbin-Watson stat1.890547    

14、;Prob(F-statistic)0.000000Inverted AR Roots      1.03Estimated AR process is nonstationary图5 AR (1)模型参数估计与检验结果Dependent Variable: LNXMethod: Least SquaresDate: 05/30/16 Time: 19:35Sample (adjusted): 1990 2009Included observations: 20 after adjustmentsConvergence achieve

15、d after 10 iterationsBackcast: 1988 1989VariableCoefficientStd. Errort-StatisticProb.  AR(1)1.0330420.003329310.32080.0000MA(2)0.0763110.2351160.3245680.7492R-squared0.994757    Mean dependent var6.582884Adjusted R-squared0.994466    S.D. dependent v

16、ar1.212894S.E. of regression0.090232    Akaike info criterion-1.878234Sum squared resid0.146551    Schwarz criterion-1.778661Log likelihood20.78234    Durbin-Watson stat1.866326Inverted AR Roots      1.03Estima

17、ted AR process is nonstationary图6 AR (1)模型参数估计与检验结果由图67模型的参数估计与检验结果对比看，可以知道，AR(1）模型中其调整后的为0.994756接近ARMA（1,2）模型中的0.994757；而AIC和SC值分别为-1.878234 , -1.778661分别大于AR(1)模型中的-2.559993，-2.453418。根据以上模型的识别与选择，我们选用了AR(1)作为最佳预测模型。估计该模型的参数及模型的相关检验结果如图7。结果表明, 模型AR(1)的参数估计值具有统计意义。其展开式为:3.4 参数的显著性检验参数的显著性检验就是要检验每一

18、个未知参数是否显著非零。这个检验的目的是为了是使模型最精简。如果某个参数不显著，即表示该参数所对应的那个自变量对因变量的影响不明显，该自变量就可以从拟合模型中删除。最终模型将由一系列参数显著非零的自变量表示。由图7模型参数估计与检验结果，可以观察到t统计量值的值均小于0.05。表明模型参数显著。3.5 预测序列走势 由预测方程及其条件方程：经预测得到2010、2011年的TIIFA值分别为396975.3亿元、484638.24亿元 ,其标准差为0.06697。预测值与真实值误差均在3%以内预测较为准确。利用此模型对20102011年TIIFA进行预测结果如表3所示：表3 20102011年模型预测值年份20102011预测值（亿元）8670.29411557.944.结 论时间序列分析的ARMA 模型预测问题, 实质上是通过对社会经济发展变化过程的分析研究, 找出其发展变化的量变规律性, 用以预测经济现象的未来。预测时不必考虑其他因素的影响, 仅从序列自身出发, 建立相应的模型进行预测, 这就从根本上避免了寻找主要因素及识别主要因素和次要因素的困难; 和回归分析相比, 可以避免了寻找因果模型中对随机扰动项的限定条件在经济实践中难以满足的矛盾。实际上这也是ARMA 模型预测与其他预测方法相比的优越性所在。本文运用时间序列的分析方法，对我国历年的国