On Thu, Oct 22, 2009 at 1:39 PM, Ben Bolker <bol...@ufl.edu> wrote: > ROBARDET Emmanuelle wrote: >> >> Dear R users, >> I'm performing some GLMMs analysis with a negative binomial link. >> I already performed such analysis some months ago with the lmer() function >> but when I tried it today I encountered this problem: >> Erreur dans famType(glmFit$family) : unknown GLM family: 'Negative >> Binomial' >> >> Does anyone know if the negative binomial family has been removed from >> this function? >> I really appreciate any response. >> Emmanuelle >> >> > > I would be extremely surprised if this worked in the past; to > the best of my knowledge the negative binomial family has > never been implemented in lmer.
I too would be extremely surprised if it had worked in the past, considering that I have never implemented it. I did exchange email with Bill Venables about it and we formulated what seems to be a reasonable approach but it hasn't made it to the top of the "To Do" list yet. Right now the big push is on code to profile the log-likelihood with respect to the parameters so we can actually get confidence intervals and, the holy grail of mixed-modeling, p-values. > One could in principle > do a glmmPQL fit with the negative binomial family > (with a fixed value of the overdispersion parameter). > glmmADMB is another option. > Can you say which version etc. you were using??? > > Follow-ups should probably be sent to r-sig-mixed-mod...@r-project.org .... > -- > View this message in context: > http://www.nabble.com/What-happen-for-Negative-binomial-link-in-Lmer-fonction--tp26013041p26015140.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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.