医学统计学：非参数检验SPSS

非参数检验的SPSS实现一、两独立样本的检验

[例1]调查某厂的铅作业工人7人和非铅作业工人10人的血铅值（μg/100g）如下，问两组工人的血铅值有无差别？

非铅作业组556791213151821铅作业组17182025344344数据准备

定义2个变量：

1.group（非铅作业组为1，铅作业组为2），2.Pb。按顺序输入数据。

统计分析

分析→非参数统计→两独立样本→Ｐb选入检验变量列表框，group选入分组变量框→定义组别，定义组1为1，组2为2→检验类型：Mann-WhitneyU→确定结果解释

结果表明，第1组的平均秩次（MeanRank）为5.95，第2组的平均秩次为13.36，Mann-Whitney

U=4.5，W=59.5，z=-2.98,精确双侧概率P=0.001，可认为铅作业组工人的血铅值高于非铅作业组。

二、多个独立样本的检验

[例2]随机抽样得右表三组人的血桨总皮质醇测定值（μg/L），试比较三组人有无差异？正常人单纯性肥胖皮质醇增多症0.40.69.81.91.210.22.22.010.62.52.413.02.83.114.03.14.114.83.75.015.63.95.915.64.67.421.67.013.624.0数据准备

定义2个变量：1.group（正常人为1，单纯性肥胖为2，皮质醇增多症为3），2.x(总皮质醇测定值)

按顺序输入数据。统计分析

分析→非参数统计→Ｋ个独立样本→x选入检验变量框，组别选入分组变量框→定义组别，最小值为1，最大值为3→确定结果解释

结果表明，1至3组的平均秩次（MeanRank）分别为9.65、11.75、25.10，χ2值（即H值）为18.1219，P=0.0001；可认为三组人的血桨总皮质醇测定值有差异，根据本例情况可看出皮质醇增多症组高于其他两组人。大鼠对别正常饲料组VitE缺乏饲料组137.225.7220.925.1331.418.8441.433.5539.834.0639.328.3736.126.2831.918.3表2正常饲料组、VitE缺乏饲料组大鼠肝中VitA含量三、配对资料秩和检验[例3]数据准备

两个变量:

1.正常饲料2.VitE

缺乏

按顺序输入数据。

统计分析

分析→非参数统计→两个相关样本检验→将正常饲料和VitE

选入检验对列框→检验类型:

选择前两种(Wilcoxon和符号检验)→确定结果解释

首先显示两变量正常饲料和VitE缺乏的例数、均数、标准差、最大值和最小值；

配对符号秩和检验（WilcoxonMatched-PairsSigned-RanksTest）结果，其平均秩分别为5.00和1.00，Z=-2.380，双侧P=0.017，可认为两组大鼠肝中VitA含量有差别，饲料中缺乏VitE会使大鼠肝中VitA含量降低；

符号检验（SignTest）的结果，双侧P=0.0703，则认为两组大鼠肝中VitA含量无差别。

在这种情况下，应取配对符号秩和检验（Wilcoxon）结果，因两法比较之下，配对符号秩和检验较为敏感，效率较高。四、随机区组设计多个样本的秩和检验

8名受试对象在相同条件下分别接受4种不同频率声音的刺激，他们的反应率（%）见下表。问4种不同频率声音的刺激的反应率（%）是否有差别？

[例4]8名受事对象对4种不同频率声音刺激的反应率（%）比较

受试号频率A频率B频率C频率D

18.49.69.811.7211.612.711.81239.49.110.49.849.88.79.91258.388.68.668.69.89.610.678.9910.611.4

87.88.28.510.8

数据准备

定义4个变量：

1.频率A2.频率B3.频率C4.频率D

统计分析

分析→非参数统计→K个相关样本检验→频率A、频率B、频率C、频率D均进入检验对列框→检验类型:选择Friedman→确定描述性统计量

N均值标准差极小值极大值频率A89.101.197.811.6频率B89.391.47812.7频率C89.901.088.511.8频率D810.861.198.612秩

秩均值频率A1.38频率B2.00频率C2.94频率D3.69检验统计量(a)N8卡方15.152df3渐近显著性0.00169表６两种疗效对肾小球肾病的疗效比较五、等级资料秩和检验疗效

消炎痛

合剂

无效

154部分缓解

69基本缓解

45完全缓解

219

[例5]等级资料的SPSS操作要点：1.定义3个变量，即疗效、疗法（或疾病）和频数。2.在进行统计分析前，需对频数加权。3.采用2个和K个独立样本的非参数检验。

[例6]某村发生一种地方病，其住户沿一条河排列，调查时对发病的住户标记为“1”，共9例，对非发病的住户标记为“0”，共17户：问病户的分布排列是呈聚集趋势，还是随机分布？01100010010000110010000101六、Runs过程（游程检验）主要功能

依时间或其他顺序排列的有序数列中，具有相同的事件或符号的连续部分称为一个游程。调用Runs过程可进行游程检验，即用于检验序列中事件发生过程的随机性分析。假定下面是由0和1组成的一个这种变量的样本：

01100010010000110010000101

其中相同的0（或相同的1）在一起称为一个游程（单独的0或1也算）。这个数据中有７个0组成的游程和７个1组成的游程。数据准备

定义1个变量为epi。

按住户顺序输入数据，发病的住户为1,非发病的住户为0。统计分析

分析→非参数统计→游程检验→epi选入检验变量列框

→选择割点→确定

1、中位数：中位数作割点，小于临界割点的为一类，大于或等于临界割点的为另一类；

2、众数：众数作割点，小于临界割点的为一类，大于或等于临界割点的为另一类；

3、均数：均数作临界割点，小于临界割点的为一类，大于或等于临界割点的为另一类；

4、定制：用户指定割点，小于定制割点的为一类，大于或等于割点的为另一类；

本例选定制，在其方框中键入1（根据需要选项，本例是0、1二分变量，故临界割点值用1），再点击OK钮即可。结果解释

检验结果可见本例游程个数为14，检验临界割点值（Testvalue）=1.00，

Z=0.325，双侧P=0.746。所以认为此地方病的病户沿河分布的情况无聚集性，而是呈随机分布。1-SampleK-S过程

调用此过程可对单样本进行Kolmogorov-SmirnovZ检验，它将一个变量的实际频数分布与正态分布（Normal）、均匀分布（Uniform）、泊松分布（Poisson）进行比较。数据准备

以SPSS自带文件cars．Sav为例，试分析马力（horse）和加速度（accel）的数据是否服从于正态分布。统计分析

分析→非参数统计→单样本

Kolmogorov-Smirnov检验→马力和加速度选入检验变量列框

→检验分布：正态分布→确定结果解释练习一

某医生欲比较某新疗法与传统疗法治疗肾综合征出血热患者的降温效果，将病人随机分为两组，分别以新疗法与传统疗法治疗，以用药开始的体温降至正常值所用的时间（小时）为疗效指标，结果见表1，试比较两种疗法的退热时间有无差别？新疗法传统疗法25363040324435483750395639594260466448195240表1新疗法与传统疗法治疗结果(小时)练习二

批次测定方法甲乙丙丁12.722.463.953.8622.322.424.313.9532.482.224.523.30表2四种方法测定甘蓝叶核黄素浓度练习三

某研究者欲研究保健食品对小鼠抗疲劳作用，将同种属的小鼠按性别和年龄相同、体重相近配成对子，共10对，并将每对中的两只小鼠随机分到保健食品两个不同的剂量组，测量小鼠负重游泳时间（min，负重5%体重），结果见表3，问不同剂量组的小鼠负重游泳时间有无差别？表3不同剂量组小鼠负重游泳时间(min)

小鼠对子号中剂量组高剂量组11415.22135.5315144176.55135.561813.5717.510810.210.2910101010.59.51113.86.8123.033.481315.25.51416.59练习四

某地监测大气中SO2的日均浓度，按不同功能区设置采样点，结果见下表，问各功能区SO2的日均浓度有无差别？对照区工业区商业区居民区1046723133830665501352307096304854080266951151851677630表4

不同功能区大气中SO2的日均浓度练习五

四种病人痰液内嗜酸性粒细胞的检查结果见表7。问四种疾病患者痰液内嗜酸粒细胞的等级有无差别?白细胞等级支气管扩张肺水肿肺癌病毒性呼吸道感染－1353＋