|
修改后可用程序如下:- data list free/Hc r1 r2 r3 r4 r5 N n1 n2 n3 n4 n5.
- begin data
- 120.269 145.16 233.98 248.56 257.31 367.49 500 100 100 100 100 100
- 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))/(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 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 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 x45=(r4-r5)**2/((N*(N+1)/12)*(1/n4+1/n5)*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 p23=1-cdf.chisq(x23,2).
- compute p24=1-cdf.chisq(x24,2).
- compute p25=1-cdf.chisq(x25,2).
- compute p34=1-cdf.chisq(x34,2).
- compute p35=1-cdf.chisq(x35,2).
- compute p45=1-cdf.chisq(x45,2).
- execute.
复制代码 |
|