我想在我的研究中引入“仪器变量”。
在许多情况下，结果（因变量； DV）是连续变量。
在连续结果（DV）的情况下，以下示例
stage1 <- lm(x2~x1+z2a, data)
data$x2hat <- fitted(stage1)
stage2 <- lm(y~x1+x2hat, data)
coeftest(stage2)
这是2SLS算法，与以下相同。
library(AER)
tsls <- ivreg(y~x1+x2 | x1+z2a, data)
coeftest(tsls)
但是，我的研究设计具有二进制结果（DV）。
R中的“ ivprobit”软件包与此分析有关。
另一方面，我对此有一些疑问。
以下是程序包中的示例代码（与程序包文档略有不同）。
library(ivprobit)
data(eco)
head(eco)
pro <- ivprobit(d2~ltass+roe+div | eqrat| ltass+roe+div+gap, eco)
summary(pro)
由于软件包文档中有一些描述，因此我很难理解这种机制。
“ gap”是工具变量，“ eqrat”是外生变量？
此外，像连续结果（DV）一样，我可以使用glm函数执行此分析吗？
这是我的意见。
model1 <- lm(eqrat~ltass+roe+div+gap, eco)
eco$hat <- fitted(model1)
model2 <- glm(d2~ltass+roe+div+hat, eco, family="binomial")
summary(model2)
谢谢。