generalized mixed model + mcmcsamp

Ben Bolker · 2009-02-11T13:20:15Z

Thomas Mang fiwi.at> writes: > > Hi, > > I have fitted a generalized linear mixed effects model using lmer > (library lme4), and the family = quasibinomial. I have tried to obtain a > MCMC sample, but on calling mcmcsamp(model1, 1000) I get the following > error which I don't understand at all: > > Error in .local(object, n, verbose, ...) : Update not yet written > > traceback() delivers: > 4: .Call(mer_MCMCsamp, ans, object) > 3: .local(object, n, verbose, ...) > 2: mc

Ben Bolker

Wed, Feb 11, 2009 5:20 AM

Thomas Mang <Thomas.Mang <at> fiwi.at> writes:

Bad news: mcmcsamp is not working at present even
for LMMs, has never worked (so far) for GLMMs, and
the author (Doug Bates) says he's not sure there's
a coherent way to deal with quasi- models in the
mcmcsamp framework (since it doesn't correspond to
any well-defined distribution).
   You may want to browse the archives of the
r-sig-mixed-models list, and/or request more
information there.

  Ben Bolker

generalized mixed model + mcmcsamp

Thread (2 messages)