[R] Error in F test on version 2.1.0

[Achim Zeileis]

On Mon, 09 May 2005 15:53:11 +0200 Uwe Ligges wrote:

> Ronaldo Reis-Jr. wrote:
> 
> > Hi,
> > 
> > I make a upgrade to R 2.1.0 and in some analysis I give an error:
> > 
> > anova(model,test="F")
> > Analysis of Deviance Table
> > 
> > Model: binomial, link: logit
> > 
> > Response: landing/total
> > 
> > Terms added sequentially (first to last)
> > 
> > 
> >     Df Deviance Resid. Df Resid. Dev    F Pr(>F)
> > NULL                    16    105.079            
> > trat  1  93.149        15    11.930 93.15      
> > Warning message:
> > NaNs produced in: pf(q, df1, df2, lower.tail, log.p) 
> > 
> > In old version it work.
> > 
> > Inte
> > Ronaldo
> 
> Have you read the posting guide?
> It asks you to specify a reproducible example ........
> Please so so (both, read the psoting guide and specify the example).

Valid point, try:

set.seed(1071)
y <- factor(rnorm(10) > 0)
x <- rnorm(10)
fm1 <- glm(y ~ 1, family = binomial)
fmx <- glm(y ~ x, family = binomial)
anova(fm1, fmx, test = "F")

and as I said in another reply already. This boils down to the problems
in pf() reported by Gordon Smyth:
  https://stat.ethz.ch/pipermail/r-devel/2005-April/032923.html
Z

> Uwe Ligges
> 
> ______________________________________________
> 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
>