%include catdata(arrests); *-- main effects model; proc logistic data=arrests; class colour year sex citizen employed; model released = colour year age sex employed citizen checks; run; *-- all two-way model, with backward elimination; proc logistic data=arrests; class colour year sex citizen employed; model released = colour | year | age | sex | employed | citizen | checks@2 / selection=backward; run;