"factor.residuals" <-
function(r, f) {
 if(is.matrix(f)) {
   rstar <- r - factor.model(f)} else {
   f <- f$loadings
    rstar <- r - factor.model(f)}
   return(rstar)}