醫(yī)學(xué)統(tǒng)計-秩和檢驗

醫(yī)學(xué)統(tǒng)計-秩和檢驗第一頁，共66頁。某實驗室觀察局部溫?zé)嶂委熜∈笠浦材[瘤的療效，以生存日數(shù)作為觀察指標，實驗結(jié)果見表1，試檢驗兩組小鼠生存日數(shù)有無差別？第二頁，共66頁。

表1兩組小鼠發(fā)癌后生存日數(shù)

第三頁，共66頁。秩和檢驗非參數(shù)統(tǒng)計（Nonparametricstatistics）不依賴于總體分布類型，也不對總體參數(shù)進行統(tǒng)計推斷的統(tǒng)計方法。參數(shù)統(tǒng)計（Parametricstatistics）要求樣本來自總體分布型是已知的（如正態(tài)分布），在這種假設(shè)的基礎(chǔ)上，對總體參數(shù)（如總體均數(shù)）進行估計和檢驗。第四頁，共66頁。秩、秩和秩

♀♂♀♀♀♂♂♀♀♂♂♂

123456789101112秩和♀=1+3+4+5+8+9=30♂=2+6+7+10+11+12=48第五頁，共66頁。例12.3某地職業(yè)病防治欲比較使用二巰基丙磺酸鈉與二巰基丁二酸鈉的驅(qū)汞效果。將22例汞中毒患者隨機分配到兩組，分別測定并計算出兩組驅(qū)汞的排汞比值，并將結(jié)果列于表12.3。試問兩藥驅(qū)汞效果有無差別？第六頁，共66頁。第七頁，共66頁。第八頁，共66頁。第九頁，共66頁。確定P值以n1及n2-n1查秩和檢驗臨界值表(附表11)。第十頁，共66頁。7984891411461510.100.050.02761540.01第十一頁，共66頁。第十二頁，共66頁。正態(tài)近似法如果n1＞10或n2-n1

＞10超出附表11的范圍，可按式（12.3）計算Z值：相同秩次較多（如超過25%）時，按式(12.4)進行校正第十三頁，共66頁。常用兩組比較的秩和檢驗檢驗分布的集中趨勢是否相同Mann-WhitneyUtestWilcoxontest當其中一組樣本含量大于10時，以上兩個統(tǒng)計量的分布均近似服從Z分布，用正態(tài)近似Z檢驗同時檢驗分布形狀及其集中趨勢thetwo-sampleKolmogorov-Smirnov第十四頁，共66頁。Wilcoxon符號秩和檢驗第十五頁，共66頁。例7.1

某市環(huán)保局對每一測定點用甲、乙兩種方法測定大氣中SO2日平均濃度(μg/m3)，測定值如表7.1。兩種方法的測定結(jié)果有無差別？第十六頁，共66頁。表7.19個測定點用甲、乙兩種方法的測定結(jié)果(μg/m3)第十七頁，共66頁。假設(shè)

第十八頁，共66頁。計算檢驗統(tǒng)計量求差值。見表7.1第(4)欄。編秩。依差值的絕對值從小到大編秩。編秩時遇差數(shù)等于零，舍去不計，同時樣本例數(shù)減1。遇絕對值相等差數(shù)，取平均秩次。秩次冠以原差值的正負號。分別計算正、負差值的秩和，任取其一為檢驗統(tǒng)計量。第十九頁，共66頁。確定P值和作出推斷結(jié)論當n50時，查附表10，T界值表。本例n=9，T=3.5，0.05>P>0.02,=0.05水準拒絕H0，差別有統(tǒng)計學(xué)意義?？梢哉J為甲、乙兩種方法測定大氣中SO2日平均濃度測定結(jié)果不同，乙法較高。第二十頁，共66頁。正態(tài)近似法若n>50超出附表10的范圍，可用Z檢驗，按（12.1)計算Z值，0.5為連續(xù)性校正數(shù)。第二十一頁，共66頁。相同“差值”（計絕對值）數(shù)多時（不包括差值為0值），改用（12.2）校正式。第二十二頁，共66頁。例12.1某研究用甲、乙兩種方法對某地方性砷中毒地區(qū)水源中砷含量(mg/L)進行測定，檢測10處，測量值如表12.1的(2)、(3)欄。問兩種方法的測定結(jié)果有無差別？第二十三頁，共66頁。第二十四頁，共66頁。第二十五頁，共66頁。第二十六頁，共66頁。注意當n<50時，根據(jù)和查附表10，T界值表(配對比較的符號秩和檢驗用)。由附表10可知，當n<5時，配對符號秩和檢驗不能得出雙側(cè)有統(tǒng)計學(xué)意義的概率，故樣本含量必須大于5。第二十七頁，共66頁。單一樣本與總體中位數(shù)比較例7.2已知某地正常人尿氟含量的中位數(shù)為2.15mmol/L。今在該地某廠隨機抽取12名工人，測得尿氟含量(mmol/L)如表7.2。問該廠工人的尿氟含量是否高于當?shù)卣Ｈ耍?/p>

第二十八頁，共66頁。第二十九頁，共66頁。H0：該廠工人尿氟含量總體中位數(shù)等于2.15mmol/LH1：該廠工人尿氟含量總體中位數(shù)大于2.15mmol/L

單側(cè)=0.05第三十頁，共66頁。(1)求差值：(2)編秩：方法同例7.1。(3)求秩和：T+=62.5，T-=3.5，取

T=3.5。第三十一頁，共66頁。查附表100.0513-530.02510-560.017-590.0055-61T=3.5第三十二頁，共66頁。710135356590.050.0250.015610.005第三十三頁，共66頁。3．確定P值，做出推斷結(jié)論

查附表10T界值表(配對比較的符號秩和檢驗用)，得單側(cè)P<0.005，按α=0.05水準拒絕H0，接受H1，差別有統(tǒng)計學(xué)意義，故可以認為該廠工人尿氟含量高于當?shù)卣Ｈ恕５谌捻?，?6頁。多個樣本比較Kruskal-Wallistest第三十五頁，共66頁。多個樣本比較Kruskal-Wallistest例

用15只小白鼠，隨機分為3組，比較小白鼠接種三種不同菌型傷寒桿菌9D、11C、DSC1后存活日數(shù)，實驗數(shù)據(jù)見下表第(1)、(3)、(5)欄，問接種不同菌型傷寒桿菌的小白鼠存活日數(shù)間有無差別？

第三十六頁，共66頁。表

小白鼠接種三種不同菌型傷寒桿菌后存活日數(shù)

9D11CDSC1存活日數(shù)

秩次

存活日數(shù)

秩次

存活日數(shù)

秩次

(1)(2)(3)(4)(5)(6)2143.579325691143.5681012.5568101114561012.51315Ri

18.5-40-61.5第三十七頁，共66頁。H0：3個總體的分布位置相同

H1：3個總體的分布位置不全相同

α=0.05計算統(tǒng)計量H：混合編秩，相同數(shù)值，取平均秩，算得各組的秩和Ri,第三十八頁，共66頁。第三十九頁，共66頁。3．確定P值，做出推斷結(jié)論

以N=15，n1=n2=n3=5查附表12，H界值表，得P<0.01，按=0.05水準拒絕H0，接受H1，差別有統(tǒng)計學(xué)意義，故可認為接種3種不同菌型傷寒桿菌的小白鼠存活日數(shù)不全相同。

第四十頁，共66頁。例12.5某醫(yī)師檢測3種卵巢功能異?；颊哐逯写冱S體素的含量(U/L)資料見表12.5第(1)、(3)、(5)欄。問3種患者血清中促黃體素的含量(U/L)是否有差別？第四十一頁，共66頁。第四十二頁，共66頁。H0：3個總體的分布位置相同

H1：3個總體的分布位置不全相同

α=0.05計算統(tǒng)計量H：混合編秩，相同數(shù)值，取平均秩，算得各組的秩和Ri,第四十三頁，共66頁。第四十四頁，共66頁。第四十五頁，共66頁。第四十六頁，共66頁。

等級資料的秩和檢驗

相同秩較多，算得的H值偏小，還需計算校正的統(tǒng)計量HC

第四十七頁，共66頁。例12.6某醫(yī)院用三種復(fù)方小葉枇杷治療老年慢性支氣管炎，數(shù)據(jù)見表12.6第(1)、(4)欄，試比較三種方劑的療效有無差異。第四十八頁，共66頁。第四十九頁，共66頁。第五十頁，共66頁。第五十一頁，共66頁。第五十二頁，共66頁。注意TheKruskal-Wallistest等同于用秩次進行單因素方差分析另外Mediantest也可用于多組的比較第五十三頁，共66頁。隨機區(qū)組設(shè)計的秩和檢驗

第五十四頁，共66頁。M檢驗(Friedman法)查表法第五十五頁，共66頁。第五十六頁，共66頁。例12.8觀察某藥不同劑量對肝功能的影響，將同種屬的28只大白鼠按窩別、性別、體重配為7個配伍組，每個區(qū)組的4只大白鼠隨機分入不同的4種劑量組，在用藥后一周測定各血清中指標DT值的變化，結(jié)果見表12.8，問此藥不同劑量對血清中指標DT值的影響有無不同？例12.8第五十七頁，共66頁。第五十八頁，共66頁。第五十九頁，共66頁。計算統(tǒng)計量M值第六十頁，共66頁。確定P值，做出推斷結(jié)論

第六十一頁，共66頁。分布近似法當處理組數(shù)k或配伍組數(shù)b超出附表13的范圍時，=k–1第六十二頁，共66頁。成組設(shè)計多個樣本資料的兩兩比較第六十三頁，共66