随机解释变量

第八章 随机解释变量

教学目的及要求：

1、理解引入随机解释变量的目的及产生的影响 2、理解估计量的渐进无偏性和一致性 3、掌握随机解释变量OLS的估计特性 4、应用工具变量法解决随机解释变量问题

第一节 随机解释变量问题

一、随机解释变量问题产生的原因

多元（k）线性回归模型：

Yi??0??1X1i??2X2i??????kXki?Ui （8-1）

其矩阵形式为：

Y?XB?U （8-2） 在多元（k）线性回归模型中，我们曾经假定，解释变量Xj是非随机的。如果Xj是随机的，则与随机扰动项Ui不相关。即：

CovXij,Ui?0 (j?1,2,???,k;i?1,2,???,n) （8-3） 许多经济现象中，这种假定是不符合实际的，因为许多经济变量是不能用控制的方法进行观测的，所以作为模型中的解释变量其取值就不可能在重复抽样中得到相同和确定的数值，其取值很难精确控制，也不易用实验方法进行精确观测，解释变量成为随机变量。又由于随机项U包含了模型中略去的解释变量，而略去的解释变量往往是同模型中相关的变量，因而就很有可能在X是随机变量的情况下与随机项U相关，这样原有的古典假设就不能满足，产生随机解释变量。

在联立方程模型以及模型中包含有滞后内生变量等情况下，如果扰动项是序列相关的，那么均有扰动项和解释变量之间的相关性的出现，模型就存在随机解释变量问题。

?? 1

例如，固定资产投资与国民收入的关系满足如下模型：

It??0??1Yt??2It?1?ut

其中，It为t期的固定资产投资，It?1为t?1期的固定资产投资，Yt为t期的国民收入，因为It?1是随机变量，故模型中存在随机解释变量。

再如，消费与收入之间的影响关系模型为 Ct??0??1Yt??2Ct?1?ut

其中，Ct为t期的消费支出，Ct?1为t?1期的消费支出，Yt是t期的收入，因为Ct?1是随机变量，故模型中存在随机解释变量。

二、随机解释变量问题的后果

模型中，在解释变量为随机变量并且与扰动项相关的情况下，应用普通最小二乘法估计参数可能会出现估计的不一致性，使得估计值产生很大的偏误，造成拟合优度检验的全面失准，F检验失效，t检验失去意义。在这种情况下，各种统计检验得到的是虚假的结果，不能作为判别估计式优劣的依据。

随机解释变量带来何种结果取决于它与随机误差项是否相关： 1）随机解释变量与随机误差项不相关

2）随机解释变量与随机误差项在小样本下相关，在大样本下渐进无关 3）随机解释变量与随机误差项高度相关 4）滞后被解释变量与随机误差项相关

第二节 随机解释变量模型的估计特性

我们讨论的估计量的性质（包括无偏性、最小方差性）都是在样本容量一定的情况下的统计性质，在数理统计上叫做小样本性质。在某些情况下，小样本时的估计量不具有某种统计性质，但是随着样本容量的增大，一个估计量在小样本时不具有的性质，大样本时就逐渐具有这种统计性质了，这种性质我们叫做大样本性质或叫做估计量的渐近统计性质。常用的渐近统计性质有渐近无偏性和一致性。

2

一、估计量的渐近无偏性

?记?(n)代表模型中参数?的估计量，其上标n表示样本容量。一般来说，n取如下的样本容量，

?(n)为一随机变量。随着样本容量n的增大，估计量??(n)构成一个估计量（随机n1?n2???nk，?变量）序列：

????＝??(n)(n1)?，?(n2)?，…，?(nk)，… （8-4）

所谓渐近理论就是讨论当n变得很大时，以上这些序列会有怎样的结果。

?(n)如果满足： 序列?n????(n)?(n)）＝? (8-5) limE（??为?的渐近无偏估计。也就是说，当样本容量越来越大，n趋于?时，?(n)?则称?的均值越

来越接近参数的真值?。

这里需要注意的是，有些估计量在小样本下是有偏的，但在大样本下是无偏的，即是渐近无偏的。例如随机变量?的样本方差

1n S??(Xi?X)

ni?12x2容易证明（在数理统计中已有证明） E(Sx)??(1?2

221) n2其中，?为总体方差。很明显，在小样本下，Sx作为?的估计量是有偏的，但随着n的无限增大，E(Sx)趋于总体的真正方差?，因此是渐近无偏的。可见，通过增加样本容量，可以改善参数估计的精度。

222

二、估计量的一致性

?如果随着样本容量的增大，估计量??或称?(n)(n)?几乎处处趋近于?真值，我们说?(n)(n)为?的一致估计量，

(n)?依概率收敛于?。如果样本容量无限增大时，??的分布收敛于?，?的方差趋于零，

?(n)就是?的一致估计量。 ? 3

?一致估计量可以记为：Plim?n???(n)?(n)??。???1或简记为Plim?式中Plim表示概率极限。

n???n????? 为简单起见，可略去上标n，记作Plim?概率极限有下列运算法则：

Plim(cX)?cPlim(X) c为常数

Plim(c1X1?c2X2)?c1?PlimX1?c2?PlimX2 c1,c2为常数

Plim(X1?X2)?Plim(X1)?Plim(X2) Plim(X1Plim(X1))?,Plim(X2)?0 X2Plim(X2)?1Plim(X?1)??Plim(X)?

这里需要弄清楚一点是，无偏性与一致性是两个截然不同的概念，无偏性可以对任何样本容量成立，而一致性则是对大样本而言的，是一种渐近性质。在大样本的条件下，一致估计量具有很高的精度，但在小样本时一致性不起作用。

?可以证明，?n??(n)为?的一致估计量，当且仅当

??(n))?? limvar?limE(?((n))?0 （8-6）

n???是渐近无偏的，且当样本容量无限增大时??的方差趋于零。 时成立。此充分必要条件说明，?上面的讨论是对随机变量而言的，对于随机向量同样有类似的结论。

三、随机解释变量模型OLS估计特性

计量经济模型中一旦出现了随机解释变量，如果仍用最小二乘法估计模型参数，不同性质的随机解释变量会出现不同的结果。

为了简单起见，我们用一元线性回归模型进行说明。给定一元线性回归模型：

Yi??0??1Xi?Ui (i?1,2,...,n) （8-7）

假设?为一随机变量，模型满足其他古典假设条件。 对式（8-7），其离差形式为：

yi??1xi?ui （8-8）

4

其中，yi?Yi?Y ， xi?Xi?X ， ui???1?Ui?U

应用普通最小二乘法，则有

?xiyi2 （8-9）

?xi把(8-8)中的yi代入（8-9），则可以得到

??ii(?1xi?ui)xiui1??xiy?x2??xi?x2??1??i?x2 i而

E(?xiuiu1xnun?x2)?E(x1i?x2??x2u2i?x2???i?x2)

i

?E(x1x2xn?x2)E(u1)?E(2)E(u2)???E()E(un) i?xi?x2i下面分三种情况讨论：

1．?和U是独立的

E(??1)??1?E????xiui???? ??x2?i?因xi和ui相互独立，并且E(ui)?0 ∴E(?xiui?x2)?0

i故有E(??)?? 2．xi与ui小样本下相关，大样本下渐近无关 小样本：E(xiui)?0

所以E(??1)??1，最小二乘法估计是有偏的。 大样本：Plim1n??(n?xiui)?0 对式（8-10）两边取概率极限可有

1Plim(??x1)??1?Plim?xiui?n?iuix2??1?Plim i1n?x2i

5

8-10）

8-11）

8-12） （ （ （

因此，在假定Plim(12xi)?0的情况下，有 ?n?)?? （8-13）Plim(?

说明最小二乘估计式也具有一致性特性。 3．xi与ui高度相关

Plim(n??1xiui)?0 ?n讨论一般情况下回归模型（8-8）式

yi??1xi?ui (i?1,2,...... n) （8-14）

假设：Var(xi)??x，Var(ui)??u，xi和ui之间的相关系数是?，如果采用普通最小二乘法估计上式，可以得到：

221xiui?xu?iin?)???Plim??1? Plim(? 1121x?iPlim?xi2nPlim???可见，如果?很高，只有当

?Cov(x,u)????u （8-15）

?Var(x)x?u（8-15）式的渐近误差才是可以忽略的。否

?x是很小的情况下，

则，最小二乘估计式将存在着很大的偏误。

第三节 随机解释变量模型的处理

如果模型中存在随机解释变量问题，则一般的随机解释变量与随机误差项之间是相关的，最小二乘估计量有偏且不一致，需要利用其他估计方法对模型参数进行估计。

一、工具变量法

工具变量（Instrument Variable, IV）法就是当随机解释变量与随机误差项相关时，寻找一个与随机解释变量高度相关，但与随机误差项不相关的变量，用该变量替代模型中的随机解释变量，进行模型的参数估计。我们称这一替代随机解释变量的变量为工具变量。

6

（一）选择工具变量的要求

作为工具变量，必须满足以下四个条件：

第一，工具变量必须是有明确经济含义的外生变量；

第二，工具变量与其替代的随机解释变量高度相关，而又与随机误差项不相关； 第三，工具变量与模型中的其他解释变量也不相关，以免出现多重共线性； 第四，模型中的多个工具变量之间不相关。

（二）工具变量的应用

工具变量对随机解释变量的替代并不是“完全的”替代，即不是用工具变量代换模型中对应的随机解释变量，而是在最小二乘法的正规方程组中用工具变量对随机解释变量进行部分替代。

对于一元线性回归模型（8-7）和（8-8）yi??1xi?ui

若x与u不相关，u满足所有的统计假定。应用OLS法，利用微分求极值的办法求出正规方程：

2???xiyi?β1?xi （8-16） ???Yi?β0?β1Xi现采用另一种方法来导出OLS正规方程。我们以xi（i?1,2,???,n)同乘以yi??1xi?ui两边，得n个式子，求和得：

?xiyi??1?xi+?xiui （8-17）

2因为x与u不相关，从而可以略去如果x与u相关，则

?xuii?0，就可以得OLS正规方程。

?xuii?0，不能用OLS法来估计参数。现在，我们要寻找一个变量Z，

用zi的离差乘以yi??1xi?ui的两边，然后求和得到一个类似于OLSZ与X高度相关而与U无关，

正规方程的方程。在这里，Z就是工具变量。

由于z与u无关，所以

?zyii ??1?zixi??ziui （8-18）

1nPlim?ziui?0 n??ni?1得：

7

?zyii ??1?zixi （8-19）

上式称为拟正规方程，从而求得

?????1????zy?zxiiii(Z???(Zi?Z)(Yi?Y)?Z)(Xi?X) （8-20）

i??Y???X?01

因此，工具变量法的基本原理在于：用工具变量代替随机解释变量X，从而利用cov(Z,U)?0克服cov(X,U)?0产生的对模型参数估计的不利影响，形成有效正规方程组并最终获得模型参数的估计量。从这一原理理解，OLS法也可以看作是一种工具变量法，即利用模型中的各解释变量作为他们自身的工具变量。

容易证明，参数工具变量估计量是有偏的、一致的估计量。 在实际经济分析中，对于工具变量的选择，一般的做法是：

对于时间序列资料，如果被解释变量Yi、随机解释变量Xi、随机误差项ui三者之间的关系有

cov(Xi,ui)?0，但cov(Xi?1,ui)?0，cov(Yi?1,ui)?0，则可用Xi?1或Yi?1作为Xi的工具变量。

（三）多元线性回归模型 对于k元线性回归模型：

Yi??0??1X1i??2X2i??????kXki?ui （8-21）

其矩阵形式为：

Y?XB?U （8-22）

假设X1i和Xki为随机解释变量，且与随机误差项ui高度相关，ui满足最小二乘法的其他假定条件，解释变量之间无多重共线性。

（1）寻找工具变量Z1i和Zki。工具变量满足以下条件：他们是有实际经济意义的变量；与其对应的随机解释变量（Z1i对应X1i，Xki对应Zki）高度相关；与随机误差项ui不相关；工具变量Z1i和Zki之间不相关；与k元线性回归模型中其他解释变量不相关。

（2）写出工具变量矩阵。除了X1t和Xkt之外，k元线性回归模型的工具变量矩阵为：

8

?1X11?1X12X???......??1X1nX21X22...X2n.....Xk1?......Xk2?? （8-23） ......??...Xkn?将X矩阵中的X1i和Xki替换为Z1i和Zki，其他外生变量和常数项均由其自身做工具变量，得Z矩阵：

?1Z11?1Z12 Z???......??1Z1nZ21Z22...Z2n.....Zk1?......Zk2?? （8-24） ......??...Zkn??，沿上述思路，用Z?同乘Y?XB?U两边，由于Z和U无关，（3）求出工具变量估计量BIV

所以有

plim1Z?U?0 n??(Z?X)?1Z?Y （8-25）BIV

?是B的一致估计量。 BIV

在EViews软件中，工具变量法是含在二阶段最小二乘法中，所以必须选择二阶段最小二乘法，在“工具变量”的提示后面，输入所有的工具变量名，即可实现工具变量法估计。

（四）工具变量法的缺陷

从理论上分析，工具变量法可以得到渐近无偏、渐近有效的参数估计量，在解释变量为随机变量并与随机误差项相关的情况下，参数估计值达到了渐近一致。但这种方法在实际应用中会遇到一定困难，主要表现在三个方面：

（1）在解释变量X与随机误差项u相关的情况下，要找寻一个既与X高度相关，又与u不相关的工作变量Z十分困难。再加上工具变量Z要具有明确的经济含义，这就更不容易。

（2）在能找到符合要求的工具变量条件下，所选择的工具变量不同，模型参数估计值也不会一致，使参数估计出现随意性。工具变量选择得当，参数估计值的质量会高一点；如果工具变量选择不得当，参数估计值就会出现较大偏误。

（3）由于使用了工具变量，有可能产生较高的标准差，不能保证参数估计值的渐近方差一定能达到最小。

9

基于上述分析，到底选择何种工具变量，还需要通过实践掌握工具变量选择的技巧。

二、含有随机解释变量的实例分析

见教材及上机操作。

本章小结

解释变量?是随机变量，叫做随机解释变量，产生随机解释变量的原因主要包括：省略了解释变量，观测误差，滞后的内生变量。随机解释变量按其与随机扰动项的关系可以分为三种不同的情况，如果随机解释变量与随机扰动项相关，我们认为它不满足最小二乘法的古典假设。在三种不同的情况下，得到的最小二乘估计的性质有所不同，随机解释变量使模型参数估计值失去一致性。工具变量法是我们选择解决随机解释变量对模型参数估计影响的方法。

10

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