You may want to try using the principal
function under psych
package.
https://personality-project.org/r/html/principal.html.
Below is an excerpt from http://www.statmethods.net/advstats/factor.html
# Principal Axis Factor Analysis
library(psych)
fit <- principal(mydata, nfactors=5, rotate="varimax")
fit # print results
mydata can be a raw data matrix or a covariance matrix.