对工业总产值的分析

1、对工业总产值的分析问题重述摘要：随着社会的发展，对工业的投入也逐渐增加。从而使工业得到迅速的发展。某地的工业也受到影响。于是有必要对工业总产值进行研究，以下是我们对该地区工业总产值的研究。分析该地区1990年1月到1997年12月工业总产值数据，建立时间序列模型分析该地区工业总产值变化特征并且用该模型预测1997年以后该地区工业总产值。关键字：时间 季节因子 差分 自相关与偏相关 ARIMA模型模型假设：某地区的工业生产总值在一段时间内保持稳步发展；符号说明：表示原时间序列；：表示时间序列的一阶对数差分；：表示序列的一阶季节差分；：表示对序列做差分计算。问题重述：下表是某地的工业总产值数据表.

2、年月199019911992199319941995199619971月1421.41757.81984.22179.12903.32996.73476.63843.842月1367.41485.71812.42408.72513.82740.32970.33181.263月1719.71893.92274.72869.434093580.93942.64404.494月1759.61969.82328.92916.73499.53746.34067.64520.185月1795.72033.72373.13022.13642.63817.94746.8994638.996月1848.1210

3、32515.83274.53871.44046.64417.2994969.937月1637.31836.322882862.933733483.93806.84146.8998月1637.61914.723212864.23463.43510.63746.34198.79月1637.62022.22441.129083663.743703.14011.14536.83910月1637.62045.12502.62911.83753.383810.74129.64783.9111月1637.62069.22608.83101.33973.1740914372.1995034.93912月163

4、7.621362823.83664.34469.024650.7994991.55545.74要求:1.根据数据分析当地工业总产值的变化特征.2.根据变化特征试建立合理的模型描绘这种特征.3.若有季节性变化,试分离出季节性变化因子,求出季节性因子.4.对残差进行白噪声检验.5.预测1998年的工业总产值.问题分析： 这是一个有关时间序列的问题，我们对数据分析得到数据有明显的增长趋势且改时间序列有季节性变化，于是需要利用Eviews软件对该时间序列进行差分变换后建立平稳的时间序列模型求解及预测。模型的的建立、求解与选择：1.时间序列特征分析： 将数据绘制成折线图，如图1所示，序列具有明显的增长趋

5、势，并包含有周期为12个月的季节波动。即有季节因子存在。 图2是序列自相关图。由图1和图2可知，改时间序列为非平稳时间序列。因此需要对其进行调整使之变成平稳系列在进行求解。图1 工业生产值折线图图2 序列自相关图为消除趋势同时减少序列的波动，即使之变成平稳时间序列。对原序列做一阶对数差分。差分后序列名为ilx，其自相关与偏向关分析图如图3所示。图3 序列ilx自相关-偏相关分析图图4 序列ilx折线图由图3，图4可见，序列的趋势基本消除，但是当k=12时，由图3知，样本的自相关系数和偏相关系数显著不为0。表明季节性还存在。因此对序列ilx做季节差分，得到新序列silx。为检验模型的预测的效果，

6、我们这将1997年的12个观测值留出，作为评价预测的精度的参照对象。建模的样本期为1990年1月至1996年12月。绘制silx自相关和偏相关分析图，如图5所示。图5 序列silx自相关-偏相关分析图由图5可知，序列样本自相关与偏相关系数很快落入随即区间，故序列趋势已基本消除，并且当k=12时，自相关与偏相关系数也明显减小。偏相关系数与0无显著差别。图5中自相关系数与0有显著性差别。我们对序列做二阶差分。查分后的得到新序列ssilx。如图6所示。图6 序列ssilx自相关-偏相关分析图由图6可见, 序列样本自相关与偏相关系数很快落入随即区间，故序列趋势已基本消除，并且当k=12时，自相关与偏相

7、关系数没有减小,反而增大。对序列进行二阶差分，序列季节性没有得到明显改善。故对该序列只需要做阶一差分即可。对系列silx进行0均值检验的结果如下： 得到该系列样本的平均数是m=-0.463，均值标准误差s=0.182，系列均值与0无显著性的差异，表明系列可以直接建立ARMA模型。 2.模型识别因为经过一阶逐期差分，序列趋势消除，故d=1；经过一阶季节差分，季节性基本消除，故D=1.所以选用ARIMAM模型。取自然对术后的工业总产值序列为ilx。观察序列silx的偏相关图，如图4所示，p=2或3比较合适；自相关图显示q=1。考虑到AR模型是线性方程估计你，相对于MA和ARMA模型的非线性估计容易

8、，参数意义也便于解释。故实际建模时用高阶AR模型替换相应的MA和ARMA模型。综上考虑，可供选择的（p,q）组合有：（2，1），（3，0），（3，1），（4，0）。由于k=12时，样本的自相关和偏相关系数都不为0，所以，P=Q=1。3.模型的建立为了方便直接对原序列x进行预测，Eviews提供了差分算子d（x，n，s）=（1-B）n（1-BS）x表明序列x做n次一阶逐期差分和一次步长为s的季节差分后的新序列。采用菜单式建立ARIMA(2,1,1) (1,1,1)12模型。Eviews软件计算的结果如下：其中，sar(s)和sma(s)分别表示季节自回归部分和季节移动平均部分变量。表1 模型参数

9、估计与相关检验的结果Dependent Variable: D(LOG(X),1,12)Method: Least SquaresSample (adjusted): 1992M04 1997M12Included observations: 69 after adjustmentsConvergence achieved after 29 iterationsMA Backcast: 1991M03 1992M03VariableCoefficientStd. Errort-StatisticProb.AR(1)0.0.0.0.6174AR(2)-0.0.-0.0.8726SAR(12)0.

10、0.1.0.2767MA(1)-0.0.-2.0.0078SMA(12)-0.0.-25.660940.0000R-squared0.Mean dependent var-0.Adjusted R-squared0.S.D. dependent var0.S.E. of regression0.Akaike info criterion-3.Sum squared resid0.Schwarz criterion-3.Log likelihood118.9088Hannan-Quinn criter.-3.Durbin-Watson stat2.Inverted AR Roots.84.73+

11、.42i.73-.42i.42-.73i.42+.73i.07+.15i.07-.15i.00+.84i-.00-.84i-.42-.73i-.42+.73i-.73-.42i-.73+.42i-.84Inverted MA Roots.99.86+.50i.86-.50i.63.50+.86i.50-.86i-.00-.99i-.00+.99i-.50-.86i-.50+.86i-.86+.50i-.86-.50i-.99由表1可见，各滞后多项式的倒数根都在单位圆内，说明这个过程既是平稳的，也是可逆的。为了检验的预测效果，现在用ARIMA(2,1,1) (1,1,1)12模型对我国1997年

12、工业总产值进行预测，预测的结果如下：图7 预测值与真实值对比图图7 预测值与真实值对比图，预测精度MAPE为6. .同理可建立ARIMA(3,1,1) (1 1 1 )模型。 计算及预测的结果如下：表2 模型参数估计与相关检验的结果Dependent Variable: D(LOG(X),1,12)Method: Least SquaresSample (adjusted): 1992M05 1997M12Included observations: 68 after adjustmentsConvergence achieved after 52 iterationsMA Backcast:

13、 1991M04 1992M04VariableCoefficientStd. Errort-StatisticProb.AR(1)-0.0.-0.0.9285AR(2)-0.0.-0.0.5778AR(3)-0.0.-0.0.3677SAR(12)0.0.0.0.4305MA(1)-0.0.-1.0.1988SMA(12)-0.0.-25.198000.0000R-squared0.Mean dependent var-0.Adjusted R-squared0.S.D. dependent var0.S.E. of regression0.Akaike info criterion-3.S

14、um squared resid0.Schwarz criterion-3.Log likelihood117.0769Hannan-Quinn criter.-3.Durbin-Watson stat2.Inverted AR Roots.82.71+.41i.71-.41i.41+.71i.41-.71i.22-.52i.22+.52i.00-.82i-.00+.82i-.41+.71i-.41-.71i-.47-.71-.41i-.71+.41i-.82Inverted MA Roots.99.86+.50i.86-.50i.50+.86i.50-.86i.46-.00-.99i-.00

15、+.99i-.50-.86i-.50+.86i-.86+.50i-.86-.50i-.99预测 图8 预测值与真实值对比图用ls命令建立ARMA(4,1,0) (1 1 1 )模型。计算结果如下表：表3 模型参数估计与相关检验的结果Dependent Variable: D(LOG(X),1,12)Method: Least SquaresIncluded observations: 67 after adjustmentsConvergence achieved after 16 iterationsMA Backcast: 1991M06 1992M05VariableCoefficien

16、tStd. Errort-StatisticProb.AR(1)-0.0.-3.0.0004AR(2)-0.0.-2.0.0154AR(3)-0.0.-2.0.0287AR(4)-0.0.-1.0.2208SAR(12)0.0.0.0.5289MA(12)-0.0.-24.006180.0000R-squared0.Mean dependent var-0.Adjusted R-squared0.S.D. dependent var0.S.E. of regression0.Akaike info criterion-3.Sum squared resid0.Schwarz criterion

17、-3.Log likelihood114.7440Hannan-Quinn criter.-3.Durbin-Watson stat2.Inverted AR Roots.81.70+.40i.70-.40i.40-.70i.40+.70i.24+.65i.24-.65i.00+.81i-.00-.81i-.40+.70i-.40-.70i-.48-.30i-.48+.30i-.70-.40i-.70+.40i-.81Inverted MA Roots.99.86-.50i.86+.50i.50-.86i.50+.86i.00+.99i-.00-.99i-.50+.86i-.50-.86i-.

18、86+.50i-.86-.50i-.99由表3可见，各滞后多项式的倒数根都在单位圆内，说明这个过程既是平稳的，也是可逆的。为了检验的预测效果，现在用ARIMA(4,1,0) (1,1,1)12模型对我国1997年工业总产值进行预测，预测的结果如下：图9 预测值与真实值对比图图9 预测值与真实值对比图，预测精度MAPE为6.同理。我们建立ARIMA(3,1,0) (1,1,1)12模型并对其进行预测。Dependent Variable: D(LOG(X),1,12)Method: Least SquaresSample (adjusted): 1992M05 1997M12Included o

19、bservations: 68 after adjustmentsConvergence achieved after 16 iterationsMA Backcast: 1991M05 1992M04VariableCoefficientStd. Errort-StatisticProb.AR(1)-0.0.-3.0.0005AR(2)-0.0.-2.0.0366AR(3)-0.0.-1.0.0779SAR(12)0.0.1.0.3195MA(12)-0.0.-23.558810.0000R-squared0.Mean dependent var-0.Adjusted R-squared0.

20、S.D. dependent var0.S.E. of regression0.Akaike info criterion-3.Sum squared resid0.Schwarz criterion-3.Log likelihood115.9755Hannan-Quinn criter.-3.Durbin-Watson stat2.Inverted AR Roots.84.72+.42i.72-.42i.42-.72i.42+.72i.07+.60i.07-.60i.00-.84i-.00+.84i-.42+.72i-.42-.72i-.60-.72-.42i-.72+.42i-.84Inv

21、erted MA Roots.99.86+.50i.86-.50i.50-.86i.50+.86i.00-.99i-.00+.99i-.50-.86i-.50+.86i-.86-.50i-.86+.50i-.99计算的结果显示ARIMA(2,1,1) (1,1,1)12模型拟合的结果明显不如其他三个模型，故不予考虑。我们利用ARIMA(3,1,0) (1,1,1)12模型进行预测。结果如下：图10 预测值与真实值对比图图10 预测值与真实值对比图，预测精度MAPE为5.90187.4.模型的选择 三个模型的参数估计和相关检验汇总列入表4和表5.表4 各种模型的参数估计（）（）-0.327-0.

22、1119-0.1420_0.0923-0.4621-0.8981（）-0.4898-0.3413-0.3062-0.1564_0.0742-0.8979（）-0.4649-0.2840-0.2216_0.11750.8937表5 各模型检验的结果（）AICSCMAPE（）0.5540-3.3-3.076.18（）0.5489-3.25-3.056.03（）0.5438-3.3-3.15.90经计算，三个模型都满足平稳条件和可逆条件，模型设定合理。比较表5中各模型检验的结果。与俩个模型相比，第三个模型的AIC和SC值较小，预测的MAPE值显示其预测的精度最高。（MAPE的取值范围在0-5之间精度极高，在10以内说明预测精度高）。虽然调整后的样本决定系数略小于前两个模型，胆预测模型的选择应力求简洁、有效，因而选择第三个模型即ARIMA(3,1,0) (1,1,1)12模型比较合适。模型的检验：残差序列自相关的LM检验：我们建立的ARIMA(3,1,0) (1,1,1)1