Draw plot using fitted probability in proc glimmix, lsmeans

BrandaLee — Wed, 20 Nov 2019 23:39:01 GMT

Dear all,

I have a GLMM model and my code runs smoothly:

proc glimmix data=production_change;
class Subject Group (ref='C');
model Response = Group / dist=binomial link=logit cl;
random intercept / subject = Subject;
lsmeans Group / adjust=Tukey plots=meanplot ilink cl;
run;

The only issues is with the "lsmeans / plots=" statement. I want to draw the graph using the fitted probability, not the estimated means (which SAS currently does for me). That is, I want to use the data in the red frame below to draw a graph, not using the "estimate" and its CL.

Anyone knows how to do that? Or do I need an ods statement, or a proc sgplot statement? I am quite new to this.

I appreciate your time and effort to help!

Re: Draw plot using fitted probability in proc glimmix, lsmeans

Rick_SAS — Fri, 22 Nov 2019 14:52:03 GMT

use the ILINK suboption to the PLOTS=MEANPLOT option:

lsmeans Group / adjust=Tukey plots=meanplot(ilink) ilink cl;

topic Draw plot using fitted probability in proc glimmix, lsmeans in Statistical Procedures

Draw plot using fitted probability in proc glimmix, lsmeans

Re: Draw plot using fitted probability in proc glimmix, lsmeans