中国能源消费影响因素分析(共21页)

1、精选优质文档-倾情为你奉上计量经济学课程论文学生姓名： 学 号： 院 部： 专 业： 班 级： 任课教师： 目录专心-专注-专业一、引言能源消费是引是指生产和生活所消耗的能源。能源消费按人平均的占有量是衡量一个国家经济发展和人民生活水平的重要标志。能源是支持经济增长的重要物质基础和生产要素。能源消费量的不断增长，是现代化建设的重要条件。我国能源工业的迅速发展和改革开放政策的实施，促使能源产品特别是石油作为一种国际性的特殊商品进入世界能源市场。随着国民经济的发展和人口的增长，我国能源的供需矛盾日益紧张。同时，煤炭、石油等常规能源的大量使用和核能的发展，又会造成环境的污染和生态平衡的破坏。可以看出

2、，它不仅是一个重大的技术、经济问题，而且以成为一个严重的政治问题。在20世纪的最后二十年里，中国国内生产总值（GDP）翻了两番，但是能源消费仅翻了一番，平均的能源消费弹性仅为0.5左右。然而自2002年进入新一轮的高速增长周期后，中国能源强度却不断上升，经济发展开始频频受到能源瓶颈问题的困扰。鉴于此，研究能源问题不仅具有必要性和紧迫性，更具有很大的现实意义。我国是一个能源大国，但是，我国人口众多，人均能源占有量不及同期发达国家的1/5。能源是任何一个国家经济发展不可缺失的物质基础。随着我国人口的继续增长，经济的快速发展，能源消费量的增加是必然的，而与年俱增的能源消费对环境造成的破坏也越来越严重

3、。因此，怎样优化能源利用结构，开发利用清洁能源，就成为我国经济发展的当务之急。这就需要我们清楚了解能源供需形势，做好影响能源消费因素分析，为能源规划及政策的制定提供科学依据，保证我国国民经济又好又快地发展。二、数据选取1、能源消费总量，在模型中用Y来表示。是指一次性能源消费总量，由煤炭、石油、天然气等组成（单位：万吨标准煤）。2、能源消费的影响因素：（1）能源生产总量，在模型中用来表示。是指一次性能源生产总量，该指标是观察全国能源生产水平、规模、构成和发展速度的总量指标（单位：万吨标准煤）。（2）全国生活能源消费总量，在模型中用来表示，是指一次性能源在在生活方面的消费量。（单位：万吨标准煤）。

4、（3）城镇居民人均可支配收入，在模型中用来表示。指城镇居民家庭人均可用于最终消费支出和其它非义务性支出以及储蓄的总和。它是家庭总收入扣除交纳的所得税、个人交纳的社会保障费以及调查户的记账补贴后的收入。（单位：元）。（4）工业能源消费总量，在模型中用来表示，是指工业方面的能源消费量。（单位：万吨标准煤）。（5）其他因素，在模型中用t表示。由于各种原因未考虑到和无法度量的因素归入随机扰动项，如能源价格变动、消费者偏好、国家的经济结构政策等。搜集到的数据见下附表1三、模型设定回归模型设定如下：Yt=0+1X1+2X2+3X3+4X4+t其中， Yt 表示能源消费总量 X1 表示能源生产总量 X2表示

5、全国生活能源消费总量 X3表示城镇居民人均可支配收入X4 表示工业能源消费总量 t 表示随机误差项0、1、2、3、4表示待定系数四、参数估计表1 回归结果Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 22:31Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C862.12012573.2160.0.7403X10.0.5.0.0000X2-0.0.-0.0.6843X

6、31.0.3.0.0042X40.0.7.0.0000R-squared0.    Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression1723.330    Akaike info criterion17.88859Sum squared resid    Schwarz criterion18.11988Log

7、 likelihood-272.2732    F-statistic14541.92Durbin-Watson stat1.    Prob(F-statistic)0.根据表1中数据，模型的估计结果为：Yt=862.1201+0.5166X1-0.1293X2+1.3019X3+0.6642X4 t = (0.3350) (5.0658) (-0.4112) (3.1375) (7.1934) R2=0. R2=0. F=14541.92 n=31五、模型检验5.1经济意义检验由回归估计结果可以看出，能源生产总

8、量、城镇居民人均可支配收入、工业能源消费总量与能源消费总量呈线性正相关，与现实经济意义理论相符。但是全国生活能源消费总量与能源消费总量呈线性负相关，与现实经济意义理论不相符。5.2统计检验1）拟合优度检验：由表1中数据可以得到R2=0.，修正的可决系数为R2=0.，这说明模型对样本的拟合很好。2）F检验：在95%的置信度下，F检验值P值小于0.05，回归方程是显著的。3）T检验：在95%的置信度下，X1，X3，X4 的t检验均值均小于0.05，表明线性作用显著，但X2的t检验均值不小于0,05，模型还需进一步完善。5.3计量经济学检验5.3.1多重共线性检验及其修正（一）相关系数检验：表2 相

9、关系数矩阵变量X1X2X3X4X1 1. 0. 0. 0.X2 0. 1. 0. 0.X3 0. 0. 1. 0.X4 0. 0. 0. 1.由相关系数矩阵可以看出，各解释变量相互之间的相关系数很高，证实确实存在严重多重共线性。（二）修正多重共线性采用逐步回归的办法，去检验和解决多重共线性问题，分别作Y对X1，X2，X3，X4的一元回归，结果如表3,表4,表5,表6所示。表3 Y 对 X1的一元回归估计结果Dependent Varia

10、ble: YMethod: Least SquaresDate: 06/08/13 Time: 22:40Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C-19219.851426.375-13.474610.0000X11.0.125.42780.0000R-squared0.    Mean dependent var.9Adjusted R-squared0.   &#

11、160;S.D. dependent var75900.69S.E. of regression3311.409    Akaike info criterion19.11048Sum squared resid3.18E+08    Schwarz criterion19.20299Log likelihood-294.2124    F-statistic15732.12Durbin-Watson stat0.    Prob(F-

12、statistic)0. 表4 Y 对X2的一元回归估计结果Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 22:41Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C-55231.908437.578-6.0.0000X211.230500.25.080140.0000R-squared0.    Mean dependent

13、var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression16206.49    Akaike info criterion22.28655Sum squared resid7.62E+09    Schwarz criterion22.37907Log likelihood-343.4416    F-statistic629.0133Durbin

14、-Watson stat0.    Prob(F-statistic)0.表5 Y对X3的一元回归估计结果Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 22:41Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C67031.063069.44621.838160.0000X313.633440.34.592720.0000R-s

15、quared0.    Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression11874.68    Akaike info criterion21.66455Sum squared resid4.09E+09    Schwarz criterion21.75706Log likelihood-333.8005&#

16、160;   F-statistic1196.656Durbin-Watson stat0.    Prob(F-statistic)0.表6 Y对X4的一元回归估计结果Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 22:42Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C5409.6261179

17、.6854.0.0001X41.0.132.96850.0000R-squared0.    Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression3123.934    Akaike info criterion18.99391Sum squared resid2.83E+08    Schwarz criteri

18、on19.08643Log likelihood-292.4057    F-statistic17680.61Durbin-Watson stat0.    Prob(F-statistic)0.整理表3,表4,表5,表6，结果如表7表7 一元回归估计结果变量X1X2X3X4参数估计值1.11.2305013.633441.t统计量125.427825.0801434.59272132.9685R20.0.0.0.R20.0.0.0.加入X4方程的R2=0.最大，所以以X4为基础，顺次加入其它变量逐步回归。表8

19、引入变量X4,X1Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 23:12Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C-6425.6122196.539-2.0.0068X40.0.6.0.0000X10.0.5.0.0000R-squared0.    Mean dependent var.9Adjusted R-squ

20、ared0.    S.D. dependent var75900.69S.E. of regression2143.538    Akaike info criterion18.27007Sum squared resid1.29E+08    Schwarz criterion18.40884Log likelihood-280.1861    F-statistic18793.04Durbin-Watson stat0. 

21、;   Prob(F-statistic)0.表9 引入X4，X2Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 23:12Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C4125.5342742.4481.0.1437X41.0.27.071540.0000X20.0.0.0.6070R-squared0.  

22、60; Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression3163.977    Akaike info criterion19.04881Sum squared resid2.80E+08    Schwarz criterion19.18759Log likelihood-292.2566    F

23、-statistic8618.092Durbin-Watson stat0.    Prob(F-statistic)0.表10 引入X4，X3Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 23:13Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C13302.412335.4925.0.0000X41.0.24.021160.

24、0000X31.0.3.0.0009R-squared0.    Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression2599.786    Akaike info criterion18.65601Sum squared resid1.89E+08    Schwarz criterion18.79478Log

25、likelihood-286.1682    F-statistic12771.20Durbin-Watson stat0.    Prob(F-statistic)0.表11 引入变量回归结果变量X1X2X3X4R2X4、X10.（5.）0.（6.）0.X4、X20.（0.）1.（27.07154）0.X4、X31.（3.）1.（24.02116）0.经比较，加入X1的方程R2=0.，改进最大，而且各参数的t检验显著，选择保留X4，再加入其它新变量逐步回归。表12 引入X4、X1、X2Dependent Vari

26、able: YMethod: Least SquaresDate: 06/08/13 Time: 23:21Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C-4641.6092169.153-2.0.0416X40.0.6.0.0000X10.0.6.0.0000X2-0.0.-2.0.0250R-squared0.    Mean dependent var.9Adjusted R-squared0.&

27、#160;   S.D. dependent var75900.69S.E. of regression1985.610    Akaike info criterion18.14515Sum squared resid1.06E+08    Schwarz criterion18.33018Log likelihood-277.2499    F-statistic14602.79Durbin-Watson stat1.  

28、  Prob(F-statistic)0. 表13 引入X4、X1、X3Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 23:22Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C1061.2542488.0340.0.6731X40.0.7.0.0000X10.0.6.0.0000X31.0.4.0.0003R-squared0. 

29、60;  Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900.69S.E. of regression1696.604    Akaike info criterion17.83056Sum squared resid    Schwarz criterion18.01559Log likelihood-272.3737    F-s

30、tatistic20004.84Durbin-Watson stat1.    Prob(F-statistic)0.表14 引入变量回归结果变量X1X2X3X4R2X4、X1、X20.(6.)-0.(-2.)0.(6.)0.X4、X1、X30.(6.)1.(4.)0.(7.)0.经比较，在X4、X1基础上加入X3后不仅使R2增大，而且t检验值也通过，所以选择保留X3，继续回归。表15 引入X4、X1、X3、X2Dependent Variable: YMethod: Least SquaresDate: 06/08/13 Time: 23:38Samp

31、le: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C862.12012573.2160.0.7403X40.0.7.0.0000X10.0.5.0.0000X31.0.3.0.0042X2-0.0.-0.0.6843R-squared0.    Mean dependent var.9Adjusted R-squared0.    S.D. dependent var75900

32、.69S.E. of regression1723.330    Akaike info criterion17.88859Sum squared resid    Schwarz criterion18.11988Log likelihood-272.2732    F-statistic14541.92Durbin-Watson stat1.    Prob(F-statistic)0.经比较，在X4、X1、X3的基础上，加入X2后

33、，不仅R2下降，而且X2参数的t检验值为负。这说明X2引起多重共线性，应予剔除。最后修正多重共线性影响的回归模型为：Yt=0+1X1+3X3+4X4+t 5.3.2异方差检验（一）图示法（1） X1、X3、X4对Y散点图图1 X1对Y散点图图2 X3对Y散点图图3 X4对Y散点图由图1，图2，图3可知，随着X1、X3、X4变量值的增加，Y值也逐渐增加，但是Y值的离散程度没有较明显的变化趋势，所以可能不存在异方差性。但是否确实存在异方差还应通过进一步的检验。（2） e2对X1、X3、X4的散点图图4 E2对x1 的散点图图5 E2 对 X3 的散点图图6 E2 对 X4 的散点图由图4,、图5、

34、图6可以看出，残差平方ei2对解释变量X1、X3、X4的散点图都分布在各个地方，大致可以看出残差平方和ei2不随X1、X3、X4的变化而变化，因此，模型很可能不存在异方差。但是否确实存在异方差还应通过进一步的检验。（二）怀特（White）检验表16 White检验结果White Heteroskedasticity Test:F-statistic2.    Probability0.Obs*R-squared16.19218    Probability0.Test Equation:Dependent Vari

35、able: RESID2Method: Least SquaresDate: 06/09/13 Time: 00:11Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C0.0.5516X1-2274.8492468.925-0.0.3673X120.0.1.0.1359X1*X30.0.1.0.3245X1*X4-0.0.-1.0.0708X3286.190311646.260.0.9806X320.0.0.0.5247X3*X4-0.0.-1.

36、0.0936X42515.1882429.1311.0.3122X420.0.2.0.0276R-squared0.    Mean dependent var.Adjusted R-squared0.    S.D. dependent var.S.E. of regression.    Akaike info criterion32.88795Sum squared resid1.83E+14    Schwarz criteri

37、on33.35053Log likelihood-499.7632    F-statistic2.Durbin-Watson stat2.    Prob(F-statistic)0.从表四可以看出，nR2=16.19218，由White检验知，在 =0.05 下，查表知0.0529=16.92>16.19，所以拒绝备择假设，不拒绝原假设，表明模型不存在异方差。5.3.3自相关检验及其修正（一）自相关检验表17 消除多重共线性后的最小二乘法估计结果Dependent Variable: YMethod: L

38、east SquaresDate: 06/13/13 Time: 01:04Sample: 1980 2010Included observations: 31VariableCoefficientStd. Errort-StatisticProb.  C1061.2572488.0320.0.6731X10.0.6.0.0000X31.0.4.0.0003X40.0.7.0.0000R-squared0.    Mean dependent var.9Adjusted R-squared0.   

39、60;S.D. dependent var75900.69S.E. of regression1696.603    Akaike info criterion17.83056Sum squared resid    Schwarz criterion18.01559Log likelihood-272.3736    F-statistic20004.87Durbin-Watson stat1.    Prob(F-statistic

40、)0.Yt=1061.257+0.4907X1+1.4008X3+0.6694X4 t = (0.4265) (6.2247) (4.2065) (7.4337)R2=0.99955 R2=0.9995 F=20004.87 DW=1.1261该回归方程可决系数较高，回归系数均显著。对样本量为31,3个解释变量的模型，5%显著水平，查DW统计表可知，dL=1.229 ,dU=1.650,模型中DW < dL ,显然模型中有自相关。（二）自相关问题处理表18 消除自相关处理结果Dependent Variable: YMethod: Least SquaresDate: 06/09/13

41、Time: 00:33Sample (adjusted): 1981 2010Included observations: 30 after adjustmentsConvergence achieved after 10 iterationsVariableCoefficientStd. Errort-StatisticProb.  C-99.334553181.945-0.0.9753X10.0.5.0.0000X30.0.2.0.0547X40.0.7.0.0000AR(1)0.0.3.0.0059R-squared0.    

42、Mean dependent var.5Adjusted R-squared0.    S.D. dependent var75587.03S.E. of regression1553.167    Akaike info criterion17.68499Sum squared resid    Schwarz criterion17.91852Log likelihood-260.2749    F-statistic17164.7

43、8Durbin-Watson stat1.    Prob(F-statistic)0.Inverted AR Roots      .56Yt=-99.3346+0.453X1+0.8781X3+0.7597X4 t = (-0.0312) (5.1551) (2.0155) (7.4222) R2=0. R2=0. F=17164.78 DW=1.7279 对样本量为31,3个解释变量的模型，5%显著水平，查DW统计表可知，dL=1.229 ,dU=1.650,模型中DW =1. > dL ,显然模型已消除自相关。六、结论及建议6.1结论1、 能源消费模型的多重共线性修正和自相关补救表明，能源消费总量与能源生产总量、城镇居民人均可支配收入及工业能源消费量存在着长期均衡的关系。2、 在多重共线性的修正过程中，可以发现，时间序列全国能源消费总量、工业能源消费量与能源消费总量具有共同变化趋势，在经济上升时期均呈现增长的趋势；在经济收缩期，又都呈现下降趋势。当这三者同时作为解释变量时，就很有可能出现多重共线性。出现多重共线性的另一原因是：抽样仅仅局限于能源消费总量影响因素的一个有限范围内。3、 在自相关的修正过程中，我们可以发现，全国生活能源消费总量、城镇居民人均可