二分类资料重复测量,分别采用glimmix和genmod程序,但结果中的“visit*group”最小二乘均值,两个程序结果相差较大,genmod程序结果更合理,接近因变量event='1'的发生率,而glimmix的结果与因变量event='1'的发生率差别较大。为什么会出现两者LSMEANS差异较大? 是否更应该采用genmod的结果?
proc glimmix data=scale ;
class Visit group (ref='2') ID ;
model HAMA_zy(event='1')=Visit group Visit*group /dist=binary link=logit ddfm=bw solution ;
random int /subject=ID type=chol;
nloptions tech=nrridg;
lsmeans Visit*group / ilink cl oddsratio diff slice=Visit;
run;
proc genmod data=scale descending;
class Visit group (ref='2') ID;
model HAMA_zy=Visit group Visit*group /dist=binomial link=logit ;
repeated subject=ID / type=exch;
lsmeans Visit*group / ilink cl diff oddsratio;
run;