topic Re: PROC GEE LSMEANS reference group in Statistical Procedures

PROC GEE LSMEANS reference group

lauragrau — Mon, 14 Dec 2020 15:49:55 GMT

Hi, I'm running the following model:

proc gee data=data descending;
class id category(ref='Non-user') fiscal_yr ;
model outcome(event='1')=category fiscal_yr /dist=bin link=logit type3;
repeated sub=id/within=fiscal_yr;
lsmeans category/ilink oddsratio diff cl;
run;

Category is a numeric variable with a format:

0= 'Non-user'

1='Varied user'

2='User'

I have very specific odds ratios I want to get from the LSMEANS. I want:

Varied vs non-user

User vs non-user

User vs varied

I have tried by changing the reference group in the class statement, which changes the reference in my parameter estimates statement, but not in my "Differences of category Least Squares Means".

Any ideas?

Thank you!

Re: PROC GEE LSMEANS reference group

StatDave — Mon, 14 Dec 2020 16:28:48 GMT

Using REF= in the CLASS statement will set the reference level in both the parameter estimates and the lsmeans difference tables. For example, the statements for a variable, A, with levels 0, 1, and 2

class a(desc ref='0');

lsmeans a / diff;

will give lsmeans differences 1 vs 0 , 2 vs 0, and 2 vs 1.

Re: PROC GEE LSMEANS reference group

lauragrau — Mon, 14 Dec 2020 17:30:55 GMT

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

Re: PROC GEE LSMEANS reference group

StatDave — Mon, 14 Dec 2020 17:49:40 GMT

The code you show (with response variable option REF='1') gives what you want in SAS release SAS 9.4 TS1M6.

Re: PROC GEE LSMEANS reference group

lauragrau — Mon, 14 Dec 2020 17:55:04 GMT

Theoretically, yes! The weird/interesting thing thing is that it does show was I want in an unadjusted model, but as soon as I add other covariates, it doesn't.