|
1统计币
在做Nemenyi秩和检验
看到有个帖子是这样的
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
13.958 32.33 22.59 12.23 16.83 35 3 11 15 6
end data.
compute H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+(r4*n4)**2/n4)/(N*(N+1))-4*(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.
跪求compute H = (12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3)+(r4*n4)**2/n4)/(N*(N+1))-4*(N+1).这个公式错在哪里,谢谢,着急!!!! |
|