scipy.special.gdtrix#

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

相对于 x 的 gdtr 的逆。

返回 p = gdtr(a, b, x) 相对于参数 x 的逆，它是伽马分布的累积分布函数。这也被称为该分布的第 p 个分位数。

参数:

aarray_like: gdtr(a, b, x) 的 a 参数值。1/a 是伽马分布的“尺度”参数。
barray_like: gdtr(a, b, x) 的 b 参数值。b 是伽马分布的“形状”参数。
parray_like: 概率值。
outndarray, 可选: 如果提供了第四个参数，它必须是一个 numpy.ndarray，其大小与 a、b 和 x 的广播结果匹配。out 则是函数返回的数组。

返回:

另请参阅

备注

CDFLIB [1] Fortran 例程 cdfgam 的包装器。

累积分布函数 p 是使用 DiDinato 和 Morris [2] 的例程计算的。x 的计算涉及搜索一个能产生所需 p 值的数值。该搜索依赖于 p 随 x 单调变化的特性。

参考文献

[1]

Barry Brown, James Lovato, 和 Kathy Russell, CDFLIB: 累积分布函数、逆和其他参数的 Fortran 例程库。

[2]

DiDinato, A. R. 和 Morris, A. H., 不完全伽马函数比及其逆的计算。ACM Trans. Math. Softw. 12 (1986), 377-393。

示例

首先评估 gdtr。

>>> from scipy.special import gdtr, gdtrix
>>> p = gdtr(1.2, 3.4, 5.6)
>>> print(p)
0.94378087442

验证逆。

>>> gdtrix(1.2, 3.4, p)
5.5999999999999996