[Rd] Standardized Pearson residuals

Brett Presnell presnell at stat.ufl.edu
Mon Mar 14 23:02:40 CET 2011

My apologies.  I guess it would help if I tried the code more than once
before posting.  That should have been:

rstandard.glm <-
           infl=influence(model, do.coef=FALSE),
           type=c("deviance", "pearson"), ...)
  type <- match.arg(type)
  res <- switch(type, pearson = infl$pear.res, infl$dev.res)
  res <- res/sqrt(1-infl$hat)
  res[is.infinite(res)] <- NaN

More information about the R-devel mailing list