bug in glm (PR#397)

Prof Brian Ripley writes:

 > > Date: Mon, 10 Jan 2000 11:53:50 +0100 (MET)
 > > From: holzer@stat.math.ethz.ch
 > > 
 > > Dear R-team
 > > 
 > > As I didn't get any answer to my bug-report last week I have taken the
 > > effort and extracted a minimal data set from my data (see below) where the
 > > following bug occurs:
 > > 

 > > > glm(SKR.ein.aus ~ ., family = binomial, data = bugdata, na.action = na.omit)

 > > Error in names<-.default(*tmp*, value = ynames) : names attribute must be the 

 > same length as the vector
Indeed, thank you (we needed an example). Alter
 > 
 >     names(w) <- ynames
 >     
 > to
 >     names(w) <- ynames[good]
 > 
 > and
 > 
 >     names(fit$effects) <-
 > 	c(xxnames[seq(fit$rank)], rep("", nobs - fit$rank))
 > 
 > to 
 >     names(fit$effects) <-
 > 	c(xxnames[seq(fit$rank)], rep("", sum(good) - fit$rank))
 > 
 > as the vector w is dropping observations with fitted probs 0 or 1.

This is actually a solution I thought of as well. However it has the
disadvantage that plot (and other functions?) gets problems:

fit.tst <- glm(SKR.ein.aus ~ ., family = binomial, data = bugdata)

Warning messages: 
1: fitted probabilities of 0 or 1 occurred in: (if (is.empty.model(mt))
glm.fit.null else glm.fit)(x = X, y = Y, 
.
.
.
7: Algorithm did not converge in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y,  

plot(fit.tst)