andre tavares correa dias
2006-Sep-25 14:42 UTC
[R] F values for glm with binomial distribution
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