From presnell at stat.ufl.edu  Mon Mar 14 23:02:40 2011
From: presnell at stat.ufl.edu (Brett Presnell)
Date: Mon, 14 Mar 2011 18:02:40 -0400
Subject: [Rd] Standardized Pearson residuals
Message-ID: <87lj0hs6m7.fsf@stat.ufl.edu>


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

rstandard.glm <-
  function(model,
           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
  res
}



[Rd] Standardized Pearson residuals
Brett Presnell