统计211

 找回密码
 立即注册

QQ登录

只需一步,快速开始

统计211 门户 查看主题

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

发布者: linxiuquan001 | 发布时间: 2010-7-25 11:13| 查看数: 84581| 评论数: 85|帖子模式

程序一(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. 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).

  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 x14=(r1-r4)**2/((N*(N+1)/12)*(1/n1+1/n4)*c).
  10. compute x15=(r1-r5)**2/((N*(N+1)/12)*(1/n1+1/n5)*c).
  11. compute x16=(r1-r6)**2/((N*(N+1)/12)*(1/n1+1/n6)*c).
  12. compute x17=(r1-r7)**2/((N*(N+1)/12)*(1/n1+1/n7)*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 x27=(r2-r7)**2/((N*(N+1)/12)*(1/n2+1/n7)*c).

  18. compute x34=(r3-r4)**2/((N*(N+1)/12)*(1/n3+1/n4)*c).
  19. compute x35=(r3-r5)**2/((N*(N+1)/12)*(1/n3+1/n5)*c).
  20. compute x36=(r3-r6)**2/((N*(N+1)/12)*(1/n3+1/n6)*c).
  21. compute x37=(r3-r7)**2/((N*(N+1)/12)*(1/n3+1/n7)*c).

  22. compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).
  23. compute x46=(r4-r6)**2/((N*(N+1)/12)*(1/n4+1/n6)*c).
  24. compute x47=(r4-r7)**2/((N*(N+1)/12)*(1/n4+1/n7)*c).

  25. compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
  26. compute x57=(r5-r7)**2/((N*(N+1)/12)*(1/n5+1/n7)*c).

  27. compute x67=(r6-r7)**2/((N*(N+1)/12)*(1/n6+1/n7)*c).

  28. compute p12=1-cdf.chisq(x12,6).
  29. compute p13=1-cdf.chisq(x13,6).
  30. compute p14=1-cdf.chisq(x14,6).
  31. compute p15=1-cdf.chisq(x15,6).
  32. compute p16=1-cdf.chisq(x16,6).
  33. compute p17=1-cdf.chisq(x17,6).

  34. compute p23=1-cdf.chisq(x23,6).
  35. compute p24=1-cdf.chisq(x24,6).
  36. compute p25=1-cdf.chisq(x25,6).
  37. compute p26=1-cdf.chisq(x26,6).
  38. compute p27=1-cdf.chisq(x27,6).

  39. compute p34=1-cdf.chisq(x34,6).
  40. compute p35=1-cdf.chisq(x35,6).
  41. compute p36=1-cdf.chisq(x36,6).
  42. compute p37=1-cdf.chisq(x37,6).

  43. compute p45=1-cdf.chisq(x45,6).
  44. compute p46=1-cdf.chisq(x46,6).
  45. compute p47=1-cdf.chisq(x47,6).

  46. compute p56=1-cdf.chisq(x56,6).
  47. compute p57=1-cdf.chisq(x57,6).

  48. compute p67=1-cdf.chisq(x67,6).

  49. 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. 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.
复制代码
参考文献:

SPSS在完全随机设计多个样本间多重比较Nemenyi秩和检验中的应用.pdf

127.31 KB, 下载次数: 2920

点评

10组的见41楼  发表于 2011-6-15 23:30

最新评论

点评 回复 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 的帖子

真的可以了啊,非常感谢啊!!真厉害,不然我的论文都不知道怎么搞了,哎呀,谢谢谢谢!!
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则


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

GMT+8, 2018-4-25 16:51 , Processed in 0.122528 second(s), 30 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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