# Leitura dos dados do arquivo "L03_CAP02_CORRAR.csv"
base <- read.table("L03_CAP02_CORRAR.csv",header=TRUE,sep=";",dec=",")
#
names(base)
head(base)
#
# Seleção de variáveis: x.vars
x.vars <- subset(base, select = ICOM:ILGE)
# Determinação do número de fatores a serem "extraídos":
library(nFactors)
av <- eigen(cor(x.vars)) # Obtem autovalores
ap <- parallel(subject=nrow(x.vars),var=ncol(x.vars),
               rep=100,cent=.05)
nS <- nScree(x=av$values, aparallel=ap$eigen$qevpea)
plotnScree(nS)
# Análise Fatorial para 6 fatores ##########
library(psych)
library(GPArotation)
# Análise para 6 fatores * sem rotação *
fit.FA <- fa(x.vars, nfactors=6, fm ="pa")
fit.FA[]
#
# Análise para 6 fatores * rotação varimax *
fit.FA.r <- fa(x.vars, nfactors=6, fm ="pa", rotate = "varimax")
fit.FA.r[]