回归模型是：Y=d+x+control+e     d是二值变量，也是要检测的变量，想看d对Y的影响
现在担心Y和d之间存在内生性问题，想用两阶段解决。
第一阶段是probit模型 d=X+control+u

我用了heckman y x c, select(d=x,c)

输出的结果并没有d的系数，但是需要报告d的系数和显著性，因为是要检测的变量。请问这种情况的内生性问题应该用什么方法检验呢？