Hi,

Sorry, I should have said the bigglm() is in the biglm package.

I have tried putting the offset in the model formula as you suggested but it doesn't make any difference. I was comparing the results of bigglm() with glm() on a small data set before I applied it to the much larger data set, this is when I noticed that bigglm() was ignoring the offset.

BW

Richard

Dennis Murphy said the following on 08/02/2011 19:55:
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

On Tue, Feb 8, 2011 at 6:16 AM, Richard Jacques <r.jacq...@sheffield.ac.uk <mailto:r.jacq...@sheffield.ac.uk>> wrote:

    Dear all,

    I have a large data set and would like to fit a logistic regression
    model using the bigglm function.  I need to include an offset in the
    model but when I do this the bigglm function seems to ignore it.

    For example, running the two models below produces the same model
    and the offset is ignored

    bigglm(y~x,offset=z,data=Test,family=binomial(link = "logit"))
    bigglm(y~x,data=Test,family=binomial(link = "logit"))

    Is it possible to fit an offset using bigglm? And if so, what am I
    doing wrong?

    Thanks,

    Richard

-- Dr Richard Jacques
    Medical Statistics Group
    School of Health and Related Research
    University of Sheffield
    Regent Court
    30 Regent Street
    SHEFFIELD
    S1 4DA

    ______________________________________________
    R-help@r-project.org <mailto: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.



______________________________________________
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.

Reply via email to