多个独立样本两两比较的Nemenyi法检验的SPSS程序（7组与3组）

linxiuquan001

程序一（7组）：

1. data list free/Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7.
   
2. begin data
   
3. * 请于此行输入Kruskal-Wallis H 检验的结果（即Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7，尽量精确）
   
4. end data.
   
5. 
   
6. compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3+(r4*n4)**2/n4+(r5*n5)**2/n5+(r6*n6)**2/n6+(r7*n7)**2/n7))/(N*(N+1))-3*(N+1).
   
7. 
   
8. compute c=H/Hc.
   
9. 
   
10. compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
    
11. compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
    
12. compute x14=(r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
    
13. compute x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c).
    
14. compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).
    
15. compute x17=(r1-r7)**2/((N*(N+1)/12)*(1/n1+1/n7)*c).
    
16. 
    
17. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
    
18. compute x24=(r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).
    
19. compute x25=(r2-r5)**2/((N*(N+1)/12)*(1/n2+1/n5)*c).
    
20. compute x26=(r2-r6)**2/((N*(N+1)/12)*(1/n2+1/n6)*c).
    
21. compute x27=(r2-r7)**2/((N*(N+1)/12)*(1/n2+1/n7)*c).
    
22. 
    
23. compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
    
24. compute x35=(r3-r5)**2/((N*(N+1)/12)*(1/n3+1/n5)*c).
    
25. compute x36=(r3-r6)**2/((N*(N+1)/12)*(1/n3+1/n6)*c).
    
26. compute x37=(r3-r7)**2/((N*(N+1)/12)*(1/n3+1/n7)*c).
    
27. 
    
28. compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).
    
29. compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).
    
30. compute x47=(r4-r7)**2/((N*(N+1)/12)*(1/n4+1/n7)*c).
    
31. 
    
32. compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
    
33. compute x57=(r5-r7)**2/((N*(N+1)/12)*(1/n5+1/n7)*c).
    
34. 
    
35. compute x67=(r6-r7)**2/((N*(N+1)/12)*(1/n6+1/n7)*c).
    
36. 
    
37. compute p12=1-cdf.chisq(x12,6).
    
38. compute p13=1-cdf.chisq(x13,6).
    
39. compute p14=1-cdf.chisq(x14,6).
    
40. compute p15=1-cdf.chisq(x15,6).
    
41. compute p16=1-cdf.chisq(x16,6).
    
42. compute p17=1-cdf.chisq(x17,6).
    
43. 
    
44. compute p23=1-cdf.chisq(x23,6).
    
45. compute p24=1-cdf.chisq(x24,6).
    
46. compute p25=1-cdf.chisq(x25,6).
    
47. compute p26=1-cdf.chisq(x26,6).
    
48. compute p27=1-cdf.chisq(x27,6).
    
49. 
    
50. compute p34=1-cdf.chisq(x34,6).
    
51. compute p35=1-cdf.chisq(x35,6).
    
52. compute p36=1-cdf.chisq(x36,6).
    
53. compute p37=1-cdf.chisq(x37,6).
    
54. 
    
55. compute p45=1-cdf.chisq(x45,6).
    
56. compute p46=1-cdf.chisq(x46,6).
    
57. compute p47=1-cdf.chisq(x47,6).
    
58. 
    
59. compute p56=1-cdf.chisq(x56,6).
    
60. compute p57=1-cdf.chisq(x57,6).
    
61. 
    
62. compute p67=1-cdf.chisq(x67,6).
    
63. 
    
64. execute.

复制代码

注：少于7组的，删掉相应的行与变量即可。
举例3组的如下：

程序二（3组）：

1. data list free/Hc r1 r2 r3 N n1 n2 n3.
   
2. begin data
   
3. * 请于此行输入Kruskal-Wallis H 检验的结果（即Hc r1 r2 r3 N n1 n2 n3，尽量精确）
   
4. end data.
   
5. 
   
6. compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).
   
7. 
   
8. compute c=H/Hc.
   
9. 
   
10. compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
    
11. compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
    
12. 
    
13. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
    
14. 
    
15. compute p12=1-cdf.chisq(x12,2).
    
16. compute p13=1-cdf.chisq(x13,2).
    
17. 
    
18. compute p23=1-cdf.chisq(x23,2).
    
19. 
    
20. execute.

复制代码

参考文献：

tyl0829

linxiuquan001 发表于 2011-1-10 22:14
回复 zf123456 的帖子

程序第一行少了一个“.”

想再问一下 我执行命令前 需要将数据进行编秩以后再执行命令 还是直接拿原数据执行命令呢

tyl0829

linxiuquan001 发表于 2011-1-10 22:14
回复 zf123456 的帖子

程序第一行少了一个“.”

好哒 看到啦 哇哇  大神就是大神啊！看到结果啦

tyl0829

zf123456 发表于 2011-1-10 17:52
**** 作者被禁止或删除 内容自动屏蔽 ****

我也是结果只显示execute,而没有数据 请问您是怎么解决的了呢

tyl0829

linxiuquan001 发表于 2011-1-3 11:21
回复 jack19850430 的帖子

并非大侠，一起交流{:1_176:}

好厉害啊大神 我照你的粘贴过来果然不提示错误了

I'm_lin.

高手，学习！

文刀火丁

高手真多啊，佩服~~~

DYING

菜鸟表示看不懂

西米

请问，我和前面有好几个童鞋的遇到了同样的问题：Hc是负值，所有的p值都是1.00  请问这个为什么呢？麻烦帮帮忙~~~ 多谢多谢！
data list free
/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
11.952 69.53 99.92 98.60 118.93 192 33 98 40 21
end data.
COMPUTE H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+(r4*n4)**2/n4)/(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 x14 = (r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
compute x23 = (r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
compute x24 = (r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).
compute x34 = (r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
compute p12 = 1-cdf.chisq(x12,2).
compute p13 = 1-cdf.chisq(x13,2).
compute p14 = 1-cdf.chisq(x14,2).
compute p23 = 1-cdf.chisq(x23,2).
compute p24 = 1-cdf.chisq(x24,2).
compute p34 = 1-cdf.chisq(x34,2).
execute.
不止这一组数据。。。。

西米

梅子 发表于 2012-8-2 20:50
请问你的这个问题解决了吗？我遇到和你一样的问题：出来的HC是负的 ，P值都是1 我之前H检验时有差异的。
...

遇到同样的问题啊，看到有人提问还以为有解决，翻了几页都没有看到、、、哪位大侠可以帮帮忙~~~？？
data list free
/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
11.952 69.53 99.92 98.60 118.93 192 33 98 40 21
end data.
COMPUTE H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+(r4*n4)**2/n4)/(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 x14 = (r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
compute x23 = (r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
compute x24 = (r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).
compute x34 = (r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
compute p12 = 1-cdf.chisq(x12,2).
compute p13 = 1-cdf.chisq(x13,2).
compute p14 = 1-cdf.chisq(x14,2).
compute p23 = 1-cdf.chisq(x23,2).
compute p24 = 1-cdf.chisq(x24,2).
compute p34 = 1-cdf.chisq(x34,2).
execute.