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

Prof Brian Ripley ripley at stats.ox.ac.uk
Sun Aug 6 11:57:44 CEST 2006


On Mon, 24 Jul 2006, Achim Zeileis wrote:

> 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.

Well, the original S code returns numeric(0): surely that is worse.

There is not enough information in the fit to recover the response 
residuals in general (the problem I see is with zero weights), so all we 
can do it throw an informative error.

> 
> Best,
> Z
> 
> ## 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 at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
> 
> 

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595



More information about the R-devel mailing list