工具变量法~[应用知识]

1、工具变量法一、工具变量法的主要思想在无限分布滞后模型中，为了估计回归系数，通常的做法是对回归系数作一些限制，从而对受限的无限分布滞后模型进行估计。在这里，考伊克模型、适应性期望模型与部分调整模型给出了很好的解决此类问题的思路。经过变换，新的模型中，随机扰动项的表达式为：考伊克模型： （ ，为衰减率） （1.1）；适应性期望模型：（ ，为期望系数）（1.2）；部分调整模型：（ ，为调整系数） （1.3）。 为原无限分布滞后模型中的扰动项，为变换后的扰动项。在原模型中的随机扰动项满足经典假设的前提下，部分调整模型也满足经典假设，但是考伊克模型与适应性期望模型的随机扰动项由于存在原随机扰动项的滞后项

2、，也就是说考伊克模型与适应性期望模型的解释变量 势必与误差项相关，因此，可能会出现上述两个模型的最小二乘估计甚至是有偏的这样严重的问题。那么，我们是否可以找到一个与高度相关但与不相关的变量来替代？在这里，一个可行的估计方法就是工具变量法。在讨论工具变量法之前，我们先来了解一下外生变量和内生变量。一般来说：一个回归模型中的解释变量有的与随机扰动项无关，我们称这样的解释变量为外生变量；而模型中有的解释变量与随机扰动项相关，我们可称这样的解释变量为内生解释变量。内生解释变量的典型情况之一就是滞后应变量为解释变量的情形，如上述考伊克模型与适应性期望模型中的。外生解释变量：回归模型中的解释变量与随机扰动

3、项无关；内生解释变量：回归模型中的解释变量与随机扰动项无关；了解了内生变量和外生变量的概念，我们接着讨论工具变量法的主要思想：工具变量法和普通最小二乘法是模型参数估计的两类重要方法，在多元线性回归模型中，如果出现解释变量与随机误差项相关（即出现内生变量）时，其回归系数的普通最小二乘估计是非一致的，这时就需要引入工具变量。工具变量，顾名思义是在模型估计过程中被作为工具使用，以替代模型中与随机误差性相关的随机解释变量（即内生变量）。满足条件：1）总体无关：工具变量与随机扰动项无关； 2）样本相关：工具变量必须与被它所代替的内生变量高度相关； 3）与模型中其他解释变量不相关，以避免出现多重共线性。做

4、了替代后，用普通最小二乘法即可得到原回归系数的一致估计量。二、工具变量法的基本原理我们分别从简单线性回归模型和多元线性回归模型两方面来具体分析工具变量法的基本原理：简单线性回归模型考虑简单线性回归模型 （2.1）其中为内生变量。则其正规方程为： （2.2）设回归模型中的解释变量与随机扰动项相关，则如前所述，普通最小二乘估计量是非一致的。现用一个工具变量来代替正规方程中的解释变量，其残差表达式不变。 （2.3）即： （2.4）解上述引入了工具变量后的正规方程可得斜率项系数的估计量为： （2.5）（2.5）式中小写字母代表相应大写字母的离差。该市所表示的估计量就是工具变量估计量，简称IV估计量，用

5、表示。易证IV估计量是一致估计量。事实上， （2.6） 若工具变量与解释变量高度相关，则表明（2.6）式中较大；若工具变量与随机扰动项渐近无关，则表明（2.6）式中随着样本容量的增加而趋向于零。故在工具变量与它相应的解释变量高度相关而与随机扰动项渐近无关的条件下，有 （2.7）样本估计总体，（2.7）表明IV估计量是一致估计量。多元线性回归模型：工具变量法可直接推广到多元线性回归模型 （2.8） 其中： 在讨论工具变量法在多元线性回归模型中的应用之前，我们先来分析工具变量的个数问题。为了一般起见，当解释变量与随机扰动项不相关时，我们把解释变量本身也作为是一个工具变量。这就是说，在我们的模型中凡

6、事预随机扰动项无关或渐近无关而与解释变量相关的变量都称为是工具变量。这样与随机扰动项无关的解释变量本身当然是与解释变量高度相关的变量，故它也是工具变量。在作了这样的约定之后对多元回归模型（2.8）来说，工具变量的个数一定不会小于解释变量（包括常数项）的个数（但可以大于解释变量的个数）。这是因为凡是与随机扰动项相关的解释变量都要有与随机扰动项无关或者渐近无关的工具变量或工具变量的线性组合，而凡与随机扰动项无关的解释变量本身就是一个工具变量（按我们上述约定）。所以工具变量的个数当然不小于解释变量（包含常数项）的个数。因此：工具变量的个数可以等于也可以大于（但不能小于）解释变量（包括常数项）的个数。

7、接下来的分析中，我们重点讨论工具变量的个数与解释变量的个数相等的情形：对于多元线性回归模型（2.8）：设解释变量的工具变量为，其中，而且某些变量可能与变量相同。若多元线性回归模型（2.8）的正规方程为： （2.9）则可通过将这个正规方程组中的解释变量换成其相应的工具变量，但残差的形式保持不变。得： (2.10)解方程组（2.10）可得回归系数的一致估计量即IV估计量。将 代入（2.10）式中，经整理得： （2.11）若令 ，,，则（2.11）式可以表示为： （2.12）若是一个满秩矩阵，则回归系数的IV估计量为： （2.13）当解释变量与随机扰动项不相关时，我们把这个解释变量本身也作为是一个工

8、具变量。把模型中凡是与随机扰动项无关或渐进无关而与解释变量相关的变量都称为解释变量。工具变量的个数一定不会小于解释变量（包括常数项）的个数。三、二阶段最小二乘法（TLS）IV估计量可以看作是两次运用最小二乘法的结果。第一阶段：求解释变量X对工具变量Z的回归，得到一个解释变量的拟合值。 （3.1）第二阶段：求应变量Y对解释变量的拟合值的回归，得到回归系数的估计值。 （3.2）注意到 为一对称的幂等矩阵，所以， （3.3）二阶段最小二乘法（TLS）简例 （3.4）、是内生变量，、是外生变量第一阶段：分别用内生解释变量对所有外生解释变量回归； （3.5） （3.6） 得、的拟合值、，也称为工具变量。

9、第二阶段：用上述拟合值代替实际值。 （3.7）我们也称，为工具变量。四、用工具变量法估计自回归模型由于考伊克模型与适应性期望模型均可化为一阶自回归模型。而自回归模型中由于含有滞后应变量作为解释变量，所以回归系数的最小二乘估计是非一致的，为了得到一致的估计可用工具变量法，但用什么变量作为工具变量也是很难作出决断的。例如对模型： （4.1）利维亚坦（Liviatan）建议用解释变量的滞后一个时期的值作为随机解释变量的工具变量。但由于大多数经济时间序列在相邻两期之间存在高度相关，从而使Liviatan的方法受到多重共线性的困扰。在这里我们看一个例子： （4.2）假定和相关。为了消除这种相关，假定我们

10、采取如下工具变量法：先求对和的回归，并从此回归得到估计值。然后做回归： （4.3）其中是从第一步回归估计处理的。在这个例子中利用两阶段最小二乘法消除了原模型中和之间的相关性。五、简述工具有效性的萨甘检验假定我们用一个工具变量来替代与误差项相关的自变量。那么工具变量又会多有效呢？也就是说我们如何知道所选的工具变量与误差项是独立的呢？萨甘提出了一个统计量（即SARG）来检验工具变量法中所使用的工具有效性。SARG所涉及的步骤如下：1、将回归方程中所包括的变量分成两组，一组是外生变量（称为），另一组是内生变量（称为）。2、选取为 1中的的变量工具，其中。3、用替代并估计原来的回归，也就是说，通过工具变量来估计原来的回归并得出残差值。4、将对一个常数项、所有的变量和所有的变量（但不包括变量）进行回归，从回归中得出。5、现在计算SARG统计量，它的定义为： （5.1）其中观测次数，原回归方程中的系数个数。在工具变量外生的虚拟假设下，萨甘证明了，SARG检验渐近服从自由度为的分布，其中指工具个数（即变量的个数），而指原回归方程中问题变量的个数。在一个具体应用中，如果计算出来的值统计显著，我们就拒绝工具的有效性。如果