Okay! That worked for the logistic. Do you know how to get that to work with a multinomial, same predictor? The outcome is 0,1,2. I want the reference group to be 1. So I want the LSMEANS differences table to have the same odds ratios described in the original question, with the outcome being 0 v 1 and 2 v 1. PROC GEE DATA=test DESCENDING ; CLASS id testcat(desc ref='0') fiscal_yr; MODEL ldl_lt100(REF='1')=testcat fiscal_yr /DIST=MULTINOMIAL LINK=GLOGIT TYPE3; REPEATED SUB=id/WITHIN=fiscal_yr; LSMEANS testcat/ILINK ODDSRATIO DIFF CL; RUN; Currently, it is showing: ldl_lt100 testcat _testcat 2 2 1 2 2 0 2 1 0 1 2 1 1 2 0 1 1 0 I want: ldl_lt100 testcat _testcat 2 2 1 2 2 0 2 1 0 0 2 1 0 2 0 0 1 0
... View more