[R] dfbetas vs dfbeta

Peter Dalgaard p.dalgaard at biostat.ku.dk
Fri Nov 13 10:28:34 CET 2009


alphaace wrote:
> Hi, I've looked around but can't find a clear answer to the difference for
> these two? Any help?
> 
> Thanks!

One is a scaled version of the other, at least for lm objects:

 > stats:::dfbetas.lm
function (model, infl = lm.influence(model, do.coef = TRUE),
     ...)
{
     xxi <- chol2inv(model$qr$qr, model$qr$rank)
     dfbeta(model, infl)/outer(infl$sigma, sqrt(diag(xxi)))
}

so it is essentially dividing by s.e.(beta), except that it uses 
infl$sigma which AFAICS is the leave-one-out residual sd.

-- 
    O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
   c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
  (*) \(*) -- University of Copenhagen   Denmark      Ph:  (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)              FAX: (+45) 35327907




More information about the R-help mailing list