####################  SÓLO ENS-PLUS

dfbetas<- function(fit, fit.s = summary(fit), fit.i = lm.influence(fit))
  {
    xxi <- diag(fit.s$cov.unscaled)
    si <- fit.i$sigma
        n.par<-length(coef(fit))
        n.datos<-fit$df+n.par
    bi <- matrix(rep(coef(fit), n.datos), byrow = T, ncol = n.par) - coef(fit.i)
    bi/(si %o% xxi^0.5)
 }