install.packages(c("psych", "Rcmdr","FactoMineR"),contriburl = "http://cran.stat.ucla.edu/bin/windows/contrib/2.6/")
library(pcaMethods)  #Change the repository to include Bioconductor (BioC Software) where this package can be found. Only if you want to estimate missing 
library(psych)


#Notes example

beer <- read.table("http://www.unt.edu/rss/class/mike/data/beer.txt",header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE)
beer2 = beer[,1:7]
cor(beer2, use="complete")  
bartlett.test(beer2)
FAmod1a=factor.pa(beer2, nfactors=3, rotate = "varimax")  #Note that SPSS won't produce this b/c as you can see, 'size' has an estimated loading over 1.  One can just interpret it as all the variance is accounted for in that variable.
cluster.plot(FAmod1a)