秩和检验最新版ppt课件VIP

.,秩和检验,ranksumtest,.,分析资料的步骤：1、确定资料的类型：分类资料、定量资料；2、选择适当的统计方法；3、作出统计推断结论。,.,定量资料的分析：,t检验Z检验F检验,分类资料的分析：,检验,.,检验的应用范围比较两个或多个样本的总体率（构成比）有无差异分析两个分类变量之间的关联性判断频数分布类型分析率的线性趋势,.,t检验的应用条件未知且n较小样本取自正态总体两样本均数比较时，两样本总体方差相等Z检验的应用条件样本取自正态总体，两样本总体方差相等未知但n较大(如n100)，n较小但已知。F检验的应用条件各样本取自正态总体各样本的总体方差相等各样本是相互独立的随机样本,.,例某实验室观察局部温热治疗小鼠移植肿瘤的疗效，以生存日数作为观察指标，实验结果见表1，试检验两组小鼠生存日数有无差别？,.,表1两组小鼠发癌后生存日数,.,内容：非参数检验的概念Wilcoxon符号秩和检验成组设计的两样本比较的秩和检验多个样本比较的秩和检验多个样本两两比较的秩和检验随机区组设计的秩和检验,.,参数统计和非参数统计,.,非参数检验适用范围,1.等级资料（有序分类资料）。如疗效按治愈、显效、有效、无效分组的资料；临床化验结果按“、+、+、+”分组的资料等。（比较构成比有无差异）2.偏态分布资料。当观察值呈偏态或极度偏态分布，而又未经变量变换或虽经变换但仍未达到正态或近似正态分布。3.分布不明的资料。如新指标分布形态不明；小样本，但不趋向正态分布资料。4.各组方差明显不齐，且不易变换达到齐性。5.组内个别观察值偏离过大的资料。6.开口分组资料。数据分组某一端或两端无明确数值的资料，只给出一个下限或上限，而没有具体数值，如60岁等。,.,非参数检验的优缺点,优点：1、不受总体分布的限制，适用范围广。2、不受资料类型的限制，可以是计量资料、等级资料以及一端或两端无确切数据的资料。,.,缺点：1、对于符合参数检验条件的资料，首先参数检验方法；如果用非参数检验，没有充分利用资料提供的信息，检验效能（power）低于参数检验。2、非参数检验一般犯第二类错误的概率比参数检验大，若要使相同，非参数检验要比参数检验需要更多的样本例数。,.,非参数检验方法:秩和检验(ranksumtest)等级相关分析符号检验游程检验,秩、秩和,秩1456811222532354045123456789101112秩和=1+3+4+5+8+9=30=2+6+7+10+11+12=48,.,秩和检验（ranksumtest）将数据从小到大排序，该序号在统计学上称为秩（rank），用数据排序的秩来代替原数据进行假设检验，这种方法称为秩和检验。,秩(rank),秩和(ranksum),秩和检验(ranksumtest),.,1,4,10,2,3,11,13,12,14,15,5,6,7,8,9,16,.,Wilcoxon符号秩和检验（Wilcoxonsignedranktest),配对设计计量差值比较单一样本与总体中位数比较,.,配对设计资料主要是对差值进行分析。通过检验配对样本的差值是否来自中位数为0的总体，来推断两个总体中位数有无差别，即推断两种处理的效应是否不同。,配对设计两样本比较,.,配对设计的两样本比较例12.1某研究用甲、乙两种方法对某地方性砷中毒地区水源中砷含量(mg/L)进行测定，检测10处，测量值如表12.1的(2)、(3)栏。问两种方法的测定结果有无差别？,.,表12.1甲、乙两种方法测定某地区10处水源中砷含量的结果(mg/L),配对差值经正态性检验得W=0.717，P=0.001,.,（1）建立检验假设，确定检验水准。H0：两种方法测定结果差值的总体中位数为0H1：两种方法测定结果差值的总体中位数不等于0=0.05,.,（2）计算检验统计量T1、求差值，见表12.1（4）2、编秩a：将差值的绝对值从小到大编秩b:当差值为0，则不计该例（n随之减少）c:如两差值相等，则取其平均秩次3、给秩添加正负号4、求秩和并确定检验统计量：分别求出差值为正与负的秩和，任取差值为正或负秩和作为统计量。本例T21.5。,.,.,（3）确定P值，并作出统计推断。A.当n50时，查附表10，T界值表。N=90.108370.055400.023420.01144检验统计量T值在上、下界值范围内，其P值大于相应的概率水平；若T值在上、下界值范围上或范围外，则P值等于或小于相应的概率水平。原则：内大外小本例T=21.5,.,以差值不等于0的数值对子数n=9查附表10,得P0.10,按照=0.05水准，不拒绝H0，差别无统计学意义。尚不能认为甲、乙两种方法测定水源中砷含量有差别。,.,由附表10可知，当n50时，如何确定P值？,注意：,.,（2）正态近似法随着n的增大，T分布逼近均数为n(n+1)/4、方差为n(n+1)(2n+1)/24的正态分布，用正态近似法(u检验），按式（12.1)计算Z值。（12.1）,.,当相同秩次（相持tie）多时（不包括差值等于0），则用（12.2）校正式：（12.2）,.,Wilcoxon配对符号秩和检验的基本思想：在配对样本中，由于随机误差的存在，其对差值的影响不可避免。假定两种处理的效应相同，则差值的总体分布为对称分布，并且差值的总体中位数为0。若此假设成立，样本差值的正秩和与负秩和应相差不大，均接近n(n+1)/4；当正负秩和相差悬殊，超出抽样误差可解释的范围时，则有理由怀疑该假设，从而拒绝H0。,.,单一样本与总体中位数比较Wilcoxon符号秩和检验目的是推断样本中位数与已知总体中位数(常为标准值或大量观察的稳定值)有无差别常用于不满足单样本检验应用条件的资料。,.,例12.2某医生从其接诊的不明原因脱发患者中随机抽取14例，测得其发铜含量(g/g)见表12.2。已知该地健康人群发铜含量的中位数为11.2g/g。问脱发患者发铜含量是否低于健康人群？,.,.,（1）建立检验假设，确定检验水准。H0：差值的总体中位数等于0，即脱发患者发铜含量与该地健康人群相同H1：差值的总体中位数小于0，即脱发患者发铜含量低于该地健康人群单侧=0.05,.,（2）计算检验统计量T求差值：编秩：方法同例12.1。求秩和：T+=4，T-=101确定检验统计量T，本例T=4或T=101。,.,（3）确定P值，并作出统计推断。查附表10T界值表(配对比较的符号秩和检验用)，本例n=14，得单侧P0.005，按=0.05水准，拒绝H0，接受H1，差别有统计学意义，故可以认为脱发患者发铜含量低于该地健康人群。,.,成组设计的两样本比较（Wilcoxonranksumtest),目的是推断连续型变量资料或有序变量资料的两个独立样本代表的总体分布位置是否有差别。,.,（1）原始数据的两样本比较例12.3某地职业病防治欲比较使用二巯基丙磺酸钠与二巯基丁二酸钠的驱汞效果。将22例汞中毒患者随机分配到两组，分别测定并计算出两组驱汞的排汞比值，并将结果列于表12.3。试问两药驱汞效果有无差别？,.,表12.3两种驱汞药物排汞效果比较,该资料为比值数据，不服从正态分布，现采用Wilcoxon秩和检验。,.,检验步骤：（1）建立检验假设，确定检验水准。H0:两种药物排汞比值的总体分布位置相同H1:两种药物排汞比值的总体分布位置不同=0.05,.,（2）计算检验统计量T将两组数据由小到大统一编秩两组中，相同数值在不同组，则必须取平均秩同一组内，可顺次编秩次。分别求出两组秩次的和取较小样本的秩和记为T，本例n2n1,T=75.5,.,（3）确定P值，并作出统计推断。查附表11T界值表(两样本比较的秩和检验），P0.01，按=0.05水准拒绝H0，可认为两种驱汞药物的排汞效果有差别。丁二酸纳组平均秩次为75.5/10=7.55，丙磺酸钠组平均秩次为177.5/12=14.79，可认为丙磺酸钠驱汞效果好于丁二酸钠。,.,若n1、n2-n1超出附表11的范围，则用正态近似法，可按式（12.3）计算Z值：(7.3)式中，N=n1+n2,0.5:连续性校正系数。,.,若持较多（如超过25），则用式（12.4）校正。（12.4）,tj为第j次相持时相同秩次的个数。,.,秩次(rank)，秩统计量是指全部观察值按某种顺序排列的位序；在一定程度上反映了等级的高低；秩和(ranksum)同组秩次之和；在一定程度上反映了等级的分布位置。对等级的分析，转化为对秩次的分析。秩和检验就是通过秩次的排列求出秩和，进行假设检验。,.,（2）等级资料的两样本比较例12.4某医生欲比较中西医疗法与西医疗法治疗急性肾盂肾炎的临床疗效，将患者随机分为两组，分别给予中西医疗法或西医疗法治疗，并观察疗效，结果见表12.4，问两种疗法疗效是否有差别？,.,表12.4两种疗法治疗急性肾盂肾炎的疗效,.,检验步骤：（1）建立检验假设，确定检验水准。H0:两种疗法治疗急性肾盂肾炎的疗效总体分布位置相同H1:两种疗法治疗急性肾盂肾炎的疗效总体分布位置不同=0.05,.,（2）计算检验统计量T求秩和：确定各等级的秩次范围，计算平均秩次；将各组段的平均秩次分别与各等级例数相乘，再求和得到T1，T2。,.,表12.4两种疗法治疗急性肾盂肾炎的疗效,.,计算Z值由于相持较多，故按式（12.3）和式（12.4）计算Zc值。,.,（3）确定P值，并作出统计推断。Zc2.1415，0.02P0.05，按=0.05水准，拒绝H0，差别有统计学意义，可认为两种疗法疗效分布不同。中西医疗法组平均秩次为6820/92=74.1，西医疗法组平均秩次为6060/68=89.1，可以认为中西医疗法治疗急性肾盂肾炎效果较好。,.,Wilcoxon秩和检验的基本思想：假设含量为n1与n2的两个样本(且n1n2)，来自同一总体或分布相同的两个总体，则n1样本的秩和T1与其理论秩和n1(N+1)/2相差不大，即T1-n1(N+1)/2仅为抽样误差所致。当二者相差悬殊，超出抽样误差可解释的范围时，则有理由怀疑该假设，从而拒绝。,.,成组设计的多个样本比较的秩和检验Kruskal-Wallistest,用于推断非正态分布定量变量或有序分类变量的多个总体分布位置有无差别。,.,（1）原始数据的多个样本比较例12.5某医师检测3种卵巢功能异常患者血清中促黄体素的含量(U/L)资料见表12.5第(1)、(3)、(5)栏。问3种患者血清中促黄体素的含量(U/L)是否有差别？,.,表12.5三种卵巢功能异常患者血清中促黄体素的含量(U/L),这是定量资料多组独立样本的比较，并且该数据属于一端无确切值的资料，不能采用方差分析，现用Kruskal-WallisH秩和检验进行分析。,.,检验步骤：（1）建立检验假设，确定检验水准。H0：3种卵巢功能异常患者血清中促黄体素的含量总体分布位置相同H1：3种卵巢功能异常患者血清中促黄体素的含量总体分布位置不全相同=0.05,.,（2）计算检验统计量H将三组数据由小到大统一编秩，遇相同数值在同一组内，可顺次编秩；当相同数值出现在不同组时，则必须求平均秩次。,表12.5三种卵巢功能异常患者血清中促黄体素的含量(U/L),.,（2）计算检验统计量H,（12.5）,本例,.,当相持较多，按式（12.5）H值偏小，可用式（12.6）进行校正得Hc值。,（12.6）,本例相同秩次较少，无需校正。,其中,tj为第j次相持时相同秩次的个数。,.,（3）确定P值，并作出统计推断。,.,Kruskal-WallisH检验的基本思想：,.,（2）等级资料的多样本比较例12.6：某医院用三种复方小叶枇杷治疗老年慢性支气管炎，数据见表12.6第(1)(4)栏，试比较三种方剂的疗效有无差异。,.,.,（1）建立检验假设，确定检验水准。H0：3种方剂疗效的总体分布位置相同H1：3种方剂疗效的总体分布位置不全相同=0.05,.,（2）计算检验统计量H编秩：与频数表两样本比较类似。先计算各等级的合计，再确定秩次范围，及平均秩次，见表第(5)、(6)、(7)栏。求各组秩和：与等级资料两样本比较类似，计算各组各等级的频数与平均秩次的乘积之和。,.,(3)计算检验统计量值H,.,由于相持出现较多，按式(12.6)计算Hc值。,.,（3）确定P值，并作出统计推断。因组数k=3，=k1=2，近似服从=2的分布，查界值表，得P0.005。按=0.05水准拒绝H0，接受H1，差别有统计学意义，故可以认为3种复方小叶枇杷方剂治疗老年慢性支气管炎的疗效有差别。若要具体回答3种复方小叶枇杷方剂治疗老年慢性支气管炎的疗效每两种之间是否有差别，还需进一步做两两比较。,.,多个样本的两两比较的秩和检验,.,（三）成组设计多个样本资料的两两比较,.,扩展的t检验,统计量t值,.,.,.,(2)列出两两比较计算表，求得t值。见表12.7。,.,.,3.确定值，作出统计推断,.,第四节随机区组设计的秩和检验,多个相关样本比较的FriedmanM检验,.,.,.,.,.,总结：秩和检验的正确应用,主要对等级资料或者分布不明确资料进行分析；秩和检验可用于任意分布(distributionfree)的资料；T检验与H检验的关系T检验H检验t检验F检验,n2,n=2,主要对等级资料或者分布不明确资料进行分析；秩和检验可用于任意分布(distributionfree)的资料；T检验与H检验的关系T检验H检验t检验F检验,.,秩和检验用于定量资料,计量资料中：极度偏态资料，或个别数值偏离过大各组离散度相差悬殊资料中含有不确定值大于5年0.0011:1024以上兼有等级和定量性质的资料,.,心肌坏死面积例,用二种食物配方饲养大白鼠，观察心肌坏死面积。分析二组间的差异。A组：(n=29,mean3.61)0,0,0,0,0,0,0,0,0,0,0.2,0.3,0.4,0.4,0.6,1,1.6,2.2,2.6,3.3,4.3,5.1