H檢驗(Kruskal-Wallis法)是用于完全隨機設計的多個樣本比較的非參數(shù)法。其具體步驟見例21.3。
例21.3 某地監(jiān)測大氣中SO2的日均濃度,按不同功能區(qū)設置采樣點,結(jié)果見表21-4“濃度”欄所示,問各功能區(qū)SO2日均濃度有無差別?
表21-4 某地1990年1月份SO2日均濃度(μg/m3)
對照區(qū) | 工業(yè)區(qū) | 商業(yè)區(qū) | 居民區(qū) | ||||
濃度(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:四個功能區(qū)SO2日均濃度總體分布相同
H1:四個功能區(qū)SO2日均濃度總體分布不同或不全相同
α=0.05
。ǘ)編秩
先將各組數(shù)據(jù)由小到大排列,再將各組數(shù)據(jù)由小到大統(tǒng)一編秩,不同組的相同數(shù)據(jù)取其平均秩次。如本例有2個630,分別在第(5)、(7)欄,其平均秩次為(13+14)/2=13.5。
(三)求各組秩和(Ri)
分別將各組秩次相加得Ri
。ㄋ)計算統(tǒng)計量H值
按式(21.4)計算。式中ni為各組觀察值個數(shù),N=Σni
公式(21.4)
本例