广义逆高斯分布#

概率密度函数由下式给出

\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).

如果 X 是 geninvgauss(p, b)，那么 1/X 的分布是 geninvgauss(-p, b)。逆高斯分布 (scipy.stats.invgauss) 是 p=-1/2 的特例。