关于多元线性回归的毕业论文

河北工程大学本科毕业设计（论文）

摘 要

许多现象往往不是简单的与某一因素有关而是要受多个因素的影响，此时就需要用两个或两个以上的影响因素作为自变量来解释因变量的变化，这就是多元回归亦称多重回归。当多个自变量与因变量之间是线性关系时，所进行的回归分析就是多元性回归。

本文的研究主要从四个部分来进行。第一章从基础内容和研究对象着手，对主要研究内容进行了简单的阐述。第二章对多元线性回归的基础进行了详细分析。第三章介绍了中国经济的现状。最后通过多元线性回归模型对我国工业生产总值进行了分析。

总的来说，本文在2007年全国各省市主要工业产品的产量与工业总产值的具体数据下，选用塑料、水泥、钢筋、平板玻璃、粗钢、盘条以及原煤等工业产品的产量作为研究对象，建立多元线性回归模型，并对模型做出参数估计.在此基础上对模型做出一定的解释，对于预测工业总产值具有一定的理论指导和现实意义。

关键词：多元线性回归模型 工业生产总值 假设检验 预测

1

河北工程大学本科毕业设计（论文）

Abstract

Many phenomena are often not simply associated with a number of factors but with

varieties. At this point we need to use two or more factors as independent variables to explain changes in the dependent variable. This is also known as multiple regression. When more than one independent variable and the dependent variable are linear relationship, the regression analysis is carried out by diversity regression.

The main research work of this thesis is divided into four parts. In the first chapter, the thesis proceed from the basic content and object of study and elaborate main content simply. In the second chapter, multiple linear regression model is analyzed detail. In the third chapter, the thesis introduces status quo of china. And at last, gross industrial production is analyzed by multiple linear regression model in this article.

Over all, this article use the specific data of the output of major industrial products and industrial output in nationwide provinces in 2007, and select the output of plastics, cement, steel, plate glass, crude steel, wire rod and raw coal as study object to establish multiple linear regression model, and then make the model parameter estimation. Based on this，we make some explanations to the model. All of these are of momentous current significance and far-reaching historical significance to the forecast of industrial production.

Key Words: Multiple linear regression model Gross industrial production Hypothetical

test Prediction

2

河北工程大学本科毕业设计（论文）

目 录

摘要.............................................................................................................................................1 Abstract.......................................................................................................................................2 1 绪论 ........................................................................................................................................ 4 2 多元线性回归分析基础 ........................................................................................................ 5

2.1 多元线性回归定义 ...................................................................................................... 5 2.2多元线性回归模型........................................................................................................6 2.2.1模型的建立及矩阵表示.....................................................................................6 2.2.2模型的假设 ........................................................................................................ 7 2.3 多元线性回归参数估计 .............................................................................................. 7

2.3.1 最小二乘估计和正规方程组 ............................................................................ 7 2.3.2 最小二乘估计的矩阵形式 ................................................................................ 8 2.4 回归拟合度评价和决定系数 ...................................................................................... 9

2.4.1 离差分解和决定系数 ........................................................................................ 9 2.4.2 决定系数的性质及修正可决系数 .................................................................. 10 2.5 统计检验 .................................................................................................................... 11

2.5.1回归参数的显著性检验（t检验） ................................................................. 11 2.5.2回归方程的显著性检验（F检验） ................................................................ 12 2.5.3 多重共线性检验 .............................................................................................. 12

2.5.4 异方差检验.......................................................................................................13 3 中国经济现状 ...................................................................................................................... 15

3.1中国经济现状 ............................................................................................................. 15 3.2 工业生产总值的概述 ................................................................................................ 15 4 工业生产总值的多因素模型分析 ...................................................................................... 15

4.1建立多因素分析模型 ................................................................................................. 16 4.2数据收集 ..................................................................................................................... 16 4.3 统计检验 .................................................................................................................... 19 4.4 计量经济学检验及模型修正 .................................................................................... 20

4.4.1 异方差检验 ...................................................................................................... 21 4.4.2 自相关检验 ...................................................................................................... 21

5 结论.......................................................................................................................................26 致 谢...................................................................................................................................27 参 考 文 献 ............................................................................................................................ 28

3

河北工程大学本科毕业设计（论文）

1绪 论

在各个方面，变量之间的关系一般来说可分为确定性的与非确定性的两种。确定性关系是指变量之间的关系可以用函数关系来表达的。另一种非确定性的即所谓的相关关系。例如人的身高与体重之间存在着关系，一般来说，人高一些，体重也要重一些，但同样高度的人，体重往往不相同。人的血压与年龄之间也存在着关系，但同年龄的人的血压往往不相同。气象中的温度与湿度之间的关系也是这样的。这是因为我们涉及的变量（如体重、血压、适度）是随机变量，上面所说的变量关系是非确定性的。此时 ，便可以用到回归分析。回归分析能帮助我们从一个变量取得的值去估计另一个变量所取的值。

工业生产总值从数值上反应一个地区的工业生产规模，是衡量一个地区的经济繁荣程度的重要指标。研究研究工业总产值与格工业产出指标之间的关系具有非常重要的现实意义，对于做好一个地区的的工业产值预测以及制定国民经济发展规划都有的非常重要的作用。

工业总产值 是指以货币表现的工业企业在一定时期内生产的已出售或可供出售的工业的产品的总量。它是反映一定时间内工业生产总规模和，总水平的重要指标，是计算工业生产发展速度和主要比例关系，计算工业产品销售率和其他经济指标的重要依据。工业总产值包括成品价值、工业性作业价值和自制半成品、在产品期末期初差额价值。工业，总产值采用“工厂法”计算，即以工业企业作为一个整体，按企业工业生产活动的最终成果来计算。但各企业之间、行业之间、地区之间存在着重复计算。其计算公式为：报告期工业总产值=报告期全部产品的成品价值+报告期工业性作业价值+（报告期自制半成品和在产品期末余额- 报告期自制半成品和在产品期初余额） 计算工业总产值采用的价格有不变价格和现行价格。

即，工业生产总值收多个因素影响，此时便需要多个影响因素来分析工业生产总值的变化。而这些变量之间的关系是线性的，这样在分析工业生产总值是用到的回归分析方法便是多元线性回归。

4

河北工程大学本科毕业设计（论文）

2 多元线性回归分析基础

2.1多元线性回归定义

在客观世界中普遍存在着变量之间的关系。变量之间的关系一般来说可分为确定性的与非确定性的两种。确定性关系是指变量之间的关系可以用函数关系来表达的。另一种非确定性的即所谓的相关关系。例如人的身高与体重之间存在着关系，一般来说，人高一些，体重也要重一些，但同样高度的人，体重往往不相同。人的血压与年龄之间也存在着关系，但同年龄的人的血压往往不相同。气象中的温度与湿度之间的关系也是这样的。这是因为我们涉及的变量（如体重、血压、适度）是随机变量，上面所说的变量关系是非确定性的。此时 ，便可以用到回归分析。回归分析能帮助我们从一个变量取得的值去估计另一个变量所取的值。

在回归分析中，如果有两个或两个以上的自变量，就称为多元回归。事实上，一种现象常常是与多个因素相联系的，由多个自变量的最优组合共同来预测或估计因变量，比只用一个自变量进行预测或估计更有效，更符合实际。因此多元线性回归比一元线性回归的实用意义更大。

在研究问题是，我们考虑一个变量受其他变量的影响时，把这变量称为因变量，记为Y，其他变量称为自变量，记为X，这时相关系数可记作

Y?f?x???, 其中f?x?为当X?x时，因变量Y的均值，即

f?x??E?Y|X?x?.

称f?x?为Y对X的回归函数，?为Y与f?x?的偏差，它是随机变量，并假定E????0。 回归函数可以是一元函数，也可以是多元函数，即

Y?f(x1,x2,?,xm)??,

其中 f(x1,x2,?,xm)?E(Y|X1?x1,X2?x2,?,Xm?xm)为m元回归函数，统称为多元回归函数。

2.2多元线性回归模型

2.2.1 模型的建立及矩阵表示 多元线性回归模型的一般形式是： Y??0??1Z1??2Z

2??Z33??Z4??Z5 ?5 ? 4 (2.1)

j?1,2,...,k）其中?（是回归系数，Y是被解释变量，z1i,z2i?，zki是k个对Y有显j著影响的解释变量（k?2），?i是 反映各种误差扰动综合影响的随机项，下标i表示第i期观察值（Yi,z1i,z2i?,zki）, i?1,2,?n。

???假设多元样本回归函数为：Yi??0??1z1i??2z2i???kzki回归残差为：?i?Yi?Yi。

5

???

河北工程大学本科毕业设计（论文）

由于有n期的观察值，这一模型实际上包含n个方程

Y1??0??1z11??1z11????kzk1??1

Y2??0??1z12??2z22????kzk2??2

?? Y2??0??1z1n??2z2n????kzkn??n

写成矩阵形式：

Y?Z???, (2.2) 其中

?Y1??1 z11 z21 ?zk1????Y2??1 z12 z22 ?zk2Y???,Z???? ? ? ?????1 z z ? z?Y?1n2nkn??n????? , ????????0???0???0???????????1???1???1?????, ?? ,????.

?? ? ?????? ????????????k??n?????k?

2.2.2 模型的假设

因为多元线性模型的建立或选择过程包含相当的主观性，所依据的理论和经验也可能不正确，因此并不能保证模型符合变量的实际关系。而如果模型本身有问题，那么分析的有效性和价值就很难有保证，为了保证所分析的变量关系符合多元线性回归分析的基本规定性，明确分析对象，保证回归分析的有效性和性质，也为了检验判断的依据，需要对多元线性回归模型作一些架设，共包括下列六条：

(1)变量Yi和X1i,X2i?,Xki,(i=1,2...n)之间，存在线性随机函数关系(2)对应每组观测数据的误差项?i，都为零均值的随机变量，即?i的数学期望(3)误差项?i的方差为常数，即Var(?i)?E??i?E(?i)??E(?i)??2 对i=1,2...n 都

22Yi??0??1X1i??2X2i????kXki??i，其中?i是随机误差项。

E(?i)=0对i=1,2...n都成立。 成立（假设（2）成立为前提）。

（4）对应不同观测数据的误差项不相关，即

Cov(?i,?j)?E(?i?E(?i))(?j?E(?j))?E(?i???j?0)对任意的 i?j都成立（假设（1）

成立为前提）。

6

河北工程大学本科毕业设计（论文）

（5）解释变量Xi(i?1,?,r)是确定性变量而非随机变量。当存在多个解释变量（r>1）时假设不同解释变量之间不存在线性关系，包括严格的线性关系和强的近似线性关系。 （6）误差项?i服从正态分布[7]。

2.3 多元线性回归参数估计

2.3.1 最小二乘估计和正规方程组

这里直接根据回归残差平方和最小的准则，推导多元线性回归模型参数的最小二乘估计量。对于多元线性回归模型Y??0??1z1????kzk??，

Y?b0?b1z1???bkzk回归残差平方和为：

?如果用b0，?，bk分别表示模型参数?0，?,?k的估计，那么样本回归方程就是

V???i2i??[Yii-(b0?b1z1i???bkzki)]2 (2.3)

当V对b0，?，bk的一阶偏导数都等于0，即下列方程组： ?2[Yi-(b0?b1z1i???bkzki)](-1)?0,

?2[Yi-(b0?b1z1i???bkzki)](-z1i)?0，

ii??

?2[Yi-(b0?b1z1i???bkzki)](-zki)?0，

同时成立时,V有最小值。对这个方程组整理，可得到如下的正规方程组：

b0?Y -(b1z1???bkzk) ,

S11b1?S12b2???S1KbK?S10,

i?

SK1b1?SK2b2???SKKbK?SK10,

其中

Sk0??(ziki-zk) (Yi-Y),k,i?1,?K，

Skj???

?(ziki-zk) (zji-zj),k,j?1,?K.

bk上述正规方程组有K+1个方程，未知数也是K+1个。只要系数矩阵非奇异即满足 的唯一的一组解，就是?0，?，?k的最小二乘估计[8]。

解释变量矩阵Z列满秩：R(Z)?k。此时，有R(Z'Z)?k，Z'Z可逆。可以解出b0，?，

2.3.2 最小二乘估计的矩阵形式

引进参数估计量，解释变量回归值和回归残差的下列向量表示：

7

河北工程大学本科毕业设计（论文）

?b0?b1?B?

????bK????Y??1?1??????，??? Y????，? ??. (2.4)

????????Y?n??n??把样本数据分别带入样本回归方程，得到回归方程组为：

Y1?b0?b1z11???bkzk1, (2.5)

? ?

Yn?b0?b1z1n???bkzkn??

写成等价的向量方程，则为：Y?ZB.

这样回归残差向量为：??Y-Y?Y-XB.

在利用向量，矩阵的运算法则，可以得到残差平方和为

V????i2????(Y?XB)(Y?XB)=YY-BXY-YXB?BXXB.

''''''''求V对b0，?，bk的偏导数，等价于V对向量B求梯度，因此最小二乘估计的正规方程

??V???b??0?组为：?BV??????2Z'?2Z'ZB?0,整理得到矩阵 形式：Z'ZB?ZY.

???V?????bn??当X?X可逆，也就是X是满秩矩阵，在上述向量方程两端左乘X?X的逆矩阵，得到：

B?(ZZ)ZY, (2.6)

'-1'这就是多元线性回归模型最小二乘估计的矩阵一般形式。

2.3.3 最小二乘估计量的性质

(1)线性性：

多元线性回归模型参数的最小二乘估计向量为：B?(Z'Z)-1Z'Y，各个参数的最小二乘估计向量为bk??(z'z)-1z'?????k?1Y，其中的?(z'z)-1z'?????k?1是矩阵(z'z)-1z'的k+1

行元素构成的行向量，上式对k=1，?，K都成立，bk正是被解释变量观测值Yi的线性组合，也就是多元线性回归参数的最小二乘估计是线性估计。

(2)无偏性：

多元线性回归的最小二乘估计也是无偏估计，即参数最小二乘估计量的数学期望都

等于相应参数的真实值，最小二乘估计向量的数学期望等于参数真实值的向量，参数真

8

河北工程大学本科毕业设计（论文）

实值是参数估计量的概率分布中心。

? ?E?(ZE(B)＝E(ZZ）ZY?E(ZZ)X(X???)'-1''-1'??'-1'?

'Z)X(X???)＝E??(ZZ)ZE(?)??'-1?

???(Z'Z)?1Z'E(?)??.

(3)最小方差性：

根据最小二乘估计公式和模型假设，可以直接导出包含各个参数估计量方差和不同

参数估计量协方差的，参数估计向量B的协方差矩阵为：

B)?Var［ZZ Var(?'?-1ZY］?Var[ZZ'?'?'-1Z(Z???)]

' ?Var[??ZZ ?ZZ?'?-1Z?]?Var[ZZ'??-1?'?-1ZVar[?][ZZ'?'?-1Z]?ZZ'?'?Z?]

-1'2'Z?I［ZZ?'?-1Z]

' ??z'z??2 (2.7)

?1

2.4 回归拟合度评价和决定系数

2.4.1 离差分解和决定系数

判断回归结果好坏基本标准，是回归直线对样本数据的逆合程度，称为“拟合度”。

回归直线的逆合度一方面取决于回归直线的选择，这就是由参数估计方法决定的，另一方面则取决于样本数据的分布。当参数估计方法固定时，主要取决于样本数据的分布。

样本数据的分布在本质上是由变量关系决定的。因此回归拟合度也是检验模型变量关系真实性，判断模型假设是否成立的重要方法。拟合度较好是对模型的支持，否则，可能意味着必须对模型进行修改。

首先需要从Y的离差中分离出由解释变量决定的部分，因变量的实际观测值与其样本均值的离差即总离差（Y?Y）可以分解为两部分：一部分是因变量的理论回归值与

??Y）其样本均值的离差（Y， 它可以看成是能够由回归直线解释的部分，称为可解释?）离差；另一部分是实际观测值与理论回归值的离差（Y?Y，它是不能由回归直线加以

解释的残差e。 对任一实际观测值Y总有：

Y?Y?(Y?Y?)?(Y??Y) (2.8)

_???Yi?Y????2????Yi?Yi????2对公式（2.8）两边平方并求和并计算，可得到：

?????_??Yi?Y???

?2.9? SST?SSR?SSE

根据最小二乘估计和回归残差的相关公式，所有Yi的离差的平方和记为

?2?2?SST =?(Yj?Yj)称为“总离差平方和”，而?(Yj?Yj)记为SSR称为“残差平方和”，

?(Yj?Yj)记为SSE称为“回归平方和”。

9

2 河北工程大学本科毕业设计（论文）

式（2.9）两边同除以SST ，得：

SSRSST?SSESST 1?

(2.10)

显而易见，各个样本观测点与样本回归直线靠的越近，SSR在SST 中所占的比重就

SSRSST越大。（2.10）式中的正是反映解释变量（或回归直线）对被解释变量决定程度的

指标，我们称它为“决定系数”（determined coefficient），通常用R2表示。计算公式为：

R2?

2.4.2 决定系数的性质及修正可决系数

决定系数是对回归模型拟合程度的综合度量，决定系数越大，模型拟合程度越高。决定系数越小，则模型对样本的拟合程度越差。决定系数R2具有如下性质：

(1) 决定系数R2具有非负性。

由决定系数的定义式可知，R2的分子分母均是不可能为负值的平方和，因此其比 值必大于零。

(2) 判定系数的取值范围为0?R2?1。

由R2的计算公式可以看出：当所有的观测值都位于回归直线上时，SSR=0，这时

2R=1，说明总离差可以完全由所估计的样本回归直线来解释；当观测值并不是全部位

SSRSST?1?SSESST.

于回归直线上时， SST>0，则SSR／SST>0，这时R2<1；当回归直线没有解释任何离差，即模型中解释变量Y与因变量Y完全无关时，Y的总离差全部归于残差平方和，即

SSR＝SST ，这时R2＝0。

(3) 判定系数是样本观测值的函数，它也是一个统计量。

判定系数R2的大小受到自变量X的个数k的影响。可以证明，增加自变量X的个数，回归平方和增大，从而使得R2增大。由于增加自变量个数引起的R2增大与拟合好坏无关，在含自变量个数k 不同的模型之间比较拟合程度时，R2就不是一个合适的指标，必须加以调整。

调整方法为：把残差平方和与总离差平方和之比的分子分母分别除以各自的自由度，变成均方差之比，以剔除自变量个数对拟合优度的影响。调整的判定系数为： R2?1-n-1n-K-1(1-R)?1-2n-1n-K-1?ε*iii2i2?(Y-Y) (2.11)

10

河北工程大学本科毕业设计（论文）

用这个调整的决定系数作为评价多元线性回归拟合度的评价标准，可以基本消除由

于解释变量数目的差异所造成的影响，更加合理和具有可比性。 2.5 统计检验

2.5.1回归参数的显著性检验（t检验）

先要找出回归系数的分布，由上述知识得知：

?j~N(?j,?cjj), (2.12)

'-1?2其中cjj为?ZZ?的第j行j列的元素。将?j标准化。一般有?未知，用?代替，得统

2??2?计量 t??j??j?~t(n?k)，以下可用t统计量来进行回归系数的假设检验。

?cjj同一元线性回归一样，要检验解释变量Zj对因变量Y的线性作用是否显著，要使用

t检验。步骤如下：

(1) 提出假设。

H0: ?j?0,j?1,2,?,k H1: ?j?0,j?1,2,?,k

(2) 在H0成立条件下，根据样本计算

?? t??j??j???j?

?cjj?cjj(3) 给定显著性水平?，查表得临界值t?/2(n?k) (4) 判断

若|t|?t?/2(n?k)，就拒绝H0，Zj对Y有显著线性作用； 若|t|?t?/2(n?k)，就接受H0，Zj对Y线性作用不显著。

2.5.2回归方程的显著性检验（F检验）

多元线性回归模型还可以进行模型总体显著性检验，也就是全体解释变量总体对被解释变量是否存在明显影响的检验，回归显著性检验的基本方法，是检验模型常数项以外所有参数同时为0的假设，使用F检验。步骤如下： (1) 提出假设。

H0:?2??3????k?0 H1:?2,?3,?,?k不全为0

(2) 选择、（根据样本）计算统计量

11

河北工程大学本科毕业设计（论文）

F?ESS/(k?1)RSS/(n?k)~F(k?1,n?k)

(3) 给定显著性水平?，查表，得F?(k?1,n?k) (4) 判断

若F?F?(k?1,n?k)，就拒绝H0，回归方程显著成立，所有自变量对Y 的影响是显著的；

若F?F?(k?1,n?k)，就接受H0，回归方程不显著，所有自变量对Y 的线性作用不显著。

2.5.3 多重共线性检验

在多元线性回归模型Y??Z中，对Z的基本假定是：矩阵的各列向量之间是线性无关的，即有：r(Z)?k(k?n),即(Z'Z)?0如果这一假定不满足，则称模型存在多重共线性。多重共线性表现为两种情况：

(1) 完全多重共线性：r(Z)?k，也就是(Z'Z)?0，(Z'Z)-1不存在。

(2) 不完全多重共线性：（实际中多为此情况）(Z'Z)?0，(Z'Z)-1对角线元素较大。而一般产生多重共线性的背景为：

（1）时间序列数据中经济变量在时间上常有共同的变动趋势； （2）经济变量之间本身具有内在联系（常在截面数据中出现）； （3）由于某种决定性因素的影响可能使各个变量向着同方向变化； （4）滞后变量引入模型，同一变量的逐次值一般都存在相互关系； 多重共线性的检验方法有：

（1）简单相关系数矩阵法（辅助手段）

此法简单易行；但要注意两变量的简单相关系数包含了其他变量的影响，并非它们真实的线性相关程度的反映；一般在0.8以上可初步判定它俩之间有线性相关。 （2）变量显著性与方程显著性综合判断；

（修正）可决系数大，F值显著大于临界值，而t值不显著；那么可认为存在多重共线性。

（3）辅助回归：

将每个解释变量对其余变量回归，若某个回归方程显著成立，则该解释变量和其余变量有多重共线性。

多重共线性的克服和处理方法有：

截面数据和时序数据结合，有时在时间序列数据中多重共线性严重的变量，在截面数据中不一定有严重的共线性。在假定截面数据估计出的参数在时间序列数据中变

12

河北工程大学本科毕业设计（论文）

化不大的前提下，可先用截面数据估计出一些变量的参数，再代入原模型估计另一些变量的参数。

变换模型形式（差分法）: 假设Z2和Z3存在高度线性相关。 设原模型为：

Yt??1??2z2. t??3z3. t??t.

将其滞后一期：

Yt-1??1??2z2. t-1??3z3. t-1??t-1.

将上述两式相减，得：

Yt-Yt-1??2(z2. t-z2. t-1)??3(z3. tz3. t-1)??t-?t-1

令?Yt?Yt-Yt-1, ?z2. t?z2. t-z2. t-1, ?z3. t?z3. t-z3. t-1, ??t??t-?t-1

则上述差分式子变成：

?Yt??2?z2. t??3?z3. t???t

差分后，?z2和?z3的共线性将明显减弱。

2.5.4 异方差检验

在回归模型的假设得到满足之后，用最小二乘法估计的模型参数具有无偏和方差在线性无偏估计方法中最小的有效性，在这些假设中，其中有一条是误差项的方差不变。如果误差项的方差随观测次数的改变而改变，或随解释变量增减而变化，则称回归模型中存在异方差。异方差可以表示为Var??i???i2或

??12?2?2?????????????????2?n????Var????E????'

其中异方差的的发现和检验方法有戈德菲尔德-夸特检验：构造统计量：

F??i2i1ei22i12?e?n?c??K?1???2??n?c???K?1???2??ei2i12i22i1?e.

如果F?F?，误差项存在明显的递增异方差性； 如果1?F?F?，误差项没有明显的异方差性。

异方差的克服和处理：如线性回归模型为Yi??0??1z1i????KzKi??i，经检验，

13

河北工程大学本科毕业设计（论文）

误差项有如下异方差性?i2?f(zji)?2，可以用

Yif(zji)??01f(zji)??1z1if(zji)f(zji)除模型各项，得到：

zkif(zji)?????K?if(zji),

新模型的误差项方差为：

Var[?if(zji)]?[1f(zji)1f(zji)]Var[?i]

2???i2

?f(zji)?

2 ?1f(zji)2 ??.

14

河北工程大学本科毕业设计（论文）

3 中国经济现状

3.1 中国经济现状

改革开放30年来，中国经济持续高速增长，相当程度上是依赖于中小企业的崛起。快速、健康和持续发展的中小企业，对经济增长的贡献有目共睹：在繁荣经济、促进增长、国际贸易、扩大就业、推动创新、提高消费能力等方面发挥着重要的作用，已成为推动我国经济社会发展的重要力量，是大企业发展的依托，是活跃市场的基本主体，也是经济活力的具体体现[1]。

回顾2008年中国：我们经历了年初的雪灾、5月的地震灾害、8月承办奥运、中国股市连连下挫，上证指数从2007年的最高点6124点一路下滑至2008年8月份的最低点2284点、半年光景约有6.7万家中小企业倒闭、国际油价的居高不下，煤、电、油、运全面紧张。

针对新局势、新变化，我国政府把防过热、防通胀的经济政策迅速调整为保发展、控通胀。确保发展和控制物价是对立的统一，既有矛盾，也可以相互促进，关键在于我们采取什么样的政策，拉动GDP的三大要素是——投资、消费和进出口，根据相关研究今年经济增长如果不超过9.4%，通货膨胀率控制在5%左右，就是一个很好的平衡点，能为明年打下一个好的基础[2]。

中国经济正处在低谷的边缘。因为我国经济面临内忧外患，内忧是通货膨胀，外患是全球经济放缓，这些都对我国经济有很大影响，我们正在经历着动荡的考验：美元走低、人民币升值、外需放缓，这对于对外依存度超过60%的中国经济，是一次巨大的挑战[3]。

对于中小企业而言在投资和出口问题上主要依赖于国家的宏观调控，就困境中的中小企业本身来说基本上是无能为力的，然而可以团结起来、集合资源，向管理要效益，向降低成本要效益，从扩大内需中要效益，那么就要进一步激励民众扩大内需、大力推动消费、刺激消费，寻找一种能够产生新的消费热情的方法上下功夫，在实现消费增值的基础上取得企业效益，从而保持企业持续健康的发展。

3.2工业生产总值的概述

工业总产值是指以货币表现的工业企业在报告期内生产的工业产品总量。工业总产值按“工厂法”计算，即以工业企业作为一个整体，按企业工业生产活动的最终成果计算[5]。企业内部不允许重复计算，不能把企业内部各个车间生产的成果相加。工业总产值包括成品价值、对外加工费收入和自制半成品、在产品期末期初差额价值[8]。

15

河北工程大学本科毕业设计（论文）

4 工业生产总值的多因素模型分析

4.1建立多因素分析模型 设

Y??0??1X1??2X2??3X3??4X4??5X5??6X6??7X7??,

其中：

Y为工业总产值（按当年价格，单位：亿元） X1为塑料制品产量（单位：万吨） X2为水泥产量（单位：万吨）

X3为平板玻璃产量（单位：万重量箱）

X4为生铁产量（单位：万吨）

X5为粗钢产量（单位：万吨） XX67为钢筋产量（单位：万吨） 为盘条产量（单位：万吨）

?为随机误差项。

4.2数据收集

研究工业生产总值与经济发展之间的关系严格来说可以对每个地区进行研究，但是具体到各个地区，根据各个地区的条件情形不同，所以对研究得出的数据也会有很大的差异，所以此次就对同一地区的工业生产总值与经济发展之间的关系作为研究，就对一个地区的数据进行收集，然后得出结论，这个结论可以反映出这个地区的工业生产值与经济发展的关系，为促进经济发展，针对某一地区的经济发展，我们可以对它进行几十年的数据收集，得出这个地区的经济发展情况。下面我们以中国2007年各省市的主要工业产品产量统计量为研究对象，通过对这些数据和工业总产值的数据建立的线性关系来预测工业总产值。

(1) 工业总产值来自中国2007年按地区分组的专用设备制造业工业企业主要经济(2) 原煤来自中国2007年按地区分组的主要工业产品产量统计(一)； (3) 由于西藏自治区一行中缺省数据太多，故删除了西藏一栏； (4) 中国2007年按地区分组的主要工业产品产量统计(六)

具体数据见下表4.1：

16

指标统计(一)；

河北工程大学本科毕业设计（论文）

表4.1 原始数据

地区 全国 北京 天津 河北 山西 内蒙古 辽宁 吉林 黑龙江 上海 江苏 浙江 安徽 福建 江西 山东 河南 湖北 湖南 广东 广西 海南 重庆 四川 贵州 云南 陕西 甘肃 青海 宁夏 新疆 工业总产值(当年价格) 10591.98 331.92 204.82 376.27 258.49 124.34 680.18 94.57 175.41 600.57 1425.76 875.07 173.27 253.94 56.09 1736.69 797.9 126.58 425.56 923.24 175.53 2.03 84.85 364.4 28.31 37.24 181.62 52.41 1.22 16.23 7.46 塑料制品 水 泥 平板玻璃 (万吨) (万吨) (万重量箱) 3305.23 32.66 40.9 137.37 13.44 6.05 159.3 21.97 14.06 72.46 327.37 768.88 111.17 106.66 12.95 315.47 121.69 53.71 23.25 754.45 21.71 2.25 22.68 83.51 7.96 16.11 5.14 10.32 0.27 2.02 39.47 136117.25 1168.6 614.79 9758.28 2780.91 2871.17 3893.2 1903.81 1645.06 959.44 11849.78 10548.51 5402.23 4500.1 5008.54 15023.89 9471.36 5638.85 5683.28 9799.57 4350.48 633.32 3000.05 6375.6 2059.06 3568.53 3175.49 1540.21 436.85 817.36 1479.28 53918.07 235.22 186.06 10031.76 971.91 1395.72 1941.21 850.76 490.21 744.53 6856.11 2917.53 472.56 1867.12 686.64 5175.61 3619.93 2178.34 1518.12 6123.03 530.34 7.68 254.18 2495.79 2.5 329.8 1175.08 591.89 92.95 61.82 113.66 生 铁 （万吨） 47651.63 780.51 1435.4 10523.01 3727.64 1260.09 4057.59 545.66 374.11 1790.36 3802.15 238.08 1517.7 477.88 1047.36 4906.67 1974.95 1679.79 1247.76 755.25 639.3 18.84 328.38 1470.73 363.19 1202.78 365.55 592.78 90.09 46.25 391.82 粗 钢 （万吨） 48928.8 810.76 1602.13 10569.29 2506.36 1040.36 4140.27 599.67 436.05 2081.58 4721.47 577.23 1663.61 588.82 1306.81 4406.91 2275.39 1778.17 1331.79 1154.03 765.67 4.54 358.36 1415.34 349.36 883.85 396.27 602.8 114.71 0.36 446.85 钢 筋 （万吨） 10275.48 296.34 185.19 1335.43 372.93 100.96 246.94 30.59 184.5 23.21 1577.49 131.63 501.04 289.47 496.39 1090.68 556.84 226.99 337.18 444.76 288.06 6.54 49.76 492.78 176.45 228.43 318.94 98.84 1.31 37.77 148.06 盘 条 (万吨) 7919.02 416.3 64.14 1171.84 465.82 76.63 453.91 69.16 33.3 78.48 949.47 183.74 257.69 222.27 332.42 638.12 660.1 178.93 224.3 284.86 207.22 0.67 102.41 195.91 118.83 183.61 53.67 140.53 7.54 5.08 142.08 注1)：数据来自中国2007年统计年鉴8w。

17

河北工程大学本科毕业设计（论文）

4.3 统计检验 由SPSS计算结果

表4.2 Variables Entered/Removed(b) Variables Variables Model Method Entered Removed 1 盘条x7, 塑料制品x1, 生铁x4, 平板0.976 Enter 玻璃x3, 水泥x2, 钢筋x6, 粗钢x5 a All requested variables entered. b Dependent Variable: 工业总产值y

表4.3 Model Summary(b) Adjusted R Std. Error of the Model R R Square Durbin-Watson Square Estimate 1 0.993(a) 0.987 0.983 245.84347 1.739 a Predictors: (Constant), 盘条x7, 塑料制品x1, 生铁x4, 平板玻璃x3, 水泥x2, 钢筋x6, 粗 钢x5

b Dependent Variable: 工业总产值y

Model 表4.4 ANOVA(b) Sum of Squares df Mean Square F Sig. 1 Regression 105370311.469 7 15052901.638 249.059 0.000(a) Residual 1390097.246 23 60439.011 Total 106760408.714 30 a Predictors: (Constant), 盘条x7, 塑料制品x1, 生铁x4, 平板玻璃x3, 水泥x2, 钢筋x6, 粗钢x5 b Dependent Variable: 工业总产值y

表4.5 Coefficients(a)

Unstandardized Standardized Collinearity Statistics Coefficients Coefficients Model t Sig. B 1 (Constant) -2.211 Std. Error 47.145 0.452 0.024 0.056 Beta 0.317 0.460 -0.281 -0.047 2.187 1.535 -0.978 0.083 0.004 1.260 -.032 0.963 0.039 0.138 0.338 0.934 0.981 0.220 0.975 Tolerance 0.027 0.006 0.007 0.001 0.001 0.005 0.004 VIF 37.156 158.745 145.392 858.829 1170.904 216.237 242.790 塑料制品x1 0.989 水泥x2 0.036 平板玻璃x3 -0.055 生铁x4 0.013 0.154 0.058 粗钢x5 0.004 0.176 0.020 钢筋x6 0.456 0.362 0.441 盘条x7 -0.016 0.499 -0.012 a Dependent Variable: 工业总产值y 18

河北工程大学本科毕业设计（论文）

即有：

??-2.211?0.989X?0.036X?0.055X?0.013X?0.004X?0.456X?0.016X Y123456722R=0.987 ,R=0.983,F=249.059

方程的决定系数较高，R2=0.987，修正可决系数R2=0.983，又回归模型拟合程度与决定系数有关，决定系数越大，模型拟合程度越高，可见本模型拟合程度较好，又F=249.059>10, 模型总体显著性检验得知模型总体显著，也就是全体解释变量总体对被解释变量存在明显影响。对t检验由上述分析结果知，除粗钢和原煤外其余预测变量都是显著的。故我们对上述模型进行计量经济学的检验，并进行修正，看是否能使模型方程得到改进。

其中我们看到平板玻璃和盘条产量的系数是负值，一般来说，平板玻璃和盘条产量等经济发展量应该与工业总产值成正比关系，但由于在研究具体某个地区的经济发展关系时，由于政策领导或其它的因为地区的特殊性的原因，造成了平板玻璃和盘条产量的系数是可能成为负值的。比如工业总产值在下降，但是因为政策或其它原因平板玻璃和盘条产量却在上升，或工业总产值在上升，但是平板玻璃和盘条产量却在下降。

19

河北工程大学本科毕业设计（论文）

4.4 计量经济学检验及模型修正

4.4.1 异方差检验

计算残差绝对值与个自变量的斯皮尔曼相关系数，结果如下表：

表4.6 Correlations

Correlations 塑料制品Correlation Spearmx1 Coefficient an's rho Sig. 1-tailed 水泥x2 N Correlation Coefficient Sig. 1-tailed N 塑料制品水 泥 平板玻璃 x1 x2 x3 1 0 31 0.735 0 31 0.738 0 31 0.615 0 31 0.685 0 31 0.592 0 31 0.711 0 31 0.590 0 31 0.735 0 31 1 0.002 31 0.855 0 31 0.547 0.001 31 0.604 0 31 0.782 0 31 0.760 0 31 0.470 0.004 31 0.738 0 31 0.855 0 31 1 0 31 0.665 0 31 0.717 0 31 0.665 0 31 0.690 0 31 0.287 0.059 31 生铁 粗 钢 钢筋 盘 条 abs_1 x4 x5 x6 x7 0.615 0.685 0.592 0.711 0 31 0 31 0 31 0 31 0.590 0 31 0.547 0.604 0.782 0.760 0.470 0.001 0 31 31 0 31 0 31 0.004 31 平板玻璃Correlation Coefficient x3 生铁x4 粗钢x5 钢筋x6 盘条x7 abs_1 Sig. 1-tailed N Correlation Coefficient Sig. 1-tailed N Correlation Coefficient Sig. 1-tailed N Correlation Coefficient Sig. 1-tailed N Correlation Coefficient Sig. 1-tailed N Correlation Coefficient Sig. 1-tailed N 0.665 0.717 0.665 0.690 0.287 0 31 1 . 31 0 31 0 31 0 31 0.059 31 0.986 0.700 0.754 0.489 0 31 0 31 0 31 0.003 31 0.986 1 0 31 0.70 0 31 . 31 0.720 0.788 0.513 0 31 0 31 0.002 31 0.720 1 0 31 0 31 0.858 0.416 0 31 0.01 31 0.547 0.001 31 0.754 0.788 0.858 1 0 31 0 31 0 31 0 31 0.489 0.513 0.416 0.547 1 0.003 0.002 0.01 31 31 31 0.001 0.001 31 31 Correlation is significant at the 0.01 level 1-tailed. 残差绝对值与各自变量的相关系数分别为0.59，0.47，0.287,0.489,0.513,0.416，0.547说明存在异方差，需采用加权最小二乘估计；

加权估计，得到最优权重为3.000，对其进行加权分析。

20

河北工程大学本科毕业设计（论文）

表4.7 Variables Entered/Removed(b,c)

Variables Variables Model Method Entered Removed 1 盘条x7, 塑料制品x1, 钢筋x6, 生铁0.536 Enter x4, 平板玻璃x3, 粗钢x5, 水泥x2 a All requested variables entered. b Dependent Variable: 工业总产值y

c Weighted Least Squares Regression - Weighted by Weight for 工业总产值y from WLS, MOD_4 abs_1** -3.000

表4.8 Model Summary(b,c)

Adjusted R Std. Error of R Square Durbin-Watson Square the Estimate 1.000 1.000 0.08592 1.581 Model R 1 1.000(a) a Predictors: (Constant), 盘条x7, 塑料制品x1, 钢筋x6, 生铁x4, 平板玻璃x3, 粗钢x5, 水泥x2 b Dependent Variable: 工业总产值y

c Weighted Least Squares Regression - Weighted by Weight for 工业总产值y from WLS, MOD_4 abs_1** -3.000

Sum Squares 2587.098 0.170 2587.268 表4.9 ANOVA(b,c) of df Mean Square 7 23 30 369.585 0.007 Model 1 Regression Residual Total F 50065.439 Sig. 0.000(a) a Predictors: (Constant), 盘条x7, 塑料制品x1, 钢筋x6, 生铁x4, 平板玻璃x3, 粗钢x5, 水泥x2 b Dependent Variable: 工业总产值y

c Weighted Least Squares Regression - Weighted by Weight for 工业总产值y from WLS, MOD_4 abs_1** -3.000

21

河北工程大学本科毕业设计（论文）

Model 1 (Constant) 塑料制品x1 水泥x2 平板玻璃x3 生铁x4 粗钢x5 钢筋x6 盘条x7 表4.10 Coefficients(a,b)

Unstandardized Standardized Coefficients Coefficients t B Std. Error Beta -13.013 0.841 0.020 -0.017 -0.008 0.080 0.405 -0.217 3.486 0.095 0.010 0.019 0.019 0.021 0.055 0.104 0.264 0.260 -0.084 -0.038 0.370 0.393 -0.163 -3.733 8.823 2.103 -0.866 -0.440 3.836 7.408 -2.090 Sig. 0.001 0.000 0.047 0.395 0.664 0.001 0.000 0.048 Collinearity Statistics Tolerance 0.003 0.000 0.000 0.000 0.000 0.001 0.000 VIF 314.686 5366.476 3333.757 2602.303 3255.656 984.360 2131.465 由上表所知，方程的决定系数较高，R2=1.000，修正可决系数R2=1.000，又回归模型拟合程度与决定系数有关，决定系数越大，模型拟合程度越高，可见本模型拟合程度较好，又F=50065.439>10, 模型总体显著性检验得知模型总体显著，也就是全体解释变量总体对被解释变量存在明显影响。但是D.W值为1.581，要比原来多元线性回归的

D.W值要小，自相关性比较严重，需要进一步讨论

4.4.2 自相关检验

自相关(Autocorrelation)是对随机扰动项之间相互独立假定的违背，指扰动项序

列相邻期之间不是随机独立而是存在相关关系，又称为序列相关。自相关主要表现在时间序列中。因此对于线性回归模型：

Y??0??1X1??2X2??3X3??4X4??5X5??6X6??7X7??.

自相关可表示为：

Cov(?i,?j)?0,i,j?1,2,??,n.

用于检验扰动项是否存在自相关的方法主要有：D-W 检验（Durbin-Watson 杜宾-

瓦特森检验），D-W 检验是Durbin 和Watson 于1951 年提出的一种自相关检验方法。它只适用于扰动项的形式为：?i???i-1??i'（?为自相关系数）的一阶自相关问题。这种方法是最常用的一阶自相关检验方法。其检验步骤如下： （1）提出假设。

H0：??0，即扰动项不存在一阶自相关； H1：??0，即扰动项存在一阶自相关。

n（2）构造统计量。定义DW 统计量为DW??(ei?2in?ei?1)2i2，对于大样本可以把D.W统计

?ei 22

河北工程大学本科毕业设计（论文）

?)，其中??为自相关系数?的估计。因为|??|?1，所以D.W检验统量写为DW?2(1??计的值域为0?DW?4。

（3）判断。根据样本容量n 和解释变量的数目p 查D.W分布表，得下临界值DL和上临界值DU，并依下列准则判断扰动项的自相关情形。

①如果0

②如果DL

③如果DU

④如果4-DU

⑤如果4-DL

由上图4.3得知Durbin-Watson Statistics=1.955873，查表得在显著水平?=0.05下，查表 n=30，k=3时，DL=1.21，DU=1.55，由于

DU=1.21

则接受零假设，扰动项不存在一阶正自相关。DW 越接近2，判断无自相关性把握越大。

由表可知，D.W值为1.739，可以看出存在正自相关性,通过迭代

??1??12D.W,??yt?1?yt??yt???xt?1?xt??xt??

??0.1305.进行迭代计算,得 来改善自相关性。?表4.11 Variables Entered/Removed(b)

Variables Variables Model Method Entered Removed 1 x71, x11, x41, x31, . x21, x61, x51(a) Enter a All requested variables entered. b Dependent Variable: y11

表4.12 Model Summary(b)

Adjusted R Std. Error of R Square Durbin-Watson Square the Estimate 0.988 0.984 243.30632 1.867 Model R 1 0.994(a) a Predictors: (Constant), x71, x11, x41, x31, x21, x61, x51 b Dependent Variable: y11

23

河北工程大学本科毕业设计（论文）

表4.13 ANOVA(b)

Model Sum of Squares 108742477.288 1361553.187 110104030.476 df 7 23 30 Mean Square 15534639.613 59197.965 F 262.418 Sig. .000(a) 1 Regression Residual Total a Predictors: (Constant), x71, x11, x41, x31, x21, x61, x51 b Dependent Variable: y11

表4.14 Coefficients(a)

Unstandardized Standardized Coefficients Coefficients t B 1.362 1.055 0.034 0.058 0.060 0.043 0.507 0.051 Std. Error 45.897 0.454 0.024 0.054 0.152 0.173 0.348 0.466 Beta 0.339 0.438 0.298 0.273 0.201 0.491 0.038 0.030 2.324 1.448 1.075 0.396 0.251 1.456 0.109 Model 1 (Constant) x11 x21 x31 x41 x51 x61 x71 Sig. 0.977 0.029 0.161 0.294 0.696 0.804 0.159 0.915 Collinearity Statistics Tolerance VIF 0.025 39.639 0.006 170.277 0.007 143.309 0.001 887.679 0.001 1198.966 0.005 211.372 0.004 222.321 a Dependent Variable: y11

方程的决定系数较高，R2=0.988，修正可决系数R2=0.984，又回归模型拟合程度与决定系数有关，决定系数越大，模型拟合程度越高，可见本模型拟合程度较好，又F=262.418>10, 模型总体显著性检验得知模型总体显著，也就是全体解释变量总体对被解释变量存在明显影响，并且D.W值为1.867，与原来的线性回归结果相比自相关性得到了进一步的优化。此时的多元线性回归方程为：

???1.362?1.055X?0.034X?0.058X?0.060X?0.043X?0.507X?0.054XY1111213141516171由表可知，D.W值为1.867，可以看出自相关性得到了进一步改善,下面再次迭代：

??1??12D.Wttt?1,??x??x???xt?1 t?t?yy??y????0.0665.进行迭代计算，得到如下结果： 来改善自相关性。?表4.15 Variables Entered/Removed(b)

Variables Variables Model Entered Removed Method 1 X72, X12, X42, X32, . Enter X22, X62, X52(a) 24

河北工程大学本科毕业设计（论文）

a All requested variables entered. b Dependent Variable: Y22

表4.16 Model Summary(b)

Model 1 R 0.994(a) R Square Adjusted R Std. Error of Square the Estimate Durbin-Watson 0.988 0.984 243.39980 2.072 a Predictors: (Constant), X72, X12, X42, X32, X22, X62, X52 b Dependent Variable: Y22

表4.17 ANOVA(b)

Model 1 Sum of Squares Regression 112177062.223 Residual 1362599.661 Total 113539661.885 df 7 23 30 Mean Square F Sig. 16025294.603 270.499 .000(a) 59243.464 a Predictors: (Constant), X72, X12, X42, X32, X22, X62, X52

b Dependent Variable: Y22

表4.18 Coefficients(a)

Unstandardized Standardized Coefficients Coefficients B Std. Error -1.459 45.400 1.120 0.453 0.032 0.024 -0.060 0.051 0.129 0.151 -0.112 0.170 0.550 0.327 -0.082 0.422 Beta t -0.032 0.361 2.472 0.414 1.365 -0.308 -1.175 0.587 0.857 -0.521 -0.659 0.534 1.682 -0.061 -0.195 1 Model (Constant) X12 X22 X32 X42 X52 X62 X72 Collinearity Statistics Sig. Tolerance VIF 0.975 0.021 0.024 40.909 0.185 0.006 176.467 0.252 0.008 131.958 0.400 0.001 898.594 0.517 0.001 1198.068 0.106 0.005 193.112 0.847 0.005 187.968 a Dependent Variable: Y22

22由上表可知方程的决定系R=0.988，修正可决系数R=0.984，与上次迭代相比基本

稳定，且均接近于1，可见本模型拟合程度较好，又F=270.499>10, 模型总体显著性检验得知模型总体显著，也就是全体解释变量总体对被解释变量存在明显影响，并且D.W值为2.072，自相关性的影响基本消除。此时的多元线性回归方程为：

???1.459?1.120YX1?0.032X2?0.060X3?0.129X4?0.112X5?0.550X6?0.082X7

5 结 论

25

河北工程大学本科毕业设计（论文）

在本文中，我们通过一系列主要工业产品产量的分析，探索通过个工业产品产量来预测工业总产值的方法。我们在本文中选取了工业总产值，塑料制品产量，水泥产量，平板玻璃产量，生铁产量，粗钢产量，钢筋产量，盘条产量。在通过对模型的异方差性和自相关性的检验以及修正之后我们得到的最终的线性回归方程为：

???1.459?1.120X?0.032X?0.060X?0.129X?0.112X?0.550X?0.082XY1234567

通过以上线性回归方程我们可以知道，工业总产值与塑料制品产量有较大关系，塑料产量每增加一个单位，工业总产值就会增加1.120个单位，水泥产量每增加1个单位总也总产值就会增加0.032个单位，生铁产量每增加一个单位，工业总产值就会增加0.129个单位，同样的钢筋每增加1个单位工业总产值就会增加0.55个单位。 但同时我们看到有些主要工业产品的产量系数是负值，一般来说，各主要工业产品的产量等经济发展量应该与工业总产值成正比关系，但由于在研究具体某个地区的经济发展关系时，由于政策领导或其它的因为地区的特殊性的原因，造成了水泥产量的系数是可能成为负值的。比如工业总产值在下降，但是因为政策或其它原因一些工业产品产量却在上升，或工业总产值在上升，但是某些工业产品的产量却在下降。

在我国，目前经济发展速度很快，工业生产总值能够保持稳定增长，我们要保持住这种增长势头，实现经济的又好又快发展，那么对于多元线性回归研究经济发展的问题就有很现实的作用。并且，在我国，经济发展还不是很平衡，对各地区进行经济发展的研究，使我们可以在宏观调控的条件下，也可以结合各地区的实际情况作出经济决策，使经济发展达到最好的效果。

由以上得到的回归方程我们可以看到变量X1和X6的系数绝对值最大，可见对Y的影响是最大，X1即机塑料制品产量，X6是钢筋产量，所以对全国各地区从目前的经济发展状况和条件看，塑料制品产量和钢筋产量对工业总产值有最大的正影响关系，也就是说，为了最高的提高全国各地区工业总产值的指标，我们可以加大对塑料制品和钢筋的生产，它对全国各地区工业总产值的指标有最明显的正向关系。这样，我们在对于资源的配置方面，可以先考虑对各地区纸塑料制品和钢筋制造的投资加大，能够达到资源的优化配置。

致 谢

26

河北工程大学本科毕业设计（论文）

在本文的体系和材料数据的收集及处理过程中及在撰写整篇论文的过程中，王小胜王导师给了我很大的帮助，在材料的选择、文章的重心、参数估计统计检验及计量经济学检验等许多方面给了我有益的指导和启迪。在此，我表示最衷心的感谢。

参考文献

27

河北工程大学本科毕业设计（论文）

[1] 谢建平.回归分析在经济中的应用[J].学术月刊,2009,11 (6):71-74. [2] 张润楚.数理统计学[M].北京:科学出版社,2010. [3] 谢宇.回归分析[M].北京:社会科学文献出版社,2010. [4] 王松桂.线性模型引论[M].北京:科学出版社,2003.

[5] 陈正昌,程炳林,陈新丰等.多变量分析方法统计软件应用[M].北京:中国税务出版社,2008.

[6] 米红 张文璋.实用现代统计分析方法与SPSS应用[M].北京:当代中国出版社,2006.

[7] 任庆玉.当代中国经济[M].北京:中国经济出版社,2005.

[8] 何东耀.中国4万亿投资解读[M].北京:中国人民大学出版社,2008.

[9] Harbert White.A Heteroskedasticity-Consistent Covariance Matrix Estimator and a DirectTest for Heteroskedasticiry[J]. The Econometric Society May,1980,17(2):817-838.

[10] Eilat P, Einav L.Determinants of international tourism: a three-dimensional panel data analysis. Applied Economics,2004,36( 12) : 1315-1327. 34

[11]王祖兰.全国各地区工业经济实力分析[J].复印报刊资料,2009,15 (7):177-181.

28

本文来源：https://www.bwwdw.com/article/8v53.html