Showing overlap and need to condense into factors
psych::VSS.scree(cor(cordat))
factanal(x = cordat, factors = 2, rotation = "varimax")
##
## Call:
## factanal(x = cordat, factors = 2, rotation = "varimax")
##
## Uniquenesses:
## unemp poverty medinc ownerocc black h.density hs
## 0.497 0.176 0.434 0.005 0.730 0.807 0.491
##
## Loadings:
## Factor1 Factor2
## unemp 0.709
## poverty 0.867 -0.268
## medinc -0.725 0.203
## ownerocc -0.454 0.888
## black 0.515
## h.density -0.438
## hs -0.710
##
## Factor1 Factor2
## SS loadings 2.756 1.104
## Proportion Var 0.394 0.158
## Cumulative Var 0.394 0.551
##
## Test of the hypothesis that 2 factors are sufficient.
## The chi square statistic is 17285.46 on 8 degrees of freedom.
## The p-value is 0
fitAfterRotation <- factanal(cordat, factors = 2, rotation = "varimax")
print(fitAfterRotation$loadings, cutoff = .40, sort = TRUE)
##
## Loadings:
## Factor1 Factor2
## unemp 0.709
## poverty 0.867
## medinc -0.725
## black 0.515
## hs -0.710
## ownerocc -0.454 0.888
## h.density -0.438
##
## Factor1 Factor2
## SS loadings 2.756 1.104
## Proportion Var 0.394 0.158
## Cumulative Var 0.394 0.551
f <- factanal(cordat, factors=2, rotation="varimax", scores="regression")
dat2 <- cbind(dat2, f$scores)
##
## Call:
## factanal(x = cordat5, factors = 2, rotation = "varimax")
##
## Uniquenesses:
## unemp poverty medinc ownerocc black h.density hs
## 0.508 0.183 0.441 0.005 0.739 0.805 0.493
##
## Loadings:
## Factor1 Factor2
## unemp 0.701
## poverty 0.863 -0.268
## medinc -0.720 0.199
## ownerocc -0.450 0.890
## black 0.503
## h.density -0.441
## hs -0.708
##
## Factor1 Factor2
## SS loadings 2.714 1.110
## Proportion Var 0.388 0.159
## Cumulative Var 0.388 0.546
##
## Test of the hypothesis that 2 factors are sufficient.
## The chi square statistic is 92285.95 on 8 degrees of freedom.
## The p-value is 0
##
## Loadings:
## Factor1 Factor2
## unemp 0.701
## poverty 0.863
## medinc -0.720
## black 0.503
## hs -0.708
## ownerocc -0.450 0.890
## h.density -0.441
##
## Factor1 Factor2
## SS loadings 2.714 1.110
## Proportion Var 0.388 0.159
## Cumulative Var 0.388 0.546