Re: [Rd] Standardized Pearson residuals

From: Brett Presnell <presnell_at_stat.ufl.edu>
Date: Mon, 14 Mar 2011 18:02:40 -0400

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
}

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Mon 14 Mar 2011 - 22:15:43 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 15 Mar 2011 - 03:50:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive