Displaying 1 result from an estimated 1 matches for "yglm".
Did you mean:
glm
2001 Oct 10
2
Pearson residuals (PR#1123)
...ponding diagonal element of the weight matrix W evaluated at the last step
of the iterative model fitting procedure (IWLS), instead of dividing by the
square root of the model variance. Here's an example of the weird effects it
produces:
> me <- c(1:100)/2
> y <- rpois(100,me)
> yglm <- glm(y~me-1, family=poisson(link="identity"))
> plot(fitted(yglm), residuals(yglm,type="response"))
plots the raw residuals... but
> plot(fitted(yglm), residuals(yglm,type="pearson"))
plots the raw residuals MULTIPLIED by sqrt(fitted(yglm))... because...