司徒頓t檢定

假說檢定的一種

司徒頓t 檢定(英語:Student's t-test)是指虛無假說成立時的任一檢定統計有司徒頓t分布統計假說檢定,屬於母數統計。司徒頓t檢定常作為檢定一群來自常態分配母體的獨立樣本期望值是否為某一實數,或是二(兩)群來自常態分配母體的獨立樣本期望值的差是否為某一實數。舉個簡單的例子,在某個學校中我們可以從某個年級中隨機抽樣一群男生,以檢定該年級男生與全校男生之身高差異程度是否如我們所假設的某個值。

由來

司徒頓t檢定是威廉·高斯特為了觀測釀酒品質於1908年所提出的,「司徒頓 (student)」則是他的筆名[1][2][3][4] 基於克勞德·健力士(Claude Guinness)聘用從牛津大學劍橋大學出來的最好的畢業生,[2]以將生物化學及統計學應用到健力士工業流程的創新政策,高斯特受僱於都柏林的健力士釀酒廠擔任統計學家。高斯特提出了t檢定以降低啤酒重量監控的成本。高斯特於1908年在《Biometrika英語Biometrika》期刊上公布t檢定,但因其老闆認為其為商業機密而被迫使用筆名,統計學論文內容也跟釀酒無關。實際上,其他統計學家是知道高斯特真實身份的。

應用

常見的應用有:

  • 單樣本檢定:檢定一個常態分布的母體的均值是否在滿足虛無假說的值之內,例如檢定一群軍校男生的身高的平均是否符合全國標準的170公分界線。
  • 獨立樣本t檢定(雙樣本):其虛無假說為兩個常態分布的母體的均值之差為某實數,例如檢定二群人之平均身高是否相等。若兩母體的變異數是相等的情況下(同質變異數),自由度為兩樣本數相加再減二;若為異質變異數(母體變異數不相等),自由度則為Welch自由度,此情況下有時被稱為Welch檢定。
  • 配對樣本t檢定(成對樣本t檢定):檢定自同一母體抽出的成對樣本間差異是否為零。例如,檢測一位病人接受治療前和治療後的腫瘤尺寸大小。若治療是有效的,我們可以推定多數病人接受治療後,腫瘤尺寸將縮小。
  • 檢定一迴歸模型的偏迴歸係數是否顯著不為零,即檢定解釋變數X是否存在對被解釋變數Y的解釋能力,其檢定統計量稱之為t-比例(t-ratio)。

前提假設

大多數的t檢定之統計量具有t = Z/s的形式,其中Zs是已知資料的函數。Z通常被設計成對於對立假說有關的形式,而s是一個比例母數使t服從於t分布。以單樣本t檢定為例, ,其中 為樣本平均數, 為樣本數, 為母體標準差。至於s在單樣本t檢定中為 ,其中 為樣本的標準差。在符合零假說的條件下,t檢定有以下前提:

  • Z 服從標準常態分布
  • (n - 1)s2 服從自由度(n - 1)的卡方分布
  • Zs互相獨立

計算

單樣本t檢定

檢定虛無假說為一群來自常態分配獨立樣本xi之母體期望值μμ0可利用以下統計量

 

其中  為樣本平均數, 為樣本標準差n樣本數。該統計量t在虛無假說:μ = μ0為真的條件下服從自由度為n − 1的t分布

配對樣本t檢定

配對樣本t檢定可視為單樣本t檢定的擴展,不過檢定的對象由一群來自常態分配獨立樣本更改為兩配對樣本之觀測值之差。

若兩配對樣本x1ix2i之差為di = x1ix2i獨立且來自常態分配,則di之母體期望值μ是否為μ0可利用以下統計量

 

其中  為配對樣本差值之平均數, 為配對樣本差值之標準差n為配對樣本數。該統計量t在虛無假說:μ = μ0為真的條件下服從自由度為n − 1的t分布

獨立雙樣本t檢定

同質變異數假設 (Homoscedasticity)、樣本數相等

若兩獨立樣本x1ix2i具有相同之樣本數n,且來自兩個母體變異數相同(同質變異數假設)的常態分配,則兩母體之期望值差μ1 - μ2是否為μ0可利用以下統計量

 

其中   為兩樣本各自的平均數, 為樣本之共同變異數。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為2n − 2的t分布

同質變異數假設 (Homoscedasticity)、樣本數不相等

若兩獨立樣本x1ix2j具有不相同之樣本數n1n2,且來自兩個母體變異數相同(同質變異數假設)的常態分配,則兩母體之期望值之差μ1 - μ2是否為μ0可利用以下統計量

 

其中 ,其中   為兩樣本各自的平均數, 為兩樣本共同之變異數。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為n1 + n2 − 2的t分布

異質變異數假設 (Heteroscedasticity)

若兩獨立樣本x1ix2j具有相同或不相同之樣本數n1n2,且兩者母體變異數不相等(異質變異數假設)的常態分配,則兩母體之期望值之差μ1 - μ2是否為μ0可利用以下統計量

 

其中 ,其中   為兩樣本各自的平均數,  分別為兩樣本之變異數。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為

 

t分布。這種方法又常稱為Welch檢定。

其它相關檢定

偏迴歸係數是否為零之檢定

以簡單線性迴歸為例

模型假設:

 

其中xii = 1, ..., n為已知,αβ為未知係數,εi殘差獨立且服從期望值0且變異數σ2未知的常態分布,yii = 1, ..., n為觀測值。我們可以檢定迴歸係數β是否相等於特定的β0,通常使β0 = 0以檢定xiyi是否存在解釋能力,在此例(簡單線性迴歸模型)即為檢定迴歸式之斜率是否為零。

  最小平方法之估計值,  為最小平方法估計值之標準誤差,則

 

在虛無假說為β = β0的情況下服從自由度為n − 2之t分布,此檢定統計量被稱作「t比率 (t-ratio)」,其中

 

由於  為殘差(即估計誤差),而   為殘差之離均平方和,我們可改寫t

 

另請參閱:F檢定

電腦軟體

大多數的試算表軟體及統計軟體,諸如QtiPlotOpenOffice.org CalcLibreOffice CalcMicrosoft ExcelSASSPSSStataDAPgretlRPython ([1]頁面存檔備份,存於網際網路檔案館))、PSPPMinitab等,都可以進行t檢定運算。

程式語言/軟體程序 函數 注釋
Microsoft Excel 2010 之前的版本 TTEST(array1, array2, tails, type) 參見 [2]
Microsoft Excel 2010 及更高版本 T.TEST(array1, array2, tails, type) 參見 [3]頁面存檔備份,存於網際網路檔案館
LibreOffice TTEST(Data1; Data2; Mode; Type) 參見 [4]頁面存檔備份,存於網際網路檔案館
Google Sheets TTEST(range1, range2, tails, type) 參見 [5]頁面存檔備份,存於網際網路檔案館
Python scipy.stats.ttest_ind(a, b, axis=0, equal_var=True) 參見 [6]頁面存檔備份,存於網際網路檔案館
Matlab ttest(data1, data2) 參見 [7]頁面存檔備份,存於網際網路檔案館
Mathematica TTest[{data1,data2}] 參見 [8]頁面存檔備份,存於網際網路檔案館
R t.test(data1, data2)
SAS PROC TTEST 參見 [9]
Java tTest(sample1, sample2) 參見 [10]頁面存檔備份,存於網際網路檔案館
Julia EqualVarianceTTest(sample1, sample2) 參見 [11]
Stata ttest data1 == data2 See [12]頁面存檔備份,存於網際網路檔案館

參見

參考文獻

  1. ^ Richard Mankiewicz, The Story of Mathematics (Princeton University Press), p.158.
  2. ^ 2.0 2.1 約翰·J·奧康納; 埃德蒙·F·羅伯遜, Gosset, MacTutor數學史檔案 (英語) 
  3. ^ Fisher Box, Joan. Guinness, Gosset, Fisher, and Small Samples. Statistical Science. 1987, 2 (1): 45–52. JSTOR 2245613. doi:10.1214/ss/1177013437. 
  4. ^ 存档副本 (PDF). [2013-08-10]. (原始內容 (PDF)存檔於2017-05-16).