分析101

工具变量原理

2021-04-25  本文已影响0人  Boye0212

在做回归时,很多时候会有\text{E}(x_t \varepsilon_t)\neq 0的情况,这也意味着不满足外生性条件\text{E}(\varepsilon|X)=0,此时的OLS估计量\hat\beta就不再满足无偏性,并且随着n的变大,它的bias也无法变小。若对此无法理解,请先掌握《小样本OLS回归梳理》中的内容。

此时该怎么办?一种解决方法是利用一些与\varepsilon无关的变量,这就是工具变量(instrumental variables,下文统称IV)。我们假设找到的IV是一些l\times 1的向量z_t,再将它排成n\times l的矩阵Z=[z_1,\cdots,z_n]'

IV需要与原来的x_t足够接近,因此Z'XXn\times k矩阵)必须满列秩。而我们寻找IV的目的,就是要让IV满足\text{E}(z_t\varepsilon_t)=0,由数据生成过程\varepsilon_t=y_t-x'_t\beta_o可知,我们要求解的就是满足\text{E}(z_t(y_t-x'_t\beta_o))=0\beta_o

我们无法知道\text{E}(z_t y_t)\text{E}(z_t x'_t),但可以用样本矩代替,即
n^{-1}\sum_{t=1}^{n}z_t(y_t-x'_t \beta_o)=Z'(y-X\beta_o)/n=0

上面的方程,若l \lt k,则有多个解,若l=kZ'X非奇异,则有唯一解\tilde \beta_n=(Z'X)^{-1}Z'y,若l \gt k,无解。在经济学理论中,往往会出现l \gt k的情形,此时尽管方程无解,但我们依旧可以寻找\beta_o,使Z'(y-X\beta_o)尽可能接近0

我们可以定义一个Z'(y-X\beta_o)0之间的二次距离:
d_n(\beta)=(Y-X\beta)'Z \hat{P}_n Z'(y-X\beta)
其中\hat{P}_n是一个l\times l的正定范数矩阵(positive definite norming matrix),它可以是随机矩阵。这里之所以选择二次距离,是因为这样在求解最优化问题时比较方便,可以直接写出一阶条件:
\dfrac{\partial d_n(\beta)}{\partial \beta} = -2X'Z\hat{P}_n Z'(y-X\beta)=0

假设X'Z\hat{P}_nZ'X非奇异,就可以得到IV估计量
\tilde \beta_n=(X'Z\hat{P}_nZ'X)^{-1}X'Z \hat{P}_nZ'y

只要选择Z\hat P_n,就可以得到各种计量经济学中的估计量。比如选择Z=X\hat P_n=(X'X/n)^{-1},那么\tilde \beta_n就变成了OLS估计量\hat \beta_n。而选择\hat P_n=(Z'Z/n)^{-1},就得到了2SLS(two-stage least squares)估计量。

IV估计量是无偏的吗?在数据生成过程y=X\beta_o+\varepsilon下,有
\begin{aligned} \tilde \beta_n=&(X'Z\hat{P}_nZ'X)^{-1}X'Z \hat{P}_nZ'y\\ =&(X'Z\hat{P}_nZ'X)^{-1}X'Z \hat{P}_nZ'(X\beta_o+\varepsilon)\\ =& \beta_o+(X'Z\hat{P}_nZ'X)^{-1}X'Z \hat{P}_nZ'\varepsilon \end{aligned}

事实上,上式的第二项我们没有理由保证它为0,哪怕有\text{E}(\varepsilon|Z)=0也无法保证。但在假设Z'\varepsilon /n \stackrel{a. s. }{\longrightarrow}0Z'X/n \stackrel{a. s. }{\longrightarrow} QQ为有限满列秩矩阵)以及\hat P_n\stackrel{a. s. }{\longrightarrow}PP为有限正定矩阵)之后,可以得到比无偏性更弱的一致性:\tilde\beta_n \stackrel{a. s. }{\longrightarrow} \beta_o

上一篇下一篇

猜你喜欢

热点阅读