I am looking for some help with factor analysis in R.
I ran a factor analysis using psych package and ended up with a four-factor model, which is fine. BUT...
I want to use the results in further analysis:
- discuss which of the factors is the most prevalent for each respondent,
- eventually, use the factor scores in regression and/or other analyses
However, the scores I get are all standardized and have a mean of 0 but min and max differ. My original data is 7-point Likert. I would like to get a similar or at least comparable range format. I tried unstandardizing them using the advice from this thread, with the following code.
unstandardized_factor_scores<-rowMeans(data)+apply(data,1,sd)*factor_scores
It did shift the distribution but still, min and max differ for each factor, the individual scores are not in the original range, and some are negative and some are, e.g. 14. Would any of you have a suggestion on what to do? All advice will be highly appreciated, thanks a lot!!!
