Factor Analysis in R

fa.alldat <- factanal(alldat[c(1,2,3)],1,scores="regression");

biplot(fa.alldat$scores,fa.alldat$loadings,cex=0.5)


factanal
using maximum likelihood factor analysis

fa
is basically using Principal Axes