[R] gee

Nirmala Ravishankar ravishan at fas.harvard.edu
Sun May 11 21:51:27 CEST 2003

I am trying to use gee() to calculate the robust standard errors for a 
logit model.  My dataset (zol) has 195019 observations; winner, racebl, 
raceas, racehi are all binary variables.   ID is saved as a vector of 
length 195019 with alternating 0's and 1's.   I get the following error 
message.  I also tried the same command with corstr set to "independence" 
and got the same error message.

> ID <- as.vector(array(0, nrow(zol)))
> k <- seq(2, nrow(zol), 2)
> ID[k] <- 1

> fm <- gee(winner ~ racebl + racehi + raceas, id = ID, data = zol, family 
= binomial(logit), corstr = "exchangeable")
[1] "Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27"
[1] "running glm to get initial regression estimate"
[1]  0.4308219 -0.1929547 -0.1741733 -0.1925523
Error in rep(0, maxclsz * maxclsz) : invalid number of copies in "rep"
In addition: Warning message: 
NAs produced by integer overflow in: maxclsz * maxclsz 

What am I doing wrong?

Nirmala Ravishankar

More information about the R-help mailing list