再次感谢lib_321老师

继续找啊。
或者找两个包，一个能做IV probit的，另一个能做weak iv检验的，比如ivpack里的anderson.rubi.ci那个函数，然后改函数代码，但需要统计高手帮忙弄清检验的原理。

Lib_321老师太给力了，感谢

#ivprobit()的源码看了下：
第一步是用内生变量对工具变量做lm回归；
第二步是用Y对工具变量和第一步的残差做probit回归，计算vcov矩阵；
第三步是用Y对外生变量和第一步的残差和拟合值做probit回归，比较第二步和第三步的系数：
第四步，用系数的差对内生变量处理，然后对工具变量做lm回归，计算vcov，与第二步的vcov相加；
最后计算系数、标准误、y对外生变量和内生变量拟合值作lm回归。
代码不复杂，调用的都是base、stats基础包里的函数，矩阵变化那里还没有想明白。只是输出结果太简单了。我感觉在第一步加上弱工具变量检验的代码是可行的。查了些弱工具变量检验的资料，没什么头绪，感觉把这个搞透的话可以写个论文啦。

lib_321老师，请教您，怎么搜索R在github里的包？谢谢

我一般在https://github.com/里右上角有个搜索栏，比如找ivmodel这个包，输入ivmodel language:R就可以
我试过：
iv probit language:R

weak iv language:R

感谢lib_321老师，你的回答给我帮了大忙

同样R控，请问你们这个问题解决了吗？

看了一下说明还是没搞明白，哪个位置放哪种变量。

Usage
ivprobit(formula, data)

Arguments
    formula y~x|y1|x2 whre y is the dichotomous l.h.s.,x is the r.h.s.   
            exogenous variables,y1 is the r.h.s. endogenous variables and
            x2 is the complete set of instruments
    data    the dataframe
然后它显示了相应的示例：

data(eco)

pro<-ivprobit(d2~ltass+roe+div|eqrat+bonus|ltass+roe+div+gap+cfa,eco)

summary(pro)
如果我符合指令的解释，

y= d2 = dichotomous l.h.s.
x= ltass+roe+div = the r.h.s. exogenous variables
y1= eqrat+bonus = the r.h.s. endogenous variables
x2= tass+roe+div+gap+cfa = the complete set of instruments