FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;

FILENAME url URL "http://192.38.117.59/~linearpredictors/datafiles/readSurgery.sas";
%include url;

DATA surgery;
SET surgery;
RUN;

PROC GENMOD DESCENDING DATA=surgery;
    CLASS surgtype;
    MODEL complication = surgtype / dist=binomial link=logit;
    estimate 'OR, abdominal vs. gynecological' surgtype 1 -1 0 / exp;
    estimate 'OR, orthopedic vs. gynecological' surgtype 0 -1 1 / exp;
    estimate 'OR, abdominal vs. orthopedic' surgtype 1 0 -1 / exp;
RUN;

/* For the Bonferroni/Sidak-correction, we adjust the significance level, alpha */

PROC GENMOD DESCENDING DATA=surgery;
    CLASS surgtype;
    MODEL complication = surgtype / dist=binomial link=logit;
    estimate 'OR, abdominal vs. gynecological' surgtype 1 -1 0 / exp alpha=0.0166666;
    estimate 'OR, orthopedic vs. gynecological' surgtype 0 -1 1 / exp alpha=0.0166666;
    estimate 'OR, abdominal vs. orthopedic' surgtype 1 0 -1 / exp alpha=0.0166666;
RUN;