topic Re: GEE GLIMMIX in SAS Procedures

GEE GLIMMIX

BonittaG — Sat, 21 Feb 2015 06:09:34 GMT

Dear all,

i had done propensity score matching. i need use proc glimmix because permit to specify beta distribution for longitudinal model

How can account for matching variable in this proc ( the strata variable into conditional logistic)

*Beta GEE;

proc glimmix data=data_;

class id time trt;

model y = trt time time trt / dist=beta s ddfm=none covb;

were i account for strata (matching variable)??? :smileyshocked:

random _residual_ / subject=id type=cs vcorr;

lsmeans diab*time /ilink cl;

run;

Best

G Bonitta

Re: GEE GLIMMIX

SteveDenham — Mon, 23 Feb 2015 19:01:28 GMT

Some things to consider:

You need a name for the matching variable. For now, I will call it pairid.

There is a variable in the lsmeans statement called diab that is not identified in either the class or model statement. What does this variable represent?

Here is what might work:

proc glimmix data=data_;

class id pairid time trt;

model y = trt time trt*time/dist=beta s ddfm=none covb;

random intercept/subject=pairid;

random time/residual subject=id type= <> vcorr;

lsmeans trt*time/ilink cl;

run;

I did not specify the structure for the R side repeated in time variable. That will depend on the spacing on the 'time' variable and what assumptions can be made. Type=CS may involve too many assumptions to be a realistic covariance structure, although it is handy for fitting unequal spacing in time.

Steve Denham

Re: GEE GLIMMIX

BonittaG — Tue, 24 Feb 2015 06:21:46 GMT

Thank You very Much

you right pardon me my mistake diab is trt ....

I'm an R user and new into SAS enviroment..

Seeing your competence i ask agian a little question .. :

how can I plot mean residual vs decile of linear predictor about the specifided model we are discuss ?

Sorry for noise ...

Thank You again

Best

Bonitta

Re: GEE GLIMMIX

BonittaG — Tue, 24 Feb 2015 07:12:01 GMT

For residuals plot that i need maby i can use this from

SAS/STAT(R) 9.2 User's Guide, Second Edition

output statement correct?

proc glimmix;

stuff code

output out=igausout pred=p resid=r;

run;

Best

Bonitta

Thank You

Re: GEE GLIMMIX

SteveDenham — Tue, 24 Feb 2015 12:39:21 GMT

Yes, the OUTPUT statement will give predicted and residual values. From there to your request for plotting is a journey through the DATA step to calculate the decile cutpoints, then the mean residual within those, and finally SGPLOT to plot the two against one another.

Steve Denham

Re: GEE GLIMMIX

BonittaG — Wed, 25 Feb 2015 12:41:55 GMT

Thank you ....can you so gentle to write the code for do that ? about the model

for plot mean residuals vs linear predictor

Thank You so much

Sorry for noise

Best

Bonitta

Re: GEE GLIMMIX

SteveDenham — Thu, 26 Feb 2015 13:09:18 GMT

I am not a coder in that way, so my attempts would be pretty horrible. Google for examples of each of the steps I mentioned, and I am sure you can find much better code than I can develop (or have time to develop).

Steve Denham