请问大家使用我使用SAS 的glimmix模块做apc模型，跑了模型之后，想做出附件1中的图形，请问大家应该怎么写命令？？？？应该先在命令中预测每个个体发生y的概率，然后再用画图命令分别做年龄、时期和世代与y预测概率的相关图吗？？？
我自己的命令如下：
        proc glimmix data=work.import maxopt=10000;
        class period cohort;
        model Y(event='1') = age x1 x2 x3 / solution cl dist=binary;
        random period cohort / solution; covtest glm/wald;
        run;




有无大神求教！！！！！