Displaying 2 results from an estimated 2 matches for "maxclsz".
Did you mean:
maxcls
2003 May 11
2
gee
...cehi + 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
2010 Jun 17
0
Modifyiing R working matrix within "gee" source code
...scale.fix <-
mm$scale.value <- mm$id <- NULL
mm[[1]] <- as.name("glm")
beta <- eval(mm, parent.frame())$coef
print(beta)
beta <- as.numeric(beta)
}
if (length(id) != length(y))
stop("Id and y not same length")
maxclsz <- as.integer(max(unlist(lapply(split(id, id), "length"))))
maxiter <- as.integer(maxiter)
silent <- as.integer(silent)
if (length(offset) <= 1)
offset <- rep(0, length(y))
if (length(offset) != length(y))
stop("offset and y not same len...