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

Thomas Lumley tlumley@u.washington.edu
Sat, 17 Nov 2001 10:40:55 -0800 (PST)


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