Draw plot using fitted probability in proc glimmix, lsmeans

[BrandaLee]

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!

[Rick_SAS]

use the ILINK suboption to the PLOTS=MEANPLOT option:

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