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
[EMAIL PROTECTED]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
