本站分享:大数据、数据分析师考试认证培训,包括:Python培训Excel培训Matlab培训SPSS培训SAS培训R语言培训Hadoop培训Amos培训Stata培训Eviews培训广告位

Kruskal-Wallis非参数检验后多重比较_kruskal wallis检验

spss培训 cdadata 2060℃

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

转载请注明:数据分析 » Kruskal-Wallis非参数检验后多重比较_kruskal wallis检验

喜欢 (2)or分享 (0)