闂傚倸鍊搁崐鎼佸磹閹间礁纾归柟闂寸绾惧綊鏌熼梻瀵割槮缁惧墽鎳撻—鍐偓锝庝簼閹癸綁鏌i鐐搭棞闂囧鏌ㄥ┑鍡欏妞ゅ繒濮风槐鎺楀焵椤掍胶绡€闁稿本顨嗛弬鈧梻浣虹帛钃辩憸鏉垮暣閸┾偓妞ゆ巻鍋撴い鏇ㄥ幘濡叉劙骞橀幇浣瑰兊闁哄鐗勯崝宀勫几閹达附鈷戦柛婵嗗閺嗘瑩鏌eΔ鈧€氫即宕洪埀顒併亜閹达絾纭剁紒鎰⒒閳ь剚顔栭崰鏇犲垝濞嗘挶鈧礁顫滈埀顒勫箖濞嗘挻鍤戦柛銊︾☉娴滈箖鏌涢…鎴濇灀闁衡偓娴犲鐓ユ繛鎴灻鈺伱归悩娆忔噽绾惧吋銇勯弴鐐村櫣闁诲骏濡囬埀顒侇問閸n噣宕戞繝鍥х畺闁冲搫鍟扮壕鍏间繆椤栫偞鏁遍悗姘偢濮婂宕掑▎鎴g獥闂佺ǹ顑呭Λ婵嗙暦閹存績妲堥柕蹇娾偓铏吅婵$偑鍊栭悧妤冪矙閹烘垟鏋嶉柣妯肩帛閸婂灚绻涢幋鐐垫嚂濞撴碍绋掔换娑㈠醇濠婂懐鐓撻梺鍝勭焿缂嶄線鐛崶顒夋晢濠㈣泛锕ラ宥呪攽閻橆偅濯伴柛銉檮閸犳岸姊洪崫鍕缂佸鍏樼瘬濞撴埃鍋撻柡灞剧洴楠炴ḿ鈧稒岣垮▓銈囩磽娓氬洤鏋熼柣鐔叉櫊閻涱噣骞掑Δ鈧崡鎶芥煟閺冨伋鐟扳枍閿熺姵鈷掗柛灞捐壘閳ь剛鍏橀幃鐐烘晜閻e瞼顔曟繝鐢靛Т濞层倗绱掗埡鍛骇闁割偅纰嶅▍鍡涙煛閸℃ḿ鐭掗柡灞诲姂閹垽宕崟鎴欏灲閺岋箓宕橀鍕亪闂佸搫鑻粔褰掑春閳╁啯濯撮柛娑橈攻椤斿鈹戦悩鎰佸晱闁搞劑浜堕獮鎰板礃瀹割喗缍庢繝鐢靛У閼归箖鏌嬮崶顒佺厽闁哄啫娲㈤埀顑藉亾闂佸搫顑嗗Λ鍐潖閾忓湱鐭欐繛鍡樺劤閸撻亶姊洪悷鐗堝暈闁诡喖鍊垮畷娲倷閸濆嫮顓洪梺鎸庢磵閸嬫捇鏌i幇顒婅含闁哄被鍔戝顒勫垂椤旇瀵栭梻浣风串缂嶄胶绮婚弽褜娼栫紓浣股戞刊鎾煕濞戞﹫榫氱痪顓涘亾闂傚倷鐒﹂崕宕囨崲閹扮増鐓€闁挎繂鎳愰弳锕€霉閸忓吋缍戦柛鎰ㄥ亾婵$偑鍊栭幐楣冨窗鎼淬劌鐒垫い鎺嗗亾婵炵》绻濆濠氭晲婢跺⿴娼婇梺缁樏崯鍧楀汲閸儲鈷戠紓浣股戠亸顏堟煕鎼达絾鏆€殿喖顭烽弫鎰緞婵犲倸绁梻浣虹帛閸旀牜绮婇幘顔肩;闁规崘顕ч幑鑸点亜閹惧鈽夊ù婊呭亾閹便劌螣閻撳簼澹曢柣搴㈢啲閹凤拷缂傚倸鍊搁崐鎼佸磹閹间礁纾归柟闂寸绾惧綊鏌熼梻瀵割槮缁炬儳缍婇弻鐔兼⒒鐎靛壊妲紒鐐劤缂嶅﹪寮婚敐澶婄闁挎繂鎲涢幘缁樼厱濠电姴鍊归崑銉╂煛鐏炶濮傜€殿喗鎸抽幃娆徝圭€n亙澹曢梺褰掓?閻掞妇鈧艾鎳橀弻锝夊棘閸喗鍊梺缁樻尭缁绘劙鍩為幋锔藉€烽梻鍫熺◥婢规洖鈹戦悙鍙夊櫤闁圭懓娲濠氬焺閸愨晛顎撶紓浣割儐椤戞瑦瀵奸崘顔解拺闁告繂瀚﹢鎵磼鐎n偄鐏撮柛鈹垮劜瀵板嫰骞囬鍌ゆ敤闂備胶绮崝妯间焊濞嗘搩鏁婇柟鐑橆殕閳锋帒霉閿濆懏鍟為柟顖氱墛娣囧﹪顢曢姀銈呭及閻庤娲橀崝鏇㈠煘閹寸姭鍋撻敐搴濇捣闁硅姤娲栭埞鎴︽倷閺夋垹浠ч梺鎼炲妼缂嶅﹪骞冮悙瑁佹椽顢旈崨顖氬箞闂備浇顫夊妯绘櫠鎼达絽顥氶悹鎭掑妿绾惧ジ鏌涘▎鎯奉亜鐣甸崱妯诲弿濠电姴鍊归幆鍫ユ煟閿濆洤鍘寸€规洖鐖奸崺鈩冩媴娓氼垰浠规繝鐢靛Х椤h棄危閸涙潙纾婚柛鏇ㄥ灠閻ゎ噣鏌熺紒銏犳灓鐎规挷绶氶弻娑㈠箛閳轰礁顥嬮梺鍝勫暙閻楀棗顔忓┑鍥ヤ簻闁哄啫娲よ闂佺粯绻冩繛濠傤潖閾忚鍏滈柛娑卞幒濮规鏌i悙瀵糕棨闁稿海鏁诲畷娲焵椤掍降浜滈柟鍝勬娴滃墽绱掗崜褑妾搁柛娆忓暣閵嗕礁鈻庨幘宕囶吅闂佺粯岣块弫鎼佹晬濠靛鈷戠紒瀣濠€浼存煟閻曞倻鐣甸挊婵嗩渻鐎n亝鎹g痪鎯с偢濡懘顢楁担鍝ョシ闂侀€炲苯澧紒璇茬墦閹即顢欑捄銊ф澑濠电偞鍨堕悷锕€鈻嶉姀銈嗏拺閻犳亽鍔屽▍鎰版煙閸戙倖瀚�闂傚倸鍊搁崐鎼佸磹閹间礁纾归柟闂寸绾惧湱鈧懓瀚崳纾嬨亹閹烘垹鍊為悷婊冪箻瀵娊鏁冮崒娑氬幈濡炪値鍘介崹鍨濠靛鐓曟繛鍡楃箳缁犳娊鏌嶈閸撴瑧绮诲澶婄?闂侇剙鍗曢崶顒夋晬婵犲﹤鎳愰悞濂告⒑閸涘﹤濮€闁哄倸鍊圭粋宥呪堪閸曗晙绨婚梺瑙勫礃濞夋盯寮告惔銊︾厽闊洢鍎崇粔顕€鏌″畝鈧崰鏍€佸璺哄耿婵炲棗绻愰悡鍌炴⒒娴e摜鏋冩い顐㈩樀瀹曞綊宕稿Δ鈧粻鏍煃閸濆嫭鍣圭紒鐘崇☉闇夐柨婵嗘噺閸熺偞绻涢崨顔炬创婵﹦绮粭鐔煎焵椤掆偓宀e潡鎳¢妶鍥╊槸闂佺硶鍓濈粙鎴犵玻濡や椒绻嗛柕鍫濇噹閺嗙喖鏌i鐑嗘Ш闁逞屽墮缁犲秹宕曢柆宓ュ洦绂掔€n亝鐎梺绋跨灱閸嬬偤鎮¢弴鐔翠簻闁归偊鍠栧瓭闂佽绻嗛弲婵堟閹烘挾鐟瑰┑鐘插閺嗐倝姊洪崫鍕効缂佹煡绠栭獮鍡涘籍閸偅鏅┑鐘诧工閹叉盯鏁愭径瀣ф嫼闂佸憡绺块崕杈ㄧ墡闂備胶绮〃鍡欏垝鎼淬劌绀嗛柟鐑樻尵缁♀偓濠殿喗锕╅崢钘夆枍濠婂牊鈷戠紒顖涙礀婢ф煡鏌i悢婵嗘处閸庡﹪鏌涢銈呮灁缂佲檧鍋撻梻浣圭湽閸ㄨ棄岣胯閻☆參姊绘担鍛靛綊鏁冮妷鈺佺畺闁稿瞼鍋戦埀顑跨窔瀵噣宕奸锝嗘珖闂備焦瀵у濠氬疾椤愶箑鍌ㄩ柣銏犳啞閳锋垿鏌涘☉姗堝伐闁诲繒鍠栭弻娑欑節閸愵亜鈷堥梺閫炲苯澧柣鏃戝墴楠炲繘鏁撻敓锟�

第三节 秩和检验

《医学统计学》在线阅读中医杂集书籍在线阅读

用秩号代替原始数据后,所得某些秩号之和,称为秩和,用秩和进行假设检验即为秩和检验。其检验假设在两组比较(成对或不成对)时,H0:F(X1)=F(X2),即两总体的分布函数相等,备择假设H1:F(X1)≠F(X2)。本法由于部份地考虑了数据的大小,故检验效力较符号检验大大提高。至于其方法、步骤,不论是查表法或计算法、也都相当简便,现举例说明如下。

相关

一、成对资料的比较

此法由Wilcoxon氏首次提出,故又称Wilcoxon氏法。

处理时可用查表法或计算法,今以例10.3分别说明如下。

查表法步骤:

1.排队,将差数按绝对值从小至大排列并标明原来的正负号,见表10.3第(5)栏,排队后与原豚鼠号已无对应关系。

2.编秩号,成对资料编秩号时较为复杂,要注意三点:

(1)按差数的绝对值自小至大排秩号,但排好后秩号要保持原差数的正负号;

(2)差数绝对值相等时,要以平均秩号表示,如表10.3中差数绝对值为4者共三人,其秩号依次应为2、3、4,现皆取平均秩号3;

(3)差数为0时,其秩号要分为正、负各半,若有一个0,因其绝对值最小,故秩号为1,分为0.5与-0.5,若有两个0,则第二个0的秩号为2,分为1与-1等等。

3.求秩号之和即将正、负秩号分别相加,本例得正秩号之和为68,负秩号之和为10,正负秩号绝对值之和应等于1/2n(n+1),可用以核对,如本例68+10=12/1(12+1)=78,证明秩号计算正确。

4.以较小一个秩号之和(R),查附表12进行判断,该表左侧为对子数,表身内部是较小秩号和,与上端纵标目之概率0.05,0.01相对应,其判断标准是

R>R0.05时P>0.05

R0.05≥R>R0.01时0.05≥P>0.01

P≤R0.01时 P≤0.01

例10.3 请以表10.1资料用秩和检验处理之。

表10.3 豚鼠给药前后灌流滴数及其秩号

豚鼠号
(1)
每分钟灌流滴数按差数绝对值排队
(5)
秩号
用药前
(2)
用药后
(3)
差数
(4)

(6)

(7)
 
1304616-2 1
2385012-4 3
34852443 
44852443 
56058-2-8 6
646641886 
726563086 
85854-4108 
946548129 
104858101610 
114436-81811 
12465483012 

68 R=10

将表中10.1中用药前后的数据求出差数,并按差数绝对值排队,结果见表10.3第(5)栏。再编秩号,为计算方便,正、负秩号分列两栏,见表10.3第(6)、(7)栏。

上例,n=12,?OR?O=10,查附表12得

R0.05=14R0.01=7

今R0.05>R>R0.01,故0.05>P>0.01,在概率0.05水平上拒绝H0,接受H1,即用药前后的相差是显著的,给药后每分钟灌流滴数比用药前增多了。

附表12中只列有n≤25时的临界值。当n值较大时亦可采用计算法。

计算法步骤:

在计算法时,对差数的排队,编秩号及求秩号之和同查表法,不同的是求得秩号之和以后的算,所用公式是:

u0.05=1.96u0.01=2.58 (10.5)

式中n为原始资料中数据的对子数,R为正秩号之和或负秩号之和,为计算方便,通常取绝对值较小的秩号之和为r 。

本例,n=12,R=-10,代入得:

U0.05<U<0.01,故0.05>P>0.01,在α=0.05水准上拒绝H0,接受H1,结论与查表法相同。

据研究,当n大于10时,上式算得的u近似正态分布,故计算法只用于n值较大时。

因本例资料接近正态分布,故曾用t检验的个别比较方法处理过,结果是:t=2.653 0.05>P>0.01,与秩和检验结论相同,但与符号检验结论不同(χ2=2.083,P>0.05),说明符号检验的检验效率比秩和与t检验都要低,比较粗糙,而秩和检验的效率与t检验较接近。

相关

二、两组资料的比较

此法又称为wilcoxon氏两样本法。

处理时也可用查表法或计算法,今以例10.4分别说明之。

查表法步骤:

1.各自排队,统一编秩号,即将两组数据分别从小到大排列,但编秩号时要两组统一进行,凡分属于两组的相等数据用平均秩号,如本例0.042共三个,取平均序号皆为8。

2.令较小样本秩号之和为r ,例数为n

3.计算R',公式为:

R'=n1(n1+n2+1)-r (10.6)

R'是同一个样本资料,当秩号倒排(即由大至小)时较小样本秩号之和。

4.以R和R'两秩号之和中较小者与附表13中R的临界值比较,以作出判断,其标准仍是:

R>R0.05时 P>0.05

R0.05≥R>R0.01时 0.05≥P>0.01

P≤R0.01时 P≤0.01

例10.4 请以表10.2资料用本法处理之。

表10.4 九名健康人与八名铅作业工人的尿铅值(mg/L)

健康人秩号铅作业工人秩号
0.00110.0428
0.00220.0428
0.01430.04810
0.02040.05011
0.03250.08214
0.03260.08615
0.04280.09216
0.054120.09817
0.06413  
n2=954n1=8R=99

先将本表10.2中两组数据各自排队并统一编秩号,结果见表10.4。

较小样本为铅作业工人组,n1=8,R=99,代入式(10.6)

R'=8(8+9+1)-99=45

R与R'两者中以R'较小,故以P'值与附表13数值比较,得R0.05=51,R0.01=45;今R'=R0.01,故P=0.01,在α=0.05水平上拒绝H0,接受H1,差别显著,故铅作业工人尿铅值比健康人高。

计算法步骤:

两组资料比较时,也可用计算法。用计算法时,对两组数据各自排队、统一编秩号同查表法,不同的是求得秩号之和以后计算,公式是:

u0.05=1.96u0.01=2.58 (10.7)

为便于计算和前后符号一致,n1作为较小样本例数,R为较小样本的秩和,n2则为较大样本的例数。

本例n1=8,R=99,n2=9代入公式得:

今?Ou?O>u0.01,故P<0.01,在α=0.01水准上拒绝H0接受H1,其结论同查表法,

据研究,当n1、n2都大于8时,算得的u近于正态分布,若例数太少,则以查表法更为精确。

本例如用t检验的团体比较处理,则t=3.169,P<0.01,二者结论一致,但与符号检验结论不同(χ2=2.930,P>0.05)同样说明符号检验较粗糙,检验效率低,而秩和检验与t检验的结论较近。

相关

三、两组等级资料的比较

等级资料又称为半计量资料,当两组等级资料比较时,用秩和检验来比较其相差是否显着比用χ2检验要恰当。两组等级资料,通常例数都较多,故一般都用计算法,其步骤与两组资料的秩和检验相似,不同的是要求各等级的平均秩号,为此,先要求得各等级的秩号范围。今举例10.5说明之。

1.求各等级的平均秩号。为此,先要求出各等级的秩号范围,如等级“-”共18+8=26例,共秩号范围自1~26。要注意的是各等级的秩号范围必须紧相联接。最后一组秩号范围的上限一定等于两组例数之和。求得各等级秩号范围后,再求其下限和上限的平均,即可算得平均秩号,如等级“一”的平均秩号为(1+26)/2=13.5。余类推。

2.求出R及其n1,为计算方便,把例数少的正常人组的秩号之和作为R其例数为n1得R=308,n1=20,n1=32

3.代入式(10.7)得u值,即可作结论。

例10.5,今有20名正常人和32名铅作业工人尿棕色素定性检查结果如下表10.5,试问其相差是否显著?

表10.5 20名正常人和32名铅作业工人尿棕色素定性检查结果

尿棕色素定性结果正常人铅作业工人合计秩号范围平均秩号例数较小组的秩和
-188261—2613.5243
+2101227—3832.565
++7739—4542.0
+++3346—4847.0
++++4449—5250.5

n1=20 n2=32 R=308

代入式(10.7)

u0.01=2.58,今u>u0.01,故P<0.01,在α=0.01水准上拒绝H0,接受H1。两组相差显著,铅作业工人尿棕色素比正常人为高。

相关 查疾病 查药材 查方剂 查穴位

四、多组资料的比较

多组资料的比较也是从排秩号开始,但不是直接用秩和进行检验,有的书籍称之为秩检验(rank test),以示与秩和检验有别,其检验假设也较复杂:在处理完全随机设计的资料时,H0:F(X)=F(X2)=F(X3)=……,即比较的各样本所对应的各总体的分布函数相等,H1:各总体的分布函数不相等或不全相等;在处理随机单位组设计的资料时,H0:P(χij=r)=1/n,即内组各秩号r之概率相等,都是1/n(r=1,2,……,n)而H1为:P=(χij=r)≠1/n。

因不同实验设计所得资料的处理也有别,故下面分别举例说明之。

(一)完全随机设计所得资料的比较

用的方法是单因素多组秩检验,称为Kruskal-Wallis氏法,或H检验。其计算步骤如下。

1.各自排队,统一编秩号。即将各组数据在本组内从小到大排队,见表10.6各含量栏,再将各组数值一起考虑编出统一秩号,见表10.6各“秩号”栏,分属不同组的相同数值用平均秩号;
2.求各组秩号之和R1以及各组数n1
3.代入下式计算H值:

(10.8)

式中N为各组例数之和,Ri和ni为各组的秩号之和以及例数:

4.查表作结论

当比较的组数多于三组,或组数虽只有三组但每组例数大于5时,H值的分布近于自由度等于组数-1的χ2分布,故可用对应的χ2值作界值。当三组比较时每组例数均不超过5时,H值与χ2值有较大偏离,此时可查附表14,直接查得H0.05和H0.01

例10.6 雄鼠20只随机分为四组,第1、2组在皮肤上涂用放射性锡(Sn113)标记的三乙基硫酸锡,涂后将皮肤暴露于空气中;第3、4组涂药后用密闭小玻璃管套使皮肤与外界空气隔开,三小时后杀死,测肝中放射物,结果如表10.6,试比较各组含量间有无显著相差?

表10.6 白鼠皮肤涂药后,肝中放射性Sn113的含量

涂干药后敞开涂湿药后敞开涂干药后密闭涂湿药后密闭

含量秩号含量秩号含量秩号含量秩号
0.0011.82110.6653.6714
0.422.52.79120.7164.4616
0.422.53.07130.7574.5118
0.5944.19150.8385.0719
0.9794.47171.49106.0220
RiR1=19R2=68 R3=36 R4=87 
nin1=5n2=5 n3=5 n4=5 

各组资料各自排队,统一编秩号,以及求各组的秩号之和Ri和例数ni见表10.6

代入式(10.8)得

本例组数为4(>3),查χ2值表,ν=4-1=3,得χ20.05,3=7.81,χ20.01,3=11.34,今H>χ20.01,3故P<0.01,在α=0.01水准上拒绝H0,接受H1,即各组肝中放射性Sn113含量差别显著。

(二)随机单位组设计所得资料的比较

用的方法是双因素多组秩检验,即Friedman氏法。

处理这种资料时可分成两步,对两个因素分别进行检验。现用例10.7说明其计算步骤:

先比较四种防护服对脉搏的影响

1.将穿四种防护服的每一受试者的脉搏数从小到大编秩号,当数值相等时用平均秩号,见表10.7各秩号栏。

2.求各防护服组秩号之和Ri

3.代入式10.9求H值

(10.9)

式中t(treatment)为处理组数,b(block)为单位组数。

4.查表作结论

当t>4或t=4且b>5或t=3且b>9时,H值的分布近于自由度ν=t-1时的χ2分布,故可查相应的χ2值与H值比较作出判断:如t、b不能满足上述条件,则所算得的H值与χ2分布有较大偏离,需查附表15作判断。

例10.7 受试者5人,每人穿四种不同的防护服时的脉搏数如表10.7,问四种防护服对脉搏的影响有无显著差别?又五个受试者的脉搏数有无显著差别?

表10.7 比较穿四种防护服时的脉搏数(次/分)

受试者防护服A防护服B防护服C防护服D
编 号脉搏秩号脉搏秩号秩号秩号脉搏秩号
1144.44143.03133.41142.82
2116.22119.24118.03110.81
3105.81114.83113.22115.84
498.01120.03104.02132.84
5103.82110.64109.83100.61
秩秩号和Ri 10 17 11 12

t=4b=5

排队、编秩号、求各比较组的Ri见表10.7所示。

将表10.7中各数代入式10.9,得

本例t=4,b=5查附表15,得H0.05=7.80,今H>H0.05,故P>0.05,在α=0.05水准上接受H0,无显著差别,故四种防护服对脉搏的影响无显著差别。

再比较五名受试者的脉搏数:

将数据列出(同表10.7),但秩号是按每种防护服中受试者脉搏的数值从小到大编定,然后求出各受试者秩号之和R,详细见表10.8

表10.8 比较五名受试者的脉搏数

受试者防护服A防护服B防护服C防护服DRi
编 号脉搏秩号脉搏秩号脉搏秩号脉搏秩号
1144.45143.05133.45142.8520
2116.24119.23118.04110.8213
3105.83114.82113.23115.8311
498.01120.04104.01132.8410
5103.82110.61109.82100.616

t=5b=4

将表10.8 所得各数据代入式10.9得

此处t>4,故查ν=5-1=4时的χ2值表,得:χ20.05,4=9.49,χ20.01,4=13.28,今χ20.05,4<H<X20.01,4,故0.05>P>0.01,在α=0.05水准上拒绝H0,接受H1差别显著;即五名受试者脉搏数相差显著,1号受试者最高,5号受试者最低。

相关 查疾病 查药材 查方剂 查穴位

五、多组资料间的两两比较

当多组间的差别显着时,则需进一步判断那些组之间的差别有显着性,这个问题的解决方法与第八章第二节中的多个均数间的两两比较很相似,在例10.6四个实验组涂放射性锡的例子中,结果为H>χ20.01,3,P<0.01,现以此为例,进一步作各组两两间比较,步骤如下:

1.将各组秩和从大到小依次排队,并求得两两间的相差,见表10.9

2.计算标准误,计算公式是:

(10.10)

式中σ为任意两个秩和之差的标准误,n为各组例数,a为处理数,此式要求各组例数相等,

3.查q值表定界限作结论

仍查方差分析时用的q值表,v→∝

各q值须与处理数相同的标准误相乘,如处理数为2的q值要乘以处理数为2时的标准误,2.77×6.77=18.75,3.64×6.77=24.64等,余类推。

例10.6资料两两间比较如下:

表10.9 每两组秩和之间的相差及其显著性

组别秩和RiRi—19Ri—36Ri—68
涂湿药后密闭8768**51**19*
涂湿药后敞开6849**32** 
涂干药后密闭3617  
涂干药后敞开19   

计算标准误:n=5,用式10.10

查q值表,得:

处理数234
q0.05,∞2.773.313.63
q0.01,∞3.644.124.40
q0.05,∞σ18.7533.1048.02
q0.01,∞σ24.6441.2058.21

两两比较后的结论见表10.9所示,结合起来看,结论是:涂湿药的比涂干药肝中放射性Sn113含量要高,涂湿药中,密闭的比敞开的含量高。

1000余本中医古籍txt电子书免费下载

下载《医学统计学》.chm chm电子书打不开?

下载所有中医书籍

《医学统计学》相关章节: