python - 如何在 Python 中集成 beta 发行版

Question

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

 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