司徒頓t分布

司徒頓t分布（Student's t-distribution），簡稱t 分布，在機率論及統計學中用於根據小樣本來估計母體呈常態分布且標準差未知的期望值。若母體標準差已知，或是樣本數足夠大時（依據中央極限定理漸進常態分布），則應使用常態分布來進行估計。其為對兩個樣本期望值差異進行顯著性測試的司徒頓t檢定之基礎。

司徒頓t 分布

機率密度函數
累積分布函數
母數	$\nu >0\!$自由度
值域	$x\in (-\infty ;+\infty )\!$
機率密度函數	${\frac {\Gamma ((\nu +1)/2)}{{\sqrt {\nu \pi }}\,\Gamma (\nu /2)\,(1+x^{2}/\nu )^{(\nu +1)/2}}}\!$
累積分布函數	${\frac {1}{2}}+{\frac {x\Gamma \left((\nu +1)/2\right)\,_{2}F_{1}\left({\frac {1}{2}},(\nu +1)/2;{\frac {3}{2}};-{\frac {x^{2}}{\nu }}\right)}{{\sqrt {\pi \nu }}\,\Gamma (\nu /2)}}$其中：$\,_{2}F_{1}$是超幾何函數
期望值	$\nu >1$時為$0$，$\nu =1$時未定義
中位數	$0$
眾數	$0$
變異數	$\nu >2$時為${\frac {\nu }{\nu -2}}\!$，否則為無窮大
偏度	$\nu >3$時為$0$
峰度	$\nu >4$時為${\frac {6}{\nu -4}}\!$
熵	${\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}}$ $\psi$: 雙Γ函數, $B$: 貝塔函數
動差母函數	未定義
特徵函數	${\frac {K_{\nu /2}({\sqrt {\nu }}|t|)({\sqrt {\nu }}|t|)^{\nu /2}}{\Gamma (\nu /2)2^{\nu /2-1}}},\;\nu >0$ $K_{\nu }(x)$: 第二類修正貝塞爾函數

司徒頓t 檢定改進了Z檢定（Z-test），因為在小樣本中，Z檢定以母體標準差已知為前提，Z檢定用在小樣本會產生很大的誤差，因此必須改用司徒頓t 檢定以求準確。但若在樣本數足夠大（普遍認為超過30個即足夠）時，可依據中央極限定理近似常態分布，以Z檢定來求得近似值，

在母體標準差數未知的情況下，不論樣本數量大或小皆可應用t檢定。在待比較的數據有三組以上時，因為誤差無法被壓低，此時可以用變異數分析（ANOVA）代替t檢定。

t 分布的推導最早由德國大地測量學家弗里德里希·羅伯特·赫爾默特於1876年提出，並由德國數學家雅各布·魯洛斯證明。

英國人威廉·高斯特於1908年再次發現並發表了t分布，當時他還在愛爾蘭都柏林的吉尼斯啤酒釀酒廠工作。酒廠雖然禁止員工發表一切與釀酒研究有關的成果，但允許他在不提到釀酒的前提下，以筆名發表t 分布的發現，所以論文使用了「司徒頓」（Student）這一筆名。之後t檢定以及相關理論經由羅納德·費雪發揚光大，為了感謝高斯特的功勞，費雪將此分布命名為司徒頓t 分布（Student's t）。

描述

假設$X$是呈常態分布的獨立的隨機變數（隨機變數的期望值為$\mu$，母體變異數為$\sigma ^{2}$但其值未知）。 令：

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

為樣本期望值，

${S_{n}}^{2}={\frac {1}{n-1}}\sum _{i=1}^{n}\left(X_{i}-{\overline {X}}_{n}\right)^{2}$

為樣本變異數，

$Z={\frac {{\overline {X}}_{n}-\mu }{\frac {\sigma }{\sqrt {n}}}}$

為呈期望值為0變異數為1的常態分布的隨機變數，

但因母體變異數$\sigma ^{2}$為未知，因此依史拉斯基定理以${S_{n}}^{2}$替換之：

$T={\frac {{\overline {X}}_{n}-\mu }{\frac {S_{n}}{\sqrt {n}}}}$

T 的機率密度函數是：

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

$\nu$ 等於n − 1。 T的分布稱為t 分布。母數$\nu$ 一般被稱為自由度。

$\Gamma$ 是伽瑪函數。 如果$\nu$是偶數,

${\frac {\Gamma ({\frac {\nu +1}{2}})}{{\sqrt {\nu \pi }}\,\Gamma ({\frac {\nu }{2}})}}={\frac {(\nu -1)(\nu -3)\cdots 5\cdot 3}{2{\sqrt {\nu }}(\nu -2)(\nu -4)\cdots 4\cdot 2\,}}\cdot$

如果$\nu$是奇數,

${\frac {\Gamma ({\frac {\nu +1}{2}})}{{\sqrt {\nu \pi }}\,\Gamma ({\frac {\nu }{2}})}}={\frac {(\nu -1)(\nu -3)\cdots 4\cdot 2}{\pi {\sqrt {\nu }}(\nu -2)(\nu -4)\cdots 5\cdot 3\,}}\cdot \!$

T 的機率密度函數的形狀類似於期望值為0變異數為1的常態分布，但更低更寬。隨著自由度$\nu$的增加，則越來越接近期望值為0變異數為1的常態分布。

t 分布密度 (紅色曲線) 在自由度為 1, 2, 3, 5, 10, 30比較於標準常態分布(藍色曲線).
前幅圖用綠色曲線表示.

1 degree of freedom	2 degrees of freedom	3 degrees of freedom
5 degrees of freedom	10 degrees of freedom	30 degrees of freedom

T分布的機率累計函數，用不完全貝塔函數I表示：

$F(t)=\int _{-\infty }^{t}f(u)\,du=1-{\tfrac {1}{2}}I_{x(t)}\left({\tfrac {\nu }{2}},{\tfrac {1}{2}}\right),$

其中

$x(t)={\frac {\nu }{t^{2}+\nu }}.$

T分布的動差為：

$E(T^{k})={\begin{cases}0&{\mbox{k odd}},0<k<\nu \\{\frac {\Gamma ({\frac {k+1}{2}})\Gamma ({\frac {n-k}{2}})^{k/2}}{{\sqrt {\pi }}\Gamma ({\frac {n}{2}})}}&{\mbox{k even}},0<k<\nu \\{\mbox{NaN}}&{\mbox{k odd}},0<\nu \leq k\\\infty &{\mbox{k even}},0<\nu \leq k\\\end{cases}}$

司徒頓t 分布信賴區間的推導

假設數量A在當T呈t-分布（T的自由度為n − 1）滿足

$\Pr(-A<T<A)=0.90\,$

這與

$\Pr(T<A)=0.95\,$是相同的

A是這個機率分布的第95個百分點

那麼

$\Pr \left(-A<{{\overline {X}}_{n}-\mu \over S_{n}/{\sqrt {n}}}<A\right)=0.9,$

等價於

$\Pr \left({\overline {X}}_{n}-A{S_{n} \over {\sqrt {n}}}<\mu <{\overline {X}}_{n}+A{S_{n} \over {\sqrt {n}}}\right)=0.9$

因此μ的90%信賴區間為：

${\overline {X}}_{n}\pm A{\frac {S_{n}}{\sqrt {n}}}$

計算

現在最方便的計算T分布的辦法是使用電子表格軟體（如Excel）或查相關線上計算網站。例如，Excel的TDIST(x,v,sides)用來計算自由度為v的T分布，如果第三個母數為1，則給出Pr(T>x)；如果第三個母數為2，則計算Pr(T>x Or T<-x).

下表列出了自由度為$\nu$的t 分布的單側和雙側區間值。例如，當樣本數量n=5時，則自由度$\nu$=4，我們就可以查找表中以4開頭的行。該行第5列值為2.132，對應的單側值為95%（雙側值為90%）。這也就是說，T小於2.132的機率為95%（即單側），記為Pr(−∞ < T < 2.132) = 0.95；同時，T值介於-2.132和2.132之間的機率為90%（即雙側），記為Pr(−2.132 < T < 2.132) = 0.9。

這是根據分布的對稱性計算得到的，

Pr(T < −2.132) = 1 − Pr(T > −2.132) = 1 − 0.95 = 0.05,

因此，

Pr(−2.132 < T < 2.132) = 1 − 2(0.05) = 0.9.

注意關於表格的最後一行的值：自由度為無限大的t-分布和常態分布等價。

單側	75%	80%	85%	90%	95%	97.5%	99%	99.5%	99.75%	99.9%	99.95%
雙側	50%	60%	70%	80%	90%	95%	98%	99%	99.5%	99.8%	99.9%
1	1.000	1.376	1.963	3.078	6.314	12.71	31.82	63.66	127.3	318.3	636.6
2	0.816	1.061	1.386	1.886	2.920	4.303	6.965	9.925	14.09	22.33	31.60
3	0.765	0.978	1.250	1.638	2.353	3.182	4.541	5.841	7.453	10.21	12.92
4	0.741	0.941	1.190	1.533	2.132	2.776	3.747	4.604	5.598	7.173	8.610
5	0.727	0.920	1.156	1.476	2.015	2.571	3.365	4.032	4.773	5.893	6.869
6	0.718	0.906	1.134	1.440	1.943	2.447	3.143	3.707	4.317	5.208	5.959
7	0.711	0.896	1.119	1.415	1.895	2.365	2.998	3.499	4.029	4.785	5.408
8	0.706	0.889	1.108	1.397	1.860	2.306	2.896	3.355	3.833	4.501	5.041
9	0.703	0.883	1.100	1.383	1.833	2.262	2.821	3.250	3.690	4.297	4.781
10	0.700	0.879	1.093	1.372	1.812	2.228	2.764	3.169	3.581	4.144	4.587
11	0.697	0.876	1.088	1.363	1.796	2.201	2.718	3.106	3.497	4.025	4.437
12	0.695	0.873	1.083	1.356	1.782	2.179	2.681	3.055	3.428	3.930	4.318
13	0.694	0.870	1.079	1.350	1.771	2.160	2.650	3.012	3.372	3.852	4.221
14	0.692	0.868	1.076	1.345	1.761	2.145	2.624	2.977	3.326	3.787	4.140
15	0.691	0.866	1.074	1.341	1.753	2.131	2.602	2.947	3.286	3.733	4.073
16	0.690	0.865	1.071	1.337	1.746	2.120	2.583	2.921	3.252	3.686	4.015
17	0.689	0.863	1.069	1.333	1.740	2.110	2.567	2.898	3.222	3.646	3.965
18	0.688	0.862	1.067	1.330	1.734	2.101	2.552	2.878	3.197	3.610	3.922
19	0.688	0.861	1.066	1.328	1.729	2.093	2.539	2.861	3.174	3.579	3.883
20	0.687	0.860	1.064	1.325	1.725	2.086	2.528	2.845	3.153	3.552	3.850
21	0.686	0.859	1.063	1.323	1.721	2.080	2.518	2.831	3.135	3.527	3.819
22	0.686	0.858	1.061	1.321	1.717	2.074	2.508	2.819	3.119	3.505	3.792
23	0.685	0.858	1.060	1.319	1.714	2.069	2.500	2.807	3.104	3.485	3.767
24	0.685	0.857	1.059	1.318	1.711	2.064	2.492	2.797	3.091	3.467	3.745
25	0.684	0.856	1.058	1.316	1.708	2.060	2.485	2.787	3.078	3.450	3.725
26	0.684	0.856	1.058	1.315	1.706	2.056	2.479	2.779	3.067	3.435	3.707
27	0.684	0.855	1.057	1.314	1.703	2.052	2.473	2.771	3.057	3.421	3.690
28	0.683	0.855	1.056	1.313	1.701	2.048	2.467	2.763	3.047	3.408	3.674
29	0.683	0.854	1.055	1.311	1.699	2.045	2.462	2.756	3.038	3.396	3.659
30	0.683	0.854	1.055	1.310	1.697	2.042	2.457	2.750	3.030	3.385	3.646
40	0.681	0.851	1.050	1.303	1.684	2.021	2.423	2.704	2.971	3.307	3.551
50	0.679	0.849	1.047	1.299	1.676	2.009	2.403	2.678	2.937	3.261	3.496
60	0.679	0.848	1.045	1.296	1.671	2.000	2.390	2.660	2.915	3.232	3.460
80	0.678	0.846	1.043	1.292	1.664	1.990	2.374	2.639	2.887	3.195	3.416
100	0.677	0.845	1.042	1.290	1.660	1.984	2.364	2.626	2.871	3.174	3.390
120	0.677	0.845	1.041	1.289	1.658	1.980	2.358	2.617	2.860	3.160	3.373
$\infty$	0.674	0.842	1.036	1.282	1.645	1.960	2.326	2.576	2.807	3.090	3.291

範例

給定一個樣本：樣本期望值和變異數分別為10和2，樣本大小為11（自由度為10）。根據公式：

${\overline {X}}_{n}\pm A{\frac {S_{n}}{\sqrt {n}}}$

可知，使用該方法統計出來的最大值，平均有90%的機率（即90%置信度/信心水準/confidence level）低於：

$10+1.37218{\frac {\sqrt {2}}{\sqrt {11}}}=10.58510.$

同理，使用該方法統計出來的最小值，平均有90%的機率（即90%置信度/信心水準/confidence level）高於：

$10-1.37218{\frac {\sqrt {2}}{\sqrt {11}}}=9.41490.$

因此，使用該方法統計出來的最大值和最小值，平均有80%的機率介於：

$10\pm 1.37218{\frac {\sqrt {2}}{\sqrt {11}}}=[9.41490,10.58510]$

兩值之間。（需注意此非代表數據的真正期望值介於這兩個值之間的機率為80%，詳情請參見信賴區間。）

參見

• 假說檢定
• 司徒頓t檢定
• 機率分布
• Levene's test

外部連結

• (en)Probability, Statistics and Estimation（頁面存檔備份，存於網際網路檔案館） 首先第112頁。