SAS script: How to plot random intercepts in logistic glimmix models?

Reply
N/A
Posts: 1

SAS script: How to plot random intercepts in logistic glimmix models?

proc glimmix;

class A B;

model Y = X/ dist=binary link=logit ddfm=bw solution;

random intercept / subject=A type = un solution ; covtest/wald;

random intercept / subject=B type=un solution ; covtest/wald;

output out=glimmixout_3 pred( blup )=PredProb/* sem ilink é em logit*/

                    pred(noblup )=PredProb_PA;

run;

and then:

ods html;

proc sgplot data=glimmixout_3;

scatter y=PredProb_PA x=X/group=A name="aaa";by A;

keylegend "dados"/title="aaa";

run;

ods html close;

I have 13 subjects (random effect). I get linear intercepts but when I plot all random effects together (group=A in the abobe example) I did not get parallel lines. I expected to see the variance of the random effects but, instead, all the lines are together. Where is my mistake?

Ask a Question
Discussion stats
  • 0 replies
  • 146 views
  • 0 likes
  • 1 in conversation