Нецентральное F-распределение#

Распределение \(\left(X_{1}/X_{2}\right)\left(\nu_{2}/\nu_{1}\right)\) if \(X_{1}\) является нецентральным хи-квадрат с \(\nu_{1}\) степеней свободы и параметра \(\lambda\), и \(X_{2}\) имеет распределение хи-квадрат с \(\nu_{2}\) степени свободы.

Есть 3 параметра формы: степени свободы \(\nu_{1}>0\) и \(\nu_{2}>0\); и \(\lambda\geq 0\).

\begin{eqnarray*} f\left(x;\lambda,\nu_{1},\nu_{2}\right) & = & \exp\left[\frac{\lambda}{2} + \frac{\left(\lambda\nu_{1}x\right)} {2\left(\nu_{1}x+\nu_{2}\right)} \right] \nu_{1}^{\nu_{1}/2}\nu_{2}^{\nu_{2}/2}x^{\nu_{1}/2-1} \\ & & \times\left(\nu_{2}+\nu_{1}x\right)^{-\left(\nu_{1}+\nu_{2}\right)/2} \frac{\Gamma\left(\frac{\nu_{1}}{2}\right) \Gamma\left(1+\frac{\nu_{2}}{2}\right) L_{\nu_{2}/2}^{\nu_{1}/2-1} \left(-\frac{\lambda\nu_{1}x} {2\left(\nu_{1}x+\nu_{2}\right)}\right)} {B\left(\frac{\nu_{1}}{2},\frac{\nu_{2}}{2}\right) \Gamma\left(\frac{\nu_{1}+\nu_{2}}{2}\right)} \end{eqnarray*}

где \(L_{\nu_{2}/2}^{\nu_{1}/2-1}(x)\) является присоединённым полиномом Лагерра.

Если \(\lambda=0\), распределение становится эквивалентным распределению Фишера с \(\nu_{1}\) и \(\nu_{2}\) степени свободы.

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