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

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

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

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

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

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

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

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

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

用15只小白鼠，隨機(jī)分為3組，比較小白鼠接種三種不同菌型傷寒桿菌9D、11C、DSC1后存活日數(shù)，實(shí)驗(yàn)數(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個(gè)總體的分布位置相同

H1：3個(gè)總體的分布位置不全相同

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

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

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

H1：3個(gè)總體的分布位置不全相同

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

等級資料的秩和檢驗(yàn)

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

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

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

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