第十章 基于秩次的非参数检验.ppt

第十章基于秩次的非参数检验 nonparametrictestoftransformationrank 要求 1 理解并能区别参数性检验与非参数性检验2 掌握配对设计和成组设计的秩和检验3 了解多个样本两两比较的秩和检验4 掌握秩和检验的SPSS统计程序操作 参数检验与非参数检验 参数检验 parametrictest 在已知总体分布类型的基础上 对总体参数进行估计和检验 非参数检验 nonparametrictest 总体分布类型未知 或非精确测量资料 或参数检验的条件得不到满足时 对资料的分布规律及特征进行统计推断 秩和检验 卡方检验 游程检验 序贯检验和Ridit分析都是非参数检验 参数检验与非参数检验的优缺点 非参数检验的主要优点是对资料的要求不象参数检验那样严 适合处理非正态 方差不齐 以及分布类型不明确等资料 其次 有些问题本身还没有适当的参数检验方法可用 而非参数检验法则能予以处理 非参数检验的不足之处是 由于它对原始数据所包含的信息利用不充分 检验效率就低了 尤其是对适宜用参数检验法处理的资料 若用非参数法分析 一般犯第二类错误的概率比参数检验大 故适宜用参数检验法的资料 应首选参数检验法 非参数检验适用对象 1 等级资料2 偏态资料3 分布类型未知4 多组资料的变异度大 方差不齐 且不易变换达到齐性5 初步分析 第一节单样本和配对设计资料的符号秩和检验 一 单样本资料的符号秩和检验 Wilcoxon法 例10 1已知某地正常人尿铅含量的中位数为2 50 mol L 今对该地随机抽取16名工人 测得尿铅含量见表10 1 问该厂工人的尿铅含量是否高于当地正常人 不是均数 表10 118名工人测得尿铅含量与2 50的比较 mol L H0 该厂工人尿铅含量与正常人基本一致H1 该厂工人尿铅含量高于正常人单侧 0 05 当差值等于0时 去掉该点 T 108T 28 1 查表法 5 n 50时 例10 1的分析 本例资料正态性检验W 0 809 P 0 0036 即尿铅含量不服从正态分布 不宜用t检验 本例n 16 T 28 单侧 当 0 05时查表9 结果 T 28 35 101 所以P 0 05 拒绝H0 可认为该厂工人的尿铅含量高于当地正常人 2 正态近似法 n 50时 式中0 5是连续性校正值 因为T值是不连续的 而Z分布是连续的 这种校正影响甚微 常可省去 当相同 差值 指绝对值 较多时 不包括差值为0者 改用校正公式 式中tj为第j个相同差值的个数 假定在差值中有2个1 5 3个6 5个13 则t1 2 t2 3 t3 5 tj3 tj 23 2 33 3 53 5 150 二 配对设计资料的符号秩和检验 主要用于配对资料的检验 H0 差值总体中位数Md 0 其方法步骤如下 1 求差值 求各对数据 xi yi 的差值di xi yi2 编秩 依差值的绝对数大小从小到大编秩 再根据差值的正 负给秩冠以正负号 编秩时 若差值为0 弃去不计 若遇有差值的绝对值相等 取其平均秩次 3 求秩和并确定检验统计量T 分别求正 负秩次之和 记作T 和T 并以绝对值较小者为T值4 确定P值和作出结论 当配对数n 50时 查配对设计T界值表 T在界内 P T在界外 P 当配对数n 50时 可用正态近似法的Z检验 例10 2 表10 2电极法和分光光度法测定氟离子浓度 例10 2对11份工业污水测定氟离子浓度 mg L 每份水样同时采用电极法和分光光度法测定 资料见表 2 3 推断不同方法测定结果有无差别 1 查表法 5 n 50时 例10 2分析H0 两法测定结果一致 0 05秩和为 T 43 5 T 11 5由于n 10 去掉0差值 查表9本例n 10 T 11 5 双侧 当 0 05时 T 11 5 8 49 所以P 0 05 不拒绝H0 认为两法测定结果没有显著性差异 2 正态近似法Z检验统计量 当配对数n 50时 式中0 5是连续性校正值 因为T值是不连续的 而Z分布是连续的 这种校正影响甚微 常可省去 当相同 差值 指绝对值 较多时 不包括差值为0者 改用校正公式 式中tj为第j个相同差值的个数 假定在差值中有3个4 5个2 2个7 则t1 3 t2 5 t3 2 tj3 tj 33 3 53 5 23 2 150 Wilcoxon配对法基本思想 本法的基本思想是假定两种处理结果的效应相同 即差值之总体分布是对称的 H0 Md 0 总体中位数为0 同理 假定某种处理无作用 则每一受试对象处理前后所得结果之差值的总体中位数亦为0 如果假设成立 则样本的正 负秩和应比较接近 若正 负秩和相差悬殊 则假设成立的可能性也小 第二节完全随机设计两个独立样本比较的秩和检验 主要用于成组比较资料的检验 其方法步骤如下 1 编秩在原假设 H0 两总体分布相同 下 将两组数据混合由小到大统一编秩 编秩时如遇有原始数据相同的 则取它们的平均秩次 2 求秩和并确定检验统计量当两样本例数不等时 以样本例数小者为n1 其秩和为统计量T 若n1 n2时 可任取一组的秩和为统计量T 3 确定P值和作出结论当n1 10 n2 n1 10时 查成组设计T界值表 附表10 T在界内 P T在界外 P 当n1 10或n2 n1 10时 可用正态近似法的Z检验 一 定量资料两组独立样本的秩和检验 例10 3用两种药物杀灭钉螺 每批用200 300只活体钉螺 用药后清点钉螺的死亡数 并计算死亡率 问两种药物杀死钉螺的效果有无差别 一 定量资料两组独立样本的秩和检验 编秩 1 查表法 当n1 10 n1 n2 10时 例10 3分析H0 两种药物杀死钉螺的效果无差别 0 05统计量T 33 5查表10本例n1 7 n1 n2 0 当 0 05时 双侧 T 33 5 36 69 所以P 0 05 拒绝H0 认为两法测定结果有显著性差异 两样本比较T界值表 2 正态近似法 当n1 10 或n1 n2 10时 式中N n1 n2当相同秩次较多时 25 改用校正公式 式中tj为第j个相同秩次的个数 二 两组有序变量资料的秩和检验 例10 4某地居民夏冬两个季节体内核黄素营养状况比较 二 两组有序变量资料的秩和检验 编秩 第三节多组独立样本比较的秩和检验 1 编秩在原假设 H0 多总体分布相同 下 将各组数据混合由小到大统一编秩 编秩时如遇有原始数据相同的 则取它们的平均秩次 2 求秩和求每组秩次总和 Ri i表示第i组 3 计算统计量H值 Kruskal Wallis 式中ni为样本量 N ni为样本总量 当各样本相同秩次较多 25 时 应用校正统计量Hc 4 确定P值和作出推断结论当组数k 3 每组例数 5时 查H界值表 附表11 得出P值 当最小样本例数大于5 则H近似服从v k 1的 2分布 由 2界值表确定P值 一 定量资料多组独立样本的秩和检验 例10 5某医院用3种不同方法治疗15例胰腺癌患者 每种方法治疗5例 治疗后生存月数见表 问3种方法对胰腺癌患者疗效有无差别 一 定量资料多组独立样本的秩和检验 编秩 求统计量 H界值表 确定P值 解 H0 3种方法治疗后患者生存月数的总体中位数相等H1 3种方法治疗后患者生存月数的总体中位数不等 0 05 本例k 3 ni 5 查表11 H0 05 5 78 所以P 0 05 可认为3种方法治疗后患者生存月数的总体中位数不等 近似法 当超出H界值表时 二 有序变量多组独立样本的秩和检验 例10 6某医院用3种方法治疗478例慢性喉炎 资料见表 问3种方法治疗慢性喉炎的疗效有无差别 二 有序变量多组独立样本的秩和检验 编秩 例10 6分析 三 多个独立样本间