司徒頓t檢定
司徒頓t 檢定(英語:)是指虛無假說成立時的任一檢定統計有司徒頓t分布的統計假說檢定,屬於母數統計。學生t檢驗常作為檢驗一群來自常態分配母體的獨立樣本之期望值是否為某一實數,或是二(两)群來自常態分配母體的獨立樣本之期望值的差是否為某一實數。舉個簡單的例子,在某個學校中我們可以從某個年級中隨機抽樣一群男生,以檢驗該年級男生與全校男生之身高差異程度是否如我們所假設的某個值。
由來
司徒頓t檢定是威廉·戈塞為了觀測釀酒品質於1908年所提出的,「司徒頓 (student)」則是他的筆名。[1][2][3][4] 基於克勞德·健力士(Claude Guinness)聘用從牛津大學和劍橋大學出來的最好的畢業生,[2]以將生物化學及統計學應用到健力士工業流程的創新政策,戈塞受雇於都柏林的健力士釀酒廠擔任統計學家。戈塞提出了t检验以降低啤酒重量监控的成本。戈塞於1908年在《Biometrika》期刊上公布t檢驗,但因其老闆認為其為商業機密而被迫使用筆名,統計學論文內容也跟釀酒無關。實際上,其他统计学家是知道戈塞真實身份的。
應用
常見的應用有:
- 单样本检验:检验一个正态分布的总体的均值是否在满足零假设的值之内,例如檢驗一群軍校男生的身高的平均是否符合全國標準的170公分界線。
- 獨立樣本t檢定(双样本):其零假设为两个正态分布的总体的均值之差為某實數,例如檢定二群人之平均身高是否相等。若两母體的變異數是相等的情况下(同質變異數),自由度為兩樣本數相加再減二;若為異質變異數(母體變異數不相等),自由度則為Welch自由度,此情況下有时被称为Welch检验。
- 配对樣本t檢定(成對樣本t檢定):檢定自同一母體抽出的成對樣本間差异是否为零。例如,檢测一位病人接受治疗前和治疗后的肿瘤尺寸大小。若治疗是有效的,我们可以推定多数病人接受治疗后,肿瘤尺寸將縮小。
- 检验一迴歸模型的偏迴歸係數是否显著不为零,即檢定解釋變數X是否存在對被解釋變數Y的解釋能力,其檢定統計量稱之為t-比例(t-ratio)。
前提假設
大多數的t檢定之統計量具有t = Z/s的形式,其中Z與s是已知資料的函數。Z通常被設計成對於對立假說有關的形式,而s是一個比例母數使t服從於t分佈。以單樣本t檢驗為例,,其中為樣本平均數,為樣本數,為总体標準差。至於s在單樣本t檢驗中為,其中為樣本的標準差。在符合零假說的條件下,t檢定有以下前提:
- Z 服從標準常態分佈
- (n - 1)s2 服從自由度(n - 1)的卡方分佈
- Z與s互相獨立
計算
單樣本t檢驗
檢驗虛無假說為一群來自常態分配獨立樣本xi之母體期望值μ為μ0可利用以下統計量
其中,為樣本平均數,為樣本標準差,n為樣本數。該統計量t在虛無假說:μ = μ0為真的條件下服從自由度為n − 1的t分佈。
配對樣本t檢驗
配對樣本t檢驗可視為單樣本t檢驗的擴展,不過檢驗的對象由一群來自常態分配獨立樣本更改為兩配對樣本之觀測值之差。
若兩配對樣本x1i與x2i之差為di = x1i − x2i獨立且來自常態分配,則di之母體期望值μ是否為μ0可利用以下統計量
其中,為配對樣本差值之平均數,為配對樣本差值之標準差,n為配對樣本數。該統計量t在虛無假說:μ = μ0為真的條件下服從自由度為n − 1的t分布。
同質變異數假設 (Homoscedasticity)、樣本數相等
若兩獨立樣本x1i與x2i具有相同之樣本數n,且來自兩個母體變異數相同(同質變異數假設)的常態分配,則兩母體之期望值差μ1 - μ2是否為μ0可利用以下統計量
其中,及為兩樣本各自的平均數,為樣本之共同方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為2n − 2的t分佈。
同質變異數假設 (Homoscedasticity)、樣本數不相等
若兩獨立樣本x1i與x2j具有不相同之樣本數n1與n2,且來自兩個母體變異數相同(同質變異數假設)的常態分配,則兩母體之期望值之差μ1 - μ2是否為μ0可利用以下統計量
其中,其中,及為兩樣本各自的平均數,為兩樣本共同之方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為n1 + n2 − 2的t分佈。
異質變異數假設 (Heteroscedasticity)
若兩獨立樣本x1i與x2j具有相同或不相同之樣本數n1與n2,且兩者母體變異數不相等(異質變異數假設)的常態分配,則兩母體之期望值之差μ1 - μ2是否為μ0可利用以下統計量
其中,其中,及為兩樣本各自的平均數,及分別為兩樣本之方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為
之t分布。這種方法又常稱為Welch檢驗。
以簡單線性迴歸為例
模型假設:
其中xi,i = 1, ..., n為已知,α與β為未知係數,εi為殘差獨立且服從期望值0且方差σ2未知的常態分佈,yi,i = 1, ..., n為觀測值。我們可以檢驗迴歸係數β是否相等於特定的β0,通常使β0 = 0以檢定xi對yi是否存在解釋能力,在此例(簡單線性迴歸模型)即為檢定迴歸式之斜率是否為零。
令與為最小平方法之估計值,與為最小平方法估計值之標準誤差,則
在虛無假說為β = β0的情況下服從自由度為n − 2之t分布,此檢定統計量被稱作「t比率 (t-ratio)」,其中
由於 為殘差(即估計誤差),而 為殘差之離均平方和,我們可改寫t為
另请参阅:F检验
電腦軟體
大多數的試算表軟體及統計軟體,諸如QtiPlot、OpenOffice.org Calc、LibreOffice Calc、Microsoft Excel、SAS、SPSS、Stata、DAP、gretl、R、Python ((页面存档备份,存于))、PSPP、Minitab等,都可以進行t檢驗運算。
编程语言/软件程序 | 函数 | 注释 |
---|---|---|
Microsoft Excel 2010 之前的版本 | TTEST(array1, array2, tails, type) | 参见 |
Microsoft Excel 2010 及更高版本 | T.TEST(array1, array2, tails, type) | 参见 (页面存档备份,存于) |
LibreOffice | TTEST(Data1; Data2; Mode; Type) | 参见 (页面存档备份,存于) |
Google Sheets | TTEST(range1, range2, tails, type) | 参见 (页面存档备份,存于) |
Python | scipy.stats.ttest_ind(a, b, axis=0, equal_var=True) | 参见 (页面存档备份,存于) |
Matlab | ttest(data1, data2) | 参见 (页面存档备份,存于) |
Mathematica | TTest[{data1,data2}] | 参见 (页面存档备份,存于) |
R | t.test(data1, data2) | |
SAS | PROC TTEST | 参见 |
Java | tTest(sample1, sample2) | 参见 (页面存档备份,存于) |
Julia | EqualVarianceTTest(sample1, sample2) | 参见 |
Stata | ttest data1 == data2 | See (页面存档备份,存于) |
參考文獻
- Richard Mankiewicz, The Story of Mathematics (Princeton University Press), p.158.
- 約翰·J·奧康納; 埃德蒙·F·羅伯遜, , (英语)
- Fisher Box, Joan. . Statistical Science. 1987, 2 (1): 45–52. JSTOR 2245613. doi:10.1214/ss/1177013437.
- (PDF). [2013-08-10]. (原始内容 (PDF)存档于2017-05-16).