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
rRecovery.r
"rRecovery" <-
function(R, loadings, diagCommunalities=FALSE) {
recoveredR <- loadings %*% t(loadings)
recovery <- list(R = R, recoveredR = recoveredR, difference = R - recoveredR)
if (diagCommunalities == FALSE) {diag(R) <- NA; diag(recoveredR) <- NA }
corr <- cor(c(R),c(recoveredR), use="pairwise.complete.obs")
recovery <- list(recovery, cor = corr)
return(recovery)
}