Обобщенное обратное гауссовское распределение

Функция плотности вероятности задается как:

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

Реализация: scipy.stats.geninvgauss