Dennis Murphy <djmuser <at> gmail.com> writes:
>
> Hi:
>
> Did you try putting the offset in the model formula, as in
>
> bigglm( y ~ offset(z) + x, ...) ?
>
> I haven't tried bigglm() personally (BTW, it's in the biglm package, which
> wasn't mentioned), but this syntax works in the standard glm() function, so
> perhaps it maps to bigglm() as well... (?)
>
> HTH,
> Dennis
A quick test suggests that bigglm does *not* respect offset
parameters in the formula, and a quick check in bigglm.function suggests
likewise:
...
tt <- terms(formula)
beta <- start
etafun <- function(x) if (is.null(beta))
rep(0, nrow(x))
else x %*% beta
...
farther down, etafun is called with the results of
model.matrix(). It should be fairly 'straightforward'
(ha!) to incorporate an offset term in this, but it might
be better to ask the package maintainer to do it rather
than to have random people like me doing the hacking ...
Ben Bolker
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
