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
}
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
