I am working on an EFA (Exploratory Factor Analysis) using the psych package. I would like to extract the scores from the result of this EFA at the observation level, so I can use them as an independent variable in a subsequent GLM (Generalized Linear Model). The idea is to model the relationship between the factor identified in the EFA and other predictors in my dataset. I can get the factor scores, but I need them to be in a vector with the same number of observations as in my dataset, not a single value repeated across the entire dataset. How can I do this?