我有一个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过程?非常感谢!!!