第九讲 基于秩次的统计方法.ppt

章 基于秩次的统计方法 ;主讲内容;我们常常遇到以下一些资料，如需比较患者和正常人的血铁蛋白、血铅值等各项指标、护理效果评分、医疗质量评估等，这类资料有如下特点： （1）资料的总体分布类型未知； （2）资料分布类型已知，但不符合正态分布； （3）某些变量可能无法精确测量如等级资料。 对这类资料可以采用非参数统计：即不考虑总体分布类型是否已知，不比较总体参数，只比较总体分布的位置是否相同的统计方法。此类资料可以采用非参数方法进行统计分析。;秩和检验概述;一、秩次和秩和;两组各有5个变量值。现在依从小到大的顺序将它们排列起来，并标明秩次，结果如下： A组 2.6 3.2 4.7 5.2 6.4 B组 1.7 2.3 2.6 3.6 3.7 秩次 1 2 3 4 5 6 7 8 9 10;原始值中有两个“2.6”，分属A、B组，它们的秩次应是3和4，然而它们的数值本来是同样大小的，哪组取“3”，哪组取“4”呢？我们计算它们的平均数（3+4）/2=3.5，作为“2.6”的秩次，称为“平均秩次”，这样才公平合理。这样两组所得的秩次及秩和如下： A组 3.5 5 8 9 10 /35.5 B组 1 2 3.5 6 7 /19.5;上面A组和B组中各有五个原始值，按顺序排列：最小值设为1，再按绝对值大小对余下的变量逐个排序，最大值为两组变量个数之和10。依次可得1，2，3.5,3.5,5,6,7,8,9,10。这10 个序号即是秩次。A组秩和就是等于3.5+5+8+9+10=39.5，B组秩和就是等于1+2+3.5+6+7=19.5。从两组的原始变量值也可以初步看出：A组偏大，B组偏小。现在得出的秩和也是A组大于B组，与由变量值所观察到的一致。;二、怎样编秩 ;1.配对资料的编秩规则;2.完全随机设计的两组或多组资料的编秩规则 ;三、怎样选用不同设计的秩和检验方法 ;若是多组独立样本资料，来自非正态总体或方差不齐的定量资料或者是有序多分类变量资料（等级资料）宜用Kruskal-Wallis秩和检验。 各种类型的秩和检验通常有查表法和正态近似法两种计算方式，分别适用于不同样本量的资料。;第一节 两独立样本差别的秩和检验;观察值X;根据数据不难检验两样本方差不齐，因此不能直接应用t检验。这里我们应用两独立样本差别的秩和检验(rank sum test);两独立样本秩和检验基本思想;A 样本;秩和检验的检验统计量T定义：;在H0 成立的条件下，T值应近n0(1+N)/2，其中N=n1+n2,n0=min(n1,n2)。如果T严重偏离n0(1+N)/2 ，提示H0可能不正确。; （1）建立检验假设H0,H1，给定检验水准α。（2）统一编秩号，分组求秩和R1,R2，计算检验统计量T。 （3）查附表，得检验界值区间。 （4）下结论：如果T位于检验界值区间内,P α,不拒绝H0；如果P≤ α，则拒绝H0 ，接受H1。;本例计算得，T=R2=47，取α =0.05，查附表得双侧检验界值区间（49，87），T位于区间外，P0.05,因此在α =0.05的水准上，拒绝H0，接受H1，认为两样本不是来自同一个总体。;;相同秩问题;例1分别用5％咪喹莫特软膏和5-氟脲嘧啶软膏治疗尖锐湿疣的随机双盲临床研究的疗效观察结果见表2第(1)、（2）列，试比较两种药物治疗尖锐湿疣的疗效。;1.建立检验假设，确定检验水准。 H0：用5％咪喹莫特软膏组与用5-氟脲 嘧啶软膏疗效相同 H1：用5％咪喹莫特软膏组与用5-氟脲 嘧啶软膏疗效不同;疗效;;本例，由公式计算得c=0.5836，代入公式，得u=0.7249，除以校正系数c，得到： uc=u/c=0.7249/0.5836≈ 1.2421.96,P0.05;;注意;第二节 配对设计资料的秩检验-Wilcoxon符号秩检验(sign rank test) ;1．H0：两处理效应相同 H1：两处理效应不相同。 设定检验水准，如α =0.05;2．编秩号。成对资料编秩号时较为复杂，要注意三点： （1）剔去差数为0的数据； （2）余下的n个差数按绝对值自小至大排秩号，但排好后秩号要保持原差数的正负号； （3）差数绝对值相等时，要以平均秩号表示。;3．求秩号和，即将正、负秩号分别相加;4. 检验统计量R取较小一个秩号和，根据R值查附表10进行判断，该表左侧为对子数，表身内部是较小秩