[R] Color-coded Biplot of Varimax rotated factors from PCA based factor analysis

I'm using R for PCA and? factor analysis. I want to create biplots of varimax rotated factors that color-code points by their classification. My research is on streams that are urban and rural. So, I want to color code them by this classification. If you just do a biplot from prcomp or princomp

[R] color-coding of biplot points for varimax rotated factors (from PCA)

ot;", length(collab)), col = c("black", "orange")) legend(x="bottomright", legend=c("Urban","Rural"), text.col=c("Blue","Green"), bg="Grey90") ### scale the plot parameters rrr<-apply(fit$scores[,1:2],2, range) (abs(r