司徒頓t檢定

来自testwiki
跳转到导航 跳转到搜索

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

由來

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

應用

常見的應用有:

  • 单样本检验:检验一个正态分布的总体的均值是否在满足零假设的值之内,例如檢驗一群軍校男生的身高的平均是否符合全國標準的170公分界線。
  • 獨立樣本t檢定(双样本):其零假设为两个正态分布的总体的均值之差為某實數,例如檢定二群人之平均身高是否相等。若两母體的變異數是相等的情况下(同質變異數),自由度為兩樣本數相加再減二;若為異質變異數(母體變異數不相等),自由度則為Welch自由度,此情況下有时被称为Welch检验。
  • 配对樣本t檢定(成對樣本t檢定):檢定自同一母體抽出的成對樣本間差异是否为零。例如,檢测一位病人接受治疗前和治疗后的肿瘤尺寸大小。若治疗是有效的,我们可以推定多数病人接受治疗后,肿瘤尺寸將縮小。
  • 检验一迴歸模型的偏迴歸係數是否显著不为零,即檢定解釋變數X是否存在對被解釋變數Y的解釋能力,其檢定統計量稱之為t-比例(t-ratio)。

前提假設

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

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

計算

單樣本t檢驗

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

t=xμ0s/n

其中i=1nx=i=1nxin為樣本平均數,s=i=1n(xix)2n1為樣本標準差n樣本數。該統計量t在虛無假說:μ = μ0為真的條件下服從自由度為n − 1的t分佈

配對樣本t檢驗

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

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

t=dμ0sd/n

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

獨立雙樣本t檢驗

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

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

t=x1x2μ02sp2/n

其中i=1nx1=(i=1nx1i)/nx2=(i=1nx2i)/n為兩樣本各自的平均數,sp2=(i=1n(x1ix1)2+i=1n(x2ix2)2)/(2n2)為樣本之共同方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為2n − 2的t分佈

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

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

t=x1x2μ0sp2/n1+sp2/n2

其中i=1n1,其中j=1n2x1=(i=1nx1i)/nx2=(i=1nx2i)/n為兩樣本各自的平均數,sp2=(i=1n(x1ix1)2+j=1n(x2jx2)2)/(n1+n22)為兩樣本共同之方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為n1 + n2 − 2的t分佈

異質變異數假設 (Heteroscedasticity)

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

t=x1x2μ0s12/n1+s22/n2

其中i=1n1,其中j=1n2x1=(i=1n1x1i)/n1x2=(j=1n2x2j)/n為兩樣本各自的平均數,s12=(i=1n(x1ix1)2)/(n11)s22=(j=1n(x2jx2)2)/(n21)分別為兩樣本之方差。該統計量t在虛無假說:μ1 - μ2 = μ0為真的條件下服從自由度為

df=(s12/n1+s22/n2)2(s12/n1)2/(n11)+(s22/n2)2/(n21)

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

其它相關檢驗

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

以簡單線性迴歸為例

Template:Main 模型假設:

yi=α+βxi+εi,

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

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

t=β^β0SEβ^𝒯n2

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

SEβ^=1n2i=1n(yiy^i)2i=1n(xix)2

由於 ε^i=yiy^i=yi(α^+β^xi)為殘差(即估計誤差),而 SSR=i=1nε^i2 為殘差之離均平方和,我們可改寫t

t=(β^β0)n2SSR/i=1n(xix)2

另请参阅:F检验

電腦軟體

大多數的試算表軟體及統計軟體,諸如QtiPlotOpenOffice.org CalcLibreOffice CalcMicrosoft ExcelSASSPSSStataDAPgretlRPython ([1]Template:Wayback)、PSPPMinitab等,都可以進行t檢驗運算。

编程语言/软件程序 函数 注释
Microsoft Excel 2010 之前的版本 TTEST(array1, array2, tails, type) 参见 [2]
Microsoft Excel 2010 及更高版本 T.TEST(array1, array2, tails, type) 参见 [3]Template:Wayback
LibreOffice TTEST(Data1; Data2; Mode; Type) 参见 [4]Template:Wayback
Google Sheets TTEST(range1, range2, tails, type) 参见 [5]Template:Wayback
Python scipy.stats.ttest_ind(a, b, axis=0, equal_var=True) 参见 [6]Template:Wayback
Matlab ttest(data1, data2) 参见 [7]Template:Wayback
Mathematica TTest[{data1,data2}] 参见 [8]Template:Wayback
R t.test(data1, data2)
SAS PROC TTEST 参见 [9]
Java tTest(sample1, sample2) 参见 [10]Template:Wayback
Julia EqualVarianceTTest(sample1, sample2) 参见 [11]
Stata ttest data1 == data2 See [12]Template:Wayback

參見

參考文獻

Template:Reflist

Template:- Template:统计学

  1. Richard Mankiewicz, The Story of Mathematics (Princeton University Press), p.158.
  2. 2.0 2.1 Template:MacTutor Biography
  3. Template:Cite journal
  4. Template:Cite web