====== F분포 (F Distribution) ======
===== 정의 =====
$U \sim \chi^{2}(\nu_{1})$, $V \sim \chi^{2}(\nu_{2})$이고 $U$와 $V$가 서로 [[독립]]이면 $F=\frac{(U/\nu_{1})}{(V/\nu_{2})}$는 [[자유도]]가 $(\nu_{1}, \ \nu_{2})$인 [[F분포]]를 따른다.
===== 표기 =====
$$ X \sim F(\nu_{1}, \ \nu_{2}) $$
* $\nu_{1}$은 분자의 [[자유도]]
* $\nu_{2}$은 분모의 [[자유도]]
===== 받침 =====
$$ x \in [ \ 0 \ , \ \infty \ ) $$
===== 확률밀도함수 =====
$$ \begin{displaymath}\begin{split} f_{\nu_{1},\nu_{2}}(x) &= \frac{\Gamma \left[[ (\nu_{1} + \nu_{2})/2 \right]] \cdot \nu_{1}^{\nu_{1}/2} \cdot \nu_{2}^{\nu_{2}/2}}{\Gamma(\nu_{1}/2) \cdot \Gamma(\nu_{2}/2)} \frac{x^{n/2 - 1}}{(\nu_{2} + \nu_{1} x)^{(\nu_{1} + \nu_{2})/2}} \\ &= \frac{\nu_{1}^{\nu_{1} / 2} \cdot \nu_{2}^{\nu_{2} / 2} \cdot x^{n/2 - 1}}{(\nu_{2} + \nu_{1} x)^{(\nu_{1} + \nu_{2})/2} \cdot B(\frac{1}{2} \nu_{1},\frac{1}{2} \nu_{2})} \end{split}\end{displaymath} $$
단, $\Gamma(x)$는 [[감마함수]], $B(x,y)$는 [[베타함수]]
set title "F Distribution PDF"
set size 1.0
set xrange [0:4]
set yrange [0:0.8]
set format x "%.1f"
set format y "%.2f"
set xlabel "x"
set ylabel "f(x)"
Binv(p,q)=exp(lgamma(p+q)-lgamma(p)-lgamma(q))
f(x,df1,df2)=Binv(0.5*df1,0.5*df2)*(df1/df2)**(0.5*df1)*x**(0.5*df1-1.0)/(1.0+df1/df2*x)**(0.5*(df1+df2))
plot f(x,3,3) title "F(3,3)", \
f(x,7,4) title "F(7,4)", \
f(x,9,6) title "F(9,6)"
===== 누적분포함수 =====
$$ \begin{displaymath}\begin{split} F_{\nu_{1},\nu_{2}} (x) &= I \left( \frac{\nu_{1} x}{\nu_{2} + \nu_{1} x} \ ; \ \frac{1}{2} \nu_{1} \ , \ \frac{1}{2} \nu_{2} \right) \\ &= 2 \nu_{1}^{(\nu_{1} - 2)/2} \cdot \left( \frac{x}{\nu_{2}} \right)^{\nu_{1} / 2} \cdot \frac{_{2}F_{1} \left( \frac{1}{2} (\nu_{1} + \nu_{2}) \ , \ \frac{1}{2} \nu_{1} \ ; \ 1 + \frac{1}{2} \nu_{1} \ ; \ - \nu_{1} x / \nu_{2}) \right) }{B(\frac{1}{2} \nu_{1},\frac{1}{2} \nu_{2})} \end{split}\end{displaymath} $$
단, $I(x;a,b)$는 ??(FIXME)함수, $_{2}F_{1} (a,b;c;z)$는 [[초기하함수]]
set title "F Distribution CDF"
set size 1.0
set xrange [0:4]
set yrange [0:1.1]
set format x "%.1f"
set format y "%.2f"
set xlabel "x"
set ylabel "F(x)"
set key left
cf(x,df1,df2)=1.0-ibeta(0.5*df2,0.5*df1,df2/(df2+df1*x))
plot cf(x,3,3) title "F(3,3)", \
cf(x,7,4) title "F(7,4)", \
cf(x,9,6) title "F(9,6)"
===== 기대값 =====
$$ E(X) = \frac{\nu_{2}}{\nu_{2} - 2} $$
===== 분산 =====
$$ Var(X) = \frac{2 \nu_{2}^{ \ 2} (\nu_{1} + \nu_{2} - 2)}{\nu_{1} (\nu_{2} - 2)^{2} (\nu_{2} - 4)} $$
===== 왜도 =====
$$ \gamma_{1} = \frac{2(2 \nu_{1} + \nu_{2} - 2)}{\nu_{2} - 6} \sqrt{\frac{2(\nu_{2} - 4)}{\nu_{1} (\nu_{1} + \nu_{2} - 2)}} $$
===== 첨도 =====
$$ \gamma_{2} = \frac{12(-16 + 20 \nu_{2} - 8 \nu_{2}^{ \ 2} + \nu_{2}^{ \ 3} + 44 \nu_{1} - 32 \nu_{1} \nu_{2} + 5 \nu_{1} \nu_{2}^{ \ 2} - 22 \nu_{1}^{ \ 2} + 5 \nu_{1}^{ \ 2} \nu_{2})}{\nu_{1} (\nu_{2} - 6)(\nu_{2} - 8)(\nu_{1} + \nu_{2} - 2)} $$
===== 특징 =====
- $F_{\alpha} (\nu_{1}, \ \nu_{2}) = \frac{1}{F_{1-\alpha} (\nu_{2}, \ \nu_{1})} $인 관계가 성립한다.
- $F = \frac{(S_{1}^{ \ 2} / \sigma_{1}^{ \ 2})}{(S_{2}^{ \ 2} / \sigma_{2}^{ \ 2})} \sim F(n_{1} - 1, \ n_{2} - 1) $인 관계가 성립한다.
===== 타 분포와의 관계 =====
* [[t분포와 F분포 관계]]