[R] F values for glm with binomial distribution

[andre tavares correa dias]

Hi Rneters,
I'm running a GLM model with a full factorial design in blocks and binomial 
error distribution. I would like to have the F values for this model but I 
got a message that "using F test with a binomial family is inappropriate in: 
anova.glm(model, test = "F")". Should I not report F statistics on this kind 
of analysis? I would appreciate any comment on this.

This is my output:

> model=glm(y~agua*fert+bloco, family=binomial)
> anova.glm(model,test="F")

Analysis of Deviance Table

Model: binomial, link: logit

Response: y

Terms added sequentially (first to last)


          Df Deviance Resid. Df Resid. Dev      F  Pr(>F)  
NULL                         43     85.018                 
agua       1    0.908        42     84.110 0.9081 0.34063  
fert       1    5.673        41     78.437 5.6727 0.01723 *
bloco      1    0.044        40     78.393 0.0444 0.83320  
agua:fert  1    0.616        39     77.777 0.6160 0.43253  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
Warning message:
using F test with a binomial family is inappropriate in: anova.glm(model, 
test = "F")

Best regards,
André





--
André Tavares Corrêa Dias
Laboratório de Ecologia Vegetal
Universidade Federal do Rio de Janeiro
CCS-IB-Departamento de Ecologia
Caixa Postal 68020
21941-970 Rio de Janeiro – RJ, Brazil
tel: +55 21 25626377
Fax: + 55 21 25626320
atcdias at biologia.ufrj.br