多个独立样本比较的秩和检验（Kruskal-WallisH）

Kruskal-Walis H检验，用于推断计量资料或等级资料的多个独立样本所来自的多个总体分布是否有差别。在理论上检验假设H0为总体分布位置相同。

例 题

比较小白鼠接种三种不同菌型伤寒杆菌9D、11C、和DSC1后存活日数，结果见下表。问小白鼠接种三种不同菌型伤寒杆菌的存活日数有无差别？

9D

2 2 2 3 4 4 4 5 7 7

5

5 6 6 6 7 8 10 12

11C

3

5 6 6 6 7 7 9 10 11 11

DSC1

SPSS操作

组别：9D=1，11C=2，DSC1=3。 1建立数据库

2SPPS操作

分析——非参数检验——K个独立样本

3结果解释

H=9.94，P=0.007,按照α=0.05的检验标准，拒绝H0，接受H1，可以认为小白鼠接种不同菌型伤寒杆菌的存活日数有差别。

两两比较

经过Kruskal-Wallis H检验可知三组的总体分布位置不同，若要

进一步推断是哪两两总体分布位置不同时，可用Nemenyi法检验。Nemenyi法在SPSS中没有直接的窗口操作，需要利用SPSS编程实现。

1编辑程序

文件——新建——语法

2程序代码

n1,n2,n3分别表示总例数和三组例数。

data list free/Hc r1 r2 r3 N n1 n2 n3 . begin data

9.94 8.4 18.78 19.27 30 10 9 11 end data. compute

H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).

HC表示总体检验的H值，r1,r2,r3分别表示三组的平均秩次，N，

compute c=H/Hc.

compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c). compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c). compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c). compute p12=1-cdf.chisq(x12,2). compute p13=1-cdf.chisq(x13,2). compute p23=1-cdf.chisq(x23,2). execute. 3结果解释 运行——全部

两两比卡方 P 较

9D和6.70 0.04 11C

9D和8.13 0.02 DSC1

11C和0.02 0.69 DSC1

9D和11C、DCSI组伤寒杆菌存活天数有统计学差异，11C组和DSC1组没有统计学差异。