Chen et al.(2006)年Journal of corporate finance 的一篇文章Ownership structure, corporate governance, and fraud中用bivariate probit 模型考察了违规行为的可能性以及违规行为被查可能性，但我们现在只能观察到被披露出来的财务违规事件，意味着我们知道z, z=Fj*Dj, 但现在需要同时得到对Fj和Dj的回归结果，并且运用bivariate probit模型（模型如下图）
请问 如何在stata里实现呢？  
假设我Fj 中对应的解释变量是firmsize, firmage, board size
假设我Dj 中对应的解释变量是roa, roe