《等级资料常用检验方法.ppt》由会员分享，可在线阅读，更多相关《等级资料常用检验方法.ppt（57页珍藏版）》请在三一办公上搜索。
1、等级资料常用检验方法,临床流行病学应用研究室 周罗晶,在医学资料中，特别是临床医学资料中，常常遇到一些定性指标，如临床疗效的评价、疾病的临床分期、症状严重程度的临床分级、中医诊断的一些临床症状等，对这些指标常采用分成若干等级然后分类计数的办法来解决它的量化问题，这样的资料我们在统计学上称为有序变量（ordered variable）或半定量资料，也称为等级资料（ranked data）。,等级资料定义：,等级资料划分的两种情况：,特点：观察结果具有等级差别。,按性质划分：如药物疗效分为痊愈、显效、好转、无效；麻醉效果分为、级等。按数量分组：数据两端不能确切测定的计量资料。如抗体滴度分为1:20
2、,1:20,1:40,1:80,1:80；年龄分为10,10,20,40,60等。,？,等级资料的分析方法是否和一般计数资料的检验方法相同呢？,等级资料的分析应该选用什么方法？,实例1 考察硝苯地平治疗老年性支气管炎的疗效，治疗组60人，用硝苯地平治疗，对照组58人，常规治疗，两组患者的性别、年龄、病程无显著性差异，治疗结果见表1。,在变量窗口“variable view”中设定变量 在数据窗口“data view”中录入数据 使用“Weight Cases”过程权重记录 SPSS软件默认一行就是一条记录，而我们是以频数格式录入数据，即相同的观测值只录入一次，另加一个频数变量（count）用于
3、记录该数值共出现的次数。因此我们使用此过程:,1、建立数据库,2、卡方检验,分析结果：,两组疗效的构成百分比差异无统计学意义。两组的疗效无差别。,结论：,（）,一般的2检验不适用于有序分类资料“等级”、“程度”、“优劣”的比较分析。因为检验只利用了两组构成比提供的信息，损失了有序指标包含的“等级”信息。,注意：,例如，假定两组的显效例数和有效例数互换，见表2。,显然，两组反映的信息是不同的，但由于两组的结构百分比无变化（仅仅是位置不同），不改变检验结果。（2=5.224，P0.05）,等级资料正确的统计分析方法：,非参数统计的秩和检验 Kendall、spearman等级相关 CMH卡方检验
4、Ridit分析 线性趋势卡方检验 有序变量的Logistic回归分析,一、非参秩和检验,由于非参数检验法不考虑数据的分布规律，检验不涉及总体参数，检验统计量多是人们在总结经验的基础上创造出来的，所以这类检验方法的特点是针对性强。但是不同设计、不同目的所用的非参数检验法是不同的。,v 单向有序行列表,在表的两个方向上的分类中，一个方向（横向）无顺序和等级概念，另一个方向（纵向）是有顺序的分类，称为单向有序行列表。,两组独立样本等级资料比较的Mann-Whitney秩和检验,以表1为例。将无效、有效、显效三个疗效等级数量化，数值用平均秩号，然后比较各组平均秩号的大小。,两组的平均秩号分别为：治疗组
5、：R1=(610.5+1940+3589)/60=65.6对照组：R2=（1410.5+2040+2489）/58=53.1 经秩和检验，u=2.169，P0.05，两组疗效差异有统计学意义，因为治疗组平均秩号大于对照组，所以治疗组疗效好。,计算两组秩号并进行秩和检验,Analyze Nonparametric Tests 2 independent Samples Test variable List：result Grouping variable：group Define groups：group1：1；group2：2，Test Type：Mann-Whitney OK,SPSS窗口操
6、作过程：,分析结果：,结论：两组疗效差异有统计学意义，且治疗组效果好于对照组。,=-2.169，P0.05,同样方法，对表2数据进行秩和检验，结果如下：,=0.731，P0.05 结论：两组疗效差异没有统计学意义。,两组配对样本等级资料比较的Wilcoxon秩和检验,多组等级资料比较的Kruskal-Wallis秩和检验,该方法对K(K2)组独立样本进行K个总体分布函数相同假设的检验，是在Wilcoxon秩和检验基础上扩展的方法，称为K-W检验。,例2 对54例牙病患者的64颗患牙的根端形态不同分为3种，X线片显示喇叭口状为A型，管壁平行状为B型，管壁由聚状为C型,比较不同根端形态患牙的疗效有
7、否差别。,1、建立数据库1.1定义变量 group:1 A型 2 B型 3 C型 result:1成功 2进步 3失败 count:例数1.2录入数据1.3权重频数,SPSS软件操作步骤：,Analyze Nonparametric Tests K independent Samples Test variable List：result Grouping variable：group Define range：minimum：1；maximum：3 Continue Test Type：Kruskal-Wallis HOK,2、分析：,3、结果,结论：按=0.05的检验水平，三组间差异有统计
8、学意义。,H=6.528，P=0.038,计算结果中显示的2值并不是2检验，只是Kruskal-Wallis Test的检验统计量H，此时近似2分布，所以按2分布的近似值来确定概率，它的自由度=组数-1。,注意：,双向有序资料行列表Kendall等级相关法和Spearman等级相关分析法,例3 某病病情与疗效的关系,两表的区别仅在于病情“极重”组和“中”组的数据进行了互换。,分别对上面两个表格数据进行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.
9、896，P=0.000,两种检验都无法表达表4和表5的差别，直观地看，表4的资料显示病情越轻者疗效越好，表5却未显示这种趋势来。,此时我们选用Kendall和Spearman等级相关分析法分别计算相关系数t和rs。,计算公式：,t=,n：总例数 m：最长对角线上的格子数 S：专用统计量,Kendall等级相关意义：当一个变量的等级为标准时，另一个变量的等级与它不一致的情况（可分析两个以及多个变量间的等级相关性）。,rs=1-,n：总例数 d：每一对值的等级差,Spearman等级相关意义：两个变量之间的等级相关性。,（只适用于分析两个变量关系）,Spearman等级相关公式：,1.建立数据库2
10、.录入数据3.权重频数4.界面操作（以表4为例）：Analyze Correlate Bivariate Row(s)：疗效 Column(s)：病情 Statistics Kendalls tau-b，Spearman OK,SPSS操作演示：,表4检验结果：,表5检验结果：,重复测量等级资料时间趋势检验 CMH卡方检验,实例4 在某药治疗闭塞性动脉炎的临床试验中,治疗26例下肢溃疡的病人溃疡改善情况见表6，评价该药有无促进溃疡愈合的作用。,本例特点：1.个体的重复测量 2.溃疡随时间的变化趋势,卡方检验,Cochran-Mantel-Hansel检验简称，包括非零相关、行平分差和一般联系3
11、种检验方法。对于这种重复测量的等级变量,行列变量均为等级变量,应该作非零相关检验。本方法仅限于检验线性趋势。,计算公式：,Qcs=,结果：2=4.7424，P=0.0294结论：溃疡改善程度随着治疗时间延长有变好的趋势。,多组等级资料的两两比较Ridit分析,Ridit是“Relative to an identified distribution integral transformation”的首个字母缩写，意指对于一个确认的分布作积分变换。,Ridit分析是一种关于等级资料进行对比组与标准组比较的假设检验方法，其基本思想是先确定一个标准组作为特定总体，求得各等级的Ridit值，标准组平均
12、Ridit值理论上可以证明等于0.5，其他各组与标准组比较，看其可信区间是否与0.5重叠，来判断组间的统计学显著性，最后得出专业解释。,Ridit分析适用范围：1.两组或两组以上等级资料的比较和分析 2.两端数据不确切的计量资料分组转换成计数资料的分析，如血清滴度等。特点：简便、直观、适用性广。,Ridit分析关键步骤：,一、确定标准组：1.利用已知的标准分布（传统方法）作为标准组。如某药物大规模的观察研究结果，计算不同疗效的R值。2.以例数最多的一组作为标准组。3.各组的例数都差不多时，可把各组观察结果合并起来作为标准组。,二、计算标准组的Ridit值：三、利用对照组计算各组的平均R值四、置
13、信区间判断五、统计检验:u检验、t检验、2检验,Ridit分析关键步骤：,Ridit分析（实例5）：,表7 三种方剂对某妇科病患者治疗效果比较,PEM3.1操作过程：,1.建立数据库 1.1 定义变量名 1.2 录入数据2.界面操作：其他统计分析 Ridit分析 分析目的：多个样本比较 分组变量：分组 类别变量：类别 频数变量：例数 确定,主要结果：,检验水准:=0.05 选用同一参照组(各样本合计为参照组)样本数 k=3 等 级 数:4,主要结果：,各样本平均R 样本名 频数合计 平均R 第1组 361 0.4819 第2组 58 0.6287 第3组 77 0.4881,主要结果：,卡方检
14、验:卡方值=13.0887 自由度 v=3 概 率 P=0.0044 结论：三种药物对妇科病的疗效差异有统计 学意义。,主要结果：,三者之间究竟有何种差异呢？,1.标准组平均R值：=248.00/496=0.5,2.标准组方差、标准差的计算：,方差=160.16 2482/496/（496-1）=0.0730标准差=0.2702,标准组方差、标准差的计算,对比组平均R值,糖衣组=0.4819 黄体酮组=0.6287 复方组=0.4881,如果等级是由劣到优的顺序排列，则平均 Ridit值越大越好，否则越小越好。但我们不能单纯以各组平均Ridit值的大小来判断各组的差异，因为Ridit值是一个点
15、估计，没有考虑抽样误差，而样本正好等于总体指标的可能性很小，因此应该结合可信区间或统计检验方法综合判断。,注意：,计算对比组可信区间,对比组标准误：,对比组95%可信区间：,计算对比组95%CI并进行判断,糖衣组95%CI：0.45340.5103黄体酮组95%CI：0.55780.6997复方组95%CI：0.42650.5497,判断：通过可信区间比较对比组与标准组的差别。看对比组95%可信区间是否包括标准组Ridit值（0.5），不包括，则P0.05，差别有统计学意义。否则反之。,对比组95%CI图,标准组,多个对比组的两两比较,统计界限：u0.05,差异无统计意义 u1.96,P0.05,差异有统计意义 u2.58,P0.01,差异有统计意义,u 1,2=3.84,u 1,3=0.18,u 2,3=2.99,小 结,
