Displaying 1 result from an estimated 1 matches for "gm1sl".
Did you mean:
gm1s
2009 Nov 20
1
different results across versions for glmer/lmer with the quasi-poisson or quasi-binomial families: the lattest version might not be accurate...
...L <- lmer(toto ~ period + (1 | herd), family = binomial)
summary(gm2sL)
#now poisson with similar data:
#essai avec jeu de données simulé:
set.seed(2)
toto<-rpois(1000,exp(period+herd.effect[herd]))
gm1s <- lmer(toto ~ period + (1 | herd), family = quasipoisson)
summary(gm1s)
gm1sL <- lmer(toto ~ period + (1 | herd), family = poisson)
summary(gm1sL)
########## here is the quasi-poisson result under R2.5.1 and Package lme4 version 0.99875-9: (all is OK)
Generalized linear mixed model fit using Laplace
Formula: toto ~ period + (1 | herd)
Family: quasipoisson(log...