What happen for Negative binomial link in Lmer fonction?
Dear R users, As I have no answer to my las reply i try again. If the binomal negative link has never been implemented in the lmer fonction, and as it was working when I did it some months ago, What kind of result did I get? This is the exemple of the answer of R when i was using lmer with negative binomal family: m1a<-lmer(atpos~ninter+saison+milieu*zone+(1|code),family=neg.bin(0.429),method="Laplace",data=manu) summary(m1a) Generalized linear mixed model fit using Laplace Formula: atpos ~ ninter + saison + milieu * zone + (1 | code) Data: manu Family: Negative Binomial(log link) AIC BIC logLik deviance 125.1 147.6 -54.57 109.1 Thank you for your answers, Best Regards Emmanuelle
E. Robardet wrote:
Thank you for your answers, I have an exemple of that i was using: m1a<-lmer(atpos~ninter+saison+milieu*zone+(1|code),family=neg.bin(0.429),method="Laplace",data=manu) summary(m1a) Generalized linear mixed model fit using Laplace Formula: atpos ~ ninter + saison + milieu * zone + (1 | code) Data: manu Family: Negative Binomial(log link) AIC BIC logLik deviance 125.1 147.6 -54.57 109.1 I think It was the version lme4 0.9975-10. Unfortunately, I have this version no more available on my computer.. I wonder if this old results are still valid.. Ben Bolker 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. 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-models at r-project.org ....
View this message in context: http://old.nabble.com/What-happen-for-Negative-binomial-link-in-Lmer-fonction--tp26013041p26439076.html Sent from the R help mailing list archive at Nabble.com.