proc freq data=dent; table trt*resp / out=os; run; proc transpose data=os out=tran; by trt; var count; run; data a; set tran; const=0; c1=log((sum(of col1-col1)+const)/(sum(of col2-col5)+const)); c2=log((sum(of col1-col2)+const)/(sum(of col3-col5)+const)); c3=log((sum(of col1-col3)+const)/(sum(of col4-col5)+const)); c4=log((sum(of col1-col4)+const)/(sum(of col5-col5)+const)); run; proc sgplot; series y=c1 x=trt; series y=c2 x=trt; series y=c3 x=trt; series y=c4 x=trt; yaxis label="Empirical Logits"; run;