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

linxiuquan001 · 发表于 2010-7-25 11:13:19

程序一（7组）：

data list free/Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7.
begin data
* 请于此行输入Kruskal-Wallis H 检验的结果（即Hc r1 r2 r3 r4 r5 r6 r7 N n1 n2 n3 n4 n5 n6 n7，尽量精确）
end data.
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).
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 x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c).
compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).
compute x17=(r1-r7)**2/((N*(N+1)/12)*(1/n1+1/n7)*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 x25=(r2-r5)**2/((N*(N+1)/12)*(1/n2+1/n5)*c).
compute x26=(r2-r6)**2/((N*(N+1)/12)*(1/n2+1/n6)*c).
compute x27=(r2-r7)**2/((N*(N+1)/12)*(1/n2+1/n7)*c).
compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
compute x35=(r3-r5)**2/((N*(N+1)/12)*(1/n3+1/n5)*c).
compute x36=(r3-r6)**2/((N*(N+1)/12)*(1/n3+1/n6)*c).
compute x37=(r3-r7)**2/((N*(N+1)/12)*(1/n3+1/n7)*c).
compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).
compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).
compute x47=(r4-r7)**2/((N*(N+1)/12)*(1/n4+1/n7)*c).
compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
compute x57=(r5-r7)**2/((N*(N+1)/12)*(1/n5+1/n7)*c).
compute x67=(r6-r7)**2/((N*(N+1)/12)*(1/n6+1/n7)*c).
compute p12=1-cdf.chisq(x12,6).
compute p13=1-cdf.chisq(x13,6).
compute p14=1-cdf.chisq(x14,6).
compute p15=1-cdf.chisq(x15,6).
compute p16=1-cdf.chisq(x16,6).
compute p17=1-cdf.chisq(x17,6).
compute p23=1-cdf.chisq(x23,6).
compute p24=1-cdf.chisq(x24,6).
compute p25=1-cdf.chisq(x25,6).
compute p26=1-cdf.chisq(x26,6).
compute p27=1-cdf.chisq(x27,6).
compute p34=1-cdf.chisq(x34,6).
compute p35=1-cdf.chisq(x35,6).
compute p36=1-cdf.chisq(x36,6).
compute p37=1-cdf.chisq(x37,6).
compute p45=1-cdf.chisq(x45,6).
compute p46=1-cdf.chisq(x46,6).
compute p47=1-cdf.chisq(x47,6).
compute p56=1-cdf.chisq(x56,6).
compute p57=1-cdf.chisq(x57,6).
compute p67=1-cdf.chisq(x67,6).
execute.

复制代码

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

程序二（3组）：

data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
* 请于此行输入Kruskal-Wallis H 检验的结果（即Hc r1 r2 r3 N n1 n2 n3，尽量精确）
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.

复制代码

参考文献：

wwy · 发表于 2010-8-18 23:15:37

领悟了！
真是高！
SPSS越接触越觉得妙不可言！

lykyoung · 发表于 2010-10-12 23:27:07

我的有十组，我改下程序去

275932488 · 发表于 2010-10-22 13:20:12

lykyoung 发表于 2010-10-12 23:27
我的有十组，我改下程序去

改完上传上来共享。{:1_168:}

grace · 发表于 2010-12-11 23:50:15

请教下高手啊，为什么我输入程序后，显示Error # 4135 in column 1. Text: 10.415
>You must not specify the column location for variables read in FREE or LIST
>format. 是不是我第三条命令输的不对啊？？
求各位高手指点下啊。。。。。。感激不尽！！！！！

275932488 · 发表于 2010-12-12 00:18:32

回复 grace 的帖子

可以把你的程序发过来看下吗？

grace · 发表于 2010-12-12 00:25:55

回复 275932488 的帖子

data list free/Hc r1 r2 r3 N n1 n2 n3
begin data
10.415 7.88 18.88 10.75 24 8 8 8
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.

275932488 · 发表于 2010-12-12 17:50:38

回复 grace 的帖子

帮你修改了一下，可以了，如下：
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
10.415 7.88 18.88 10.75 24 8 8 8
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.

275932488 · 发表于 2010-12-12 18:11:16

回复 grace 的帖子

程序第1行与第5行各少了一个“.”，这个“.”是用来结束语句的，不能省略。

grace · 发表于 2010-12-12 20:43:27

回复 275932488 的帖子

真的可以了啊，非常感谢啊！！真厉害，不然我的论文都不知道怎么搞了，哎呀，谢谢谢谢！！

帐号		自动登录	找回密码
密码			立即注册

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

点评