我现在需要做一个关于公司财务舞弊行为（0/1变量）的预测变量模型，根据附件这篇已发表的文献，现在这类研究需要采用bivariate probit regression中的partial observability model来处理。

关于partial observability model的处理，我已经能够在Stata中实现（简要说明参见附件Stata中biprobit的说明或下面这个页面）。无奈的是，正如Stata说明文档所说，这个函数经常会遭遇难以收敛的问题。
https://www.stata.com/support/faqs/statistics/bivariate-probit-with-partial-observability/

我了解到现在R语言计算功能强大，或许可以克服这个问题。我也在网上找到了声称可以解决bivariate probit regression的package（见附件），但是因为本人对于这块的统计原理和编程知识都很缺乏，所以看了package的说明书之后不知道如何迁移运用，实现原本已经用Stata可以实现的函数估计。

不知哪位高手有这方面专长，可否指点一二？无论能否成功，我都愿意以现金酬谢。