平衡不完全配伍组试验的秩和检验医学统计学教研室(226001)陈峰,陆守曾提要在进行随机 区组试验设计时,若试验的处理数大于每一配伍组中能容纳的观察单位数时,需用平衡不完全配伍 组设计。本文介绍了当资料不服从正态或各处理组方差不齐时所用的秩和检验方法,当每个配伍组 中实验单位数等于处理组数时,该法与Fridman法等价。关键词平衡不完全配伍组试验;秩 和检验在进行随机区组试验设计时,若试验的处理数大于每一配伍组中能容纳的观察单位数时,需 用平衡不完全配伍组(BIB)设计。笔者曾详细讨论了平衡不完全配伍组试验的线性模型解法, 该法基于方差分析中的回归方法,对各处理总和进行调整,以消除区组的干扰[1]。该法要求数 据之分布近似正态且各处理组方差达到齐性。当条件不满足时,不能直接用方差分析。此时可用秩 和检验。1检验步骤与参数方法相对应,以v表处理组数;b表配伍组数;r表每个处理组内的重 复数;k表每个配伍组内的例数。检验步骤如下。1.1编秩对每个配伍组中的k个观察值依大小 顺序从1到k编秩,遇到等值时编以等秩。1.2计算v个处理组的秩合计Rj(j=1,2,… ,v);1.3建方检验假设H0:各处理之效应相同;Hl:各处理之效应不同或不完全相同。 1.4按式(1)计算统计量D。此时D近似服从自由度为v-1的X2分布;1.5各处理效应 的两两比较。在拒绝H。时,各处理之效应Rj的两两比较可用推广了的r检验,统计量t由式( 2)、(4)计算,其中SRj-Rj为秩和差的标准误,处理组与对照组比较时,t值由式(3 )计算。统计量t近似服从自由度为bk-v-b+1的t分布。2D的原分布如各处理之效应相 同,则每个区组中秩次的排列是随机的,共有K种,则b个区组共有(K)b种排列方法,每种排 列之出现为等可能的,其概率为1/(K)b。如每种处理的重复数r足够大,则该处理组之秩和 Rj之分布将趋近于正态分布,从而其中:由此可推得统计量D近似服从自由度式(1)与式(8 )是等价的。特别当k=v,r=b时,BIB设计即为完全随机区组设计,统计量D即为完全随 机区组设计的Friedman秩和检验统计量。在出现等秩时,方差为:3实例分析以《中国医 学百科全书·医学统计学》第57页"平衡不完全配伍组试验"条目中资料为例[3].由于该资 料中9个处理组之方差不齐(X2=84.38,p<0.001),12个区组之方差亦不齐。 现用秩和检验。首先将各区组中3个观察值从小到大依次排序,并求出各处理组之秩和Rj(见表 1)。本组v=9,k=3,b=12,r=4。H0:9种处理的皮质酮含量相同;H1:9种 处理的皮质酮含量不同或不全相同。α=0.05。由式(1)求得:表1大白鼠以500R的x 线照射前后不同时间血浆皮质酮含量的秩次分布按自由度(9-1)查X2界值表即知:P<0. 05。按α=0.05水准拒绝H0,接受H1,故认为照射前后不同时间血浆皮质酮含量不同。 将各实验组之效应与对照组之效应作比较。由式(4)得:SRj-Rj=2.2913,按式( 3)计算t值,列于表2,此时自由度为16,查t界值表,结果见表2。表2各实验组与对照组 比较P<0.05结论:大白鼠经X线照射后,血浆皮质酮含量在5小时左右升至高峰,随后逐渐 降低。照射后12小时基本上接近照射前水平。4几点说明4.1当k=v,r=b时,BTB设 计即为完全随机区组设计,此时,统计量D与Fridman的统计量等价,说明该方法的连续性 。4.2当需分析区组效应时,只要将处理与区组的位置互换,按类似方法分析。4.3当资料服 从正态分布,且各处理组方差达到齐性时,应该用方差分析的方法,此时的秩和检验效能只有方差 分析的95.5%。反之,秩和检验的效能将大于方差分析者。参考文献||1.陈峰.平衡不完 全配伍组试验的简便分析法.中国卫生统计1993;10(1):16.2.ConoverW J.Practicalnonparametricstatistics.2ed.USA: JohnWiley&Sons,1980:309~315.3.杨树勤主编.中国医学百科全 书.医学统计学.第1版.上海:科学技术出版社,1985:57~62.THERANKSU MTESTOFTHEBALANCEDINCOMPLETEBLOCKTRIAL¥Chen Feng;LuShouzeng(DepartmentofMedicalStatisti cs,NantongMedicalCollege,Jiangsu)Abstract:I ntherandomizedcompleteblockdesigneverytreat mentisappliedineveryblock.However,itissomet imesimpracticalorimpossibleforallofthetreat menttobeappliedtoeachblock.Inthesesituation s,thebalancedincompleteblock(BIB)designisem ployed.Thispaperpresentedaranksumtestthatma ybeusedifthenormalityorhomogeneityofvarianc eassumptionsarenotmet.Theranksumtestreduces totheFridman'stestifthenumberoftreatmentseq ualsthenumberofexperimentalunitsperblock.Th emethodisillustratedwithanexample.Keywords: thebalancedincompleteblock(BIB)trial;ranksu mtest平衡不完全配伍组试验的秩和检验@陈峰,陆守曾平衡不完全配伍组试验;秩和检验在 进行随机区组试验设计时,若试验的处理数大于每一配伍组中能容纳的观察单位数时,需用平衡不 完全配伍组设计。本文介绍了当资料不服从正态或各处理组方差不齐时所用的秩和检验方法,当每 个配伍组中实验单位数等于处理组数时,该法与Fridman法等价。1.陈峰.平衡不完全配 伍组试验的简便分析法.中国卫生统计1993;10(1):16.2.ConoverWJ. Practicalnonparametricstatistics.2ed.USA:Jo hnWiley&Sons,1980:309~315.3.杨树勤主编.中国医学百科全书. 医学统计学.第1版.上海:科学技术出版社,1985:57~62.block.Inthe sesituations,thebalancedincompleteblock(BIB )designisemployed.Thispaperpresentedaranksu mtestthatmaybeusedifthenormalityorhomogenei tyofvarianceassumptionsarenotmet.Theranksum testreducestotheFridman'stestifthenumberoft reatmentsequalsthenumberofexperimentalunits perblock.Themethodisillustratedwithanexampl e.Keywords:thebalancedincompleteblock(BIB)t rial;ranksumtest平衡不完全配伍组试验的秩和检验@陈峰,陆守曾平衡不完全 配伍组试验;秩和检验在进行随机区组试验设计时,若试验的处理数大于每一配伍组中能容纳的观 察单位数时,需用平衡不完全配伍组设计。本文介绍了当资料不服从正态或各处理组方差不齐时所 用的秩和检验方法,当每个配伍组中实验单位数等于处理组数时,该法与Fridman法等价。 1.陈峰.平衡不完全配伍组试验的简便分析法.中国卫生统计1993;10(1):16.2 .ConoverWJ.Practicalnonparametricstatistics.2ed.USA:JohnWiley&Sons,1980:309~315.3.杨树勤主编.中国医学百科全书.医学统计学.第1版.上海:科学技术出版社,1985:57~62.
More abstracts about the 平衡不完全配伍组试验的秩和检验