如何在 Python 中集成 beta 发行版

Question

在 R 中，以下用于计算参数为 10 和 20 的 beta 分布上点 0 和 0.5 之间的积分：

 integrate(function(p) dbeta(p,10,20),0,0.5)

结果是：

0.9692858 absolute error < 6.6e-08

这如何在 Python 中完成？

Answer 1

您可以使用 的.cdf属性scipy.stats.beta。对于适当的间隔使用差异，例如

betacdf = scipy.stats.beta(10,20).cdf
betacdf(0.5)-betacdf(0.2)
# 0.9200223098258666