scipy.special.gdtrix#

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

逆 gdtr vs x。

返回关于参数 x 的逆，其中 p = gdtr(a, b, x) 是伽马分布的累积分布函数。这称为分布的 pth 分位数。

参数:

a类似数组: a 是 gdtr(a, b, x) 的参数值。 1/a 是伽马分布的“比例”参数。
b类似数组: b 是 gdtr(a, b, x) 的参数值。 b 是伽马分布的“形状”参数。
p类似数组: 概率值。
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