原创投稿微信公众号（敬请关注）
开启辅助访问 统计211使用手册

统计211

 找回密码
 立即注册

QQ登录

只需一步，快速开始

统计211 › 门户 › 查看主题

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

发布者: linxiuquan001 | 发布时间: 2010-7-25 11:13| 查看数: 103251| 评论数: 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.

  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).

  8. compute c=H/Hc.

  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).

  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).

  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).

  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).

  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).

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

  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).

  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).

  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).

  55. compute p45=1-cdf.chisq(x45,6).
  56. compute p46=1-cdf.chisq(x46,6).
  57. compute p47=1-cdf.chisq(x47,6).

  59. compute p56=1-cdf.chisq(x56,6).
  60. compute p57=1-cdf.chisq(x57,6).

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

  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.

  6. compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3))/(N*(N+1))-3*(N+1).

  8. compute c=H/Hc.

  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).

  13. compute x23=(r2-r3)**2/((N*(N+1)/12)*(1/n2+1/n3)*c).

  15. compute p12=1-cdf.chisq(x12,2).
  16. compute p13=1-cdf.chisq(x13,2).

  18. compute p23=1-cdf.chisq(x23,2).

  20. execute.
复制代码
参考文献：

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

127.31 KB, 下载次数: 5917

点评

linxiuquan001
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 的帖子

真的可以了啊，非常感谢啊！！真厉害，不然我的论文都不知道怎么搞了，哎呀，谢谢谢谢！！
123456789下一页
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

浏览过的版块


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

GMT+8, 2025-11-1 02:02 , Processed in 0.094686 second(s), 29 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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