Determining number of factors in exploratory factor analysis I'm using R's factanal function for factor analysis.  I know from reading that there are various ways of picking how many factors to use in the analysis.  I don't know which to choose, or how to do any of them.
Here's the data I have so far from factanal.  I don't understand what SS loadings are, or why degrees_of_freedom is not = min(#rows,#columns) - #factors.
Just judging from the Cumulative Var, which I think I understand, I would guess that 2 is the right number of factors, but am I right?  And if so, how can I convince others that this is the right number of factors?
factanal(x = charges[3:8], factors = 1)

#                Factor1
# SS loadings      4.779
# Proportion Var   0.797

Test of the hypothesis that 1 factor is sufficient.
The chi square statistic is 279.13 on 9 degrees of freedom.
The p-value is 6.9e-55 
factanal(x = charges[3:8], factors = 2, scores = "regression")

#                Factor1 Factor2
# SS loadings      2.817   2.544
# Proportion Var   0.470   0.424
# Cumulative Var   0.470   0.894

Test of the hypothesis that 2 factors are sufficient.
The chi square statistic is 77.1 on 4 degrees of freedom.
The p-value is 7.15e-16 
factanal(x = charges[3:8], factors = 3)

#                Factor1 Factor2 Factor3
# SS loadings      2.769   2.618   0.063
# Proportion Var   0.461   0.436   0.010
# Cumulative Var   0.461   0.898   0.908

The degrees of freedom for the model is 0 and the fit was 0.1047 
 A: There are several approaches to determining the number of factors to extract for exploratory factor analysis (EFA). However, practically all of them boil down to be either visual, or analytical.
Visual approaches are mostly based on visual representation of factors' eigenvalues (so called scree plot - see this page and this page), depending on extracted factor number. A scree plot allows to determine the number of factors to extract by detecting an area, where the curve makes relatively sharp drop (called "elbow"). Note that the scree plot term is also applicable to principal component analysis (PCA) - for a basic example, see this page.
Analytical approaches are based on various criteria and heuristics, including Kaiser criterion (eigenvalue greater than one), variance explained criterion (this heuristic's cut-off values vary from 0.8-0.9 to as low as 0.5, depending on researcher's specific goals), parallel analysis, Very Simple Structure (VSS) criterion, Velicer's MAP test and other techniques (see more details here and here, as well as via links within).
While I've tried to answer your question briefly and give mainly an overview of the topic, there are many nice answers to similar or related questions on Cross Validated, which I highly recommend you to review. For example, for PCA vs. EFA basic arguments, see this discussion. For much more advanced arguments on this topic, see this and this discussion. For applying VSS criterion, using R, see this discussion. For parallel analysis, see this discussion. 
A: I think that I cant help you specifically with R, but a common step on verifying the number of factors is to either look at the loadings and retaining as many factors as those ones with a value of 1 or more. This is the most straight way to pick the factors. If you want, you can manage to get the "scree plot", wich gives you a more visual insight on the number of factors but is more subjective. In that plot, when you get a "fall" in the curve plotted might indicate the number of factors. This approach is good when you want to reduce the number of variables of a big dataset with many variables highly correlated. So two is the number of factors to extract in your case.
Instead, if you want to really cluster variables and this factor has a meaning on people or economic behavior, you might want to select only the two or three highest loadings, because you will be able to access the Component Plot in Rotated Space, i.e. the graphical representation of the variables in a two-dimensional (or three-dimensional) space, wich gives you clues on interpreting the factor if this is your goal. 
Since you put the theoretical references for your studies at the end of your work, you will be telling others that you are not doing something sparse - it´s based on proofs stated by others. Statistically speaking, you could run a correlation analysis between the scores obtained from the pca with the response variable, for example. If the PCA was ok and the original variables are also correlated with this response variable, your factor will also be and you´ll have a very solid and understandable number to convince people.
Hope this helped you in some way.If I have said anything wrong, anyone please correct me.
