有个GEE的问题需要请教一下:
用proc genmod做binary模型,param=GLM, 有三个自变量,都是categorical的。
最后生产的结果有type III GEE analysis, 同时我还用contrast去对某一变量的几个选项进行两两比较。
如果type III analysis的结果中,某一变量并不显著,是否可以说在控制其他变量的情况下,这个变量对因变量无影响。
Contrast的结果,有些变量之间的两两比较呈现显著,这个结果是基于其他变量在reference level的情况下吗?
contrast和type III analysis的结果相反,应该以哪个为准?我对sas不熟悉,之前以为这两个类似anova和Tukey's HSD。
Estimate是其他变量在reference level的结果吗?
代码如下:
proc genmod data=data DESCENDING;
class enrollment gender prim school Z/ PARAM=GLM ;
model enrollment =school school*gender gender prim prim*school prim*gender /dist=binomial link=logit type3;
contrast "A vs H" school 1 -1 0 0 0/WALD;
estimate "A" intercept 1 school 1 0 0 0 0;
repeated subject=Z/modelse type=cs;
run;
谢谢