[R] glm binomial with no successes

From: juli pausas <pausas_at_gmail.com>
Date: Wed, 27 Feb 2008 10:49:03 +0100

Dear all,
I have a question on glm, family binomial. I do not see significant differences between the levels of a factor (treatment) if all data for a level is 0; and replacing a 0 for a 1 (in fact reducing the
difference), then I detect the significant difference that I expected. Is there a way to overcome this problem? or this is an expected behaviour ? Here is an example:

s <- c(2,4,4,5,0,0,0,0)
f <- c(31,28,28,28,32,37,34,35)
tr <- gl(2, 4)
sf <- cbind(s,f) # numbers of successes and failures summary(glm(sf ~ tr, family="binomial")) # tr ns

sf[8,1] <- 1
summary(glm(sf ~ tr, family="binomial")) # tr significative **

Thanks for any suggestion



R-help_at_r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Wed 27 Feb 2008 - 10:09:34 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 27 Feb 2008 - 12:30:16 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive