非参数检验SPSS.ppt

非参数检验的SPSS实现1一、两独立样本的检验 例13-1调查某厂的铅作业工人7人和非铅作业工人10人的血铅值（g/100g）如下，问两组工人的血铅值有无差别？非铅作业组非铅作业组5 5 6 7 9 12 13 15 18 215 5 6 7 9 12 13 15 18 21铅作业组铅作业组 17 18 20 25 34 43 44 17 18 20 25 34 43 442数据准备定义2个变量：group（非铅作业组为1，铅作业组为2），Pb。按顺序输入数据。3统计分析分析非参数统计两独立样本b选入检验变量列表框，group选入分组变量框定义组别，定义组1为1，组2为2 检验类型：Mann-Whitney U确定4结果解释 结果表明，第1组的平均秩次（Mean Rank）为5.95，第2组的平均秩次为13.36，Mann-WhitneyU=4.5，W=59.5，z=-2.98,精确双侧概率P=0.001，可认为铅作业组工人的血铅值高于非铅作业组。5二、多个独立样本的检验例13-2随机抽样得以下三组人的血桨总皮质醇测定值（g/L），试比较有无差异？正常人正常人单纯性肥胖单纯性肥胖皮质醇增多症皮质醇增多症0.40.40.60.69.89.81.91.91.21.210.210.22.22.22.02.010.610.62.52.52.42.413.013.02.82.83.13.114.014.03.13.14.14.114.814.83.73.75.05.015.615.63.93.95.95.915.615.64.64.67.47.421.621.67.07.013.613.624.024.06数据准备定义分组变量为 group（正常人为1，单纯性肥胖为2，皮质醇增多症为3），总皮质醇测定值为x。按顺序输入数据。7统计分析分析非参数统计个独立样本x选入检验变量框，组别选入分组变量框定义组别，最小值为1，最大值为3 确定8结果解释 结果表明，1至3组的平均秩次（Mean Rank）分别为9.65、11.75、25.10，2 值（即H值）为 18.1219，P=0.0001；可认为三组人的血桨总皮质醇测定值有差异，根据本例情况可看出皮质醇增多症组高于其他两组人。9三、配对资料秩和检验 例13-3研究饲料中缺乏Vit E对大鼠肝中Vit A含量的关系，将大鼠按性别相同、体重相近的原则配成8对，并将每对大鼠随机分为2组（正常饲料组、Vit E缺乏饲料组），一定时间后杀死大鼠，测定肝中Vit A含量，结果如下表，问：饲料中缺乏Vit E对大鼠肝中Vit A含量有无影响？10大鼠对别大鼠对别正常饲料组正常饲料组Vit E Vit E 缺乏饲料组缺乏饲料组1 137.237.225.725.72 220.920.925.125.13 331.431.418.818.84 441.441.433.533.55 539.839.834.034.06 639.339.328.328.37 736.136.126.226.28 831.931.918.318.311数据准备 两个变量:正常饲料;Vit E 缺乏,按顺序输入数据。12统计分析分析非参数统计两个相关样本检验将正常饲料和Vit E 选入检验对列框检验类型:选择前两种(Wilcoxon和符号检验)确定13结果解释首先显示两变量正常饲料和Vit E 缺乏的例数、均数、标准差、最大值和最小值；配对符号秩和检验（Wilcoxon Matched-Pairs Signed-Ranks Test）结果，其平均秩分别为5.00 和1.00，Z=-2.380，双侧P=0.017，可认为两组大鼠肝中Vit A含量有差别，饲料中缺乏Vit E会使大鼠肝中Vit A含量降低；符号检验（Sign Test）的结果，双侧P=0.0703，则认为两组大鼠肝中Vit A含量无差别。在这种情况下，应取配对符号秩和检验（Wilcoxon）结果，因两法比较之下，配对符号秩和检验较为敏感，效率较高。14四、配伍组资料秩和检验例13-4用某药治疗血吸虫病患者，在治疗前和治疗后一周、二周和四周各测定7名患者血清SGPT值的变化，以观察该药对肝功能的影响，结果如下表，问：患者四个阶段的血清SGPT值有无不同 15患者编号患者编号治疗前治疗前治治 疗疗 后后一周一周二周二周四周四周1 1636318818813813854542 290902382382202201441443 35454300300838392924 445451401402132131001005 5545417517515015036366 6727230030016316390907 764642072071851858787用某药治疗血吸虫病患者治疗前后血清SGPT值16数据准备定义4个变量：治疗前为before、治疗后一周为w1、二周为w2、四周为w4，按顺序输入各组SGPT数据。17统计分析分析非参数统计K个相关样本检验before、w1、w2和w4均进入检验对列框检验类型:选择Friedman 确定18结果解释首先显示的是四个变量before、w1、w2、w4的例数、均数、标准差、最大值和最小值。接着显示检验结果：秩均值分别1.29、3.86、3.00和1.86，2=16.714，P=0.001，可认为患者四个阶段的血清SGPT值有差别。19 例13-5 某研究者欲比较消炎痛与消炎痛加皮质激素制剂（简称合剂）治疗肾小球肾病的疗效，将64例肾小球肾病患者随机分为两组，分别用消炎痛与合剂治疗，全程用药后病情分为完全缓解、基本缓解、部分缓解与无效四个等级，数据见下表，试比较两种药物治疗肾小球肾病的疗效有无不同？20表 两种疗效对肾小球肾病的疗效比较疗效疗效消炎痛消炎痛合剂合剂完全缓解完全缓解219基本缓解基本缓解45部分缓解部分缓解69无无 效效154212024/5/22 周三22五、等级资料秩和检验等级资料的SPSS操作要点：1.定义3个变量，即疗效、疗法（或疾病）和频数。2.在进行统计分析前，需对频数加权。3.采用2个和K个独立样本的非参数检验。23例13-6某村发生一种地方病，其住户沿一条河排列，调查时对发病的住户标记为“1”，共9例，对非发病的住户标记为“0”，共17户：问病户的分布排列是呈聚集趋势，还是随机分布？0 1 1 0 0 0 1 0 0 1 0 0 0 0 1 1 0 0 1 0 0 0 0 1 0 124主要功能依时间或其他顺序排列的数据序列中，具有相同的事件或符号的连续部分称为一个游程。调用Runs过程可进行游程检验，即用于检验序列中两类事件发生过程是否随机。六、Runs过程（游程检验）25假定下面是由0和1组成的一个这种变量的样本：0 1 1 0 0 0 1 0 0 1 0 0 0 0 1 1 0 0 1 0 0 0 0 1 0 1 其中相同的0（或相同的1）在一起称为一个游程（单独的0或1也算）。这个数据中有个0组成的游程和个1组成的游程。26数据准备定义住户变量为epi。按住户顺序输入数据，发病的住户为1,非发病的住户为0。27统计分析分析非参数统计游程检验epi选入检验变量列框 选择割点确定281、中位数：中位数作割点，其值小于割点为一类，大于或等于割点的为另一类；2、众数：众数作割点，其值小于割点为一类，大于或等于临界割点的为另一类；3、均数：均数作割点，其值小于割点为一类，大于或等于临界割点的为另一类；4、定制：用户指定割点，其值在临界割点之下的为一类，大于或等于割点的为另一类；本例选定制，在其方框中键入1（根据需要选项，本例是0、1二分变量，故临界割点值用1），再点击OK钮即可。29结果解释检验结果可见本例游程个数为14，检验临界割点值（Test value）=1.00，Z=0.325，双侧 P=0.746。所以认为此地方病的病户沿河分布的情况无聚集性，而是呈随机分布。30七、1-Sample K-S过程 调用此过程可对单样本进行Kolmogorov-Smirnov 检验，它将一个变量的实际频数分布与正态分布（Normal）、均匀分布（Uniform）、泊松分布（Poisson）进行比较。31数据准备以salesSav为例32统计分析分析非参数统计单样本Kolmogorov-Smirnov检验revenue选入检验变量列框 检验分布：正态分布确定33结果解释K-S正态性检验的结果显示，Z值=0.750，双侧P值=0.627，可认为该厂的销售收益服从正态分布。34练习一某医生欲比较某新疗法与传统疗法治疗肾综合征出血热患者的降温效果，将病人随机分为两组，分别以新疗法与传统疗法治疗，以用药开始的体温降至正常值所用的时间（小时）为疗效指标，结果见表1，试比较两种疗法的退热时间有无差别？新疗法新疗法传统疗法传统疗法25363040324435483750395639594260466448195240表1 新疗法与传统疗法治疗结果(小时)35练习二 批批 次次测测 定定 方方 法法甲甲乙乙丙丙丁丁12.722.46 3.95 3.86 22.32 2.42 4.31 3.95 32.48 2.22 4.52 3.30 表2 四种方法测定甘蓝叶核黄素浓度36练习三 某研究者欲研究保健食品对小鼠抗疲劳作用，将同种属的小鼠按性别和年龄相同、体重相近配成对子共10对，并将每对中的两只小鼠随机分到保健食品两个不同的剂量组，测量小鼠负重游泳时间（min，负重5%体重），结果见表3，问不同剂量组的小鼠负重游泳时间有无差别？37表3 不同剂量组小鼠负重游泳时间(min)小鼠对子号小鼠对子号中剂量组中剂量组高剂量组高剂量组11415.22135.5315144176.55135.561813.5717.510810.210.2910101010.59.51113.86.8123.033.481315.25.51416.5938练习四某地监测大气中SO2的日均浓度，按不同功能区设置采样点，结果见下表，问各功能区SO2的日均浓度有无差别？对照区对照区工业区工业区商业区商业区居民区居民区1046723133830665501352307096304854080266951151851677630表4 不同功能区大气中SO2的日均浓度39练习五 四种病人痰液内嗜酸性粒细胞的检查结果见表7。问四种疾病患者痰液内嗜酸粒细胞的等级有无差别?白细胞等级白细胞等级支气管扩张支气管扩张肺水肿肺水肿肺癌肺癌病毒性呼吸病毒性呼吸道感染道感染1353257595336221表5 四种疾病痰液嗜酸粒细胞等级比较40练习六从某装瓶机出来得30盒化妆品的重量（单位：克）如下，问该瓶装机工作是否正常？71.671.071.870.370.572.971.071.070.171.871.970.370.969.371.267.367.667.767.668.168.067.569.867.569.770.069.170.471.069.941练习七以SPSS自带文件carsSav为例，试分析马力（horse）和加速度（accel）的数据是否服从于正态分布。422024/5/22 周三43