inst/efas/functions/printLoadings.R

In PsyChiLin/EFAshiny: EFAshiny : An User-friendly Application for Exploratory Factor Analysis

printLoadings <- function (x, digits = 2, cutoff = 0, sort = T,decreasing = T, ...) {
        Lambda <- unclass(x)
        p <- nrow(Lambda)
        factors <- ncol(Lambda)
        if (sort) {
                #mx <- max.col(abs(Lambda))
                #ind <- cbind(1L:p, mx)
                #mx[abs(Lambda[ind]) < 0.5] <- factors + 1
                #Lambda <- Lambda[order(mx, 1L:p), ]
                L <- list()
                mx <- max.col(abs(Lambda))
                for (i in 1:factors){
                        Ltemp <- Lambda[mx == i,]
                        Ltemp <- Ltemp[order(Ltemp[,i],decreasing = decreasing),]
                        L[[i]] <- Ltemp
                }
                Lambda <- do.call(rbind,L)
        }
        cat("\nLoadings:\n")
        fx <- format(round(Lambda, digits))
        names(fx) <- NULL
        nc <- nchar(fx[1L], type = "c")
        fx[abs(Lambda) < cutoff] <- paste(rep(" ", nc), collapse = "")
        newx <- print(fx, quote = FALSE, ...) # I assigned this to a variable
        vx <- colSums(x^2)
        varex <- rbind(`SS loadings` = vx)
        if (is.null(attr(x, "covariance"))) {
                varex <- rbind(varex, `Proportion Var` = vx/p)
                if (factors > 1) 
                        varex <- rbind(varex, `Cumulative Var` = cumsum(vx/p))
        }
        cat("\n")
        #print(round(varex, digits))
        invisible(newx) #previously returned x
}