各位大神：
我参考一篇文献用Heckman两阶段模型做稳健性检验，解释变量是ESG评级分歧。文献中采用Heckman两阶段模型的原因是“鉴于ESG评级机构发布的ESG绩效评级结果并非包含中国A股的所有上市公司，对于未被纳入ESG评级的公司将不在回归样本中”，即将考虑到的是解释变量的缺失问题，将解释变量ESG分歧转换为虚拟变量进行probit回归。
文献如图。

我找到的相关Heckman两阶段模型的资料都是针对被解释变量y的缺失，将y转换为虚拟变量进行probit分析。
我自己尝试着写了一下stata命令，但是第二阶段回归的结果中没有x对y的系数。解释变量：ESG_Dis，生成的虚拟变量：ESG_Dis_dummy
被解释变量：COD，工具变量：ESG_range
控制变量：ESG_Ave  Agency_num $Controls
我首先尝试的是分步手工进行。

生成的结果如下：

我又尝试twostep的heckman模型
结果如下：


就是两种方法跑出来的结果都是没有主要解释变量ESG_Dis对于被解释变量COD的系数。
求求大神解答如何！！！谢谢！！！