Help using Base SAS procedures

GEE GLIMMIX

Reply
New Contributor
Posts: 4

GEE GLIMMIX

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

Respected Advisor
Posts: 2,655

Re: GEE GLIMMIX

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

New Contributor
Posts: 4

Re: GEE GLIMMIX

Posted in reply to SteveDenham

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

New Contributor
Posts: 4

Re: GEE GLIMMIX

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

Respected Advisor
Posts: 2,655

Re: GEE GLIMMIX

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

New Contributor
Posts: 4

Re: GEE GLIMMIX

Posted in reply to SteveDenham

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

Respected Advisor
Posts: 2,655

Re: GEE GLIMMIX

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

Ask a Question
Discussion stats
  • 6 replies
  • 443 views
  • 0 likes
  • 2 in conversation