|
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
* 请于此行输入Kruskal-Wallis H 检验的结果(即Hc r1 r2 r3 r4 N n1 n2 n3 n4,尽量精确)
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. |
|