Re: [R] naive variance in GEE

2008-09-09 Thread Prof Brian Ripley
On Mon, 8 Sep 2008, Qiong Yang wrote: 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? They are computed quantities from iterations with

Re: [R] naive variance in GEE

2008-09-09 Thread Thomas Lumley
On Mon, 8 Sep 2008, Qiong Yang wrote: Hi, The standard error from logistic regression is slightly different from the naive SE from GEE under independence working correlation structure. Yes Shouldn't they be identical? Anyone has insight about this? No, they shouldn't. They are different

Re: [R] naive variance in GEE

2008-09-09 Thread Thomas Lumley
Sorry, I misread your message. Prof Ripley is right, as usual -- the estimates use different stopping criteria and so are just numerically different. -thomas On Tue, 9 Sep 2008, Thomas Lumley wrote: On Mon, 8 Sep 2008, Qiong Yang wrote: Hi, The standard error from logistic

[R] naive variance in GEE

2008-09-08 Thread 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)