scipy.special.gdtrib#
- scipy.special.gdtrib(a, p, x, out=None) = <ufunc 'gdtrib'>#
gdtr
的逆 vs b。返回相对于参数 b 的
p = gdtr(a, b, x)
的逆函数,即伽玛分布的累积分布函数。- 参数:
- aarray_like
a 参数值为 gdtr(a, b, x)。1/a 是伽马分布的“尺度”参数。
- parray_like
概率值。
- xarray_like
非负实数值,来自伽马分布的定义域。
- 出ndarray,可选
如果给出了第四个参数,它必须是一个 numpy.ndarray,其大小与 a, b 和 x 的广播结果匹配。out 则是函数返回的数组。
- 返回:
- b标量或ndarray
b 参数的值使得 p = gdtr(a, b, x)。b 是伽玛分布的“形状”参数。
注释
CDFLIB [1] Fortran 例程 cdfgam 的包装器。
累积分布函数 p 是使用 DiDinato 和 Morris [2] 的例程计算的。b 的计算涉及搜索一个值,该值产生所需的 p 值。搜索依赖于 p 随 b 的单调性。
参考文献
[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, gdtrib >>> p = gdtr(1.2, 3.4, 5.6) >>> print(p) 0.94378087442
验证逆命题。
>>> gdtrib(1.2, p, 5.6) 3.3999999999723882