data(ohio)
fit.ex <- geeglm(resp ~ age + smoke + age:smoke, id=id, data=ohio,
family=binomial, corstr="exch", scale.fix=TRUE)
QIC(fit.ex)
fit.ar <- geeglm(resp ~ age + smoke + age:smoke, id=id, data=ohio,
family=binomial, corstr="ar1", scale.fix=TRUE)
QIC(fit.ex)
Run the code above in your browser using DataLab