I'm finding it hard to create a biplot in ggplot2 from a principal components analysis done in psych.
library(psych)
data("Thurstone")
pca1<-principal(Thurstone, nfactors=2, rotate="varimax")
#install.packages("ggfortify")
library(ggfortify)
autoplot(pca1)