factor.stats(r, f, phi=NULL, n.obs = NA,alpha=.1)
VSS
, ICLUST
, and principal
for this fit statistic.factanal
(which seems to be Bartlett's test) :
$\chi^2 = (n.obs - 1 - (2 * p + 5)/6 - (2 * factors)/3)) * f$
Note that this is different from the chi square reported by the sem package which seems to use
$\chi^2 = (n.obs - 1 - (2 * p + 5)/6 - (2 * factors)/3)) * f$fa
with fm="pa" for principal axis factor analysis, fa
with fm="minres" for minimum residual factor analysis (default). factor.pa
also does principal axis factor analysis, but is deprecated, as is factor.minres
for minimum residual factor analysis. See principal
for principal components.v9 <- sim.hierarchical()
f3 <- fa(v9,3)
factor.stats(v9,f3,n.obs=500)
f3o <- fa(v9,3,fm="pa",rotate="Promax")
factor.stats(v9,f3o,n.obs=500)
Run the code above in your browser using DataLab