Kruskal-Wallis非参数检验后多重比较
kruskal wallis检验、kruskal wallis h检验、kruskal–wallis 检验
SPSS的秩和检验没有两两比较,“Mann-Whitney U检验”是用于两独立样本秩和检验(教科书上写的是:Wilcoxon Rank Sum test,威尔克逊秩和检验,两者等价),不是用于两两比较的。 应选用”Kruskal-Wallis H检验“,若有统计学意义,可进一步用”Nemenyi检验”作两两比较,只是需手工计算,或使用其他软件。 也有人利用SPSS的特点,可进行两两比较的SPSS秩和检验方法。 步骤1 编秩: Transform>>Rank Cases…, DATA>>Variable(s),单击OK。 运行后在原数据中产生新变量:RDATA,即按照变量DATA值编秩后的新变量。 步骤2 Analyxe>>Compare Means>>One-Way ANOVA…, RDATA>>Dependent List, GROUP>>Factor。 单击“Post Hoc…”,选择SNK法,单击Continue, 单击OK。
******STATA实现方法:
. list
+—————–+
| g1 g2 x |
|—————–|
1. | 1 1 74.77 |
2. | 2 4 79.76 |
3. | 2 1 80 |
4. | 2 3 81.93 |
5. | 1 4 82.95 |
|—————–|
6. | 1 5 83.03 |
7. | 2 5 83.36 |
8. | 1 3 85.33 |
9. | 2 2 94.12 |
10. | 3 4 96.46 |
|—————–|
11. | 1 2 97.37 |
12. | 3 3 98.21 |
13. | 3 5 98.52 |
14. | 3 1 100 |
15. | 3 2 100 |
+—————–+
. egen rx=rank(x)
. list
+————————+
| g1 g2 x rx |
|————————|
1. | 1 1 74.77 1 |
2. | 2 4 79.76 2 |
3. | 2 1 80 3 |
4. | 2 3 81.93 4 |
5. | 1 4 82.95 5 |
|————————|
6. | 1 5 83.03 6 |
7. | 2 5 83.36 7 |
8. | 1 3 85.33 8 |
9. | 2 2 94.12 9 |
10. | 3 4 96.46 10 |
|————————|
11. | 1 2 97.37 11 |
12. | 3 3 98.21 12 |
13. | 3 5 98.52 13 |
14. | 3 1 100 14.5 |
15. | 3 2 100 14.5 |
+————————+
. oneway rx g1,bon
Analysis of Variance
Source SS df MS F Prob > F
————————————————————————
Between groups 176.4 2 88.2 10.27 0.0025
Within groups 103.1 12 8.59166667
————————————————————————
Total 279.5 14 19.9642857
Bartlett’s test for equal variances: chi2(2) = 1.5154 Prob>chi2 = 0.469
Comparison of rank of (x) by g1
(Bonferroni)
Row Mean-|
Col Mean | 1 2
———+———————-
2 | -1.2
| 1.000
|
3 | 6.6 7.8
| 0.012 0.004
. kwcomp x g1
Kruskal-Wallis Multiple Compare (秩和检验后多重比较)
Kruskal-Wallis equality-of-populations rank test
+———————+
| g1 | Obs | Rank Sum |
|—-+—–+———-|
| 1 | 5 | 31.00 |
| 2 | 5 | 25.00 |
| 3 | 5 | 64.00 |
+———————+
chi-squared = 8.820 with 2 d.f.
probability = 0.0122
chi-squared with ties = 8.836 with 2 d.f.
probability = 0.0121
Kruskal-Wallis H 检验后的多重比较
grp vs grp 平 均 秩 次 平均秩次差 t-value p(双测) p(单侧)
_________________________________________________________________________________
1 vs 2 6.2000 5.0000 1.2000 0.6467 0.5300 0.2650
1 vs 3 6.2000 12.8000 6.6000 3.5570 0.0039 0.0020
2 vs 3 5.0000 12.8000 7.8000 4.2038 0.0012 0.0006