广义逆高斯分布

概率密度函数由下式给出:

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

where \(x > 0\) is a real number and the parameters \(p, b\) satisfy \(b > 0\). \(K_v\) is the modified Bessel function of second kind of order \(v\) (scipy.special.kv).

如果 Xgeninvgauss(p, b) ,然后分配 1/Xgeninvgauss(-p, b) 。逆高斯分布 (scipy.stats.invgauss )是p=-1/2的特例。

实施: scipy.stats.geninvgauss