T分布

スチューデントのt分布

確率密度関数
累積分布関数
母数	${\displaystyle \nu >0}$ 自由度
台	${\displaystyle (-\infty ,+\infty )}$
確率密度関数	${\displaystyle {\frac {\Gamma ({\frac {\nu +1}{2}})}{{\sqrt {\nu \pi }}\,\Gamma ({\frac {\nu }{2}})}}\left(1+{\frac {t^{2}}{\nu }}\right)^{-({\frac {\nu +1}{2}})}}$
累積分布関数	${\displaystyle {\begin{matrix}{\frac {1}{2}}+t\Gamma \left({\frac {\nu +1}{2}}\right)\cdot \\[0.5em]{\frac {{}_{2}\!F_{1}\left({\frac {1}{2}},{\frac {\nu +1}{2}};{\frac {3}{2}};-{\frac {t^{2}}{\nu }}\right)}{{\sqrt {\pi \nu }}\,\Gamma ({\frac {\nu }{2}})}}\end{matrix}}}$ ここで、2F1 は 超幾何関数
期待値	${\displaystyle 0}$（ただし ${\displaystyle \nu >1}$）
中央値	${\displaystyle 0}$
最頻値	${\displaystyle 0}$
分散	${\displaystyle \nu >2}$ の場合、${\displaystyle {\frac {\nu }{\nu -2}}}$ ${\displaystyle 1<\nu \leq 2}$ の場合、${\displaystyle \infty }$（無限大）
歪度	${\displaystyle 0}$（ただし ${\displaystyle \nu >3}$ の場合）
尖度	${\displaystyle {\frac {6}{\nu -4}}}$（ただし ${\displaystyle \nu >4}$ の場合）
エントロピー	${\displaystyle {\begin{matrix}{\frac {\nu +1}{2}}\left[\psi ({\frac {1+\nu }{2}})-\psi ({\frac {\nu }{2}})\right]\\[0.5em]+\log \left[{\sqrt {\nu }}B({\frac {\nu }{2}},{\frac {1}{2}})\right]\end{matrix}}}$ ${\displaystyle \psi }$：ディガンマ関数, ${\displaystyle B}$：ベータ関数
モーメント母関数	なし
特性関数	${\displaystyle {\frac {K_{\nu /2}({\sqrt {\nu }}|t|)({\sqrt {\nu }}|t|)^{\nu /2}}{\Gamma (\nu /2)2^{\nu /2-1}}}}$ ただし、${\displaystyle \nu >0}$ の場合。 ${\displaystyle K_{\nu }(x)}$ はベッセル関数[1]

統計学および確率論において、t分布（ティーぶんぷ、またはスチューデントのt分布、英: Student's t-distribution）は、連続確率分布の一つであり、正規分布する母集団の平均と分散が未知で標本サイズが小さい場合に平均を推定する問題に利用される。また、2つの平均値の差の統計的有意性を検討するt検定で利用される。t分布は、一般化双曲型分布の特別なケースである。

t分布は1908年にウィリアム・シーリー・ゴセットにより発表された。当時の彼はビール醸造会社であるギネスビールに雇用されており、ギネスビールでは秘密保持のため従業員による科学論文の公表を禁止していたので、彼はこの問題を回避するため「スチューデント」というペンネームを使用して論文を発表した[2]。

その後、ロナルド・フィッシャーがこの論文の重要性を見抜きスチューデントのt分布と呼んだため、このように呼ばれるようになった。

導出

X₁, …, X_n を平均 μ、分散 σ^2 の正規分布に従う独立な確率変数とする。また標本平均を

${\displaystyle {\overline {X}}={\frac {X_{1}+\cdots +X_{n}}{n}}}$

とし、不偏分散を

${\displaystyle S^{2}={\frac {1}{n-1}}\sum _{i=1}^{n}(X_{i}-{\overline {X}})^{2}}$

とする。ここで次の変数

${\displaystyle t={\frac {{\overline {X}}-\mu }{S/{\sqrt {n}}}}}$

を考えると、これは

${\displaystyle f(t)={\frac {\Gamma ((\nu +1)/2)}{{\sqrt {\nu \pi \,}}\,\Gamma (\nu /2)}}(1+t^{2}/\nu )^{-(\nu +1)/2}}$

（ただし ν = n − 1, Γ はガンマ関数）という確率密度関数に従うことが、ゴセットによって示された。ここで t の従う分布をt 分布（またはスチューデント分布）と呼ぶ。ν は自由度と呼ばれる。この分布は ν によるが、元の正規分布の母標準偏差σ にはよらないという重要な性質を持っている。

この確率密度関数は、元の正規分布の母数であるμおよびσが既知と仮定しているので、厳密には条件付確率密度関数${\displaystyle f(t\mid \mu ,\sigma )}$と書くべきものである。μおよびσを確率変数と考え、その確率密度関数を適当に仮定し(例えばテーブル状の一様分布関数)、ベイズの定理を適用することによって、標本平均${\displaystyle {\overline {X}}}$ および不偏標準偏差${\displaystyle S}$が既知の場合の条件付確率密度関数${\displaystyle f(t|{\overline {X}},S)}$を計算することができる(もう少し正確に言えば、まず条件付確率密度関数${\displaystyle f({\overline {X}},S\mid \mu ,\sigma )}$を求め、これにベイズの定理を適用して${\displaystyle f(\mu ,\sigma \mid {\overline {X}},S)}$を求め、さらにσについて積分して${\displaystyle f(t\mid {\overline {X}},S)}$を求める)。実はこの関数は${\displaystyle f(t\mid \mu ,\sigma )}$と全く同じ形をしている。つまり、

${\displaystyle {\begin{aligned}f(t\mid {\overline {X}},S)&={\frac {\Gamma ((\nu +1)/2)}{{\sqrt {\nu \pi \,}}\,\Gamma (\nu /2)}}(1+t^{2}/\nu )^{-(\nu +1)/2}\\t&={\frac {{\overline {X}}-\mu }{S/{\sqrt {n}}}}\end{aligned}}}$

である。これが、t分布が母標準偏差σ にはよらないという性質の反映である。不偏標準偏差${\displaystyle S}$は既知であるから、tの確率分布から母平均値μの確率分布を求めることができ、これを用いてμの区間推定や、仮説検定を行うことができる。

t分布を用いた母集団の平均値μの区間推定では、t=0について対称な区間で、その区間に亘る確率密度の積分値が95%となる区間（95%信頼区間）を考え、これに対応するμの区間を信頼区間 (CI) とする方法が広く用いられている（99%信頼区間を用いる場合も有る）。

t分布を用いた母集団の平均値μの仮説検定では、tの値が予め定めたα水準の下での信頼区間（95%あるいは99%）に含まれるか否かを判定基準とし、含まれる場合は母集団の平均値がμであるという仮説（帰無仮説）は棄却されず、区間からはみ出す場合は仮説を棄却する。

累積分布関数

累積分布関数は、正則不完全ベータ関数を用いて以下のように表される。

${\displaystyle \int _{-\infty }^{t}f(u)\,du=I_{x}\left({\frac {\nu }{2}},{\frac {\nu }{2}}\right)={\frac {B\left(x;{\frac {\nu }{2}},{\frac {\nu }{2}}\right)}{B\left({\frac {\nu }{2}},{\frac {\nu }{2}}\right)}}}$

ただし、

${\displaystyle x={\frac {t+{\sqrt {t^{2}+\nu }}}{2{\sqrt {t^{2}+\nu }}}}.}$

モーメント

t分布のモーメントは以下の式で表される。

• k が奇数の場合

${\displaystyle E(t^{k})={\begin{cases}0,&\quad 0<k<\nu \\{\mbox{undefined}},&\quad 0<\nu \leq k\end{cases}}}$

• k が偶数の場合

${\displaystyle E(t^{k})={\begin{cases}{\frac {\Gamma ({\frac {k+1}{2}})\Gamma ({\frac {\nu -k}{2}})\nu ^{k/2}}{{\sqrt {\pi }}\Gamma ({\frac {\nu }{2}})}},&\quad 0<k<\nu \\\infty ,&\quad 0<\nu \leq k\end{cases}}}$

特別なケース

ν の値により、簡単な形となる。

ν = 1 の場合

コーシー分布と一致する。

累積分布関数：

${\displaystyle F(t)={\frac {1}{2}}+{\frac {1}{\pi }}\arctan(t).}$

確率密度関数：

${\displaystyle f(t)={\frac {1}{\pi (1+t^{2})}}.}$

ν = 2 の場合

累積分布関数：

${\displaystyle F(t)={\frac {1}{2}}\left[1+{\frac {t}{\sqrt {2+t^{2}}}}\right].}$

確率密度関数：

${\displaystyle f(t)={\frac {1}{\left(2+t^{2}\right)^{3/2}}}.}$

ν → ∞ の場合

自由度 ν が ∞（無限大）に近づくにつれ、t分布は正規分布に近づく。

出典

1. Hurst, Simon, The Characteristic Function of the Student-t Distribution, Financial Mathematics Research Report No. FMRR006-95, Statistics Research Report No. SRR044-95
2. Walpole, Ronald; Myers, Raymond; Ye, Keying. Probability and Statistics for Engineers and Scientists. Pearson Education, 2002, 7th edition, pg. 237

参考文献

• 栗原伸一 (2011年7月25日). 入門統計学. オーム社. ISBN 978-4-274-06855-3 

関連項目

• 確率分布
• 非心t分布
• コーシー分布
• 正規分布