等级资料常用检验方法

等级资料常用检验方法演示文稿第一页，共52页。优选等级资料常用检验方法第二页，共52页。

等级资料划分的两种情况：特点：观察结果具有等级差别。按性质划分：如药物疗效分为痊愈、显效、好转、无效；麻醉效果分为Ⅰ、Ⅱ、Ⅲ、Ⅳ级等。按数量分组：数据两端不能确切测定的计量资料。如抗体滴度分为>1:20,1:20,1:40,1:80,<1:80；年龄分为<10,10~,20~,40~,≥60等。第三页，共52页。？等级资料的分析方法是否和一般计数资料的检验方法相同呢？

等级资料的分析应该选用什么方法？第四页，共52页。实例1

考察硝苯地平治疗老年性支气管炎的疗效，治疗组60人，用硝苯地平治疗，对照组58人，常规治疗，两组患者的性别、年龄、病程无显著性差异，治疗结果见表1。

第五页，共52页。

在变量窗口“variableview”中设定变量在数据窗口“dataview”中录入数据使用“WeightCases”过程权重记录

SPSS软件默认一行就是一条记录，而我们是以频数格式录入数据，即相同的观测值只录入一次，另加一个频数变量（count）用于记录该数值共出现的次数。因此我们使用此过程:1、建立数据库第六页，共52页。2、卡方检验分析结果：

第七页，共52页。

两组疗效的构成百分比差异无统计学意义。两组的疗效无差别。结论：（×）第八页，共52页。

一般的χ2检验不适用于有序分类资料——“等级”、“程度”、“优劣”的比较分析。因为检验只利用了两组构成比提供的信息，损失了有序指标包含的“等级”信息。注意：第九页，共52页。例如，假定两组的显效例数和有效例数互换，见表2。显然，两组反映的信息是不同的，但由于两组的结构百分比无变化（仅仅是位置不同），不改变检验结果。（χ2=5.224，P>0.05）第十页，共52页。等级资料正确的统计分析方法：

非参数统计的秩和检验

Kendall、spearman等级相关

CMH卡方检验

Ridit分析线性趋势卡方检验有序变量的Logistic回归分析第十一页，共52页。一、非参秩和检验

由于非参数检验法不考虑数据的分布规律，检验不涉及总体参数，检验统计量多是人们在总结经验的基础上创造出来的，所以这类检验方法的特点是针对性强。但是不同设计、不同目的所用的非参数检验法是不同的。

第十二页，共52页。v

单向有序行列表

在表的两个方向上的分类中，一个方向（横向）无顺序和等级概念，另一个方向（纵向）是有顺序的分类，称为单向有序行列表。

第十三页，共52页。两组独立样本等级资料比较的Mann-Whitney秩和检验

以表1为例。将无效、有效、显效三个疗效等级数量化，数值用平均秩号，然后比较各组平均秩号的大小。

第十四页，共52页。两组的平均秩号分别为：治疗组：R1=(6×10.5+19×40+35×89)/60=65.6对照组：R2=（14×10.5+20×40+24×89）/58=53.1

经秩和检验，u=2.169，P<0.05，两组疗效差异有统计学意义，因为治疗组平均秩号大于对照组，所以治疗组疗效好。

计算两组秩号并进行秩和检验第十五页，共52页。SPSS窗口操作过程：第十六页，共52页。分析结果：第十七页，共52页。结论：两组疗效差异有统计学意义，且治疗组效果好于对照组。μ=-2.169，P<0.05

第十八页，共52页。同样方法，对表2数据进行秩和检验，结果如下：

第十九页，共52页。μ=0.731，P>0.05

结论：两组疗效差异没有统计学意义。

第二十页，共52页。两组配对样本等级资料比较的Wilcoxon秩和检验多组等级资料比较的Kruskal-Wallis秩和检验

该方法对K(K>2)组独立样本进行K个总体分布函数相同假设的检验，是在Wilcoxon秩和检验基础上扩展的方法，称为K-W检验。

第二十一页，共52页。例2

对54例牙病患者的64颗患牙的根端形态不同分为3种，X线片显示喇叭口状为A型，管壁平行状为B型，管壁由聚状为C型比较不同根端形态患牙的疗效有否差别。第二十二页，共52页。1、建立数据库1.1定义变量group:1A型2B型3C型

result:1成功2进步3失败

count:例数1.2录入数据1.3权重频数SPSS软件操作步骤：第二十三页，共52页。2、分析：第二十四页，共52页。3、结果第二十五页，共52页。结论：按α=0.05的检验水平，三组间差异有统计学意义。H=6.528，P=0.038第二十六页，共52页。注意：第二十七页，共52页。

双向有序资料行列表——Kendall等级相关法和Spearman等级相关分析法

例3

某病病情与疗效的关系

第二十八页，共52页。

两表的区别仅在于病情“极重”组和“中”组的数据进行了互换。第二十九页，共52页。

分别对上面两个表格数据进行χ2检验和多组等级资料比较的Kruskal-Wallis秩和检验。对于表4：χ2=40.000，P=0.000

H=24.896，P=0.000对于表5：χ2=40.000，P=0.000

H=24.896，P=0.000两种检验都无法表达表4和表5的差别，直观地看，表4的资料显示病情越轻者疗效越好，表5却未显示这种趋势来。第三十页，共52页。此时我们选用Kendall和Spearman等级相关分析法分别计算相关系数t和rs。计算公式：t=n：总例数m：最长对角线上的格子数S：专用统计量Kendall等级相关意义：当一个变量的等级为标准时，另一个变量的等级与它不一致的情况（可分析两个以及多个变量间的等级相关性）。

第三十一页，共52页。rs=1-n：总例数d：每一对值的等级差Spearman等级相关意义：两个变量之间的等级相关性。（只适用于分析两个变量关系）Spearman等级相关公式：第三十二页，共52页。1.建立数据库2.录入数据3.权重频数4.界面操作（以表4为例）：Analyze——Correlate——Bivariate——Row(s)：疗效

——Column(s)：病情——Statistics——Kendall’stau-b，Spearman——OKSPSS操作演示：第三十三页，共52页。表4检验结果：第三十四页，共52页。表5检验结果：第三十五页，共52页。

重复测量等级资料时间趋势检验——CMH卡方检验

实例4

在某药治疗闭塞性动脉炎的临床试验中,治疗26例下肢溃疡的病人溃疡改善情况见表6，评价该药有无促进溃疡愈合的作用。

第三十六页，共52页。本例特点：

1.个体的重复测量

2.溃疡随时间的变化趋势

ＣＭＨ卡方检验,Cochran-Mantel-Hansel检验简称，包括非零相关、行平分差和一般联系3种检验方法。对于这种重复测量的等级变量,行列变量均为等级变量,应该作非零相关检验。本方法仅限于检验线性趋势。第三十七页，共52页。计算公式：

Qcs=结果：χ2=4.7424，P=0.0294结论：溃疡改善程度随着治疗时间延长有变好的趋势。

第三十八页，共52页。多组等级资料的两两比较——Ridit分析

Ridit是“Relativetoanidentifieddistributionintegraltransformation”的首个字母缩写，意指对于一个确认的分布作积分变换。Ridit分析是一种关于等级资料进行对比组与标准组比较的假设检验方法，其基本思想是先确定一个标准组作为特定总体，求得各等级的Ridit值，标准组平均Ridit值理论上可以证明等于0.5，其他各组与标准组比较，看其可信区间是否与0.5重叠，来判断组间的统计学显著性，最后得出专业解释。第三十九页，共52页。Ridit分析适用范围：

1.两组或两组以上等级资料的比较和分析

2.两端数据不确切的计量资料分组转换成计数资料的分析，如血清滴度等。特点：简便、直观、适用性广。第四十页，共52页。Ridit分析关键步骤：

一、确定标准组：

1.利用已知的标准分布（传统方法）作为标准组。如某药物大规模的观察研究结果，计算不同疗效的R值。

2.以例数最多的一组作为标准组。

3.各组的例数都差不多时，可把各组观察结果合并起来作为标准组。

第四十一页，共52页。二、计算标准组的Ridit值：三、利用对照组计算各组的平均R值四、置信区间判断五、统计检验:

u检验、t检验、χ2检验Ridit分析关键步骤：

第四十二页，共52页。Ridit分析（实例5）：

表7三种方剂对某妇科病患者治疗效果比较第四十三页，共52页。PEM3.1操作过程：

1.建立数据库

1.1定义变量名

1.2录入数据2.界面操作：其他统计分析

——Ridit分析

——分析目的：多个样本比较

——分组变量：分组

——类别变量

：类别

——频数变量：例数

——确定

第四十四页，共52页。主要结果：检验水准:α=0.05选用同一参照组(各样本合计为参照组)

样本数k=3等级数:4第四十五页，共52页。主要结果：第四十六页，共52页。各样本平均R

───────────────────

样本名频数合计平均R

───────────────────

第1组3610.4819

第2组580.6287

第3组770.4881

───────────────────

主要结果：第四十七页，共52页。卡方检验:

卡方值=13.0887

自由度v=3

概率P=0.0044

结论：三种药物对妇科病的疗效差异有统计学意义。主要结果：第四十八页，共52页。三者之间究