scipy.special.betaincc#

scipy.special.betaincc(a, b, x, out=None) = <ufunc 'betaincc'>#

正则不完全 beta 函数的补集。

计算正则不完全 beta 函数的补集,定义为 [1]

\[\bar{I}_x(a, b) = 1 - I_x(a, b) = 1 - \frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)} \int_0^x t^{a-1}(1-t)^{b-1}dt,\]

对于 \(0 \leq x \leq 1\)

参数:
a, barray_like

正实值参数

xarray_like

实值,使得 \(0 \leq x \leq 1\),积分的上限

outndarray,可选

函数值的可选输出数组

返回值:
标量或 ndarray

正则不完全 beta 函数的值

另请参阅

betainc

正则不完全 beta 函数

betaincinv

正则不完全 beta 函数的逆

betainccinv

正则不完全 beta 函数补集的逆

beta

beta 函数

scipy.stats.beta

beta 分布

备注

在版本 1.11.0 中新增。

参考

[1]

NIST 数学函数数字图书馆 https://dlmf.nist.gov/8.17

示例

>>> from scipy.special import betaincc, betainc

betainc(a, b, x) 的值接近 1 时,朴素的计算 1 - betainc(a, b, x) 会损失精度

>>> 1 - betainc(0.5, 8, [0.9, 0.99, 0.999])
array([2.0574632e-09, 0.0000000e+00, 0.0000000e+00])

使用 betaincc,我们可得到正确值

>>> betaincc(0.5, 8, [0.9, 0.99, 0.999])
array([2.05746321e-09, 1.97259354e-17, 1.96467954e-25])