scipy.special.gdtrix#
- scipy.special.gdtrix(a, b, p, out=None) = <ufunc 'gdtrix'>#
相對於 x 的 gdtr 反函數。
傳回相對於參數 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
機率值。
- outndarray, optional
如果給定第四個參數,則它必須是一個 numpy.ndarray,其大小與 a、b 和 x 的廣播結果相符。out 接著會是此函數傳回的陣列。
- 傳回:
- x純量或 ndarray
參數 x 的值,使得 p = gdtr(a, b, x)。
註解
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