[R] naive variance in GEE

[Qiong Yang]

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"))