医学统计学精品教学(汕头大学)第十章-非参数检验

2020 4 6 1 第十章非参数秩和检验 吴库生汕头大学医学院预防医学教研室 2020 4 6 秩和检验 2 参数统计与非参数统计 1 参数统计 parametricstatistics 样本所来自的总体分布具有某个已知的函数形式 如正态分布 而其中有的参数是未知的 统计分析的目的是对这些未知的参数进行估计或检验 2020 4 6 秩和检验 3 但实际上有些医学研究资料并不符合上述条件 即使用变量变换的方法也难以达到应用参数方法的要求 这就需要一种不依赖于总体分布的具体形式的统计方法来解决此类问题 由于这类方法不受总体参数的影响 故称非参数检验法 它检验的是分布 而非参数 2020 4 6 秩和检验 4 2 非参数统计 non parametricstatistics 不考虑研究对象总体分布的具体形式 也不对总体参数进行统计推断 而是通过检验样本所代表的总体分布形式是否一致来得出统计结论 由于这类方法不受总体参数的限制 故称非参数统计法 有时也称为任意分布检验 distribution freestatistics 2020 4 6 秩和检验 5 3 非参数统计的优点 适用范围广 不论样本所代表的总体分布形式如何 甚至是未知的 都可适用 资料可用 等级 符号 表示 收集方便 多数非参数检验方法简便 易于理解和掌握 2020 4 6 秩和检验 6 4 非参数统计的缺点 对符合用参数检验的资料 如用非参数统计 会丢失信息 导致检验效率下降 犯第II类错误的可能性比参数检验大 即 值 非参数检验的有些问题的计算仍显繁冗 2020 4 6 秩和检验 7 何时使用非参数检验 等级资料 半定量资料 当指标只能用严重程度 轻 中 重 优劣等级 好 中 差 治疗转归 治愈 好转 无效 等形式表达时 原始资料并非定量的 这些数据比 定量 粗 而比一般的 定性 细 是相继的若干阶梯 但毗邻的阶梯之间并非等距离 亦不能度量 偏态分布资料 2020 4 6 秩和检验 8 何时使用非参数检验 个别数据偏离过大资料这里指的是随机的偏离而不属于 过失误差 所致 还有一种情况就是数据的某一端无确定数值 如 50mg 等等 只给一个上限或下限 而没有具体数值 各组离散程度相差悬殊 方差不齐 的资料如经变量变换达到方差齐性 则可用参数检验 否则的话就要用非参数法处理 2020 4 6 秩和检验 9 何时使用非参数检验 分布形态不明的资料小样本 但又不趋向于正态分布 不稳定 例数少 分布不明确的资料 符合参数检验条件 首选参数检验不符合条件 非参数检验 2020 4 6 秩和检验 10 秩和检验 秩次 rank 就是将观察值按顺序由小到大排列 排名次 并用序号 如1 2 3 代替变量值本身 秩次即通常意义上的序号 秩和 即秩次的和 秩转换 将数值变量转换成秩次 再计算检验统计量的过程 秩和检验 就是通过秩次的排列求秩和进行假设检验的方法 2020 4 6 秩和检验 11 Contents 第一节配对资料的符号秩和检验 Wilcoxon配对法 第二节两独立样本比较的秩和检验 Wilcoxon两样本法 第三节完全随机设计多个样本比较的秩和检验 Kruskal allis检验 2020 4 6 秩和检验 12 第一节两配对样本差值的符号秩和检验 Wilcoxonsignedranktest 符号秩和检验由Wilcoxon于1945年提出 作为配对t检验的替代方法 检验配对资料的差值是否来自中位数为0的总体 在数据满足配对t检验的要求时 符号秩和检验的功效是配对t检验效能的95 左右 2020 4 6 秩和检验 13 符号秩和检验的基本思想 在H0成立的条件下 两配对处理效应相同 两配对样本的差数的正负是随机的 则正差数的秩和与负差数的秩和应该相差不会太大 即总体中位数为0 反之 若两秩和相差太悬殊 则认为H0成立的可能性受到怀疑 从而拒绝H0 接受H1 2020 4 6 秩和检验 14 例1 分别用氰化高铁法与光电比色法测得10个血样的血红蛋白含量 问两法测得结果有无差别 2020 4 6 医学统计学 秩和检验 15 表1两种方法测定血清血红蛋白含量血样号氰化高铁法光电比色法115 012 5212 812 9313 513 0411 312 8512 315 4612 112 9713 713 0812 914 1912 912 11013 012 9 配对差值经检验不符合正态分布 不满足t检验条件 2020 4 6 秩和检验 16 符号秩和检验的基本步骤 建立检验假设 确定检验水准H0 Md 0 两处理效应相同 H1 Md 0 两处理效应不相同 0 05 2020 4 6 秩和检验 17 2 求两组数据的差数d 对差数的绝对值编秩号 差数为0的数据忽略不计 对子数减1 余下的n i个差数按绝对值由小到大排秩号 但排好后的秩号要保持原差数的正负号 差数绝对值相等时 则取其平均秩次 2020 4 6 医学统计学 秩和检验 18 表1两种方法测定血清血红蛋白含量血样号氰化高光电比差值负秩正秩铁法色法 1 2 3 4 2 3 5 6 115 012 52 59212 812 9 0 1 1 5313 513 00 53411 312 8 1 5 8512 315 4 3 1 10612 112 9 0 8 5 5713 713 00 74812 914 1 1 2 7912 912 10 85 51013 012 90 11 5T 32T 23验算 T T n n 1 2 55 统计量T 对差数的绝对值编秩 2020 4 6 秩和检验 19 3 分别求正 负差数的秩和 以绝对值较小者为T验算 T T n n 1 24 确定P值和作出推断结论 小样本时 n 50 n为d 0的对子数 通过查p226附表8 T界值表 得到P值 原则是 若统计量T处于某T界值范围之内 其P值大于或等于相应概率 反之 则小于相应概率 当T值恰好等于附表中的界值时 其确切概率常小于相应的概率水平 2020 4 6 秩和检验 20 本例中 T 23对子数n 10 查p226附表8得23 1045 故P 0 10下结论 在 0 05水准上不拒绝H0 故不能认为两法测得血红蛋白含量有差别 2020 4 6 秩和检验 21 大样本时 对子数n 50 通过公式进行u转化 采用正态近似检验 式10 1 式中0 5是连续性校正数 因为T值是不连续的 而u分布是连续的 秩和标准差 平均秩和 2020 4 6 秩和检验 22 当存在同一秩较多时 按 式10 1 计算的u值要比实际的u值小 应采用 公式10 2 对u加以校正 其中 tj是相同秩号的个数 校正系数 2020 4 6 秩和检验 23 符号秩和检验的SPSS实现SPSS StatisticalPackageforSocialScience 社会科学统计软件包 一直是国际上最流行且最具权威性的统计分析软件之一 SPSS最显著的特点是菜单和对话框操作方式 无需编程 且图形功能强大 易于掌握和操作 2020 4 6 秩和检验 24 2020 4 6 秩和检验 25 符号秩和检验的SPSS实现 1定义变量 2020 4 6 秩和检验 26 符号秩和检验的SPSS实现 2输入数据 2020 4 6 秩和检验 27 2020 4 6 秩和检验 28 符号秩和检验的SPSS实现 2020 4 6 秩和检验 29 符号秩和检验的SPSS实现 2020 4 6 秩和检验 30 WilcoxonSignedRanksTest 结果 2020 4 6 医学统计学 秩和检验 31 第二节两独立样本秩和检验 Wilcoxon两样本比较法 观测值XA71422364048639841 00 29 81B356101718203914 75 11 73方差齐性检验 P 0 05两样本方差不齐 不能应用t检验 2020 4 6 医学统计学 秩和检验 32 采用Wilcoxon两独立样本秩和检验 基本思想 假定原假设H0成立 样本来自两个相同总体 合并两样本 将观察值由小到大排列 编写秩次 在H0成立的条件下 T值应接近平均秩和n0 1 N 2 其中N n1 n2 n0 min n1 n2 如果T值严重偏离平均秩和 则提示H0可能是不正确的 1 2 3 16 2 68n0 1 N 2 8 1 16 68 2020 4 6 医学统计学 秩和检验 33 采用Wilcoxon秩和检验 一 建立检验假设 确定检验水准H0 样本来自两个相同总体 样本的每个观察值来自两总体的概率均为0 5 H1 样本来自两个不同总体 样本的每个观察值来自两总体的概率不等 0 05 2020 4 6 秩和检验 34 二 编秩将两组数据由小到大统一编秩 如遇相同原始数据 若在同一组 仍按数据出现顺序编秩 若在不同组 则取其平均秩次 本例编秩结果见表2 2020 4 6 医学统计学 秩和检验 35 表2两种疗法的退热时间新疗法旧疗法退热时间秩号退热时间秩号25136530240932344113544813 53765015397 55616397 5591742106018461264194813 51952024021n1 10T1 66 5n2 11T2 164 5 两样本混合后统一编秩求得秩号 分别求秩和 2020 4 6 医学统计学 秩和检验 36 三 分组求秩和并确定检验统计量当样本例数不等时 取较小例数组的秩和为统计量T 如n1 n2 则取任一组秩和为统计量T 本例中 n1 10 n2 11 T min T1 T2 T1 66 5 2020 4 6 医学统计学 秩和检验 37 四 查p227附表9 n1 n2均较小 确定检验界值的区间 得到P值并做出推断结论原则 找到n1和n2 n1对应的四行界值 如果T位于检验界值区间内 P 不拒绝H0 如果T值恰好等于界值 P 如果T位于检验界值区间外 P 拒绝H0 接受H1 2020 4 6 医学统计学 秩和检验 38 在本例中 T 66 5 取 0 05 根据附表9双侧检验界值0 01区间为 73 147 T位于区间外 所以P 0 01 P 0 01 因此在 0 05的水准上 拒绝H0 接受H1 认为新疗法与旧疗法的退热时间不同 2020 4 6 医学统计学 秩和检验 39 大样本情况下统计量的计算 当样本例数较大 超出附表9的范围时 可以对T值进行变换 记为u 式10 3 可利用u检验法对假设做出检验 平均秩和 2020 4 6 医学统计学 秩和检验 40 在相同秩次较多时需要按 式10 4 计算校正的检验统计量 式10 4 式中tj为第j个相同秩次的个数 在无相同秩次时 C 1 即为10 3式 2020 4 6 秩和检验 41 Wilcoxon秩和检验的基本步骤总结 建立检验假设 确定检验水准 混合编秩 分组求秩和 计算检验统计量T值 查附表9 得检验界值区间 确定P值 作出结论 2020 4 6 医学统计学 秩和检验 42 混合编秩的注意事项 如果有相同的数值 出现相同秩次 ties 时可按以下方法处理 相同数据在同一组 仍按顺序编秩 相同数据在不同组 取其平均秩次 样本较小时 如果相同秩次较多 检验结果会存在偏性 应提高测量精度 避免出现较多的同秩 大样本时 要对公式 10 3 进行适当校正 以减少相同秩次带来的偏秩 2020 4 6 医学统计学 秩和检验 43 P173例11 3两独立样本秩和检验SPSS实现 1定义变量 定义组别 2020 4 6 医学统计学 秩和检验 44 两独立样本秩和检验SPSS实现 2数据输入 2020 4 6 医学统计学 秩和检验 45 两独立样本秩和检验SPSS实现 2020 4 6 医学统计学 秩和检验 46 两独立样本秩和检验SPSS实现 检验变量 分组变量 2020 4 6 医学统计学 秩和检验 47 定义分组 2020 4 6 医学统计学 秩和检验 48 结果 2020 4 6 医学统计学 秩和检验 49 二两样本等级资料的比较 例3就表3中 1 2 3 列的资料 比较巴非林与阿司匹林两种药物解热镇痛的疗效有无差别 2020 4 6 医学统计学 秩和检验 50 表3巴非林与阿司匹林疗效比较疗效频数巴非林阿司匹林 1 2 3 痊愈4512显效2717有效2217无效64合计10050单向有序列联表资料 用检验比较疗效将损失信息 应采用秩和检验 单向有序变量 2020 4 6 医学统计学 秩和检验 51 建立检验假设 确定检验水准 H0 两组疗效相同H1 两组疗效不相同 0 05 2020 4 6 医学统计学 秩和检验 52 混合编秩 求各组秩和表3巴非林与阿司匹林疗效比较计算表疗效频数合计累计秩次平均秩次和巴非林阿司匹林频数范围秩次巴非林阿司匹林 1 2 3 4 5 6 7 8 2 7 9 3 7 痊愈451257571 57291305348显效27174410158 10179 52146 51351 5有效221739140102 14012126622057无效6410150141 150145 5873582合计100501506986 54338 5 T T1T2 有序变量 疗效合并 混合编秩 求和 1 57 2 2020 4 6 医学统计学 秩和检验 53 计算检验统计量 确定P值由于本例样本量较大 且存在同秩 故采用校正后的统计量uc先计算检验统计量T n1 100 n2 50 T 例数较少组对应的秩和 T2 4338 5 T 2020 4 6 医学统计学 秩和检验 54 根据公式 10 3 计算u根据C 根据t界值表 P 0 001 2020 4 6 医学统计学 秩和检验 55 根据P值 下结论按 0 05的检验水准 拒绝H0 接受H1 故可以认为两种药物的疗效不同 注意 Wilcoxon秩和检验与一般文献中使用的Mann Whitney的U检验是独立提出的 但是两者的检验结果完全等价 前者用T统计量 后者用U统计量 可任意选取 2020 4 6 医学统计学 秩和检验 56 Wilcoxon两样本比较法基本思想 假定原假设H0成立 合并两样本 将观察值由小到大排列 编写秩次 在H0成立的条件下 T值应接近平均秩和n0 1 N 2 其中N n1 n2 n0 min n1 n2 如果T值严重偏离平均秩和 则提示H0可能是不正确的 2020 4 6 医学统计学 秩和检验 57 两样本等级资料秩和检验的SPSS实现 定义变量 2020 4 6 医学统计学 秩和检验 58 2输入数据 2020 4 6 医学统计学 秩和检验 59 3频数加权 1 2 3 2020 4 6 医学统计学 秩和检验 60 4开始分析 2020 4 6 医学统计学 秩和检验 61 检验变量 分组变量 定义分组 2020 4 6 医学统计学 秩和检验 62 结果 2020 4 6 秩和检验 63 第三节完全随机设计多组差别的秩和检验 Kruskal Wallis检验 Kruskal Wallis检验的基本思想 如果各组均来自同一总体 则各组的平均秩和近似相等 2020 4 6 秩和检验 64 Kruskal Wallis检验的基本步骤 建立检验假设 确定检验水准 混合编秩 分组求秩和Ri i 1 2 k k是处理的个数 计算检验统计量H 确定P值 作出推断结论 2020 4 6 秩和检验 65 例3试检验表4中三组人的血浆总皮质醇含量有无差别 2020 4 6 医学统计学 秩和检验 66 表3三组人血浆总皮质醇测定值 g L 正常人单纯性肥胖皮质醇增多症测定值秩次测定值秩次测定值秩次0 1110 1722 79200 5240 3332 81210 6160 5552 92220 6980 6673 59230 7790 8610 53 86250 8610 51 13144 08261 02121 38164 30271 08131 63174 30281 27152 04195 96291 92183 75246 6230秩和96 5117 5251例数101010 2020 4 6 秩和检验 67 建立检验假设 确定检验水准 H0 三组人的血浆总皮质醇含量总体分布相同H1 三组人的血浆总皮质醇含量总体分布不全相同 0 05 2020 4 6 秩和检验 68 混合编秩 分组求秩和Ri i 1 2 k k是处理的个数 将各组数据混合 由小到大编秩 遇有原始数据相同时 若在同一组 则按顺序编秩 若分属不同组 取其平均秩次 再将各组秩号分别相加 求出各组的秩和Ri 2020 4 6 医学统计学 秩和检验 69 表3三组人血浆总皮质醇测定值 g L 正常人单纯性肥胖皮质醇增多症测定值秩次测定值秩次测定值秩次0 1110 1722 79200 5240 3332 81210 6160 5552 92220 6980 6673 59230 7790 8610 53 86250 8610 51 13144 08261 02121 38164 30271 08131 63174 30281 27152 04195 96291 92183 75246 6230秩和Ri96 5117 5251例数ni101010 统一编秩 分组求秩和 2020 4 6 秩和检验 70 计算检验统计量H检验统计量H按式 10 5 计算 式10 5 式中 Ri为各组的秩和 ni为各组样本的含量 N为总的样本量 2020 4 6 秩和检验 71 当各组相同秩次较多时 计算所得H值偏