scipy.special.gdtrix#

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

gdtr 的逆函数与 x 的关系。

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

参数:

aarray_like: a 参数值为 gdtr(a, b, x)。1/a 是伽马分布的“尺度”参数。
barray_like: b 参数的 gdtr(a, b, x) 的值。b 是伽马分布的“形状”参数。
parray_like: 概率值。
出ndarray，可选: 如果给出了第四个参数，它必须是一个大小与 a, b 和 x 的广播结果匹配的 numpy.ndarray。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