[R] why does lmer not give p valuses for quasibinomial family?

[Spencer Graves]

	  Please provide a replicatable example.  Consider the following:

 > library(lme4)
 > library(mlmRev)
 > fit1q <- lmer(use~1+(1|district), data=Contraception, 
family=quasibinomial)
 > fit2q <- lmer(use~age+(1|district), data=Contraception, 
family=quasibinomial)
 > anova(fit1q, fit2q)
Data: Contraception
Models:
fit1q: use ~ 1 + (1 | district)
fit2q: use ~ age + (1 | district)
       Df     AIC     BIC  logLik  Chisq Chi Df Pr(>Chisq)
fit1q  2  2538.5  2549.6 -1267.2
fit2q  3  2537.9  2554.7 -1266.0 2.5254      1     0.1120
 >
	  hope this helps,
	  spencer graves

Szentirmai Istvan wrote:

> Dear All,
> 
> I'm running a binomial model using the lmer() function, and I get p values 
> for the parameter estimates only with family=binomial, but not with 
> quasibinomial? Why is that so? I wanted to use quasibinomial family, because 
> my data were overdispersed.
> 
> Thanks,
> Istvan
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html