[Rd] residuals.glm for glm(..., y = FALSE)

From: Achim Zeileis <Achim.Zeileis_at_wu-wien.ac.at>
Date: Mon 24 Jul 2006 - 19:16:45 GMT

I just ran into a problem with residuals.glm() for GLMs fitted with   glm(..., y = FALSE)
For such fits, the residuals of type "deviance", "pearson" and "response" can't be computed (see example below). The reason is that residuals.glm() uses
  y <- object$y
which obviously can't work for such fits.

I've checked ?glm and ?residuals.glm and did not find it mentioned explicitely...looks like a bug to me.


## example for poisson GLM from ?glm

  d.AD <- data.frame(treatment = gl(3,3), outcome = gl(3,1,9),     counts = c(18,17,15,20,10,20,25,13,12))   glm.D93 <- glm(counts ~ outcome + treatment, family = poisson,     data = d.AD, y = FALSE)

## residuals cannot be computed

  residuals(glm.D93, type = "deviance")
  residuals(glm.D93, type = "pearson")
  residuals(glm.D93, type = "response")

R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue Jul 25 05:19:10 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sun 06 Aug 2006 - 10:27:06 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.