我们热爱生命科学!-生物行

第三节 多组资料的比较

时间:2006-07-01 21:49来源:大众医药网 作者:admin 点击: 104次

第三节 多组资料的比较

  H检验(Kruskal-Wallis法)是用于完全随机设计的多个样本比较的非参数法。其具体步骤见例21.3。

  例21.3 某地监测大气中SO2的日均浓度,按不同功能区设置采样点,结果见表21-4“浓度”栏所示,问各功能区SO2日均浓度有无差别?

表21-4 某地1990年1月份SO2日均浓度(μg/m3

对照区 工业区 商业区 居民区 浓度(1) 秩次(2) 浓度(3) 秩次(4) 浓度(5) 秩次(6) 浓度(7) 秩次(8) 10 1 467 9 231 6 338 7 30 2 665 15 501 11 352 8 30 3 709 18 630 13.5 485 10 40 4 802 19 669 16 511 12 51 5 851 20 677 17 630 13.5 Ri 15   81   63.5   50.5 ni 5   5   5   5

  (一)建立假设

  H0:四个功能区SO2日均浓度总体分布相同

  H1:四个功能区SO2日均浓度总体分布不同或不全相同

  α=0.05

  (二)编秩

  先将各组数据由小到大排列,再将各组数据由小到大统一编秩,不同组的相同数据取其平均秩次。如本例有2个630,分别在第(5)、(7)栏,其平均秩次为(13+14)/2=13.5。

  (三)求各组秩和(Ri

  分别将各组秩次相加得Ri

  (四)计算统计量H值

  按式(21.4)计算。式中ni为各组观察值个数,N=Σni

  公式(21.4)

  本例

  (五)确定P值,作出推论

  若组数K=3,每组例数≤5,可查附表21-3“秩和检验用H界值表”得出P值;若超出附表21-3的范围,可按v=k-1查x2界值表得出P值。本例k=4,超出附表21-3范围,按v=4-1=3查x2界值表,x20.01(3)=11.34,P<0.01,按α=0.05检验水准拒绝H0,可认为四种功能区SO2日均浓度有差别。

附表21-1 符号秩和检验临界值表

对子数n T0.05 T0.02 T0.01 对子数n T0.05 T0.02 T0.01 6 0 — — 16 29 23 19 7 2 0 — 17 34 27 23 8 3 1 0 18 40 32 27 9 5 3 1 19 46 37 32 10 8 5 3 20 52 43 37 11 10 7 5 21 58 49 42 12 13 9 7 22 65 55 48 13 17 12 9 23 73 62 54 14 21 15 12 24 81 69 61 15 25 19 15 25 89 76 68

附表21-2 等级总和数临界值(双侧检验)

n2较大n P

n1=    较 小  n 2 3 4 5 6 7 8 9 10 11 12 13 14 15 4 0.05     10                       0.01     —                       5 0.05   6 11 17                     0.01   — — 15                     6 0.05   7 12 18 26                   0.01   — 10 16 23                   7 0.05   7 13 20 27 36                 0.01   — 10 17 24 32                 8 0.05 3 8 14 21 29 38 49               0.01 — — 11 17 25 34 43               9 0.05 3 8 15 22 31 40 51 63             0.01 — 6 11 18 26 35 45 56             10 0.05 3 9 15 23 32 42 53 65 78           0.01 — 6 12 19 27 37 47 58 71           11 0.05 4 9 16 24 34 44 55 68 81 96         0.01 — 6 12 20 28 38 49 61 74 87         12 0.05 4 10 17 26 35 46 58 71 85 99 115       0.01 — 7 13 21 30 40 51 63 76 90 106       13 0.05 4 10 18 27 37 48 60 73 88 103 119 137     0.01 — 7 14 22 31 41 53 65 79 93 109 125     14 0.05 4 11 19 28 38 50 63 76 91 106 123 141 160   0.01 — 7 14 22 32 43 54 67 81 96 112 129 147   15 0.05 4 11 20 29 40 52 65 79 94 110 127 145 164 185 0.01 — 8 15 23 33 44 56 70 84 99 115 133 151 171 16 0.05 4 12 21 31 42 54 67 82 97 114 131 150 169   0.01 — 8 15 24 34 46 58 72 86 102 119 137 155   17 0.05 5 12 21 32 43 56 70 84 100 117 135 154     0.01 — 8 16 25 36 47 60 74 89 105 122 140     18 0.05 5 13 22 33 45 58 72 87 103 121 139       0.01 — 8 16 26 37 49 62 76 92 108 125       19 0.05 5 13 23 34 46 60 74 90 107 124         0.01 3 9 17 27 38 50 64 78 94 111         20 0.05 5 14 24 35 48 62 77 93 110           0.01 3 9 18 28 39 52 66 81 97           21 0.05 6 14 25 37 50 64 79 95             0.01 3 9 18 29 40 53 68 83             22 0.05 6 15 26 38 51 66 82               0.01 3 10 19 29 42 55 70               23 0.05 6 15 27 39 53 68                 0.01 3 10 19 30 43 57                 24 0.05 6 16 28 40 55                   0.01 3 10 20 31 44                   25 0.05 6 16 28 42                     0.01 3 11 20 32                     26 0.05 7 17 29                       0.01 3 11 21                       27 0.05 7 17                         0.01 4 11                         28 0.05 7                           0.01 4                          

附表21-3 秩和检验用H界值表

n n1 n2 n3 P 0.05 0.01 7 3 2 2 4.71   3 3 1 5.14   8 3 3 2 5.36   4 2 2 5.33   4 3 1 5.21   5 2 1 5.00   9 3 3 3 5.60 7.20 4 3 2 5.44 6.44 4 4 1 4.97 6.67 5 2 2 5.16 6.53 5 3 1 4.96   10 4 3 3 5.73 6.75 4 4 2 5.45 7.04 5 3 2 5.25 6.82 5 4 1 4.99 6.95 11 4 4 3 5.60 7.14 5 3 3 5.65 7.08 5 4 2 5.27 7.12 5 5 1 5.13 7.31 12 4 4 4 5.69 7.65 5 4 3 5.63 7.44 5 5 2 5.34 7.27 13 5 4 4 5.26 7.76 5 5 3 5.71 7.54 14 5 5 4 5.64 7.79 15 5 5 5 5.78 7.98 (责任编辑:泉水)
顶一下
(10)
100%
踩一下
(0)
0%
------分隔线----------------------------
发表评论
请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
评价:
表情:
用户名: 验证码:点击我更换图片