我想根据以下模型绘制工作性别差异的预测年龄(包括两个年龄变化:age和age2)、时期和队列变化,我应该如何添加SAS代码?
proc glimmix data=model noitprint noclprint;
class period cohort job;
model job =age age2 sex syz edu_2 edu_3 party hukou mar diqu_2 diqu_3 diqu_4 fedu_2 fedu_3
/ solution CL
dist = mult link = glogit DDFM=BW ;
random intercept sex/sub=period group=job;
random intercept sex/sub=cohort group=job;
covtest GLM / WALD;
NLOPTIONSTECHNIQUE=NRRIDG;
run;
谢谢!