之前用的是Stata 软件，完全找不到解决问题的方向，有朋友建议用 Matlab，所以过来求救！还不知道是不是问对了地方，可以解决问题的话，愿意赠送我那为数不多的论坛币。问题如下：

劳动力市场Y中一共有J个部分Yj ，在任何一个部分Yj 中个体i的工资为yij。
由于个体的选择问题，观测结果存在可能选择性偏差。

  （1）
工资和个体选择的结果可以表示为（y, ys）。（u, us）服从二元正态分布，且相关系数为 ρ。
使用贝叶斯公式，可以计算得到j部分的工资分布如下：

(2）
其中φ和 Φ分别代表标准正态分布的密度和累积密度函数。
让P(i∈Yj)=πj 代表个体i属于j部分的可能性，那么整个劳动力市场的工资分布可以表示为：

（3）
请问：在使用probit得到个体选择的方程（1）之后，如何求方程（2）和（3）？？？是要编程才能计算得出来么？？？

想了很久了，没有出路，急用，，，请教各位，一点提示也可以，谢谢！！！！
请加QQ：157056527 注明人大经济论坛