[有新版本了！！]Kruskal-Wallis及其事后的Nemenyi检验【软件】

abc886y365hxg

    本检验15元一次！！（少抽一包烟，少去和一圈）   
    用matlab开发了一个统计检验的程序：Kruskal-Wallis及其事后的Nemenyi检验。主要是在以前的版本上升了级，有需要的联系我哈，也希望大家多多发帖，多多捧场本论坛板块！！
这个检验的简介如下：（很多人会错误的把有一类适合这个检验的做成卡方检验哦，核心期刊上见惯不怪了！！再此好心提醒）
      当遇到的资料不符合正态分布，或为评分资料、等级频数资料、开口资料等情况时，通常就需要考虑采用进行非参数检验。此种情况下的多组完全随机资料，检验方法通常采用Kruskal-Wallis检验，如果处理结果显示总体存在统计学差异，就需要进行Nemenyi检验，以确定是哪些组别之间存在差异！
新版本界面：
【计量资料】

【计数资料】

老版本界面及操作展示：

abc886y365hxg

喜羊羊 发表于 2012-5-11 16:04
非常谢谢您噶，这么热心，我从昨天一直弄到现在，我最后找到了stata 的 编程程序。刚刚计算出来结果，和您 ...

it's my pleasure

喜羊羊

abc886y365hxg 发表于 2012-5-11 14:01
统计量H=56.6716>临界值9.4880；拒绝H0,有差异性

i            j         卡方       显著性

就在给您上传数据后，我发现了了一篇文献，教怎么用stata do 命令来做nemenyi test。 所以我试了下，刚刚做出来。所以请您不要误会了噶，因为马上就要毕业了，但我才开始做统计，所以有点着急。

喜羊羊

abc886y365hxg 发表于 2012-5-11 14:01
统计量H=56.6716>临界值9.4880；拒绝H0,有差异性

i            j         卡方       显著性

非常谢谢您噶，这么热心，我从昨天一直弄到现在，我最后找到了stata 的 编程程序。刚刚计算出来结果，和您帮我做的一模一样。非常非常谢谢您

abc886y365hxg

统计量H=56.6716>临界值9.4880；拒绝H0,有差异性

i            j         卡方       显著性
1           2       2.605……无显著差异  
1           3       2.325……无显著差异  
1           4       0.428……无显著差异  
1            5       43.679……有显著差异
2           3       0.005……无显著差异  
2           4       2.807……无显著差异  
2            5       36.595……有显著差异
3           4       2.597……无显著差异  
3            5       35.709……有显著差异
4            5       14.573……有显著差异

abc886y365hxg

喜羊羊 发表于 2012-5-11 09:07
您好版主，请问下您能帮我做下吗？非常谢谢，因为我不知道怎么编程不是统计专业的。非常谢谢， ...

中午回来给你回复

喜羊羊

您好版主，请问下您能帮我做下吗？非常谢谢，因为我不知道怎么编程不是统计专业的。非常谢谢，我会向周围朋友推广的。我把数据传上。
研究变量            分组
67.94        0
60.99        0
72.98        0
59.05        0
70.04        0
55.14        0
89.52        0
74.02        0
21.98        0
65.76        0
65.03        0
94.99        0
72.20        0
79.27        0
88.71        0
77.34        0
75.98        0
73.33        0
82.20        0
71.54        0
72.70        0
77.22        0
69.15        0
75.00        0
68.22        0
48.40        0
64.04        0
83.15        0
68.64        0
71.77        0
41.40        0
65.09        0
72.29        0
27.68        0
76.96        0
82.28        0
63.45        0
76.70        0
84.59        0
63.50        0
64.35        0
76.88        0
75.22        0
19.38        0
60.42        0
78.08        0
90.35        0
69.82        0
73.02        0
58.17        0
69.74        0
59.82        0
65.73        0
47.24        0
67.60        0
55.30        0
84.91        0
66.38        0
72.53        0
86.15        0
60.57        0
72.33        0
64.94        0
66.42        0
68.11        0
76.59        0
85.07        0
83.17        0
67.34        0
71.57        0
81.17        0
51.71        0
26.91        0
71.07        0
82.67        0
70.16        0
60.00        0
68.21        0
59.76        0
71.48        0
79.46        0
64.62        0
71.07        0
31.01        0
74.10        0
71.49        1
81.13        1
79.63        1
82.82        1
75.64        1
83.42        1
54.91        1
62.18        1
72.15        1
78.53        1
70.50        1
72.21        1
77.39        1
90.32        1
80.56        2
71.38        2
91.37        2
96.02        2
87.80        2
80.31        2
66.54        2
47.81        2
82.08        2
73.73        2
41.49        2
78.13        2
55.24        2
83.33        2
51.69        3
45.20        3
62.09        3
82.27        3
71.72        3
80.39        3
44.90        3
66.12        3
33.73        3
79.34        3
76.55        3
2.17                4
0.00          4
0.00         4
0.00         4
2.23         4
0.71         4
1.04         4
2.12         4
2.77         4
0.63         4
0.39         4
0.00         4
3.50         4
0.00         4
0.83         4
0.63         4
5.05         4
0.41         4
0.34         4
0.00         4

abc886y365hxg

艺颂雨 发表于 2012-4-12 22:45
谁能帮忙看一下，为什么运行后，没有任何结果出来啊？
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin da ...

你能像27楼那样把每组数据和分组情况列出来吗？并告诉我数据资料的类型：计数或者计量？

艺颂雨

谁能帮忙看一下，为什么运行后，没有任何结果出来啊？
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
24.662 116.73 76.86 118.23 200 20 85 95
end data.
compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).
compute c=H/Hc.
compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
compute p12=1-cdf.chisq(x12,2).
compute p13=1-cdf.chisq(x13,2).
compute p23=1-cdf.chisq(x23,2).
execute.

abc886y365hxg

riskeen 发表于 2011-11-26 18:29
麻烦帮忙看下，谢谢

你这个用kruskal-wallis检验出来有差异，但用Nemenyi检验进行事后检验的时候，却未显示那两组之间有差异，可能事后检验应该选用其他方法吧（kruskal-wallis检验的事后检验不止Nemenyi检验一种，各有适用范围）
具体处理结果：
统计量H=26.6461拒绝H0,有差异性
组别       组别      卡方      差异性
1           2       7.226……无显著差异  
1           3       0.409……无显著差异  
1           4       0.014……无显著差异  
1           5       9.724……无显著差异  
1           6       0.001……无显著差异  
1           7        0.74……无显著差异  
1           8       1.068……无显著差异  
1           9       0.258……无显著差异  
2            3       11.071……无显著差异
2           4       7.883……无显著差异  
2           5       0.185……无显著差异  
2           6       7.388……无显著差异  
2           7       3.341……无显著差异  
2           8       2.738……无显著差异  
2           9       4.754……无显著差异  
3           4        0.27……无显著差异  
3            5       14.118……无显著差异
3           6       0.371……无显著差异  
3           7       2.248……无显著差异  
3           8       2.798……无显著差异  
3           9       1.315……无显著差异  
4            5       10.483……无显著差异
4           6       0.008……无显著差异  
4           7        0.96……无显著差异  
4           8       1.329……无显著差异  
4           9       0.393……无显著差异  
5           6       9.911……无显著差异  
5           7       5.099……无显著差异  
5           8       4.346……无显著差异  
5           9       6.815……无显著差异  
6           7       0.792……无显著差异  
6           8       1.131……无显著差异  
6           9       0.289……无显著差异  
7           8        0.03……无显著差异  
7           9       0.124……无显著差异  
8           9       0.276……无显著差异  

riskeen

麻烦帮忙看下，谢谢