秩和检验综合概述(ppt 42页).ppt

1 本资料来源 2 非参数检验 山东大学公共卫生学院 3 参数统计 parametricstatistics 以样本来自已知分布 如正态分布 总体为假设基础 对未知总体参数进行估计或检验的方法 如t检验 u检验和方差分析等 非参数统计 nonparametricstatistics 总体分布不易确定 或分布呈明显偏态 方差不齐又无适当的变量转换方法以满足参数统计条件 非参数统计比较分布而不是比较参数 4 非参数统计的优点与缺点 优点 适用范围广 不受总体分布的限制 对数据的要求不严 如某些指标难以准确测定 只能以严重程度 优劣等级 先后次序等表示的资料也可应用 方法简便 易于理解和掌握 缺点 如果对符合参数检验的资料应用非参数检验 因不能充分利用资料提供的信息 会使检验效能低于参数检验 若要使检验效能相同往往需要更大的样本含量 5 主要内容 配对设计的符号秩和检验成组设计两样本比较的秩和检验成组设计多样本比较的秩和检验随机区组设计资料的秩和检验 6 第一节配对设计的符号秩和检验 Wilcoxon符号秩和检验 Wilcoxonsignedranktest 是推断其差值是否来自中位数为零的总体的方法 可用于 配对设计的两样本的比较单一样本与总体中位数的比较 7 一 配对设计的两样本的比较 例7 1某医院用中草药 抗苯一号 治疗9例苯中毒患者 其治疗前后的白细胞总数如表7 1 问该药是否对患者的白细胞总数有影响 8 查表法 n 25 1 建立检验假设 确定检验水准H0 Md 0H1 Md 0 0 052 计算检验统计量T值 1 求各对的差值 2 编秩按差值的绝对值由小到大编秩 依差值的正负给秩次冠以正负号 编秩时 若差值为0 舍去不计 若差值的绝对值相等 则取其平均秩次 配对设计两样本比较 9 3 求秩和并确定统计量T分别求出正 负秩次之和 正秩和以T 表示 负秩和以T 表示 可任取正秩和或负秩和为统计量T T T n n 1 2 n为不等于0的对子数 3 确定P值 做出推断结论当n 25时 以T值查附表10T界值表 配对设计用 若检验统计量T值在T界值范围内 则P值大于相应的概率水平 若T值在T界值范围外或等于界值 则P值小于或等于相应的概率水平 配对设计两样本比较 10 本例n 9 T 15 5 查附表10T界值表 得双侧P 0 05 故按 0 05检验水准 不拒绝H0 尚不能认为患者治疗前后的白细胞总数的差别有统计学意义 即尚不能认为中草药 抗苯一号 对苯中毒患者的白细胞总数有影响 当n 5时 应用秩和检验不能得出双侧有统计学意义的概率 故n必须大于5 配对设计两样本比较 11 正态近似法 若n 25 超出附表的范围 用正态近似法作u检验 若相同秩次较多时 不包括差值为0者 计算校正uc 其中tj为第j j 1 2 个相同秩次的个数 配对设计两样本比较 12 基本思想 假定从一总体中 随机抽取n 5的样本 按上述步骤 可以求出T T 当重复抽取所有可能组合的样本 秩和T 的分布见下表 13 14 T分布以均数为中心 均数处频数最多 左右对称 向两侧逐渐减少 当H0成立时 从总体随机抽取n 5的一个样本 所得T值在均数附近的概率最大 而T值远离均数的概率较小 随着n增大 T分布逐渐逼近正态分布 其均数为n n 1 4 方差为n n 1 2n 1 24 当n 25时 T分布较好的近似正态分布 配对设计符号秩和检验 n必须大于5 T分布 15 二 单一样本与总体中位数的比较 例7 2已知某地正常人尿氟含量的中位数为2 15mmol L 今在该地某厂随机抽取12名工人 测得尿氟含量如表 问该厂工人的尿氟含量是否高于当地正常人 16 1 建立检验假设 确定检验水准H0 M 2 15H1 M 2 15单侧 0 052 计算检验统计量 1 求差值差值为各观察值与已知总体中位数之差 2 编秩次方法同配对设计 3 求秩和并确定统计量本例T 62 5 T 3 5 取T 3 5 3 确定P值 做出推断结论本例n 11 T 3 5 查T界值表 得P 0 005 按 0 05检验水准 拒绝H0 接受H1 可认为 单一样本与总体中位数比较 17 第二节成组设计两样本比较的秩和检验 成组设计两样本比较的秩和检验 Wilcoxon两样本比较法 研究目的是推断两样本分别代表的总体分布是否相同 原始数据的两样本比较频数表资料 或等级资料 的两样本比较 18 一 原始数据的两样本比较 例7 3某医师为研究血铁蛋白与肺炎的关系 随机抽查了10名肺炎患者和16名正常志愿者 测得血铁蛋白 g L 见表7 3 问肺炎患者与正常人血铁蛋白含量有无差别 1 建立检验假设 确定检验水准H0 肺炎患者与正常人的血清铁蛋白总体分布相同H1 肺炎患者与正常人的血清铁蛋白总体分布不同 0 05 成组设计两样本比较的秩和检验 19 20 2 计算统计量T值 1 编秩次将2组数据由小到大统一编秩 遇有相同数据且不同组时取平均秩次 2 求秩和并确定统计量T2组秩次分别相加 其对应的秩和分别为183 5和167 5 若2组例数相等 则任取一组的秩和为统计量 若2组例数不等 以样本例数较小者为n1 其对应的秩和为统计量T 3 确定P值 做出推断结论以n1 n2 n1 和检验统计量T值查附表11T界值表确定P值 若T值在界值范围内 则P值大于相应的概率 若T值在界值范围外或等于界值 则P值小于或等于相应概率 本例 成组设计两样本比较的秩和检验原始数据的两样本比较 21 若n1或n2 n1超出了附表11T界值表的范围 可用正态近似法 式中N n1 n2 0 5为连续性校正数 若相同秩次较多 如超过25 时应进行校正 其中tj为第j个相同秩次的个数 成组设计两样本比较的秩和检验原始数据的两样本比较 正态近似法 22 二 频数表资料 或等级资料 的两样本比较 例7 4用某药治疗不同病型的老年慢性支气管炎病人 疗效见表7 4 试比较该药对2种病型的疗效 23 1 建立检验假设 确定检验水准H0 两种病型的病人疗效总体分布相同H1 两种病型的病人疗效总体分布不同 0 052 计算检验统计量 1 编秩本例为等级资料 先计算各等级的合计人数 再确定各等级的秩次范围 并计算其平均秩次 2 求秩和并计算检验统计量以各等级的平均秩次分别乘以两组相应等级的例数 再分别求秩和得到T1 T2 成组设计两样本比较的秩和检验频数表或等级资料的两样本比较 24 25 本例T1 8780 5 T2 12955 5 n1 82 n2 126 故检验统计量T 8780 5 因n1 82 需用u检验 又因该等级资料的相同秩次过多 需进行校正 3 确定P值 做出推断结论uc 0 54130 05 按 0 05检验水准 不拒绝H0 尚不能认为该药对两种病型支气管病人的疗效分布不同 成组设计两样本比较的秩和检验频数表或等级资料的两样本比较 26 第三节成组设计多样本比较的秩和检验 Kruskal Wallis检验 该法是由Kruskal和Wallis在Wilcoxon秩和检验的基础上扩展的方法 亦称为K W检验或H检验 研究目的是推断各样本分别代表的总体分布有无差别 原始数据的比较频数表资料的比较 27 一 原始数据的比较 例7 5某院外科用3种手术方法治疗肝癌患者15例 每组5例 每例术后生存月数如表7 5 试问3种不同手术方法治疗肝癌的效果有无不同 28 1 建立检验假设 确定检验水准H0 3个总体的分布位置相同H1 3个总体的分布位置不同或不全相同 0 052 计算统计量H值 1 编秩将各组数据统一由小到大编秩 如遇有相等数值且不在同组者取平均秩次 2 求秩和并计算检验统计量分别将各组秩次相加 得出各组的秩和Ti 计算H值 Ti为各组的秩和 ni为各组例数 N ni 成组设计多样本比较的秩和检验原始数据的比较 29 3 确定P值 做出推断结论以N 15 n1 n2 n3 5查附表12 1H界值表 得P 0 05 按 0 05检验水准 拒绝H0 接受H1 可认为3种手术方法术后生存月数不同 当组数k 2 每组例数ni 10 或k 3 ni 5时 H分布近似服从于 k 1的x2分布 查x2界值表 确定P值 当相同秩次较多时 如超过25 应计算校正值HC HC H C 式中C 1 tj3 tj N3 N 成组设计多样本比较的秩和检验原始数据的比较 30 二 频数表资料的比较 例7 6某医院在研究胎盘过早剥离者的出血情况时 将妊娠时间分为早 中 晚3个阶段 用来分析不同妊娠阶段时胎盘过早剥离的失血量 资料见表7 6 问妊娠妇女在不同阶段时胎盘过早剥离的出血量有无差别 31 1 建立检验假设 确定检验水准H0 3个妊娠阶段胎盘过早剥离时失血量总体分布相同 H1 3个阶段总体的分布位置不同或不全相同 0 052 计算统计量H 1 编秩同 频数表资料的两样本比较 2 求秩和并计算检验统计量求秩和的方法同 频数表资料的两样本比较 计算检验统计量H值 成组设计多样本比较的秩和检验频数表资料的比较 32 因相同秩次较多 须进行校正 3 确定P值 做出推断结论因本例组数k 3 k 1 2 近似服从 2的x2分布 查附表9x2界值表 得P 0 05 按 0 05检验水准 不拒绝H0 尚不能认为妊娠妇女在不同阶段时胎盘过早剥离的出血量的差别有统计学意义 成组设计多样本比较的秩和检验频数表资料的比较 33 三 多个样本间两两比较的秩和检验 当经多样本的H检验拒绝H0时 只能得出各组的总体分布位置不同或不全相同的结论 常需进一步作组间的多重比较 以推断哪两个总体分布不同 可用扩展了的t检验 各组例数相等或不相等时均适用 式中 和nA nB分别为任2个对比组 A B 的平均秩和与样本含量 k为处理组数 N为各处理组的总例数 H为秩和检验中算得的统计量H值或Hc值 成组设计多样本比较的秩和检验 34 例7 7对例7 5资料作两两比较 成组设计多样本比较的秩和检验多个样本间两两比较 35 1 建立检验假设 确定检验水准H0 3个处理组中任两个总体的分布相同H1 任2个总体的分布不同 0 052 计算统计量t值如1组与2组比较3 确定P值 做出推断结论以 N k和各对比组的t值 查t界值表得P值 成组设计多样本比较的秩和检验多个样本间两两比较 36 37 第四节随机区组设计资料的秩和检验 Friedman检验 该法由M Friedman在符号检验的基础上提出来的 也称M检验 用于配伍组设计资料的比较 原始数据的比较多个样本间两两比较的秩和检验 38 一 原始数据的比较 例7 8不同种系雌性大白鼠注射不同剂量雌激素后子宫重量 g 如下表 试比较4组雌性大白鼠子宫重量的差别有无统计学意义 39 1 建立检验假设 确定检验水准H0 3种剂量雌激素注射后子宫重量总体分布相同H1 3个总体分布不同或不全相同 0 052 计算统计量M值 1 编秩先将各配伍组内数据由小到大编秩 相同数值者取平均秩次 2 求秩和并计算检验统计量计算各处理组的秩和Ti k为处理组数本例 随机区组设计资料的秩和检验原始数据的比较 40 3 确定P值 做出推断结论以处理组数k和配伍组数b查附表12 2M界值表 得出P值 本例 k 3 b 4 查M界值表 M0 05 3 4 26 M 180 05 按