PROC GENMOD DATA=donner; CLASS Status Sex; MODEL Status = Age|Sex / DIST=BINOMIAL; ESTIMATE 'Female 20' Intercept 1 Sex 1 Age 20 Sex*Age 20; RUN; PROC LOGISTIC DATA=donner; CLASS Status Sex; MODEL Status = Age|Sex; CONTRAST 'Female 20' Intercept 1 Sex 1 Age 20 Sex*Age 20 / ESTIMATE=PROB; RUN;