|
data list free/Hc r1 r2 r3 r4 N n1 n2 n3 n4.
begin data
20.751 222.34 245.41 273.19 302.56 521 129 131 130 131
end data.
compute H=(12*((r1*n1)**2/n1+(r2*n2)**2/n2+(r3*n3)**2/n3+(r4*n4)**2/n4))-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.
还有想请问下因为之前不知道这个编程,曾经自己手动做过两组之间的比较,数据跟这个结果是不同的,这个差异?
这个方法对各组数据的分布情况有要求吗? |
|