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
--
You look tired.
--
|> // | \\ [***********************************]
| ( ?? ?? ) [Ronaldo Reis J??nior ]
|> V [UFV/DBA-Entomologia ]
| / \ [36570-000 Vi??osa - MG ]
|> /(.''`.)\ [Fone: 31-3899-4007 ]
| /(: :' :)\ [chrysopa at insecta.ufv.br ]
|>/ (`. `'` ) \[ICQ#: 5692561 | LinuxUser#: 205366 ]
| ( `- ) [***********************************]
|>> _/ \_Powered by GNU/Debian Woody/Sarge
On Mon, 9 May 2005 10:34:21 -0300 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.I also stumbled over this, the other day. The source of the problem was already discussed in this thread: https://stat.ethz.ch/pipermail/r-devel/2005-April/032923.html Z> Inte > Ronaldo > -- > You look tired. > -- > |> // | \\ [***********************************] > | ( ?? ?? ) [Ronaldo Reis J??nior ] > |> V [UFV/DBA-Entomologia ] > | / \ [36570-000 Vi??osa - MG ] > |> /(.''`.)\ [Fone: 31-3899-4007 ] > | /(: :' :)\ [chrysopa at insecta.ufv.br ] > |>/ (`. `'` ) \[ICQ#: 5692561 | LinuxUser#: 205366 ] > | ( `- ) [***********************************] > |>> _/ \_Powered by GNU/Debian Woody/Sarge > > ______________________________________________ > 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 >
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 >
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 > RonaldoHave 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). Uwe Ligges
On Mon, 9 May 2005, Achim Zeileis wrote:> 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.htmlHmm, 1) It is not documented that df=Inf is valid for pf(), so the problem is rather using an "F" test when it is inappropriate. If you specify an F test, you almost certainly meant a quasibinomial model. 2) This does work in the current R-patched. -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595