醫(yī)學(xué)統(tǒng)計(jì)學(xué)第十一章-非參數(shù)檢驗(yàn)(課堂)-PPT幻燈片

參數(shù)檢驗(yàn)parametrictest

(1)總體分布類(lèi)型已知(2)由樣本統(tǒng)計(jì)量推斷未知總體參數(shù)。

t檢驗(yàn)：

F檢驗(yàn):

非參數(shù)檢驗(yàn)(nonparametrictest)對(duì)數(shù)據(jù)的總體分布類(lèi)型不作嚴(yán)格假定，又稱(chēng)任意分布檢驗(yàn)(distribution-freetest)，它直接對(duì)總體分布的位置作假設(shè)檢驗(yàn)。

應(yīng)用非參數(shù)檢驗(yàn)的情況1.不滿足正態(tài)和方差齊性條件的小樣本資料；2.總體分布類(lèi)型不明的小樣本資料；3.一端或二端是不確定數(shù)值（如＜0.002、＞65等）的資料；4.單向有序列聯(lián)表資料；第一節(jié)配對(duì)樣本比較的Wilcoxon符號(hào)秩檢驗(yàn)一、配對(duì)樣本差值的中位數(shù)和0比較例8-1對(duì)12份血清分別用原方法和新方法測(cè)谷－丙轉(zhuǎn)氨酶，結(jié)果見(jiàn)表，問(wèn)兩法所得結(jié)果有無(wú)差別？1．H0：差值的總體中位數(shù)Md=0 H1：差值的總體中位數(shù)Md≠0

α＝0.052．求差值；依其絕對(duì)值從小到大編秩次（i）絕對(duì)值相等者取平均秩次；（ii）將差值的正負(fù)標(biāo)在秩次之前；（iii）零差值不參與編秩3．分別求正負(fù)秩次之和，任取正秩和或負(fù)秩和為T(mén)值4．根據(jù)統(tǒng)計(jì)量T確定對(duì)應(yīng)的P值（i）n<=50小樣本時(shí)，查表（附表9）(ii)n>50大樣本時(shí)，正態(tài)近似編號(hào)原法新法差值正秩負(fù)秩160761682142152105319524348114808221.55242240-21.56220220071902051578253813691982434591038446411236190-4610129510053合計(jì)54.511.5本例n=11,T=11.5,查附表9，得10～56P>0.05,不拒絕H0,兩法所得結(jié)果無(wú)差別正態(tài)近似二．單個(gè)樣本中位數(shù)和總體中位數(shù)比較例8.2已知正常人尿氟含量的中位數(shù)為45.30μmol/L,今在該地某廠隨機(jī)抽取12名工人，問(wèn)測(cè)得尿氟含量見(jiàn)表，問(wèn)該廠工人的尿氟含量是否高于當(dāng)?shù)卣Ｈ说哪蚍浚縃0：該廠工人的尿氟含量的總體中位數(shù)Md=45.3H1：該廠工人的尿氟含量總體中位數(shù)Md>45.3

α＝0.05第二節(jié)兩個(gè)獨(dú)立樣本比較Wilcoxon秩和檢驗(yàn)

Wilcoxonranksumtest1．原始數(shù)據(jù)的兩樣本比較

2．頻數(shù)表資料和有序（等級(jí)）數(shù)據(jù)的兩樣本比較1．原始數(shù)據(jù)的兩樣本比較

例8-3對(duì)10名肺癌病人和12名矽肺0期工人用X片測(cè)量肺門(mén)橫徑右側(cè)距RD值，問(wèn)肺癌病人的RD值是否高于矽肺0期工人？表8-5肺癌病人和矽肺0期工人RD值比較

(同一組可直接寫(xiě)秩號(hào))H0：肺癌病人的RD值和矽肺0期總體分布位置相同H1

：肺癌病人的RD值高于矽肺0期工人α＝0.05查附表10(樣本含量較?。┎楦奖?0，n1=10,n2-n1=2得89~141,P<0.05,按拒絕H0,接受H1，肺癌病人的RD值高于矽肺0期工人

2.正態(tài)近似性檢驗(yàn)（n1>10或n2-n1>10)2．有序（等級(jí)）數(shù)據(jù)的兩樣本比較例8-439名吸煙工人和40名不吸煙工人的碳氧血紅蛋白含量見(jiàn)表，問(wèn)吸煙工人的碳氧血紅蛋白含量是否高于不吸煙工人的含量？

表吸煙與不吸煙工人的HbCO(%)含量的比較

①先確定各等級(jí)的合計(jì)人數(shù)、秩范圍和平均秩，再計(jì)算兩樣本各等級(jí)的秩和②本例T=1917；

③計(jì)算u值查附表2（t界值表1.64）得單側(cè)p<0.05，按水準(zhǔn)拒絕H0，接受H1，可認(rèn)為吸煙工人的HbCO（%）含量高于不吸煙工人的HbCO（%）含量。第三節(jié)K個(gè)獨(dú)立樣本檢驗(yàn)完全隨機(jī)設(shè)計(jì)多個(gè)樣本比較的Kruskal-WallisH檢驗(yàn)原始數(shù)據(jù)的多個(gè)樣本比較

Kruskal-WallisH檢驗(yàn)H0:多個(gè)總體分布位置相同

H1:多個(gè)總體分布位置不全相同

例8-5用三種藥物殺滅釘螺，每批用200只活釘螺，用藥后清點(diǎn)每批釘螺的死亡數(shù)，再計(jì)算死亡率，結(jié)果見(jiàn)表，問(wèn)三種藥物殺滅釘螺的效果有無(wú)差別？本例N=15，n1=n2=n3=5，查附表11得H0.05=5.78，按水準(zhǔn)拒絕H0，接受H1，可認(rèn)為三種藥物殺滅釘螺的效果不同。例8.6比較小白鼠接種三種不同菌型傷寒桿菌9D、11C和DSC1后的存活日數(shù)，問(wèn)小白鼠接種三種不同菌型傷寒桿菌存活日數(shù)有無(wú)差別？2．有序（等級(jí)）數(shù)據(jù)的多個(gè)樣本比較

例8-7四種疾病患者痰液內(nèi)的嗜酸性粒細(xì)胞的檢查結(jié)果，問(wèn)四種疾病患者痰液內(nèi)的嗜酸性粒細(xì)胞有無(wú)差別？

有相同秩兩獨(dú)立樣本比較，若n1和n2較大，可用Wilcoxon秩和檢驗(yàn)u檢驗(yàn)公式，也可以用Kruskal-WallisH檢驗(yàn)公式。兩者的關(guān)系是：H(或Hc)=u2二、多個(gè)獨(dú)立樣本作兩兩比較的Nemenyi法檢驗(yàn)

為校正系數(shù)據(jù)查表得P<0.05,可認(rèn)為小白鼠接種11C的存活日數(shù)高于接種9D的存活日數(shù)

據(jù)查表得P<0.05,可認(rèn)為小白鼠接種11DSC1的存活日數(shù)高于接種9D的存活日數(shù)據(jù)查表得P>0.05,不可認(rèn)為小白鼠接種11C的存活日數(shù)與接種DSC1的存活日數(shù)有差別第四節(jié)隨機(jī)區(qū)組設(shè)計(jì)多個(gè)樣本比較的FriedmanM檢驗(yàn)

一、多個(gè)相關(guān)樣本比較的FriedmanM檢驗(yàn)

例8-9

8名受試對(duì)象在相同實(shí)驗(yàn)條件下分別接受4種不同頻率的聲音刺激，他們的反映率資料見(jiàn)表，問(wèn)4種不同頻率的聲音刺激的反映率有無(wú)差別？在每個(gè)配伍組內(nèi)編秩次H0:4種頻率的聲音刺激的反映率總體分布位置相同H1:4種頻率的聲音刺激的反映率總體分布位置不全相同(1)將每個(gè)區(qū)組數(shù)據(jù)編秩(2)計(jì)算各樣本的秩和Ri,平均秩和1.查表法：當(dāng)n<=15和g<=15,查附表12，M界值表本例n=8,g=4查表得105，p<0.05,拒絕H0,4種頻率的聲音刺激的反映率總體分布位置不全相同2.近似法(g>4或g=4,n>5或g=3,n>9)二、兩兩比較的q檢驗(yàn)

H0:任兩反映率總體分布位置相同H1:任兩反映率總體分布位置不同本例n=8,g=4,qvap頻率A與頻率D比較7.87214<0.01頻率A與頻率C比較5.32213<0.01頻率A與頻率B比較2.13212>0.05頻率B與頻率D比較5.75213<0.01頻率B與頻率C比較3.192120.01~0.05頻率C與頻率D比較2.55212>0.05已知總體分布類(lèi)型，對(duì)未知參數(shù)進(jìn)行統(tǒng)計(jì)推斷依賴于特定分布類(lèi)型，比較的是參數(shù)

參數(shù)檢驗(yàn)（parametrictest）

非參數(shù)檢驗(yàn)（nonparametricte