数据类型为
分析四组（药物组、贴组、药物加贴组、对照组）之间的疗效（临床控制，显效，好转，无效）的差异，非参数等级资料比较的Kruskal-Wallis H秩和检验 差异有统计学意义，现在想做两两比较，不知道怎么做，求高手指点，书上有Nemenyi方法，在网上找到代码，可结果出不来，求指点：四组：
data list free
/Hc r1 r2 r3 r4N n1 n2 n3 n4.
begin data
27.625 44 33.4623.29 13.25 56 14 14 14 14
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.
前面的数据改后结果出来都是1 可能不对 应该怎么弄？？？？？
再用Ordinal 回归做行不行，求说明？？？？
高手在哪里 救救小弟吧
扫码加好友，拉您进群
全部版块
我的主页
收藏
