实用回归分析小论文

1、影响成品钢材量的多元回归分析理论上成品钢材的需求量的影响因素主要有经济发展水平、收入水平、产业发展、人民生活水平提高、能源转换技术等因素。在此，收集的数据选择与其相关的四个因素：原油产量、生铁产量、原煤产量、发电量，19801997的有关数据如下表。原始数据（中国统计年鉴）年份y1980 2716.2010595.003802.406.203006.201981 2670.1010122.003416.606.203092.701982 2902.0010212.003551.006.663277.001983 3072.0010607.003738.007.153514.001984 337

2、2.0011461.304001.007.893770.001985 3693.0012489.504384.008.724107.001986 4058.0013068.805064.008.944495.001987 4356.0013414.005503.009.284973.001988 4689.0013704.605704.009.805452.001989 4859.0013764.105820.0010.545848.001990 5153.0013830.606238.0010.806212.001991 5638.0014009.206765.0010.876775.001

3、992 6697.0014209.707589.0011.167539.001993 7716.0014523.008739.0011.518395.001994 8482.0014608.209741.0012.409281.001995 8979.8015004.9410529.2713.6110070.301996 9338.0215733.3910722.5013.9710813.101997 9978.9316074.1411511.4113.7311355.53将中国成品钢材的需求量设为y，作为被解释变量，而原油产量、生铁产量、原煤产量、发电量作为解释变量，通过建立这些经济变量的线

4、性模型来研究影响成品钢材需求量的原因。一、 模型的设定设因变量y与自变量、的一般线性回归模型为:y = +是随机变量，通常满足；Var()=二 参数估计系数a模型非标准化系数标准系数tSig.B标准 误差试用版1(常量)170.287494.572.344.736原油（万吨)-.041.090-.031-.457.655生铁（万吨）.554.170.6083.267.006原煤（万吨）-17.818115.468-.018-.154.880发电量（亿千瓦时）.389.199.4381.952.073a. 因变量: 成品钢材（万吨）再用spss做回归线性，根据系数表得出回归方程为： 再做回归预测

5、，得出如下截图：故当原油产量为16225.86万吨，生铁产量为12044.54万吨，原煤产量为13.87万吨以及发电量为12334.89亿千瓦时时，成品钢材量预测值为10727.33875万吨；当原油产量为17453万吨，生铁产量为12445.96万吨，原煤产量为14.54万吨以及发电量为13457亿千瓦时时，成品钢材量预测值为10727.33875万吨。三 回归方程检验描述性统计量均值标准 偏差N成品钢材（万吨）5465.00282460.3492618原油（万吨)13190.63721875.7887318生铁（万吨）6489.95442700.7967618原煤（万吨）9.96832.5

6、401818发电量（亿千瓦时）6220.87942768.1119118相关性成品钢材（万吨）原油（万吨)生铁（万吨）原煤（万吨）发电量（亿千瓦时）Pearson 相关性成品钢材（万吨）1.000.909.998.961.997原油（万吨).9091.000.912.973.920生铁（万吨）.998.9121.000.962.997原煤（万吨）.961.973.9621.000.971发电量（亿千瓦时）.997.920.997.9711.000Sig. （单侧）成品钢材（万吨）.000.000.000.000原油（万吨).000.000.000.000生铁（万吨）.000.000.000.0

7、00原煤（万吨）.000.000.000.000发电量（亿千瓦时）.000.000.000.000.N成品钢材（万吨）1818181818原油（万吨)1818181818生铁（万吨）1818181818原煤（万吨）1818181818发电量（亿千瓦时）1818181818由相关系数表看出，因变量与各个自变量的相关系数都很高，都在0.9 以上，说明变量间的线性相关程度很高，适合做多元线性回归模型 。模型汇总b模型RR 方调整 R 方标准 估计的误差1.999a.997.997140.71641a. 预测变量: (常量), 发电量（亿千瓦时）, 原油（万吨), 原煤（万吨）, 生铁（万吨）。b.

8、因变量: 成品钢材（万吨）由=0.997以及调整之后的=0.997知，模型对样本观测数据的拟合度很好。Anovab模型平方和df均方FSig.1回归1.026E842.566E71296.001.000a残差257414.4041319801.108总计1.029E817a. 预测变量: (常量), 发电量（亿千瓦时）, 原油（万吨), 原煤（万吨）, 生铁（万吨）。b. 因变量: 成品钢材（万吨）由F=1296.001 P值=0.000<=0.05，故拒绝原假设，认为自变量联合起来对因变量有显著影响，通过F检验。系数a模型非标准化系数标准系数tSig.B标准 误差试用版1(常量)170

9、.287494.572.344.736原油（万吨)-.041.090-.031-.457.655生铁（万吨）.554.170.6083.267.006原煤（万吨）-17.818115.468-.018-.154.880发电量（亿千瓦时）.389.199.4381.952.073a. 因变量: 成品钢材（万吨） 因为=-0.154 P值=0.880远远大于=0.05，故接受原假设，认为对因变量y没有显著影响，故应剔除。用后退法剔除变量后，再做回归线性，得如下表：系数a模型非标准化系数标准系数tSig.B标准 误差试用版1(常量)170.287494.572.344.736原油（万吨)-.041.

10、090-.031-.457.655生铁（万吨）.554.170.6083.267.006原煤（万吨）-17.818115.468-.018-.154.880发电量（亿千瓦时）.389.199.4381.952.0732(常量)197.734445.099.444.664原油（万吨)-.053.045-.041-1.172.261生铁（万吨）.564.150.6203.760.002发电量（亿千瓦时）.371.153.4172.417.0303(常量)-309.403105.079-2.944.010生铁（万吨）.591.150.6493.937.001发电量（亿千瓦时）.311.147.350

11、2.125.051a. 因变量: 成品钢材（万吨）Anovad模型平方和df均方FSig.1回归1.026E842.566E71296.001.000a残差257414.4041319801.108总计1.029E8172回归1.026E833.422E71857.513.000b残差257885.8841418420.420总计1.029E8173回归1.026E825.131E72718.023.000c残差283174.3241518878.288总计1.029E817a. 预测变量: (常量), 发电量（亿千瓦时）, 原油（万吨), 原煤（万吨）, 生铁（万吨）。b. 预测变量: (常

12、量), 发电量（亿千瓦时）, 原油（万吨), 生铁（万吨）。c. 预测变量: (常量), 发电量（亿千瓦时）, 生铁（万吨）。d. 因变量: 成品钢材（万吨）最后剔除、两个自变量，得出新的回归方程为：F=2718.023 P值=0.000 故拒绝原假设，通过F检验。 四 基本假定检验1、异方差检验等级相关系数检验做abs(e)与x的等级相关系数，得出表如下相关系数原油（万吨)生铁（万吨）原煤（万吨）发电量（亿千瓦时）abseSpearman 的 rho原油（万吨)相关系数1.000.998*.994*.994*.263Sig.（双侧）.000.000.000.291N1818181818生铁（

13、万吨）相关系数.998*1.000.989*.988*.294Sig.（双侧）.000.000.000.236N1818181818原煤（万吨）相关系数.994*.989*1.000.997*.227Sig.（双侧）.000.000.000.365N1818181818发电量（亿千瓦时）相关系数.994*.988*.997*1.000.212Sig.（双侧）.000.000.000.399N1818181818abse相关系数.263.294.227.2121.000Sig.（双侧）.291.236.365.399.N1818181818*. 在置信度（双测）为 0.01 时，相关性是显著的。

14、由表中P值全大于0.01，故接受原假设，模型中不存在异方差。2、自相关检验模型汇总b模型RR 方调整 R 方标准 估计的误差Durbin-Watson1.999a.997.997140.71641.922a. 预测变量: (常量), 发电量（亿千瓦时）, 原油（万吨), 原煤（万吨）, 生铁（万吨）。b. 因变量: 成品钢材（万吨）DW=0.992 n=18 k=5,查表得出，故DW落入无法确定的领域。自相关性不明显，由此也看出DW检验的局限性。3、多重共线性系数a模型非标准化系数标准系数tSig.共线性统计量B标准 误差试用版容差VIF1(常量)170.287494.572.344.736原

15、油（万吨)-.041.090-.031-.457.655.04124.672生铁（万吨）.554.170.6083.267.006.006180.105原煤（万吨）-17.818115.468-.018-.154.880.01473.861发电量（亿千瓦时）.389.199.4381.952.073.004261.480a. 因变量: 成品钢材（万吨）共线性诊断a模型维数特征值条件索引方差比例(常量)原油（万吨)生铁（万吨）原煤（万吨）发电量（亿千瓦时）114.8781.000.00.00.00.00.002.1186.435.02.00.00.00.003.00437.106.36.03.0

16、3.09.004.00188.718.57.46.28.21.165.000137.100.05.51.68.70.84a. 因变量: 成品钢材（万吨）因为VIF中有两个远远的大于10，故模型存在严重的多重共线性。由共线性诊断表中数据得出，、之间存在多重共线性。故先剔除变量，再做线性回归，得出表如下：系数a模型非标准化系数标准系数tSig.共线性统计量B标准 误差试用版容差VIF1(常量)59.520538.378.111.914原油（万吨)-.108.092-.082-1.178.258.04721.142生铁（万吨）.871.054.95616.144.000.06615.176原煤（万吨

17、）118.233100.887.1221.172.261.02146.956a. 因变量: 成品钢材（万吨）共线性诊断a模型维数特征值条件索引方差比例(常量)原油（万吨)生铁（万吨）原煤（万吨）113.9171.000.00.00.00.002.0797.042.02.00.06.003.00333.686.33.04.66.134.00187.743.64.96.28.87a. 因变量: 成品钢材（万吨）由于模型中仍然存在多重共线性，故继续剔除VIF最大的变量，再做线性回归，得出表如下系数a模型非标准化系数标准系数tSig.共线性统计量B标准 误差试用版容差VIF1(常量)-282.1314

18、58.218-.616.547原油（万吨)-.017.049-.013-.345.735.1685.967生铁（万吨）.920.0341.01026.871.000.1685.967a. 因变量: 成品钢材（万吨）共线性诊断a模型维数特征值条件索引方差比例(常量)原油（万吨)生铁（万吨）112.9221.000.00.00.002.0766.204.03.00.193.00238.315.971.00.81a. 因变量: 成品钢材（万吨）系数a模型非标准化系数标准系数tSig.共线性统计量B标准 误差试用版容差VIF1(常量)-282.131458.218-.616.547原油（万吨)-.017.049-.013-.345.735.