如何在 Python 中集成 beta 发行版

Posted 2023-04-15

【中文标题】如何在 Python 中集成 beta 发行版

【英文标题】：How to integrate beta distribution in Python

【发布时间】：2019-12-13 06:23:27

【问题描述】：

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

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

结果是：

0.9692858 absolute error < 6.6e-08

如何在 Python 中做到这一点？

【问题讨论】：

scipy.stats.beta(10,20).cdf(0.5)

@PaulPanzer 感谢您的回答。似乎 cdf 只接受一个参数（积分的上限）。如果我想像这个 scipy.stats.beta(10,20).cdf(0.2, 0.5) 一样将 0.2 到 0.5 积分，它会引发错误。另外，请回答问题，以便我确认。

【参考方案1】：

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

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