%include catdata(berkeley); title h=1.8 'Berkeley Admissions Data'; goptions vsize=6.2 in hsize=7 in htext=1.6; data berkeley; set berkeley; dept1AG = (gender='F') * (dept=1); format dept dept.; proc catmod order=data data=berkeley; format dept dept. ; population dept gender; weight freq; response / out=predict; model admit = dept / ml noiter noprofile title="Model (AD, DG)"; run; title 'Model: logit(Admit) = Dept' a=-90 'Probability (Admitted)' h=3 a=-90 ' '; legend1 position=(bottom inside left) offset=(4,3) mode=share cborder=blue across=1 shape=symbol(6,1.5) label=('Gender') value=(c=black 'Female' c=red 'Male'); axis1 order=(-3 to 2) offset=(4) label=(a=90 'Log Odds (Admitted)'); axis2 label=('Department') offset=(4); %pscale(lo=-4, hi=3, anno=pscale, prob=%str(0.05,.1,.25,.5,.75,.9)); %catplot(data=predict, class=gender, x=dept, z=1.96, type=FUNCTION, legend=legend1, anno=pscale); *gskip; proc catmod order=data data=berkeley; format dept dept. ; weight freq; response / out=predict; model admit = dept gender / ml noiter noprofile title="Model (AD, AG, DG)" ; run; %pscale(lo=-4, hi=3, anno=pscale, prob=%str(0.05,.1,.25,.5,.75,.9)); title 'Model: logit(Admit) = Dept Gender' a=-90 'Probability (Admitted)' h=3 a=-90 ' '; legend1 position=(bottom inside left) offset=(4,3) mode=share cborder=blue across=1 shape=symbol(6,1.5) label=('Gender') value=(c=black 'Female' c=red 'Male'); axis1 order=(-3 to 2) offset=(4) label=(a=90 'Log Odds (Admitted)'); axis2 label=('Department') offset=(4); symbol1 i=none v=circle h=1.8 c=black; symbol2 i=none v=dot h=1.8 c=red ; %catplot(data=predict, class=gender, x=dept, type=FUNCTION, z=1.96, anno=pscale, legend=legend1);