醫(yī)學統(tǒng)計學-t檢驗

t檢驗醫(yī)學統(tǒng)計學

第五講t

檢驗

計量資料兩個樣本均數(shù)比較的假設(shè)檢驗

總體標準差已知，或總體標準差未知但樣本含量n較大（n≥60

）——

z檢驗（z-test）

總體標準差未知，樣本量較小——t檢驗（t-test，又稱Student’st-test）

主要內(nèi)容單樣本均數(shù)的t檢驗配對樣本均數(shù)的t檢驗兩獨立樣本均數(shù)的t檢驗正態(tài)性檢驗兩樣本的方差的齊性檢驗兩總體方差不等時的均數(shù)比較的t’檢驗第七章兩樣本均數(shù)比較的假設(shè)檢驗（P113~125）主要內(nèi)容第一節(jié)單樣本均數(shù)t檢驗

用于樣本均數(shù)與已知總體均數(shù)μ0（理論值、標準值或經(jīng)過大量觀察所得到的穩(wěn)定值）的比較檢驗所代表的未知總體均數(shù)μ是否與已知總體均數(shù)μ0有差別應(yīng)用條件：資料服從正態(tài)分布單樣本均數(shù)的t檢驗（onesamplet-test）

（7-1）例7-1

通過以往大量資料得知某地20歲男子平均身高為168cm，今隨機測量當?shù)?6名20歲男子，得其身高均數(shù)為172cm，標準差為14cm。問當?shù)噩F(xiàn)在20歲男子的平均身高是否比以往高？根據(jù)經(jīng)驗身高服從正態(tài)分布，可用單樣本均數(shù)的t

檢驗。

（1）建立檢驗假設(shè)，確定檢驗水準

H0:μ=μ0，現(xiàn)在該地20歲男子的平均身高與以往20歲男子的平均身高相等

H1:μ>μ0，現(xiàn)在20歲男子的平均身高高于以往20歲男子的平均身高單側(cè)α=0.05（3）確定P值，作出推斷結(jié)論查t界值表，得單側(cè)概率0.10<P<0.20，按

=0.05水準，不拒絕H0，尚不能認為該地20歲男子平均身高比以往要高。

樣本均數(shù)與已知總體均數(shù)的比較，n≥60時，亦可近似地采用z檢驗（2）計算檢驗統(tǒng)計量第二節(jié)配對樣本均數(shù)的t檢驗配對樣本均數(shù)的t

檢驗

（paired/matchedt-test）

又稱配對t檢驗，適用于配對設(shè)計的計量資料兩相關(guān)樣本均數(shù)的比較。配對設(shè)計（paireddesign）有兩種：

同源配對異源配對同源配對：同一受試對象分別接受兩種不同處理如：分別用微生物法和放免法測定同一份血標本的血漿葉酸含量分別用A、B兩種儀器測定同一研究對象的Hb含量異源配對：為消除混雜因素的影響，將實驗對象按某些重要特征（對結(jié)果變量有影響但非本次研究的處理因素）相同或相近的原則配成若干對子，再隨機地將每對中的兩個研究對象分別分入兩組，施以相應(yīng)的處理如：同窩、同性別的兩只動物配成一對，然后隨機地分入甲、乙兩組，給予甲飼料和乙飼料

配對設(shè)計的數(shù)據(jù)具有一一對應(yīng)的特征分析時首先求出各對數(shù)據(jù)間的差值d，將d作為變量值計算均數(shù)。若兩處理因素的效應(yīng)無差別，理論上差值d的總體均數(shù)μd應(yīng)為0，故可將該檢驗理解為樣本均數(shù)所對應(yīng)的總體μd與總體均數(shù)0的比較配對t檢驗的應(yīng)用條件：差值d服從正態(tài)分布檢驗統(tǒng)計量按前述公式（7-1）構(gòu)造如下：（7-2）例7-2某醫(yī)院用A、B兩種血紅蛋白測定儀器檢測了16名健康男青年的血紅蛋白含量（g/L），檢測結(jié)果見表7-1。問：兩種血紅蛋白測定儀器的檢測結(jié)果是否有差別。表7-1兩種儀器檢測16名健康男青年血紅蛋白（g/L）的結(jié)果編號儀器A儀器B差值d（1）（2）（3）（4）=（3）-（2）1113140272125150253126138124130120-105150140-1061451450713513508105115109128135710135130-511100120201213013331311014737141151251015120114-61615516510同源配對

對差值d進行正態(tài)性檢驗（Shapiro-Wilk統(tǒng)計量W=0.949，P=0.470），認為其滿足正態(tài)性，可用配對樣本均數(shù)的t

檢驗。（1）建立檢驗假設(shè)，確定檢驗水準

H0：μd=0，A、B兩種儀器檢測的血紅蛋白總體平均差異為0H1：μd≠0，A、B兩種儀器檢測的血紅蛋白總體平均差異不為0α=0.05（2）計算檢驗統(tǒng)計量（3）確定P值，做出推斷結(jié)論查t界值表，t0.05,15=2.131，t0.02,15=2.602，得0.02＜P＜0.05，按

=0.05水準，拒絕H0，接受H1，差異有統(tǒng)計學意義，可認為A、B兩種血紅蛋白測定儀器檢測結(jié)果的平均差異不為0，儀器B檢測的血紅蛋白較高。

例7-3

某研究者采用配對設(shè)計進行實驗，比較兩種抗癌藥物對小白鼠肉瘤抑瘤效果。先將10只染有肉瘤的小白鼠按體重大小配成5個對子，每個對子內(nèi)2只小白鼠隨機接受兩種抗癌藥物，以肉瘤的重量為指標，實驗結(jié)果見表7-2。問：兩種藥物的抑瘤效果有無差別？表7-2不同藥物作用后小白鼠肉瘤重量（g）編號A藥B藥差值d（1）（2）（3）（4）=（2）-（3）10.820.650.1720.730.540.1930.430.340.0940.410.210.2050.680.430.25異源配對第三節(jié)

兩獨立樣本均數(shù)的t檢驗

又稱成組t檢驗，適用于完全隨機設(shè)計兩獨立樣本均數(shù)的比較，檢驗兩獨立樣本均數(shù)所代表的未知總體均數(shù)是否有差別。應(yīng)用條件：獨立性兩組數(shù)據(jù)均服從正態(tài)分布兩樣本對應(yīng)的兩總體方差相等兩獨立樣本均數(shù)的t

檢驗

（two-samplet-testforindependentsamples）

當兩總體方差相等，即時,檢驗統(tǒng)計量可按公式（7-1）在H0：μ=μ1-μ2=0條件下構(gòu)造如下：

（7-3）（7-4）

若兩總體方差不等，即時，可采用數(shù)據(jù)變換或近似t檢驗（t'檢驗）或基于秩次的非參數(shù)檢驗。

例7-4

為了解內(nèi)毒素對肌酐的影響，將20只雌性中年大鼠隨機分為甲組和乙組。甲組的每只大鼠不給予內(nèi)毒素，乙組的每只大鼠則給予3mg/kg的內(nèi)毒素。測得兩組大鼠的肌酐（mg/L）結(jié)果如表7-3。問：內(nèi)毒素是否對肌酐有影響？表7-3兩組大鼠的肌酐（mg/L）結(jié)果

甲組乙組6.28.53.76.85.811.32.79.43.99.36.17.36.75.67.87.93.87.26.98.2

本例為完全隨機設(shè)計，正態(tài)性檢驗結(jié)果顯示兩組資料均滿足正態(tài)性，兩樣本方差齊性檢驗結(jié)果顯示滿足方差齊性→可用兩獨立樣本均數(shù)的t檢驗（1）建立檢驗假設(shè)，確定檢驗水準

H0：

1=

2，即內(nèi)毒素對肌酐無影響

H1：

1≠

2，即內(nèi)毒素對肌酐有影響

=0.05（2）計算檢驗統(tǒng)計量（3）確定P值，作出推斷結(jié)論查t界值表，得0.001<P<0.002，按

=0.05水準，拒絕H0，接受H1，差別有統(tǒng)計學意義，可認為內(nèi)毒素對肌酐有影響，具有升高肌酐的作用。兩獨立樣本均數(shù)比較，當兩組例數(shù)基本相等且n1+n2≥60時，亦可近似地采用z檢驗。

兩樣本均數(shù)t檢驗的應(yīng)用條件①獨立性（independence）②正態(tài)性（normality）③方差齊性（homogeneity）①專業(yè)知識②專業(yè)知識，必要時可對各組資料進行正態(tài)性檢驗③方差齊性檢驗判斷依據(jù)第四節(jié)正態(tài)性檢驗

許多統(tǒng)計方法僅適于正態(tài)分布或近似正態(tài)分布資料如：用均數(shù)和標準差描述計量資料的分布特征

t、z檢驗、方差分析等，均要求樣本資料服從正態(tài)分布因此，選定統(tǒng)計方法時，首先檢驗資料是否服從正態(tài)分布。正態(tài)性檢驗（NormalityTest）正態(tài)分布資料應(yīng)同時具有正態(tài)對稱性（偏度，skewness）和正態(tài)峰（峰度，kurtosis）兩個特征。正態(tài)性檢的方法

圖示法：頻率-頻率圖（P-Pplot）和分位數(shù)-分位數(shù)圖（Q-Qplot）

計算法：1）矩法：偏度系數(shù)Υ1、

g1

，σg1

，峰度系數(shù)Υ2、

g2，σg2→z檢驗

2）計算綜合統(tǒng)計量：夏皮羅-威爾克

Shapiro-Wilk法（W檢驗）達戈斯提諾D′Agostino法（D檢驗）

Shapiro-Francia法（W′檢驗）

3）正態(tài)分布的擬合優(yōu)度檢驗：皮爾遜x2

檢驗、對數(shù)似然比檢驗

柯爾莫哥洛夫Kolmogorov-Smirov法負偏態(tài)圖7-3正態(tài)分布與偏態(tài)分布及峰態(tài)分布曲線示意圖圖7-1血清總膽固醇正態(tài)分布P-PPlot圖7-2血清總膽固醇正態(tài)分布Q-QPlotSPSS:

樣本含量3≤n≤5000時，以Shapiro-Wilk法為準

n>5000時，以Kolmogorov-Smirnov法為準

SAS：n≤2000時，以Shapiro-Wilk法為準

n>2000時，以Kolmogorov-Smirnov法為準正態(tài)性檢驗方法選擇的SPSS和SAS規(guī)則：第五節(jié)

兩樣本的方差的齊性檢驗

方差齊性檢驗

（Testofhomogeneityofvariance）

檢驗方法：F檢驗、Bartlett檢驗、Levene檢驗

F檢驗、Bartlett檢驗要求資料服從正態(tài)分布

F檢驗只用于兩樣本，Bartlett和Levene可用于兩樣本或多樣本

Levene不依賴總體分布形式，更為穩(wěn)?。⊿PSS采用）Bartlett法檢驗和Levene檢驗：第八章兩樣本方差齊性的F檢驗：檢驗統(tǒng)計量F（7-13）F值為兩個樣本方差之比，如僅是抽樣誤差的影響，它一般不會偏離1太遠。求得F值后，查F界值表可得P值（F值愈大，P值愈?。?，按檢驗水準

作出推斷結(jié)論。例7-6

對例7-4，用F檢驗判斷兩總體肌酐的方差是否不等。

n1=n2=10，S1=1.6985（mg/L），S2=1.5967（mg/L）（1）建立檢驗假設(shè)，確定檢驗水準

H0：σ12=σ22

，即兩組肌酐的總體方差相等

H1：σ12≠σ22，即兩組肌酐的總體方差不等

=0.10（取稍大的

可減少II型錯誤，提高把握度1-β）（2）計算檢驗統(tǒng)計量（3）確定P值，作出推斷結(jié)論以

1=9、

2=9查F界值表，F(xiàn)0.10,(9,9)=3.18，得P>0.10，按

=0.10水準，不拒絕H0，差別無統(tǒng)計學意義，尚不能認為兩組肌酐的總體方差不等。不易接受H0，結(jié)論較保守第六節(jié)

兩總體方差不等時均數(shù)比較的t’檢驗

亦稱近似t檢驗（separatevarianceestimationt-test）包括：

t’檢驗校正自由度：Satterthwaite法近似t檢驗（統(tǒng)計軟件普遍使用）

Welch法近似t檢驗校正臨界值：Cochran&Cox法近似t檢驗t檢驗的SPSS操作過程one-sample

t-test

某工廠用自動打包機打包，每包標準質(zhì)量為100kg。某日開工后在試機中共打了9個包，測得9包的質(zhì)量（kg）為：99.3、98.7、100.5、101.2、98.3、99.7、99.5、102.1、100.5。請判斷打包機工作是否正常。（數(shù)據(jù)見張文彤，Part3：onesamplet.sav）

例7-1沒有原始數(shù)據(jù)的單樣本t檢驗語法datalistfree/nmeansmu.begindata.1617214168enddata.computet=abs(mean-mu)/(s/sqrt(n)).COMPUTEp=1-CDF.T(t,n-1).EXECUTE.formatt(f10.5)p(f10.5).list.了解內(nèi)容不要求掌握例題：八年制教材第二版-SPSS數(shù)據(jù)CH07-02.savPaired