3) Calcular a matriz de correlação, removendo os valores em falta em listwise, e efetuar a análise fatorial, com extração de 5 fatores e rotação varimax.
Entretanto, tentei fazer outras rotações para além da Varimax e deparei com um problema: os resultados do R não coincidem com os do SPSS.
A rotação varimax está no pacote base e faz a normalização de Kaiser por omissão. As outras rotações estão noutros pacotes e tem que se pedir explicitamente para fazer a normalização de Kaiser. Para isso deve passar-se o parâmetro normalize=TRUE, quando se chama a função “principal”. Mas ainda assim só consegui acertar a rotação oblimin. A promax e a quartimax não coincidem com o SPSS. [1] [2]
Pesquisei mais e descobri que há uma função kaiser que pode ser invocada com dois parâmetros: a solução não rodada, e o tipo de rotação a fazer. Produz uma solução rodada com a normalização de Kaiser.
Consegui resolver todas as rotações exceto a promax, que vem no pacote psych e funciona de forma diferente.