[R] what does the it when there is a zero events in the Logistic Regression with glm?

S. Sun shuguang at gmail.com
Tue Nov 22 04:31:32 CET 2005


Dear all,

I have a question about the glm. When the events of an observation is 0, 
the logit function on it is Inf. I wonder how the glm solve it.

An example:
Treat Events Trials
A     0      50
B     7      50
C     10     50
D     15     50
E     17     50

Program:

treat <- factor(c("A", "B", "C", "D", "E"))
events <- c(0, 7, 10, 15, 17)
trials <- rep(50, 5)
glm(cbind(events, trials-events)~treat, family=binomial)

What's wrong with it? And are there better ideas?

--
Sh. Sun




More information about the R-help mailing list