Friday, May 26, 2006

無聊的SAS程式碼

options nodate nonotes ps=60;
data All;
input Z X Y count @@;
cards;
0 0 0 2
0 0 1 5
0 1 0 8
0 1 1 13
1 0 0 5
1 0 1 5
1 1 0 11
1 1 1 34
;
proc freq;
weight count;
tables X*Y/ nocol norow nopercent chisq cmh;
tables Z*Y/ nocol norow nopercent chisq cmh;
tables X*Z/ nocol norow nopercent chisq cmh;
tables X*Z*Y/ nocol norow nopercent chisq cmh;
run;

data All;
input Z X Y count @@;
cards;
0 0 0 2
0 0 1 5
0 1 0 8
0 1 1 13
1 0 0 5
1 0 1 5
1 1 0 11
1 1 1 34
;
proc catmod order=data;
weight count;
population X Z;
model Y=;
proc catmod order=data;
weight count;
population X Z ;
model Y=X;
proc catmod order=data;
weight count;
population X Z;
model Y=Z;
proc catmod order=data;
weight count;
population X Z ;
model Y=X Z;
proc catmod order=data;
weight count;
population X Z;
model Y=X|Z;
run;

No comments: