[R] re fit with binomial model (lme4)

[Douglas Bates]

Thanks for sending that, Ben.  I have incorporated a slightly modified
version, with attribution, in the lme4 sources.

On Tue, Apr 28, 2009 at 8:03 AM, Ben Bolker <bolker at ufl.edu> wrote:

> I think the following works, but use with caution -- I just
> wrote it last night.
>
> setMethod("refit", signature(object = "mer", newresp = "matrix"),
>          function(object, newresp, ...)
>      {
>        ## newresp <- as.double(newresp[!is.na(newresp)])
>        wts <- rowSums(newresp)
>        newresp <- newresp[,1]/wts
>        stopifnot(length(newresp) == object at dims["n"])
>        object at y <- newresp
>        object at pWt <- wts
>        lme4:::mer_finalize(object)
>      })