Re: [R] Error in F test on version 2.1.0

From: Ronaldo Reis-Jr. <chrysopa_at_gmail.com>
Date: Tue 10 May 2005 - 00:21:00 EST

Em Seg 09 Mai 2005 10:53, Uwe Ligges escreveu:
> 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).
>
> Uwe Ligges

> tot <- c(16,20,23,18,10,18,25,20,24,24,29,12,23,23,21,23,11)
> trat <- as.factor(rep(c("A","B"),c(8,9)))
> dep <- c(0,1,1,0,1,0,1,1,15,13,12,4,12,10,9,10,5)
> tot <- c(16,20,23,18,10,18,25,20,24,24,29,12,23,23,21,23,11)
> m <- glm(dep/tot~trat,family=binomial,weights=tot)
> anova(m,test="F")

Analysis of Deviance Table

Model: binomial, link: logit

Response: dep/tot

Terms added sequentially (first to last)

     Df Deviance Resid. Df Resid. Dev      F Pr(>F)
NULL                    16    105.685              
trat  1   96.117        15      9.568 96.117       
Mensagem de aviso:
NaNs produzidos in: pf(q, df1, df2, lower.tail, log.p)

INte
Ronaldo

-- 
It may or may not be worthwhile, but it still has to be done.
--
|>   // | \\   [***********************************]
|   ( §   § )  [Ronaldo Reis J˙nior                ]
|>      V      [UFV/DBA-Entomologia                ]
|    /     \   [36570-000 Višosa - MG              ]
|>  /(.''`.)\  [Fone: 31-3899-4007                 ]
|  /(: :'  :)\ [chrysopa@insecta.ufv.br            ]
|>/ (`. `'` ) \[ICQ#: 5692561 | LinuxUser#: 205366 ]
|    ( `-  )   [***********************************]
|>>  _/   \_Powered by GNU/Debian Woody/Sarge

______________________________________________
R-help@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
Received on Tue May 10 00:31:48 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:40 EST