司徒頓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) \,_2F_1 \left ( \frac{1}{2},(\nu+1)/2;\frac{3}{2};-\frac{x^2}{\nu} \right)} {\sqrt{\pi\nu}\,\Gamma (\nu/2)}$其中：$\,_2F_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} <pre> \frac{\nu+1}{2}\left[ \psi(\frac{1+\nu}{2}) - \psi(\frac{\nu}{2}) \right] \\[0.5em] </pre> + \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分布（Student's t-distribution），簡稱t 分布，在機率論及統計學中用於根據小樣本來估計母體呈常態分布且標準差未知的期望值。若母體標準差已知，或是樣本數足夠大時（依據中央極限定理漸進常態分布），則應使用常態分布來進行估計。其為對兩個樣本期望值差異進行顯著性測試的司徒頓t檢定之基礎。

司徒頓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分布的機率累計函數，用不完全貝塔函數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 − 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\pm1.37218 \frac{\sqrt{2}}{\sqrt{11}}=[9.41490,10.58510]$$

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

## 外部連結

- (en) Probability, Statistics and Estimation 首先第112頁。