[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?
Thanks,
Nirmala Ravishankar
More information about the R-help
mailing list