scipy.special.gdtrib#
- scipy.special.gdtrib(a, p, x, out=None) = <ufunc 'gdtrib'>#
gdtr
相对于 b 的逆。返加马分布的累积分布函数
p = gdtr(a, b, x)
的 b 参数的逆值。- 参数:
- a类数组
a 的参数值 gdtr(a, b, x)。 1/a 是伽玛分布的“规模”参数。
- p类数组
概率值。
- x类数组
非负实值,来自伽马分布的域。
- outndarray,可选
如果给出第四个参数,它必须是一个 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