[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