Hi,

The standard error from logistic regression is slightly different from the naive SE from GEE under independence working correlation structure.

Shouldn't they be identical? Anyone has insight about this?

Thanks,
Qiong

a<-rbinom(1000,1)
b<-rbinom(1000,2,0.1)
c<-rbinom(1000,10,0.5)
summary(gee(a~b, id=c,family="binomial",corstr="independence"))$coef
summary(glm(a~b,family="binomial"))

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to