"fitted.statfacmodBurSt" <- function (object) { fun.copyright <- "Placed in the public domain in 2006 by Burns Statistics" fun.version <- "fitted.statfacmodBurSt 002" ans <- object$loadings %*% t(object$loadings) diag(ans) <- 1 ans <- t(object$sdev * ans) * object$sdev ans }