data aids; input race $ azt $ yes no @@; cases = yes + no; datalines; white y 14 93 white n 32 81 black y 11 52 black n 12 43 ; proc genmod order=data; class azt race; model yes/cases = azt race / dist=bin link=logit obstats type3; run; ***************************************************************; data aids_freq; input race $ azt $ symptom $ count @@; datalines; White Imm Yes 14 White Imm No 93 White Del Yes 32 White Del No 81 Black Imm Yes 11 Black Imm No 52 Black Del Yes 12 Black Del No 43 ; proc freq order=data; weight count; tables race*azt*symptom/cmh; run;