在 MATLAB 中，binofit 返回二项式概率分布和置信区间成功的最大似然估计值。
statsmodels.stats.proportion.proportion_confint 也返回置信区间，但找不到二项式概率分布的最大似然估计函数。在 MATLAB for Python 中，有什么函数可以建议作为binofit函数吗？
我认为您建议的功能已经足够好了。我运行了一些比较Matlab binofit和Python statsmodels.stats.proportion.proportion_confint的测试。该测试像测试一样，像 100K 实验一样，[phat,pci] = binofit(x,n,alpha)与min_conf,max_conf = proportion_confint(x,n,alpha=alpha,method='beta')。
对于 x 和 n 的值介于 0 和 10000 之间，来自 Matlab 和 Python 的 RMSE 十个置信区间限制低于 5e-6。使用 alpha=0.05 和 0.01 进行测试。
我知道这不是严格的演示，但对于我的项目，我决定将置信区间的两个估计值视为等效。
尝试使用以下两个库之一：statsmodels 或 scipy。
我不知道它是否正是您要找的，但我希望您仍然觉得它有用。