统计211

标题: 如何解决SPSS中Nemenyi法检验的编程问题，为何最后P值全为1？ [打印本页]
作者: 梅子    时间: 2012-8-2 20:09
标题: 如何解决SPSS中Nemenyi法检验的编程问题，为何最后P值全为1？
各路高手请帮忙看一下~
data list free
/Hc r1 r2 r3 r4 r5 N n1 n2 n3 n4 n5.
begin data
120.269 145.16 233.98 248.56 257.31 367.49 500 100 100 100 100 100
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)/(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 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 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 x45 = (r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*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 p23 = 1-cdf.chisq(x23,2).
compute p24 = 1-cdf.chisq(x24,2).
compute p25 = 1-cdf.chisq(x25,2).
compute p34 = 1-cdf.chisq(x34,2).
compute p35 = 1-cdf.chisq(x35,2).
compute p45 = 1-cdf.chisq(x45,2).
execute.
为什么最后得的结果P值全为1？“COMPUTE H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+ (r4*n4)**2/n4+(r5*n5)**2/n5)/(N*(N+1))-3*(N+1).”中的12是固定的吗？以及“3*(N+1).”中的3也是固定的吗？？？
作者: 275932488    时间: 2012-8-5 23:00
修改后可用程序如下：
  1. data list free/Hc r1 r2 r3 r4 r5 N n1 n2 n3 n4 n5.
  2. begin data
  3. 120.269 145.16 233.98 248.56 257.31 367.49 500 100 100 100 100 100
  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))/(N*(N+1))-3*(N+1).
  6. 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).

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

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

  20. compute x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*c).

  22. compute p12=1-cdf.chisq(x12,2).
  23. compute p13=1-cdf.chisq(x13,2).
  24. compute p14=1-cdf.chisq(x14,2).
  25. compute p15=1-cdf.chisq(x15,2).

  27. compute p23=1-cdf.chisq(x23,2).
  28. compute p24=1-cdf.chisq(x24,2).
  29. compute p25=1-cdf.chisq(x25,2).

  31. compute p34=1-cdf.chisq(x34,2).
  32. compute p35=1-cdf.chisq(x35,2).

  34. compute p45=1-cdf.chisq(x45,2).

  36. execute.
复制代码
作者: 梅子    时间: 2012-8-7 17:51
275932488 发表于 2012-8-5 23:00
修改后可用程序如下：

谢谢啦！用你改后的程序P值就不一样了。但是我怎么觉得你的和我的基本上是一样的啊，除了分行那里之外。我的哪里出错啦？？？
作者: 梅子    时间: 2012-8-7 18:11
275932488 发表于 2012-8-5 23:00
修改后可用程序如下：

请问编程时和输入法状态有关吗？
作者: abc886y365hxg    时间: 2012-8-7 22:07
梅子 发表于 2012-8-7 18:11
请问编程时和输入法状态有关吗？

原数据放上来，我再给你试试界面的
作者: 莫强求    时间: 2016-4-6 14:08
请问楼主您begin data 下边的数据是什么？H检验的结果不就是有无显著性么，没有数字啊？
作者: 275932488    时间: 2016-4-6 21:41
莫强求 发表于 2016-4-6 14:08
请问楼主您begin data 下边的数据是什么？H检验的结果不就是有无显著性么，没有数字啊？

http://www.tj211.com/thread-911-1-1.html参考下这个帖子



欢迎光临 统计211 (http://www.tj211.com/) Powered by Discuz! X3.2