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

linxiuquan001

云破月 发表于 2012-1-9 19:24
我的数据分为4组，每组5个样本，自己写的Nemenyi检验的程序如下，但老是出不了结果，求高人指点！跪谢！
d ...

修改后的程序如下：

1. 
   
2. data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
   
3. begin data
   
4. 12.406 3.00 14.60 10.20 14.20 20 5 5 5 5
   
5. end data.
   
6. Compute
   
7. H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3+(r4*n4)**2/n4))/(N*(N+1))-3*(N+1).
   
8. compute c=H/Hc.
   
9. compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
   
10. compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
    
11. compute x14=(r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
    
12. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
    
13. compute x24=(r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).
    
14. compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
    
15. compute p12=1-cdf.chisq(x12,2).
    
16. compute p13=1-cdf.chisq(x13,2).
    
17. compute p14=1-cdf.chisq(x14,2).
    
18. compute p23=1-cdf.chisq(x23,2).
    
19. compute p24=1-cdf.chisq(x24,2).
    
20. compute p34=1-cdf.chisq(x34,2).
    
21. execute.

复制代码

yixuelau

您好，我的程序语言编好后出现如下：_SLINE 0 2.
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
20.77 46.91 40.08 24.06 67 13 17 20 17
end date.
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 x12=(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).
excute.
_SLINE OFF.
而我的程序语言是：
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
20.77 46.91 40.08 24.06 67 13 17 20 17
end date.
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 x12=(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).
excute.
请您帮我看看，我是哪出错了，非常感谢

275932488

yixuelau 发表于 2012-2-21 02:34
您好，我的程序语言编好后出现如下：_SLINE 0 2.
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin ...

错误太多，请重新编辑：
1 end data 而不是date
2 execute 而不是excute
3 x12重复了
。。。。

yixuelau

275932488 发表于 2012-2-21 13:26
错误太多，请重新编辑：
1 end data 而不是date
2 execute 而不是excute

谢谢，我知道了，修改后出结果了，太感谢您的帮助

275932488

yixuelau 发表于 2012-2-22 00:00
谢谢，我知道了，修改后出结果了，太感谢您的帮助

{:1_176:}

Miki^.⌒γ

还有没有人呀

pyudx1

data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
16.808 20.30 21.60 7.0 31 10 10 11
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.
这是我们的程序，但是总是运行error
data list free/Hc r1 r2 r3 N n1 n2 n3.
execute.

>Error # 119.  Command name: execute
>The in-line data must appear before a procedure can be executed and before the
>data can be modified in the data editor.
>This command not executed.
compute p23=1-cdf.chisq(x23,2).

>Error # 4285 in column 25.  Text: x23
>Incorrect variable name: either the name is more than 64 characters, or it is
>not defined by a previous command.
>This command not executed.
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.
大侠麻烦看看

linxiuquan001

pyudx1 发表于 2012-4-13 10:24
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
16.808 20.30 21.60 7.0 31 10 10 11

你的程序可正常运行，如下：

1. data list free/Hc r1 r2 r3 N n1 n2 n3.
   
2. begin data
   
3. 16.808 20.30 21.60 7.0 31 10 10 11
   
4. end data.
   
5. compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).
   
6. compute c=H/Hc.
   
7. compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
   
8. compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
   
9. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
   
10. compute p12=1-cdf.chisq(x12,2).
    
11. compute p13=1-cdf.chisq(x13,2).
    
12. compute p23=1-cdf.chisq(x23,2).
    
13. execute.
    

复制代码

运行结果见附件：
结果.sav (1.5 KB, 下载次数: 298)

2012-4-13 11:01 上传

点击文件名下载附件

pyudx1

linxiuquan001 发表于 2012-4-13 11:01
你的程序可正常运行，如下：运行结果见附件：

大侠，谢谢你了，但是不知道为什么我总是下载不了你发的附件了。呵呵。难道是我的等级太低了？

linxiuquan001

pyudx1 发表于 2012-4-13 11:49
大侠，谢谢你了，但是不知道为什么我总是下载不了你发的附件了。呵呵。难道是我的等级太低了？

你没有统计币了，下载一个附件需要1个统计币。