第二章经典线性回归模型.ppt

《第二章经典线性回归模型.ppt》由会员分享，可在线阅读，更多相关《第二章经典线性回归模型.ppt（115页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1第二章经典线性回归第二章经典线性回归模型模型现在学习的是第1页，共115页第一节线性回归模型的概念第二节线性回归模型的估计第三节拟合优度第四节非线性关系的处理第五节假设检验第六节预测第七节虚拟变量现在学习的是第2页，共115页第一节 线性回归模型的概念 一.双变量线性回归模型 我们在上一章给出的需求函数的例子 Q=+P+u （2.1）是一个双变量线性回归模型，模型中只有两个变量，一个因变量，一个解释变量，由解释变量的变动来解释因变量的变动，或者说用因变量对解释变量进行线性回归，因而称为双变量线性回归模型双变量线性回归模型，亦称简单线性回归模型简单线性回归模型。让我们再看一个例子。C=+D+u

2、 (2.2)这是凯恩斯消费函数，其中C为消费支出，D为个人可支配收入，u为扰动项（或误差项）。现在学习的是第3页，共115页 此模型中，方程左端的消费支出（C）为因变量（或被解释变量），方程右端的个人可支配收入（D）为解释变量（或自变量）。和是未知参数，由于双变量线性回归模型的图形是一条直线，因而和习惯上又分别称为截距和斜率。这里斜率的含义是解释变量增加一个单位所引起的因变量的变动。例如在(2.2)式中，的含义是个人可支配收入增加一个单位所引起的消费的增加量，经济学中称之为边际消费倾向（MPC）。截距的含义是解释变量为0时的值。截距有时有经济含义，但大多数情况下没有，因此，在计量经济分析中，通

3、常不大关注的取值如何。现在学习的是第4页，共115页 在教学中，我们习惯上采用Y表示因变量，X表示解释变量，双变量线性回归模型的一般形式为：Y=+X+u 在实践中，此模型被应用于因变量和解释变量的一组具体观测值 和 （t=1,2,n），因而模型表示为：=+ut t=1,2,n (2.3)它表明，对于n个时期t=1,2,n，该模型成立。更一般的形式为：=+ui,i=1,2,.,n (2.4)即模型对X和Y的n对观测值（i=1,2,n）成立。(2.3)式一般用于观测值为时间序列的情形，在横截面数据的情形，通常采用(2.4)式。现在学习的是第5页，共115页二、多元线性回归模型 在许多实际问题中，我

4、们所研究的因变量的变动可能不仅与一个解释变量有关。因此，有必要考虑线性模型的更一般形式，即多元线性回归模型：t=1,2,n 在这个模型中，Y由X1、X2、X3、XK所解释，有K+1个未知参数0、1、2、K。这里，“斜率”j的含义是其其它它变变量量不不变变的的情情况况下下，Xj改变一个单位对因变量所产生的影响。现在学习的是第6页，共115页 例例2.2 2.2 食品需求方程食品需求方程 其中，Y=在食品上的总支出 X=个人可支配收入 P=食品价格指数 用美国1959-1983年的数据，得到如下回归结果（括号中数字为标准误差）：Y和X的计量单位为10亿美元(按1972不变价格计算).现在学习的是第

5、7页，共115页多元线性回归模型中斜率系数的含义上例中斜率系数的含义说明如下：价格不变的情况下，个人可支配收入每上升10亿美元（1个billion），食品消费支出增加1.12亿元（0.112个 billion）。收入不变的情况下，价格指数每上升一个点，食品消费支出减少7.39亿元（0.739个billion）现在学习的是第8页，共115页回到一般模型 t=1,2,，n即对于n组观测值，有现在学习的是第9页，共115页其矩阵形式为：其中 现在学习的是第10页，共115页第二节 线性回归模型的估计 一经典一经典线性回归模型的统计假设（1）E(ut)=0,t=1,2,n 即各期扰动项的均值（期望值）

6、均为0。均值为0的假设反映了这样一个事实：扰动项被假定为对因变量的那些不能列为模型主要部分的微小影响。没有理由相信这样一些影响会以一种系统的方式使因变量增加或减小。因此扰动项均值为0的假设是合理的。现在学习的是第11页，共115页（2）E(ui uj)=0,ij 即各期扰动项互不相关。也就是假定它们之间无自相关或无序列相关。实际上该假设等同于：cov(ui,uj)=0,ij这是因为：cov(ui,uj)=Eui-E(ui)uj-E(uj)=E(uiuj)根据假设（1）（3）E(ut2)=2,t=1,2,n 即各期扰动项的方差是一常数，也就是假定各扰动项具有同方差性。这是因为：Var(ut)=E

7、ut-E(ut)2=E(ut2)根据假设（1）现在学习的是第12页，共115页（4）Xjt是非随机量，j=1,2,k t=1,2,n（5）（K+1）n;即观测值的数目要大于待估计的参数的个数 （要有足够数量的数据来拟合回归线）。（6）各解释变量之间不存在严格的线性关系。上述假设条件可用矩阵表示为以下四个条件：现在学习的是第13页，共115页A1.E(u)=0 A2.由于 显然，仅当 E(ui uj)=0,ij E(ut2)=2,t=1,2,n 这两个条件成立时才成立，因此，此条件相当前面条件(2),(3)两条，即各期扰动项互不相关，并具有常数方差。现在学习的是第14页，共115页 A3.X 是

8、一个非随机元素矩阵。A4.Rank(X)=(K+1)n.-相当于前面(5)(6)两 条 即矩阵X的秩=（K+1)n 满足条件（A1）（A4）的线性回归模型称为经经典典线线性性回归模型回归模型或古典线性回归模型（古典线性回归模型（CLR模型）模型）。当然，为了后面区间估计和假设检验的需要，还要加上一条：A5.各期扰动项服从正态分布。，t=1,2,n现在学习的是第15页，共115页二、最小二乘估计1.最小二乘原理 为了便于理解最小二乘法的原理，我们用双变量线性回归模型作出说明。对于双变量线性回归模型Y=+X+u，我们的任务是，在给定X和Y的一组观测值 (X1,Y1),(X2,Y2),.,(Xn,Y

9、n)的情况下,如何求出 Yt=+Xt+ut 中 和 的估计值 和 ,使得拟合的直线为“最佳”。直观上看，也就是要求在X和Y的散点图上穿过各观测点画出一条“最佳”直线，如下图所示。现在学习的是第16页，共115页 *et *YXXt 图 2.2 Yt现在学习的是第17页，共115页残差残差 拟合的直线 称为拟合的回归线.对于任何数据点(Xt,Yt),此直线将Yt 的总值 分成两部分。第一部分是Yt的拟合值或预测值 ：,t=1,2,n第二部分，et，代表观测点对于回归线的误差，称为拟合或预测的残残差差（residuals）：）：t=1,2,n 即 t=1,2,n现在学习的是第18页，共115页残差

10、平方和 我们的目标是使拟合出来的直线在某种意义上是最佳的，直观地看，也就是要求估计直线尽可能地靠近各观测点，这意味着应使残差总体上尽可能地小。要做到这一点，就必须用某种方法将每个点相应的残差加在一起，使其达到最小。理想的测度是残差平方和，即现在学习的是第19页，共115页最小二乘法 最小二乘法就是选择一条直线，使其残差平方和达到最小值的方法。即选择 和 ，使得达到最小值。现在学习的是第20页，共115页 运用微积分知识，使上式达到最小值的必要条件为：即现在学习的是第21页，共115页整理，得：此二式称为正规方程。解此二方程，得：其中：样本均值 离差现在学习的是第22页，共115页2多元线性回归

11、模型的最小二乘估计多元线性回归模型的最小二乘估计在多元线性回归模型的情况下，在多元线性回归模型的情况下，我们的模型是：问题是选择 ，使得残差平方和最小。残差为：现在学习的是第23页，共115页要使残差平方和 为最小，则应有：我们得到如下K+1个方程（即正规方程）：现在学习的是第24页，共115页按矩阵形式，上述方程组可表示为：现在学习的是第25页，共115页=即现在学习的是第26页，共115页三.最小二乘估计量 的性质 我们的模型为 估计式为 1 的均值现在学习的是第27页，共115页（由假设3）(由假设1)即 这表明，OLS估计量 是无偏估计量。现在学习的是第28页，共115页2 的方差为求

12、Var()，我们考虑 现在学习的是第29页，共115页不难看出，这是 的方差-协方差矩阵，它是一个(K+1)(K+1)矩阵，其主对角线上元素为各系数估计量的方差，非主对角线上元素为各系数估计量的协方差。现在学习的是第30页，共115页由上一段的(2.19)式，我们有因此现在学习的是第31页，共115页 请注意，我们得到的实际上不仅是 的方差，而且是一个方差-协方差矩阵，为了反映这一事实，我们用下面的符号表示之：为方便起见，我们也常用 表示 的方差-协方差矩阵，因此上式亦可写作：需要注意的是，这里 不表示方差向量，而是方差-协方差矩阵。现在学习的是第32页，共115页4 高斯-马尔科夫定理对于

13、以及标准假设条件A1A4,普通最小二乘估计量(OLS估计量）是最佳线性无偏估计量（BLUE）。我们已在上一段中证明了无偏性，下面证明线性和最小方差性。现在学习的是第33页，共115页由OLS估计量 的公式 可知,可表示为一个矩阵和因变量观测值向量 的乘积：其中 是一个(K+1)*n 非随机元素矩阵。因而 是线性估计量。现在学习的是第34页，共115页 现设 为 的任意一个线性无偏估计量，即其中 是一个(K+1)*n非随机元素矩阵。则 显然，若要 为无偏估计量，即 ，只有 ，为（K+1）阶单位矩阵。现在学习的是第35页，共115页 的方差为：我们可将 写成 从而将 的任意线性无偏估计量 与OLS

14、估计量 联系起来。现在学习的是第36页，共115页由 可推出：即 因而有 由 从而 ，因此上式中间两项为0，我们有现在学习的是第37页，共115页 因此 最后的不等号成立是因为 为半正定矩阵。这就证明了OLS估计量 是 的所有线性无偏估计量中方差最小的。至此，我们证明了高斯-马尔科夫定理。现在学习的是第38页，共115页4 的分布的分布 我们在前面列出的假设条件（A5）表明，N(0,),t=1,2,n即各期扰动项服从均值为0、方差为 的正态分布。考虑到假设条件（A3），即 是一个非随机元素矩阵，则由前面(2.20)式：我们有：现在学习的是第39页，共115页这表明，是N个正态分布变量 的线性函

15、数，因而亦为正态分布变量，即 （2.22）由此可知，系数估计量向量的每个元素都是正态分布的，即 j0,1,k （2.23）其中cjj为矩阵中 的（j1,j1）元素（主对角线上第j1个元素）。现在学习的是第40页，共115页第三节 拟合优度一决定系数R2 在估计了线性回归模型之后，一个很自然的问题是，估计出的回归线与观测值拟合得好不好？这就是拟合优度要解决的问题。拟合优度的一个通行的测度是因变量Y的（样本）变差被模型所解释的比例，也就是因变量Y的变差被诸解释变量所解释的比例。这个统计量称为决决定定系系数数（coefficient of determination），记做 ，定义为：现在学习的是第

16、41页，共115页其中，=残差平方和 ESS为Explained Sum of Squares的缩写；RSS为Residual Sum of Squares的缩写；TSS为Total Sum of Squares的缩写。现在学习的是第42页，共115页 决定系数 R2 计量了Y的总变差中可以归因于X和Y之间关系的比例，或者说Y的变动中可以由X的变动来解释的比例。它是回归线对各观测点拟合紧密程度的测度。我们有：完全拟合，：X与Y完全不存在线性关系，现在学习的是第43页，共115页 的值越高，拟合得越好。但什么是高？并没有绝对的标准，要根据具体问题而定。此外，回归中使用时间序列数据还是横截面数据也

17、有不同的标准。对时间序列数据来说，的值在0.8、0.9 以上是很常见的事,而在横截面数据的情况下，0.4、0.5的 值也不能算低。为方便计算，我们也可以用矩阵形式表示 。现在学习的是第44页，共115页 我们有：残差 其中，残差平方和：现在学习的是第45页，共115页而 将上述结果代入R2的公式，得到：这就是决定系数R2 的矩阵形式。现在学习的是第46页，共115页二修正决定系数：残差平方和的一个特点是，每当模型增加一个解释变量，并用改变后的模型重新进行估计，残差平方和的值会减小。由此可以推论，决定系数是一个与解释变量的个数有关的量：解释变量个数增加 减小 R2 增大也就是说，人们总是可以通过

18、增加模型中解释变量的方法来增大 R2 的值。因此，用 R2 来作为拟合优度的测度，不是十分令人满意的。现在学习的是第47页，共115页为此，我们定义修正决定系数 （Adjusted )如下：现在学习的是第48页，共115页 是经过自由度调整的决定系数，称为修正决定系数。我们有：（1）（2）仅当K=0时，等号成立。即 （3）当K增大时，二者的差异也 随之增大。（4）可能出现负值。现在学习的是第49页，共115页三例子三例子 下面我们给出两个简单的数值例子，以帮助理解这两节的内容.例2.3 Yt=1+2X2 t+3X3 t+u t 设观测数据为：Y：3 1 8 3 5 X2：3 1 5 2 4 X

19、3：5 4 6 4 6 试求各参数的OLS估计值，以及 。解：我们有现在学习的是第50页，共115页现在学习的是第51页，共115页现在学习的是第52页，共115页现在学习的是第53页，共115页现在学习的是第54页，共115页 例例2.4 设 n=20,k=3,R2=0.70，求 。解：下面改变n的值，看一看 的值如何变化。我们有 若n=10，则 =0.55 若n=5，则 =-0.20 由本例可看出，有可能为负值。这与R2不同 （）。现在学习的是第55页，共115页 第四节 非线性关系的处理 迄今为止，我们已解决了线性模型的估计问题。但在实际问题中，变量间的关系并非总是线性关系，经济变量间的

20、非线性关系比比皆是。如大家所熟悉的柯布-道格拉斯生产函数:就是一例。在这样一些非线性关系中，有些可以通过代数变换变为线性关系处理，另一些则不能。下面我们通过一些例子来讨论这个问题。现在学习的是第56页，共115页一.线性模型的含义 线性模型的基本形式是:其特点是可以写成每一个解释变量和一个系数相乘的形式。线性模型的线性包含两重含义：（1）变量的线性 变量以其原型出现在模型之中，而不是以X2或X之类的函数形式出现在模型中。（2）参数的线性 因变量Y是各参数的线性函数。现在学习的是第57页，共115页二线性化方法 对于线性回归分析，只有第二种类型的线性才是重要的，因为变量的非线性可通过适当的重新定

21、义来解决。例如，对于 此方程的变量和参数都是线性的。现在学习的是第58页，共115页 参数的非线性是一个严重得多的问题，因为它不能仅凭重定义来处理。可是，如果模型的右端由一系列的X或eX项相乘，并且扰动项也是乘积形式的，则该模型可通过两边取对数线性化。例如，需求函数 其中，Y=对某商品的需求 X=收入 P=相对价格指数 =扰动项可转换为：现在学习的是第59页，共115页 用X,Y,P的数据，我们可得到logY,logX和logP,从而可以用OLS法估计上式。logX的系数是的估计值，经济含义是需求的收入弹性，logP的系数将是的估计值，即需求的价格弹性。弹性（elasticity）是一变量变动

22、1%所引起的另一变量变动的百分比。其定义为 本例中，需求的收入弹性是收入变化1%，价格不变时所引起的商品需求量变动的百分比。需求的价格弹性是价格变化1%，收入不变时所引起的商品需求量变动的百分比。现在学习的是第60页，共115页三例子例2.5 需求函数 本章1中，我们曾给出一个食品支出为因变量，个人可支配收入和食品价格指数为解释变量的线性回归模型例子（例2.2）。现用这三个变量的对数重新估计（采用同样的数据），得到如下结果（括号内数字为标准误差）：回归结果表明，需求的收入弹性是0.64,需求的价格弹性是-0.48，这两个系数都显著异于0。现在学习的是第61页，共115页 例2.6 柯布-道格拉

23、斯生产函数 用柯布和道格拉斯最初使用的数据（美国1899-1922年制造业数据）估计经过线性化变换的模型得到如下结果（括号内数字为标准误差）：从上述结果可以看出，产出的资本弹性是0.23，产出的劳动弹性为0.81。现在学习的是第62页，共115页例2.7 货币需求量与利率之间的关系 M=a(r-2)b这里，变量非线性和参数非线性并存。对此方程采用对数变换 logM=loga+blog(r-2)令Y=logM,X=log(r-2),1=loga,2=b 则变换后的模型为：Yt=1+2Xt+ut 现在学习的是第63页，共115页 将OLS法应用于此模型，可求得1和2的估计值 ，从而可通过下列两式求

24、出a和b估计值：应当指出，在这种情况下，线性模型估计量的性质（如BLUE,正态性等）只适用于变换后的参数估计量 ，而不一定适用于原模型参数的估计量 和 。现在学习的是第64页，共115页 例2.8 上例在确定货币需求量的关系式时，我们实际上给模型加进了一个结束条件。根据理论假设，在某一利率水平上，货币需求量在理论上是无穷大。我们假定这个利率水平为2%。假如不给这一约束条件，而是从给定的数据中估计该利率水平的值，则模型变为：M=a(r-c)b 式中a,b,c均为参数。仍采用对数变换，得到 log(Mt)=loga+blog(rt-c)+ut t=1,2,n 我们无法将log(rt-c)定义为一个

25、可观测的变量X,因为这里有一个未知量c。也就是说，此模型无法线性化。在这种情况下，只能用估计非线性模型参数值的方法。现在学习的是第65页，共115页四非线性回归 模型 Y=a(X-c)b是一个非线性模型，a、b和c是要估计的参数。此模型无法用取对数的方法线性化，只能用非线性回归技术进行估计，如非线性最小二乘法（NLS）。该方法的原则仍然是残差平方和最小。计量经济软件包通常提供这类方法，本书第五章将对非线性回归方法作较深入的介绍，这里仅给出有关非线性最小二乘法的大致步骤如下：现在学习的是第66页，共115页非线性回归方法的步骤1 首先给出各参数的初始估计值（合理猜测值）;2用这些参数值和X观测值

26、数据计算Y的各期预测 值（拟合 值）;3计算各期残差，然后计算残差平方和e2;4对一个或多个参数的估计值作微小变动；5计算新的Y预测值 、残差平方和e2；6若新的e2小于老的e2，说明新参数估计值 优于老估计值，则以它们作为新起点；7重复步骤4，5，6，直至无法减小e2为止。8最后的参数估计值即为最小二乘估计值。现在学习的是第67页，共115页第五节 假设检验 本节讨论经典线性回归模型的区间估计和假设检验问题。我们的模型是：在第二节中我们证明了在扰动项服从正态分布的假设（A5）下，j0,1,k 其中 cjj 为矩阵 中的（j1,j1）元素（主对角线上第j1个元素）。这一结果为基于OLS估计量的

27、假设检验提供了坚实的基础。现在学习的是第68页，共115页一、的置信区间 我们可构造一个检验统计量 该变量服从均值为0、标准差为1的标准正态分布。与估计量相联系的概率分布的标准差，通常称为标准标准误差误差（standard error），用Se表示。的标准误差为：现在学习的是第69页，共115页 如果为已知，则由于检验统计量z服从标准正态分布，因而我们可以立即给出总体参数 的95%的置信区间为：但实际上，我们一般无法知道扰动项分布的方差 ，而必须根据观测值数据估计出 ，然后再来考虑 的置信区间的计算问题。现在学习的是第70页，共115页 1 2 的估计的估计可以证明，2的无偏估计量是 式中 是

28、残差平方和，分母是 的自由度，这是因为我们在估计 的过程中，失去了（K+1）个自由度。2.的置信区间的置信区间我们重新定义 的标准误差为：现在学习的是第71页，共115页则检验统计量 不再服从标准正态分布，而是服从自由度为（n-k-1）的t分布，即这里n和k分别为观测值和解释变量的数目。故 的（1）置信区间为：其中为显著性水平，通常取0.05。现在学习的是第72页，共115页例2.9 回到食品需求的例子（例2.2）：其中，Y=在食品上的总支出，X=个人可支配收入，P=食品价格指数 用美国1959-1983年的数据，得到如下回归结果（括号中数字为标准误差）：求 的95置信区间。现在学习的是第73

29、页，共115页由回归结果可知，我们不难得到 的95置信区间为：即为0.10580.1182。现在学习的是第74页，共115页二、假设检验的逻辑和步骤二、假设检验的逻辑和步骤 假设检验始于一个给定的假设，即所谓“原假设”，亦称“零假设”，然后计算检验统计量，这个检验统计量在原假设成立的假定下的概率分布是已知的。下一步是判断计算出的检验统计量的值是否不大可能来自此分布，如果判断是不大可能，则表明原假设不大可能成立。我们用一个例子来说明上述有关假设检验的思路。设有一个原假设规定 的值为 ，这里 是研究人员选择的一个值，如果这个原假设（H0：）成立，我们知道统计量 现在学习的是第75页，共115页应服

30、从自由度为(n-k-1)的t分布，即如果原假设不成立，则备择假设H1：成立。用于计算t的所有的量都是已知的，可以用估计值 及其标准误差Se()算出t的值，因此t可作为检验统计量用于假设检验，如果算出的t值绝对值过大，落入t分布的尾部，意味着原假设不大可能成立，因为在原假设成立的情况下，得到这样一个t值的概率很小。现在学习的是第76页，共115页 由上面的说明不难看出，假设检验可以说就是检验是否出现了小概率事件，如果出现小概率事件，则拒绝原来关于总体参数的假设；如果检验表明得到的样本值并不属于小概率事件，即若我们的假设成立，得到该样本值的概率不算小，则我们不能拒绝原来的假设，或者说，我们“接受”

31、原假设。问题是，我们上面提到的概率究竟应该小到什么程度才算小。一般说来，这取决于我们愿意承担的拒绝一个正确的假设和接受一个错误的假设这两方面的风险。在实践中，一般习惯于取5%作为拒绝假设的临界水平，称为5%的显著性水平。现在学习的是第77页，共115页假设检验的具体步骤是：（1）建立关于总体参数的原假设和备择假设；（2）计算检验统计量，检验原假设（是否出现小概率事件）；（3）得出关于原假设是否合理的结论。例例2.10 仍用食品需求的例子（例2.2）试检验原假设：。现在学习的是第78页，共115页原假设：H0：1=0.12备择假设：H1：10.12我们有：用=n-k-1=25-2-1=22查t表

32、，截断两侧5%面积的t临界值 tc =2.074 故拒绝原假设H0:。现在学习的是第79页，共115页三、系数的显著性检验 在假设检验中，有关斜率系数 是否为0 的假设检验特别重要。如果通过检验，接受 的原假设，则表明Xj和Y没有关系，即Xj对Y的变动没有影响。在这种情况下，可考虑从模型中剔除Xj。这类检验称为系数的显著性检验系数的显著性检验。1单个系数显著性检验 目的是检验某个解释变量的系数j是否为0，即该解释变量是否对因变量有影响。原假设 H0：j=0 备择假设 H1：j0现在学习的是第80页，共115页单个系数显著性检验的检验统计量是自由度为 n-k-1 的 t 统计量：t(n-k-1)

33、其中，为矩阵 主对角线上第 j+1个元素。而现在学习的是第81页，共115页 例例2.11 仍用食品需求的例子（例2.2），回归结果如下（括号中数字为标准误差）：试检验价格的系数的显著性。解：原假设 H0：备择假设 H1：查t表，故拒绝原假设H0。结论：显著异于0，P对Y有影响。现在学习的是第82页，共115页2若干个系数的显著性检验（联合假设检验）有时需要同时检验若干个系数是否为0，这可以通过建立单一的原假设来进行。设要检验g个系数是否为0，即与之相对应的g个解释变量对因变量是否有影响。不失一般性，可设原假设和备择假设为：H0:1=2=g=0 H1:H0不成立 (即X1,Xg中某些变量对Y有

34、 影响)现在学习的是第83页，共115页分析：这实际上相当于检验g个约束条件 1=0，2=0，g=0 是否同时成立。若H0为真，则正确的模型是：据此进行回归（有约束回归），得到残差平方和 SR是H0为真时的残差平方和。现在学习的是第84页，共115页若H1为真，正确的模型即原模型：据此进行无约束回归（全回归），得到残差平方和S是H1为真时的残差平方和。现在学习的是第85页，共115页 如果H0为真，则不管X1,Xg这g个变量是否包括在模型中，所得到的结果不会有显著差别，因此应该有：S SR如果H1为真，则由上一节中所讨论的残差平方和e2的特点，无约束回归增加了变量的个数，应有 S SR 通过检

35、验二者差异是否显著地大，就能检验原假设是否成立。现在学习的是第86页，共115页所使用的检验统计量是：F(g,n-k-1)其中，g为分子自由度，n-k-1为分母自由度。使用 的作用是消除具体问题中度量单位的影响，使计算出的 F 值是一个与度量单位无关的量。现在学习的是第87页，共115页例2.12 给定20组Y,X1,X2,X3的观测值，试检验模型 中X1和X3对Y是否有影响？解：（1）全回归 估计 得到：S=e2=25 （2）有约束回归 估计 得到：SR=e2=30现在学习的是第88页，共115页 原假设 H0:1=3=0 备择假设 H1:H0不成立 我们有：n=20,g=2,k=3 用自由

36、度（2，16）查F分布表，5%显著性水平下，F=1.6 FC=3.63,故接受H0。结论：X1和X3对Y无显著影响现在学习的是第89页，共115页3全部斜率系数为0的检验 上一段结果的一个特例是所有斜率系数均为0的检验，即回归方程的显著性检验：H0：1=2=K K=0 也就是说，所有解释变量对Y均无影响。注意到 g=K，则该检验的检验统计量为：现在学习的是第90页，共115页 分子分母均除以 ，有 从上式不难看出，全部斜率为0的检验实际是检验R2的值是否显著异于0，如果接受原假设，则表明因变量的行为完全归因于随机变化。若拒绝原假设，则表明所选择模型对因变量的行为能够提供某种程度的解释。现在学习

37、的是第91页，共115页四四检验其他形式的系数约束条件 上面所介绍的检验若干个系数显著性的方法，也可以应用于检验施加于系数的其他形式的约束条件，如 检验的方法仍是分别进行有约束回归和无约束回归，求出各自的残差平方和 SR 和 S，然后用 F 统计量进行检验。当然，单个系数的假设检验，如 H0：3=1.0，亦可用t检验统计量进行检验。现在学习的是第92页，共115页例2.13 Cobb-Douglas生产函数 Y=AKL 试根据美国制造业1899-1922年数据检验规模效益不变的约束：+=1解：（1）全回归 现在学习的是第93页，共115页（2）有约束回归：将约束条件代入，要回归的模型变为：Y=

38、AKL1-为避免回归系数的不一致问题，两边除以L，模型变换为：Y/L=A(K/L)回归，得：现在学习的是第94页，共115页 由回归结果得到的约束回归和全回归的残差平方和分别为 SR=0.0716 S=0.0710 （3）检验 原假设 H0:+1 备择假设 H1:+1 本例中，g=1,K=2,n=24 现在学习的是第95页，共115页 用自由度（1，21）查F表，5%显著性水平下，Fc=4.32 F=0.18 Fc=4.32 故接受原假设H0:+1 （4）结论 我们的数据支持规模收益不变的假设。现在学习的是第96页，共115页五、回归结果的提供和分析五、回归结果的提供和分析1.回归结果提供的格

39、式回归结果提供的格式 在论文、专著或报告中提供回归分析结果时一般应采用简洁而通行的格式，以便于交流。通行的格式有以下两种：（1）这里116.7、0.112和0.739分别为常数项和两个斜率系数的估计值，括号中提供的是 的标准误差。现在学习的是第97页，共115页（2）括号中数字分别是原假设 、和 成立时的t值。由此可见，这两种格式的唯一区别就在于括号中数字的含义不同。正因为如此，人们在论文或著作中提供回归结果时，必须在适当地方说明括号中数字是标准误差还是t值。需要说明的是，提供回归结果的标准格式中一般还包括检验一阶自相关的DW检验值，我们将在下一章“自相关”一节中介绍。现在学习的是第98页，共

40、115页2.回归结果的分析回归结果的分析结果的分析主要包括以下内容：（1）系数估计值。首先是分析系数的符号是否正确，系数值的大小是否恰当，是否符合理论预期和常识。上一段例中斜率系数一正一负，符合经济理论，数值大小也大致合理。（2）拟合情况。例中 很高，拟合较理想。（3）系数的显著性。例中斜率系数的t值分别为37.33和6.48，表明这些系数显著异于0，X和P对Y有影响。（4）根据DW检验值说明是否存在扰动项的自相关。如何说明，将在下一章中介绍。现在学习的是第99页，共115页第六节 预测 我们用OLS法对多元回归模型的参数进行了估计之后，如果结果理想，则可用估计好的模型进行预测。预测指的是对诸

41、自变量的某一组具体值 来预测与之相对应的因变量值 。当然，要进行预测，有一个假设前提应当满足，即拟合的模型在预测期也成立。现在学习的是第100页，共115页 点预测值由与给定的诸X值对应的回归值给出，即 而预测期的实际Y值由下式给出：其中u0是从预测期的扰动项分布中所取的值。现在学习的是第101页，共115页预测误差可定义为：两边取期望值，得因此，OLS预测量是一个无偏预测量。现在学习的是第102页，共115页 预测误差的方差为：从 e0 的定义可看出，e0 为正态变量的线性函数，因此，它本身也服从正态分布。故现在学习的是第103页，共115页由于 为未知，我们用其估计值代替它，有 则 的95

42、%置信区间为：即 现在学习的是第104页，共115页例例2.14 用例2.4的数据，预测X2=10，X3=10的Y值。解：由例2.4我们已得到：现在学习的是第105页，共115页因此 的95%置信区间为：或 3.66至24.34之间.现在学习的是第106页，共115页 第七节 虚拟变量（Dummy variables）一虚拟变量的概念 在回归分析中，常常碰到这样一种情况，即因变量的波动不仅依赖于那种能够很容易按某种尺度定量化的变量（如收入、产出、价格、身高、体重等），而且依赖于某些定性的变量（如性别、地区、季节等）。在经济系统中，许多变动是不能定量的。如政府的更迭（工党-保守党）、经济体制的改

43、革、固定汇率变为浮动汇率、从战时经济转为和平时期经济等。现在学习的是第107页，共115页 这样一些变动都可以用0-1变量来表示，用1表示具有某一“品质”或属性，用0表示不具有该“品质”或属性。这种变量在计量经济学中称为“虚虚拟拟变变量量”。虚拟变量使得我们可以将那些无法定量化的变量引入回归模型中。下面给出几个可以引入虚拟变量的例子。例1：你在研究学历和收入之间的关系，在你的样本中，既有女性又有男性，你打算研究在此关系中，性别是否会导致差别。现在学习的是第108页，共115页例2：你在研究某省家庭收入和支出的关系，采集的样本中既包括农村家庭，又包括城镇家庭，你打算研究二者的差别。例3：你在研究

44、通货膨胀的决定因素，在你的观测期中，有些年份政府实行了一项收入政策。你想检验该政策是 否对通货膨胀产生影响。上述各例都可以用两种方法来解决，一种解决方法是分别进行两类情况的回归，然后检验参数是否不同。另一种方法是用全部观测值作单一回归，将定性因素的影响用虚拟变量引入模型。现在学习的是第109页，共115页二虚拟变量的使用方法1 截距变动 设Y表示消费，X表示收入，我们有：假定不变。对于5年战争和5年和平时期的数据，我们可分别估计上述两个模型，一般将给出 的不同值。现引入虚拟变量D,将两式并为一式：其中，现在学习的是第110页，共115页 此式等价于下列两式：截距变动，斜率不变 在包含虚拟变量的

45、模型中，D的数据为0，0，0，0，0，1，1，1，1，1。估计结果如下图所示：应用t检验，2是否显著 可以表明截距项在两个时 期是否有变化。现在学习的是第111页，共115页2 斜率变动 如果我们认为战时和平时的消费函数中，截距项不变，而斜率不同，即变动，则可用下面的模型来研究两个时期边际消费倾向的差异：其中，D=不难看出，上式相当于下列两式：同样，包括虚拟变量的模型中，2是否显著可以表明斜率在两个时期是否变化。现在学习的是第112页，共115页3斜率和截距都变动在这种情况下，模型可设为：其中，D=此式等价于下列两个单独的回归式：引进了虚拟变量的回归模型对于检验两个时期中是否 发生结构性变化很

46、方便。如上例中，相当于检验 H0:2=4=0现在学习的是第113页，共115页4季节虚拟变量的使用 许多变量展示出季节性的变异(如商品零售额、电和天然气的消费等)，我们在建立模型时应考虑这一点，这有两种方法：（1）在估计前对数据进行季节调整；（2）采用虚拟变量将季节性差异反映在模型 中。现在学习的是第114页，共115页例2.15 设Y=购买汽车的实际支出额 X=实际总消费支出 用美国1973(1)-1980(2)的季度数据（按1975年价格计算），得回归结果如下：这一结果很不理想，低R2值，低t值，X的符号也不对。考虑到可能是季节性变异的问题，我们建立下面的模型：现在学习的是第115页，共115页