我在胡良平老师的《现代统计学与SAS应用》看到了类似的问题。我把他的程序改了一下。
DATA abc;
DO x=1,2,3,4;
DO y=1,2,3,4;
INPUT f @@;
OUTPUT;
END;
END;
CARDS;
70 22 4 2
27 24 9 3
16 23 13 7
9 20 15 14
;
PROC FREQ data=abc;
WEIGHT f;
TABLES y*x / CHISQ;
RUN;
PROC REG data=abc;
FREQ f;
MODEL y=x;
RUN;
quit;
DATA a;
r=4;
c=4;
kf=71.4325;
b=0.40534;
sb=0.04479;
dff=(r-1)*(c-1);
dfr=1;
dfb=dff-dfr;
vb=sb*sb;
kr=ROUND(b*b/vb, 0.001);
kb=kf-kr;
pr=1-PROBCHI(kr,dfr);
pr=ROUND(pr, 0.0001);
IF pr<0.0001 THEN pr=0.0001;
pb=1-PROBCHI(kb,dfb);
pb=ROUND(pb, 0.0001);
IF pb<0.0001 THEN pb=0.0001;
FILE PRINT;
PUT #2 @10 'CHISQ-R' @20 'P-R'
#2 @30 'CHISQ-B' @40 'P-B'
#3 @10 kr @20 pr
#3 @30 kb @40 pb;
RUN;


雷达卡



京公网安备 11010802022788号







