多组等级资料比较的秩和检验及组间两两比较的SAS实现
郭志武
【摘 要】目的 探讨多组等级资料的秩和检验及组间两两比较的SAS实现.方法 结合实例介绍多组等级资料Kruskal-Wallis H检验和两两比较Nemenyi检验的计算方法和步骤,通过编制SAS程序一次性完成Kruskal-Wallis H检验和Nemenyi检验.结果 运行SAS程序得到可靠结果.结论 编制SAS程序可以有效实现多组等级资料的秩和检验及组间两两比较.
【期刊名称】《中国医院统计》
【年(卷),期】2018(025)003
【总页数】2页(P233-234)
【关键词】等级资料;秩和检验;两两比较;SAS
【作 者】郭志武
金融高管排毒养颜胶囊的作用与功效【作者单位】518033 深圳市中医院
【正文语种】中 文
多组独立样本等级资料是临床常见的资料类型。针对多组独立样本等级资料比较的假设检验通常采用Kruskal-Wallis H检验，若检验结果具有统计学意义，则可以进一步采用Nemenyi检验进行组间两两比较[1]。Kruskal-Wallis H检验利用通用的统计软件可以完成，但Nemenyi检验却不能直接从统计软件包的集成功能或模块中完成，需要另外编写自定义程序。刘伟等采用SPSS菜单操作结合编程实现了Nemenyi检验[2]，本文通过编制SAS程序一次性完成Kruskal-Wallis H检验和Nemenyi检验。
1 方法与步骤
参考文献[1]已详细介绍了Kruskal-Wallis H检验和Nemenyi检验的原理、方法和步骤。本文只针对多组独立样本的等级资料结合实例介绍计算方法和步骤。
文献[1]第131页实例如下：
4种疾病患者痰液内嗜酸性粒细胞的检查结果见表1。问4种疾病患者痰液内嗜酸性粒细胞有无差别？
表1 4种疾病患者痰液内嗜酸性粒细胞比较嗜酸性粒细胞(1)支气管扩张(2)肺水肿(3)肺癌(4)病毒性呼吸道感染(5)合计(6)秩范围(7)平均秩(8)-0353111～116 +25751912～3021 ++95332031～5040.5+++62201051～6055.5Ri739.5 436.5 409.5 244.5 ni17 15 17 11 60Ri 43.5029.1024.0922.23
(1)首先进行Kruskal-Wallis H检验
检验统计量H按以下公式计算：
(1)
式中,ni为各样本例数，N为样本总例数(N=∑ni)，Ri为各样本秩和。
由于等级资料存在较多相同秩，需按式(2)进行校正，计算HC值：
(2)
(3)
式中,C为校正系数，tj为各等级相同秩次的个数(表1第6列)。雾灵山风景区
点钞视频
统计量HC近似服从自由度为(组数-1)的χ2分布，据此可求得P值并得出结论。
(2)若Kruskal-Wallis H检验具有统计学意义，则进一步进行组间两两比较的Nemenyi检验
组间比较的统计量为χ2，按以下公式计算：
(4)
式中,和分别为参与比较的两组的平均秩，ni和nj分别为两组的样本例数，N为样本总例数，C为按式(3)计算的校正系数。
统计量χ2服从自由度为(组数-1)的χ2分布，据此可求得P值并得出结论。
2 SAS程序实现
以表1数据为例，编制SAS程序如下：
data temp;
input gn en; /*gn—组数，en—等级分类数*/
佩服英语
do g=1 to gn;
input gp $ @@;
do effect=1 to en;
input f@@; /*各组各等级的频数，等级按从小到大排列*/安昌古镇
output;
end;
end;
datalines;
4 4
A 0 2 9 6
B 3 5 5 2
C 5 7 3 2
D 3 5 3 0
;
run;
/*Kruskal-Wallis检验*/
proc npar1way data=temp wilcoxon;
class gp;
var effect;
freq f;
run;
/*Nemenyi检验*/高怡平
data a; /*改写数据*/
t temp;
do i=1 to f;
瞻怎么组词output;
end;
run;
proc rank data=a out=b; /*生成秩*/
var effect;
ranks r;
run;
proc sql;
/*计算各组秩和、平均秩、例数*/
create table c as
