python 工具变量回归_stata工具变量法一例:使用2SLS进行ivreg2估计及其检验

7e944b4db355d1e10b4c7fc4433681c4.png

作为OLS回归不符合假定的问题,还包括解释变量与随机扰动项不相关。如果出现了违反该假设的问题,就需要找一个和解释变量高度相关的、同时和随机扰动项不相关的变量,作为工具变量进行回归。工具变量通常采用二阶段最小二乘法(2SLS)进行回归,当随机扰动项存在异方差或自相关的问题,2SLS就不是有效率的,就需要用GMM等方法进行估计,除此之外还需要对工具变量的弱工具性和内生性进行检验。

1

数据与要求

以stata自带的auto.dta数据为例,在stata输入如下命令,即可得到:

sysuse auto

数据展示如下:

3b7b9169d3179f8fce2992e978a758a0.png

数据为美国 1978 年汽车数据,包括产地、车名、行使里程、重量等变量

2

构造结构方程与2SLS估计

构造如下工具变量结构方程:

  • 5
    点赞
  • 23
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
假设你要估计的模型是: $y_{it} = \beta x_{it} + \alpha_i + \delta_t + \epsilon_{it}$ 其中,$y_{it}$ 是因变量,$x_{it}$ 是感兴趣的内生变量,$\alpha_i$ 和 $\delta_t$ 分别是个体固定效应和时间固定效应,$\epsilon_{it}$ 是误差项。 使用核心解释变量的滞后一期作为工具变量进行估计,首先需要创建工具变量。以下是Stata代码: ``` // 生成工具变量 gen x_lag1 = L.x gen iv = L.x_lag1 ``` 其中,L.x 表示 $x_{it}$ 的一期滞后值,即 $x_{i,t-1}$,x_lag1 是生成的变量名称;L.x_lag1 表示 $x_{it}$ 的两期滞后值,即 $x_{i,t-2}$,iv 是生成的工具变量名称。 接下来,进行2SLS估计。以下是Stata代码: ``` // 估计第一阶段回归,得到x_lag1对iv的影响 reg x_lag1 iv // 得到x_lag1的预测值,即x_lag1的影响 predict x_lag1_hat, xb // 估计第二阶段回归使用x_lag1_hat代替x_lag1 ivreg y x (iv = x_lag1_hat), vce(cluster id) ``` 其中,第一阶段回归使用iv作为内生变量,x_lag1作为工具变量进行估计,得到iv对x_lag1的影响;第二阶段回归使用x_lag1_hat代替x_lag1,使用iv作为工具变量进行估计,得到感兴趣的内生变量x对因变量y的影响。 上述代码中,vce(cluster id) 表示进行聚类稳健标准误估计,其中id为个体或地理区域的唯一标识符。 如果你需要使用空间滞后模型进行估计,可以将空间滞后项加入到模型中。以下是修改后的代码: ``` // 生成工具变量 gen x_lag1 = L.x gen iv = L.x_lag1 // 生成空间滞后项 spmat W = r(spatialweightmatrix) gen w_y_lag1 = W*y_lag1 gen w_x_lag1 = W*x_lag1 // 估计第一阶段回归,得到x_lag1对iv的影响 reg x_lag1 iv w_x_lag1 // 得到x_lag1的预测值,即x_lag1的影响 predict x_lag1_hat, xb // 估计第二阶段回归使用x_lag1_hat代替x_lag1 ivreg y x w_y_lag1, (iv = x_lag1_hat), vce(cluster id) ``` 其中,W 是空间权重矩阵,y_lag1 是因变量的一期滞后值,w_y_lag1 是 y_lag1 的空间滞后项,w_x_lag1 是 x_lag1 的空间滞后项。 注意,在使用空间滞后模型进行估计时,需要进行聚类稳健标准误估计,以消除空间自相关带来的偏误。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值