廣義逆高斯分佈#

機率密度函數由下式給出

\begin{eqnarray*} f(x; p, b) = x^{p-1} \exp(-b(x + 1/x)/2) / (2 K_p(b)), \end{eqnarray*}

其中 \(x > 0\) 是一個實數,參數 \(p, b\) 滿足 \(b > 0\)\(K_v\) 是第二類修正貝索函數,階數為 \(v\) (scipy.special.kv)。

如果 Xgeninvgauss(p, b),則 1/X 的分佈是 geninvgauss(-p, b)。逆高斯分佈 (scipy.stats.invgauss) 是一個特殊情況,其中 p=-1/2。

實作:scipy.stats.geninvgauss