Main page: https://www.rdocumentation.org/packages/datasets/versions/3.6.2/topics/Harman23.cor
Harman23.cor## $cov
## height arm.span forearm lower.leg weight bitro.diameter
## height 1.000 0.846 0.805 0.859 0.473 0.398
## arm.span 0.846 1.000 0.881 0.826 0.376 0.326
## forearm 0.805 0.881 1.000 0.801 0.380 0.319
## lower.leg 0.859 0.826 0.801 1.000 0.436 0.329
## weight 0.473 0.376 0.380 0.436 1.000 0.762
## bitro.diameter 0.398 0.326 0.319 0.329 0.762 1.000
## chest.girth 0.301 0.277 0.237 0.327 0.730 0.583
## chest.width 0.382 0.415 0.345 0.365 0.629 0.577
## chest.girth chest.width
## height 0.301 0.382
## arm.span 0.277 0.415
## forearm 0.237 0.345
## lower.leg 0.327 0.365
## weight 0.730 0.629
## bitro.diameter 0.583 0.577
## chest.girth 1.000 0.539
## chest.width 0.539 1.000
##
## $center
## [1] 0 0 0 0 0 0 0 0
##
## $n.obs
## [1] 305require(stats)
(Harman23.FA <- factanal(factors = 1, covmat = Harman23.cor))##
## Call:
## factanal(factors = 1, covmat = Harman23.cor)
##
## Uniquenesses:
## height arm.span forearm lower.leg weight
## 0.158 0.135 0.190 0.187 0.760
## bitro.diameter chest.girth chest.width
## 0.829 0.877 0.801
##
## Loadings:
## Factor1
## height 0.918
## arm.span 0.930
## forearm 0.900
## lower.leg 0.902
## weight 0.490
## bitro.diameter 0.413
## chest.girth 0.351
## chest.width 0.446
##
## Factor1
## SS loadings 4.064
## Proportion Var 0.508
##
## Test of the hypothesis that 1 factor is sufficient.
## The chi square statistic is 611.44 on 20 degrees of freedom.
## The p-value is 1.12e-116for(factors in 2:4) print(update(Harman23.FA, factors = factors))##
## Call:
## factanal(factors = factors, covmat = Harman23.cor)
##
## Uniquenesses:
## height arm.span forearm lower.leg weight
## 0.170 0.107 0.166 0.199 0.089
## bitro.diameter chest.girth chest.width
## 0.364 0.416 0.537
##
## Loadings:
## Factor1 Factor2
## height 0.865 0.287
## arm.span 0.927 0.181
## forearm 0.895 0.179
## lower.leg 0.859 0.252
## weight 0.233 0.925
## bitro.diameter 0.194 0.774
## chest.girth 0.134 0.752
## chest.width 0.278 0.621
##
## Factor1 Factor2
## SS loadings 3.335 2.617
## Proportion Var 0.417 0.327
## Cumulative Var 0.417 0.744
##
## Test of the hypothesis that 2 factors are sufficient.
## The chi square statistic is 75.74 on 13 degrees of freedom.
## The p-value is 6.94e-11
##
## Call:
## factanal(factors = factors, covmat = Harman23.cor)
##
## Uniquenesses:
## height arm.span forearm lower.leg weight
## 0.127 0.005 0.193 0.157 0.090
## bitro.diameter chest.girth chest.width
## 0.359 0.411 0.490
##
## Loadings:
## Factor1 Factor2 Factor3
## height 0.886 0.267 -0.130
## arm.span 0.937 0.195 0.280
## forearm 0.874 0.188
## lower.leg 0.877 0.230 -0.145
## weight 0.242 0.916 -0.106
## bitro.diameter 0.193 0.777
## chest.girth 0.137 0.755
## chest.width 0.261 0.646 0.159
##
## Factor1 Factor2 Factor3
## SS loadings 3.379 2.628 0.162
## Proportion Var 0.422 0.329 0.020
## Cumulative Var 0.422 0.751 0.771
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 22.81 on 7 degrees of freedom.
## The p-value is 0.00184
##
## Call:
## factanal(factors = factors, covmat = Harman23.cor)
##
## Uniquenesses:
## height arm.span forearm lower.leg weight
## 0.137 0.005 0.191 0.116 0.138
## bitro.diameter chest.girth chest.width
## 0.283 0.178 0.488
##
## Loadings:
## Factor1 Factor2 Factor3 Factor4
## height 0.879 0.277 -0.115
## arm.span 0.937 0.194 0.277
## forearm 0.875 0.191
## lower.leg 0.887 0.209 0.135 -0.188
## weight 0.246 0.882 0.111 -0.109
## bitro.diameter 0.187 0.822
## chest.girth 0.117 0.729 0.526
## chest.width 0.263 0.644 0.141
##
## Factor1 Factor2 Factor3 Factor4
## SS loadings 3.382 2.595 0.323 0.165
## Proportion Var 0.423 0.324 0.040 0.021
## Cumulative Var 0.423 0.747 0.787 0.808
##
## Test of the hypothesis that 4 factors are sufficient.
## The chi square statistic is 4.63 on 2 degrees of freedom.
## The p-value is 0.0988