如何在R中编写删失版EASI代码？先估计probit再设定EASI需求模型

• R语言中的easi包无法在EASI设定中估计probit阶段的概率密度函数（PDF）与累积分布函数（CDF），这会打乱输出构建与后处理的全流程。
• 已有可实现该功能的Stata代码，但无授权许可使用。当前需判断两种方案的成本效益：
  • 方案一：花费时间将Stata代码转换为R代码，但缺乏可对比的基准结果验证转换正确性
  • 方案二：直接购买Stata授权许可
• 尝试用ChatGPT转换核心的130行Stata程序/函数，结果不可靠且不稳定；且该代码依赖Stata底层Mata模块，与R不存在直接映射关系。

内容的提问来源于stack exchange，提问作者MarioF