When should we use Factanal and when should we go for varimax rotation in R? Recently, I had an assignment where we were asked to fetch factors from 9 to 10 variables. So, my approach was using principal function with varimax rotation. But, we can go for Factanal function as well. If we need to run multiple regression using the factors, which should be the better one to use Factanal or Principal with Varimax rotation?