《预防医学》 第三节 多组资料的比较

H检验(Kruskal-Wallis法)是用于完全随机设计的多个样本比较的非参数法 。其具体步骤见例21.3 。
例21.3 某地监测大气中SO[XB]2[/XB]的日均浓度 , 按不同功能区设置采样点 , 结果见表21-4“浓度”栏所示 , 问各功能区SO[XB]2[/XB]日均浓度有无差别?
表21-4 某地1990年1月份SO[XB]2[/XB]日均浓度(μg/m[SB]3[/SB])
对照区工业区商业区居民区浓度(1)秩次(2)浓度(3)秩次(4)浓度(5)秩次(6)浓度(7)秩次(8)101467923163387302665155011135283037091863013.548510404802196691651112515851206771763013.5R[XB]i[/XB]15 81 63.5 50.5n[XB]i[/XB]5 5 5 5
(一)建立假设
H[XB]0[/XB]:四个功能区SO[XB]2[/XB]日均浓度总体分布相同
H[XB]1[/XB]:四个功能区SO[XB]2[/XB]日均浓度总体分布不同或不全相同
α=0.05
(二)编秩
先将各组数据由小到大排列 , 再将各组数据由小到大统一编秩 , 不同组的相同数据取其平均秩次 。如本例有2个630 , 分别在第(5)、(7)栏 , 其平均秩次为(13+14)/2=13.5 。
(三)求各组秩和(R[XB]i[/XB])
分别将各组秩次相加得R[XB]i[/XB]
(四)计算统计量H值
按式(21.4)计算 。式中ni为各组观察值个数 , N=Σni
[img]yufangyixue120.jpg[alt][/alt][/img]公式(21.4)
本例[imgz]yufangyixue121.jpg[alt][/alt][/img]
(五)确定P值 , 作出推论
若组数K=3 , 每组例数≤5 , 可查附表21-3“秩和检验用H界值表”得出P值;若超出附表21-3的范围 , 可按v=k-1查x[SB]2[/SB]界值表得出P值 。本例k=4 , 超出附表21-3范围 , 按v=4-1=3查x[SB]2[/SB]界值表 , x[SB]2[/SB][XB]0.01(3)[/XB]=11.34,P<0.01,按α=0.05检验水准拒绝H[XB]0[/XB] , 可认为四种功能区SO[XB]2[/XB]日均浓度有差别 。

    推荐阅读