最近也遇到同样的问题,我用Control Function Method解决的。建议了解一下Control Function Method,该方法针对内生变量是离散情形,可以实现第一阶段是probit, logit, tobit回归,第二阶段是OLS。 推荐文献:Control Function Methods in Applied Econometrics, Jeffrey M. Wooldridge ,2015.
此外,针对内生变量是离散情形,直接用2sls也是能得到一致估计的,只是有效性没有Control Function Method那么好。也就是说内生变量是离散情形也是可以用2sls进行估计。推荐文献:Instrumental Variables and the Search for Identification: From Supply and Demand to Natural Experiments,Angrist & Krueger,2001,第12页有提到。