https://github.com/cran/nFactors
Tip revision: 923d0cc1f43c36debbea1f1fb06e4de448065380 authored by Gilles Raiche on 31 August 2019, 09:11:55 UTC
version 2.3.3.1
version 2.3.3.1
Tip revision: 923d0cc
componentAxis.r
"componentAxis" <-
function(R, nFactors=2) {
nVar <- dim(R)[2]
acp <- principalComponents(R)
values <- acp$values[(1:nFactors)]
varExplained <- round((values/nVar)*100, 2)
cumVarExplained <- round(cumsum(varExplained), 2)
loadings <- acp$vectors[,(1:nFactors)] %*% diag(values^0.5) # F1 * diag(E)
communalities <- apply(loadings*loadings,1,sum)
apa <- list(values = values,
varExplained = varExplained,
cumVarExplained = cumVarExplained,
loadings = loadings,
communalities = communalities)
return(apa)
}