第二讲非参数统计检验资料.docVIP

第二讲 非参数检验 1. 实验目的 1.了解非参数假设检验基本思想； 2.会用SAS软件中的proc npar1way过程进行非参数假设检验和proc freq过程进行列联表的独立性检验。 2. 实验要求 1.会用SAS软件建立数据集，并进行统计分析； 2.掌握proc npar1way过程进行非参数假设检验的基本步骤； 3.掌握proc freq过程进行列联表的独立性检验的基本步骤。 3. 实验基本原理 3.1 符号检验 两种方法的处理效果无显著性差异 令 统计量 表示新方法的处理效果优于对照方法的配对组总数。若新方法的处理效果显著的优于对照方法，则的值应明显偏大。因此，若对给定的置信水平，有 ，则拒绝。 为真时，（1）服从二项分布。拒绝域为： （2）由中心极限定理可知，当的零分布趋于标准正态分布。拒绝域为： 3.2 Wilcoxon秩和检验 （1）单边假设检验 　两种方法的处理效果无显著性差异 as ：新方法优于对照方法。 用于检验的统计量为： 若对给定的置信水平，有 ，则拒绝。且的分布列为： 根据观测结果计算的观测值，计算检验的p值： 然后将值与显著水平作比较，若，则拒绝，否则接受。 （2）双边假设检验 给定的显著水平应该满足： 　　　　 仅由上式还不能唯一确定，当我们对两种方法谁优谁劣不得而知时，通常取 若利用p值进行检验，设,计算概率值 　 由对称性可知，检验的p值为上述两概率中小于1/2的那一个的2倍。例如则。求出p值后，若pa，拒绝,否则接受。 （3）列联表的独立性检验 方法的处理效果无显著性差异 ：表示格子概率，表示三维列联表中事件发生的理论频数。将概率用相应的频率频率去估计。 令 其中： 将样本数据代入统计量进行检验。然后将P与显著水平作比较，若，拒绝，否则接受。 4. 实验相关SAS知识 （1）独立样本的秩检验——proc npar1way过程 proc npar1way过程的基本语句形式为 proc npar1way [options]; class variables;（proc npar1way过程不可缺少的语句） exact;（求出检验的精确p值） var variables; 其中“options”可包含以下选项的部分或全部： ①DATA=数据集名：指定要分析的数据集。 ②ANOVA：对原始数据执行标准的单因素方差分析。 ③WILCOXON：进行wilcoxon型秩和检验。当有两种处理方法时，进行的是wilcoxon秩和检验；当有多种处理方法时，进行Kruskall-Wallis检验。 ④EDF：进行基于样本经验分布函数的非参数检验，包括Smirnov检验。 若省略这些选项，SAS系统将给出所有基于秩以及经验分布函数的非参数检验方法的分析结果。 （2）列联表的独立性检验 proc freq过程的基本语句形式为 proc freq [options]; tables variable1*variable2*……/options; weight variable; 其中“options”可包含以下选项的部分或全部： ①DATA=数据集名：指定要分析的数据集。 ②chisq:要求对生成的每个二维列联表的独立性作检验，并计算依赖于统计量的关联度。 ③cellchi2:要求输出每个格子对总统计量的贡献。 ④expected:在独立性假定下输出各格子的期望频数。 ⑤deviation:要求输出每个格子上的频数与期望频数之差。 ⑥nocol:不输出二维列联表各格子的列百分数。 ⑦norow: 不输出二维列联表各格子的行百分数。 ⑧nofreq:不输出格子频数。 ⑨nopercent:不输出各格子的百分数。 ⑩noprint:不输出列联表，但允许输出各分析结果。 5. 实验举例 5.1 Wilcoxon秩和检验(单边和双边假设检验) 例1 为了解一种新的术后护理方法和原护理方法相比是否可以显著缩短病人手术后的恢复时间,随机的将做完某种手术的18位病人分为两组,每组9人,按不同方法护理,观测到他们的恢复时间(单位：天)如下： 原方法：20,21,24,30,32,36,40,48,54 新方法：19,22,25,26,28,29,34,37,38 在下检验新方法是否显著的缩短了病人手术后的恢复时间。 Wilcoxon秩和单边假设检验SAS程序如下： data a1; input method $ time@@; cards; a 20 a 21 a 24 a 30 a 32 a 36 a 40 a 48 a 54 b 19 b 22 b 25 b 26 b 28 b 29 b 34 b