dfbetas<-function(fit,lms=summary(fit),lmi=lm.influence(fit)){
xxi<-diag(lms$cov.unscaled)
si<-lmi$sigma
bi<-t(coef(fit)-t(coef(lmi)))
bi/(si %o% xxi^0.5)
}
dffits<-function(fit,lmi=lm.influence(fit)){
h<-lmi$hat
e<-residuals(fit)
si<-lmi$sigma
h^0.5*e/(si*(1-h))
}
stanres<-function(fit,lms=summary(fit),lmi=lm.influence(fit)){
h<-lmi$hat
e<-residuals(fit)
s<-lms$sigma
si<-lmi$sigma
e/(s*(1-h)^.5)
}
studres<-function(fit,lmi=lm.influence(fit)){
h<-lmi$hat
e<-residuals(fit)
si<-lmi$sigma
e/(si*(1-h)^.5)
}