11/5/2014

data death;
input vict $ def $ yes no;
tol=yes+no;
datalines;
w w 53 414
w b 11 37
b w 0 16
b b 4 139
;
run;

proc print data=death;
run;

proc genmod data=death;
class vict def;
model yes/tol=vict def/ dist=bin link=logit;
run;

proc logistic data=death;
class vict def;
model yes/tol=vict def/lackfit;
run;

proc logistic data=death;
class vict def/param=ref;
model yes/tol=vict def/lackfit;
run;


data aids;
input race $ azt $ y n;
total=y+n;
datalines;
white yes 14 93
white no 32 81
black yes 11 52
black no 12 43
;
*Asian yes 6 29
*Asian no 7 33;
run;

proc print data=aids;
run;

proc genmod data=aids;
class race azt;
model y/total=race azt/ dist=bin link=logit;
run;

proc genmod data=aids;
class race azt;
model y/total=race azt race*azt/ dist=bin link=logit;
run;


data one;
set aids;
outcome="yes"; count=y; output;
outcome="no"; count=n; output;
drop y n total;
run;

proc print data=one;
run;

proc freq data=one;
weight count;
tables race*azt*outcome /chm1;
run;

proc logistic data=aids;
class race (ref="white") azt(ref="no") / param=ref;
model y/total =azt race;
run;


proc logistic data=aids;
class race (ref="black") azt(ref="no") / param=ref;
model y/total =azt race;

run;