如何在 Python 中集成 Beta 分布

How to integrate beta distribution in Python

在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 中完成?

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

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