求助：双重样本选择模型如何用biprobit在stata中实现 ，有两个选择过程，怎么分别计算逆米尔斯比？公示如下图
参考文献为国有部门和非国有部门中的性别工资差异—— 基于双重样本选择模型的经验研究
我编程了下，例如：
biprobit dwage jiuye education1 age married1 diqu1  mingzhu1 hangye1 zhiye1
predict c1,xb1

. predict c2,xb2
gen b=(1-0.87^2)^(1/2) 其中相关系数0.87已知

gen a1=(c2-0.87*c1)/b

. gen t1=normalden(c1)*normal(a1)*invF(c1,c2,0.87)
(4130 missing values generated) ————————这里有很多缺失值，不知道怎么回事，为什么F的逆函数里计算出来有很多缺失值，到了最后两个修正项还直接一个值都没有了。。求大仙指点，是不是我在预测的时候代码错了啊，就是这个predict c1,xb1

. predict c2,xb2？？在二元probit里应该怎么预测因变量dwage， jiuye的值呢