From tlumley@u.washington.edu  Sat Nov 17 18:40:55 2001
From: tlumley@u.washington.edu (Thomas Lumley)
Date: Sat, 17 Nov 2001 10:40:55 -0800 (PST)
Subject: [Rd] pearson residuals in glm for binomial response (PR#1175)
In-Reply-To: <200111161035.LAA27495@pubhealth.ku.dk>
Message-ID: <Pine.A41.4.33.0111171039050.137786-100000@homer34.u.washington.edu>

On Fri, 16 Nov 2001 J.T.Kent@leeds.ac.uk wrote:

> R version 1.3.0
> OS: SunOS 5.7, but I think the same problem occurs with Windows
>
> An incorrect formula seems to be used to calculate the pearson residuals
> for a generalized linear model with a binomial response.  Here is a
> simple program which gives (a) the pearson residuals calculated directly,
> (b) the pearson residuals from glm, and (c) the deviance residuals from
> glm.  The first and last columns are quite close.  The middle differs from
> the first by a factor of n (n=10 here).  This problem did not occur with
> an earlier version of R I tried, version 0.64.2.
>

Yes, it is a bug, resulting from an incorrect fix to an earlier bug that
had the sign wrong for decreasing link functions like the reciprocal.
It's already been reported and fixed.

	-thomas

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


[Rd] pearson residuals in glm for binomial response (PR#1175)
Thomas Lumley