back to first page..

Communalities

​http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg00284.html

cor(data)
library(mva)
d.pca<-princomp(data,cor=T)
d.pca
summary(d.pca)
loadings(d.pca)
d.pca$scores
d.FA <- factanal(factors = 3, covmat=cov(data))
d.FA
d.FA <- factanal(data, factors=3))
100*(1 - d.FA$uniquenesses)