У меня есть модель GLMM, и мой код proc glimmix работает гладко:
proc glimmix data=production;
class Subject Group (ref='Control');
model Response = Group / dist=binomial link=logit cl;
random intercept / subject = Subject;
lsmeans Group / adjust=Tukey plots=meanplot ilink cl;
run;
Для lsmeans / plots = meanplot SAS предоставляет график, используя оценку среднего значения и его CL,Однако мне нужен график, использующий подобранную вероятность и ее CL - то есть данные в красной рамке ниже:
Кто-нибудь знает, как это сделать? Это легко сделать в операторе "plots =", или мне нужно использовать другую процедуру proc? Большое спасибо !!!