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

chandler

为什么我的结果做出来，卡方都是负数，而且每一个P值都等于1呢？原数据经Kruskal-Wallis H法秩和检验提示组间炎症活动度差异有显著性，P值小于0.001，卡方值等于77.767。H值是负的。。。
data list free/Hc r1 r2 r3 r4 r5 r6 r7 r8 N n1 n2 n3 n4 n5 n6 n7 n8.
begin data
77.767 166.95 254.9 251.2 265.43 274.85 230.12 134 134.86 415 102 42 35 50 13 101 4 68
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+(r7*n7)**2/n7)+(r8*n8)**2/n8)/(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 x17=(r1-r7)**2/((N*(N+1)/12)*(1/n1+1/n7)*c).
compute x18=(r1-r8)**2/((N*(N+1)/12)*(1/n1+1/n8)*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 x27=(r2-r7)**2/((N*(N+1)/12)*(1/n2+1/n7)*c).
compute x28=(r2-r8)**2/((N*(N+1)/12)*(1/n2+1/n8)*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 x37=(r3-r7)**2/((N*(N+1)/12)*(1/n3+1/n7)*c).
compute x38=(r3-r8)**2/((N*(N+1)/12)*(1/n3+1/n8)*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 x47=(r4-r7)**2/((N*(N+1)/12)*(1/n4+1/n7)*c).
compute x48=(r4-r8)**2/((N*(N+1)/12)*(1/n4+1/n8)*c).
compute x56=(r5-r6)**2/((N*(N+1)/12)*(1/n5+1/n6)*c).
compute x57=(r5-r7)**2/((N*(N+1)/12)*(1/n5+1/n7)*c).
compute x58=(r5-r8)**2/((N*(N+1)/12)*(1/n5+1/n8)*c).
compute x67=(r6-r7)**2/((N*(N+1)/12)*(1/n6+1/n7)*c).
compute x68=(r6-r8)**2/((N*(N+1)/12)*(1/n6+1/n8)*c).
compute x78=(r7-r8)**2/((N*(N+1)/12)*(1/n7+1/n8)*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 p17=1-cdf.chisq(x17,2).
compute p18=1-cdf.chisq(x18,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 p27=1-cdf.chisq(x27,2).
compute p28=1-cdf.chisq(x28,2).
compute p34=1-cdf.chisq(x34,2).
compute p35=1-cdf.chisq(x35,2).
compute p36=1-cdf.chisq(x36,2).
compute p37=1-cdf.chisq(x37,2).
compute p38=1-cdf.chisq(x38,2).
compute p45=1-cdf.chisq(x45,2).
compute p46=1-cdf.chisq(x46,2).
compute p47=1-cdf.chisq(x47,2).
compute p48=1-cdf.chisq(x48,2).
compute p56=1-cdf.chisq(x56,2).
compute p57=1-cdf.chisq(x57,2).
compute p58=1-cdf.chisq(x58,2).
compute p67=1-cdf.chisq(x67,2).
compute p68=1-cdf.chisq(x68,2).
compute p78=1-cdf.chisq(x78,2).
execute.

tongji

回复 chandler 的帖子

H值是负的？？那Hc也是负的吗？

ahfyhp

我用的是SPSS13.0，运行之后只有输入的第三行数字，没有结果。并且出现下列提示：
>Error # 4285 in column 27.  Text: nl
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 11.  Text: H
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.

C




>Error # 4285 in column 50.  Text: c
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x12
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x13
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x14
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x23
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x24
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.


>Error # 4285 in column 25.  Text: x34
>Undefined error #4285 - Cannot open text file "C:\Program Files\SPSS\en\windows\spss.err": No such file or directory
>This command not executed.
我的编程是：
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
10.76 452.70 373.28 442.03 427.76 850 105 120 140 485
end data.
compute H=(12*((r1*n1)**2/nl+(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.

ahfyhp

恳请高手给解答！{:1_176:}

275932488

回复 ahfyhp 的帖子

“compute H=(12*((r1*n1)**2/nl+(r2*n2)**2/n2+(r3*n3)**2/n3+(r4*n4)**2/n4))/(N* (N+1))-3*(N+1).”这条程序中的红色字体部分错了，改成n1。
帮你修改了下，可以了。程序如下：
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
10.76 452.70 373.28 442.03 427.76 850 105 120 140 485
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.

ahfyhp

谢谢谢谢~~你果然是个高手！{:1_168:}

ahfyhp

我想请教一个目前很纠结的问题：
我做的3个病例组和1个对照组的病例对照研究，需要看一个单核苷酸多态性对3个病例组的发病是否有影响，影响有多大。用多分类的Logistic回归分析得出的相关基因型的OR值比较大，4左右，并且OR值的可信区间也跨度很大，1~20，但是P值部室很小，0.03左右。这个结果有没有问题？用多分类的Logistic回归进行多组间的比较还需要进行Bonfferroni校正吗？
恳请高手能指点指点我。{:1_176:}

ahfyhp

请问这个Nemenyi法可以用于多个实验组与一个对照组之间的比较吗？

275932488

回复 ahfyhp 的帖子

不是{:1_165:}

275932488

回复 ahfyhp 的帖子

个人觉得可以。还可以参考以下链接：
http://gw.seu.edu.cn/%BF%CE%BC%F ... %A7/Tjwlkc04/J4.htm