[R] What happen for Negative binomial link in Lmer fonction?
Emmanuelle Robardet
e.robardet at AFSSA.FR
Fri Nov 20 12:24:27 CET 2009
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.
More information about the R-help
mailing list