统计211

 找回密码
 立即注册

QQ登录

只需一步,快速开始

楼主: linxiuquan001
打印 上一主题 下一主题

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

  [复制链接]
51
 楼主| 发表于 2011-12-1 22:13:37 | 只看该作者
riskeen 发表于 2011-11-26 11:07
为什么我的程序提示:
>Warning # 522
>Undefined error #522 - Cannot open text file "D:\spss\en\wind ...

少了其中一组的样本量
52
发表于 2012-1-1 12:01:11 | 只看该作者
我的怎么到了execute就停止了,没有结果呢?下面是我的程序,请大侠帮我看看
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
21.435 24.29 14.29 15.43 4 28 7 7 7 7
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.


53
发表于 2012-1-2 10:39:17 | 只看该作者
hyxxjpp 发表于 2012-1-1 12:01
我的怎么到了execute就停止了,没有结果呢?下面是我的程序,请大侠帮我看看
data list free/Hc r1 r2 r3  ...

结果会新生成一个数据表。
54
发表于 2012-1-9 19:24:56 | 只看该作者
我的数据分为4组,每组5个样本,自己写的Nemenyi检验的程序如下,但老是出不了结果,求高人指点!跪谢!
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4 n5.
begin data
12.406 3.00 14.60 10.20 14.20 20 5 5 5 5  
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.
55
发表于 2012-2-2 13:32:55 | 只看该作者
请问我有6组数据两两比较,但我的程序不能执行,请楼主看看,谢谢!
data list free/Hc r1 r2 r3 r4 r5 r6 N n1 n2 n3 n4 n5 n6.6 E/ W0 Z7 l+ R; r% ~0 t* C2 n
begin data
59.512 192.01 145.28 162.96 180.10 47.85 185.72
end data.# M- t8 E5 x& D. l$ W

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))/(N*(N+1))-3*(N+1).

compute c=H/Hc.
# B% B2 o! g7 R( O' d
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).+ u0 ^0 z- x. z- U
compute x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c)." {8 m6 X5 N* B; A" u
compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).& [.
m1 t; d# W" f

compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).* Z% K, \5 n# a- T9 s# A4 z
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).. u2
D9 p. F: i+ G: l$ v* }7 H4 p& xs" D( Q7 L1 n8 _

compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).! K1 }, t  S; B. e4 }0 n
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).

7 k6 y! Q) `1 y2 u7 U5 Z
compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).- E0 g1 a  a6 d, ?1 p, ?$ q
compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).


compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).2 N% m% M, M  v% P/ T. f
/ Y4 {' P1 o
' z( s, p2 y6 L
compute p12=1-cdf.chisq(x12,2).
compute p13=1-cdf.chisq(x13,2).0 A5 L* d6 F3 d& n/ k
compute p14=1-cdf.chisq(x14,2).
compute p15=1-cdf.chisq(x15,2).+ s# u, K  x! x, C8 _. T( a
compute p16=1-cdf.chisq(x16,2).
9 T5 I' L  X) ^: @+ i
6 y. f# k$ g. D& ]7 S) U" ^
compute p23=1-cdf.chisq(x23,2).! @* F0 A$ i: ]; d5 ^
compute p24=1-cdf.chisq(x24,2).
compute p25=1-cdf.chisq(x25,2).
compute p26=1-cdf.chisq(x26,2).

) |; M  t$ J# w6 b, [& N. i5 N
compute p34=1-cdf.chisq(x34,2).0 p) N; F$ p! w. }1 b, o
compute p35=1-cdf.chisq(x35,2).
compute p36=1-cdf.chisq(x36,2).  i/ b/ c1 o1 K( ~! P5 v: z
~7 F3 A; H# B% @! H3 a1 s8 y
0 w9 }$ g) W5 s
compute p45=1-cdf.chisq(x45,2).5 M( p' w, K7 X$ `3 L6 m
compute p46=1-cdf.chisq(x46,2).


compute p56=1-cdf.chisq(x56,2).

5 B: ~( }) q# i- w
execute.



56
发表于 2012-2-2 17:40:00 | 只看该作者
发错了,应该是data list free/Hc r1 r2 r3 r4 r5 r6 N n1 n2 n3 n4 n5 n6.
begin data
59.512 192.01 145.28 162.96 180.10 47.85 185.72 320 80 62 70 60 30 18
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))/(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 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 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 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 x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
compute p12=1-cdf.chisq(x12,2).
compute p13=1-cdf.chisq(x13,2).
compute p14=1-cdf.chisq(x14,2).
compute p15=1-cdf.chisq(x15,2).
compute p16=1-cdf.chisq(x16,2).
compute p23=1-cdf.chisq(x23,2).
compute p24=1-cdf.chisq(x24,2).
compute p25=1-cdf.chisq(x25,2).
compute p26=1-cdf.chisq(x26,2).
compute p34=1-cdf.chisq(x34,2).
compute p35=1-cdf.chisq(x35,2).
compute p36=1-cdf.chisq(x36,2).
compute p45=1-cdf.chisq(x45,2).
compute p46=1-cdf.chisq(x46,2).
compute p56=1-cdf.chisq(x56,2)
execute.



57
发表于 2012-2-2 17:42:48 | 只看该作者
仍不能正确执行,请问这是为何?
58
发表于 2012-2-13 16:03:22 | 只看该作者
请帮忙看看,我的对吗?怎么没法运行呢?
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data
10.974 9.29 14.50 3.80 18 7 6 5.
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.
59
 楼主| 发表于 2012-2-15 08:48:34 | 只看该作者
流沙浪 发表于 2012-2-13 16:03
请帮忙看看,我的对吗?怎么没法运行呢?
data list free/Hc r1 r2 r3 N n1 n2 n3.
begin data

数据最后多了一个“.”,修改后的程序如下:
  1. data list free/Hc r1 r2 r3 N n1 n2 n3.
  2. begin data
  3. 10.974 9.29 14.50 3.80 18 7 6 5
  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.
复制代码
60
 楼主| 发表于 2012-2-15 08:53:51 | 只看该作者
昆仑过客 发表于 2012-2-2 17:42
仍不能正确执行,请问这是为何?

倒数第二行少了一个结束“.”,修改后的程序如下:
  1. data list free/Hc r1 r2 r3 r4 r5 r6 N n1 n2 n3 n4 n5 n6.
  2. begin data
  3. 59.512 192.01 145.28 162.96 180.10 47.85 185.72 320 80 62 70 60 30 18
  4. end data.
  5. Compute
  6. 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))/(N*(N+1))-3*(N+1).
  7. compute c=H/Hc.
  8. compute x12=(r1-r2)**2/((N*(N+1)/12)*(1/n1+1/n2)*c).
  9. compute x13=(r1-r3)**2/((N*(N+1)/12)*(1/n1+1/n3)*c).
  10. compute x14=(r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
  11. compute x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c).
  12. compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).
  13. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).
  14. compute x24=(r2-r4)**2/((N*(N+1)/12)*(1/n2+1/n4)*c).
  15. compute x25=(r2-r5)**2/((N*(N+1)/12)*(1/n2+1/n5)*c).
  16. compute x26=(r2-r6)**2/((N*(N+1)/12)*(1/n2+1/n6)*c).
  17. compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
  18. compute x35=(r3-r5)**2/((N*(N+1)/12)*(1/n3+1/n5)*c).
  19. compute x36=(r3-r6)**2/((N*(N+1)/12)*(1/n3+1/n6)*c).
  20. compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).
  21. compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).
  22. compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
  23. compute p12=1-cdf.chisq(x12,2).
  24. compute p13=1-cdf.chisq(x13,2).
  25. compute p14=1-cdf.chisq(x14,2).
  26. compute p15=1-cdf.chisq(x15,2).
  27. compute p16=1-cdf.chisq(x16,2).
  28. compute p23=1-cdf.chisq(x23,2).
  29. compute p24=1-cdf.chisq(x24,2).
  30. compute p25=1-cdf.chisq(x25,2).
  31. compute p26=1-cdf.chisq(x26,2).
  32. compute p34=1-cdf.chisq(x34,2).
  33. compute p35=1-cdf.chisq(x35,2).
  34. compute p36=1-cdf.chisq(x36,2).
  35. compute p45=1-cdf.chisq(x45,2).
  36. compute p46=1-cdf.chisq(x46,2).
  37. compute p56=1-cdf.chisq(x56,2).
  38. execute.
复制代码
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则


免责声明|关于我们|小黑屋|联系我们|赞助我们|统计211 ( 闽ICP备09019626号  

GMT+8, 2025-3-31 23:56 , Processed in 0.090310 second(s), 19 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表