内生性与工具变量估计方法

  • 格式:doc
  • 大小:302.00 KB
  • 文档页数:8

内生性与工具变量估计方法
一 一元模型的IV 估计
采用MROZ 数据,进行练习。

估计教育对工资收入的回报:
01log()wage educ ββμ
=++
为了便于比较首先得到OLS 估计结果,在命令窗口输入
smpl 1 428
equation eq01.ls log(wage) c educ
教育的系数估计值表明,每多接受一年教育可得到月11%的回报。

接下来,我们用父亲的受教育程度(fatheduc )作为educ 的工具变量。

我们必须认为fatheduc 与u 不相关;第二个要求是educ 与fatheduc 相关。

为了验证第二点,作一个educ 对fatheduc 的回归。

equation eq02.ls educ c fatheduc
可以看出,educ 与fatheduc 之间存在统计显著的正相关。

采用fatheduc 作为educ 的工具变量,进行工具变量回归。

equation eq03.tsls log(wage) c educ @ fatheduc
IV 估计量的标准误是OLS 标准误的2.5倍,这在我们的意料之中。

二 多元模型的IV 估计 采用card 数据,进行练习。

估计教育对工资收入的回报:
012log()var wage educ Control iables βββμ
=+++
为了便于对照,先做OLS 回归 Smpl 1 3010
Equation eq01.ls log(wage) c educ exper expersq black smsa south smsa66 reg662 reg663 reg664 reg665 reg666 reg667 reg668 reg669
在这个例子中,受教育程度的工具变量是标志着一个人是否在一所四年制大学附近成长的虚拟变量(nearc4)。

为了验证受教育程度与该虚拟变量的偏相关性,先做educ对nearc4以及其他所有外生变量的回归:
Equation eq02.ls educ c nearc4 exper expersq black smsa south smsa66 reg662 reg663 reg664 reg665 reg666 reg667 reg668 reg669
Nearc4的系数估计值意味着,在其他因素固定的情况下,曾住在大学附近的人所受的教育比不在大学附近长大的人平均多出约1/3年。

我们接下来进行工具变量回归,以nearc4作为educ的IV。

Equation eq03.tsls log(wage) c educ exper expersq black smsa south smsa66 reg662 reg663 reg664 reg665 reg666 reg667 reg668 reg669 @ nearc4 exper expersq black smsa south smsa66 reg662 reg663 reg664 reg665 reg666 reg667 reg668 reg669
三 多个工具变量的回归方法 重新回到MROZ 数据,进行练习。

估计教育对工资收入的回报:
012log()var wage educ control iables βββμ
=+++
我们认为,受教育程度educ 具有内生性,并且采用父亲和母亲的受教育程度fatheduc 、motheduc 来作为它的工具变量。

为了便于比较首先得到OLS 估计结果,在命令窗口输入
smpl 1 428
equation eq04.ls log(wage) c educ exper expersq
我们检验educ与fatheduc、motheduc存在偏相关关系,做educ对fatheduc、motheduc以及所有其他外生变量的回归:
smpl 1 428
Equation eq05.ls educ c exper expersq fatheduc motheduc
结果显示,educ与fatheduc、motheduc高度相关。

随后,我们进行IV估计:
(1)直接进行工具变量回归
采用父亲和母亲的受教育程度fatheduc、motheduc来作为educ的工具变量。

smpl 1 428
Equation eq06.tsls log(wage) c educ exper expersq @ fatheduc motheduc exper expersq
(2)两步回归法
先用educ对fatheduc、motheduc以及所有其他外生变量OLS回归,得到educ的拟合值;然后将log(wage)对educ的拟合值以及其他外生变量进行OLS回归。

在eq05中,我们已经做了第一步的回归,因此只需要得到educ的拟合值:
eq05.fit educf
第二步如下
Equation eq07.ls log(wage) c educf exper expersq
(3)另一种方法
先用educ对fatheduc、motheduc以及所有其他外生变量OLS回归,得到educ的拟合值;然后将educ的拟合值作为educ的IV,进行工具变量回归。

第一步的工作前面已经完成,只需要做第二步的工作
Equation eq08.tsls log(wage) c educ exper expersq @ educf exper expersq
对比三种操作方法的结果,可以看出,三种方法所得到的回归系数估计值都完全相同;但是,第二种方法所得到的估计系数标准误与其他两种方法的结果不同,这个计算得到的标准误和t检验统计量是不正确的。

因此,我们应该避免采用第2种方法,而采用其他两种方法。

四检验解释变量的内生性
MROZ数据,我们检验educ的内生性。

第一步,将educ对所有外生变量回归,这个工作已经在eq05中完成,然后我们生成该方程的残差。

eq05.makeresid resid_eq05
第二步,在原结构方程中加入上述残差,用OLS回归检验该残差的显著性。

如果残差系数统计显著异于0,则断定educ是内生的。

Equation eq09.ls log(wage) c educ exper expersq resid_eq05
五检验工具变量的有效性(与扰动项不相关)
过度识别约束的检验
第一步,用2sls估计结构方程,得到残差。

第二步,将残差对所有外省变量回归,获得R2。

第三步,在所有IV都与扰动项不相关的零假设之下nR2服从自由度为q的卡方分布,其中q为模型之外的工具变量个数减去内生解释变量个数。

如果nR2足够大以至于超过临界值,则拒绝原假设,认为至少部分IV不是外生的。

仍然以教育程度的工资回报mroz为例子。

第一步,工具变量回归,这个工作已经在eq06中完成,然后我们生成该方程的残差。

eq06.makeresid resid_eq06
第二步,用上述残差对所有外生变量回归:
Equation eq10.ls resid_eq06 c exper expersq fatheduc motheduc
第三步,计算卡方统计量的具体值:
scalar overidentificationtest=428*eq10.@r2
计算结果为0.378,这在2(1)
中是一个非常小的值,因此不拒绝原假设。