[R] save scores from sem

I'm not aware of any routine that those the job, although I think that it could be relatively easily done by multiplication the manifest variable vector with the estimates for the specific effect. To make an example: v1; v2; v3; v4 are manifest variables that loads on one y latent variablein

Re: [R] save scores from sem

Dear Marko, I can't quite tell whether this is an original question (as suggested by no Re: in the subject field and no text from an earlier message) or an answer to a question already posed (as suggested by the phrasing of the message); if the latter, I apologize for missing the original

Re: [R] save scores from sem

Dear Joris, thanks for your reply - it is the sem case which interests me. Suppose for example I use sem to construct a CFA for a set of variables, with a single latent variable, then this could be equivalent to a PCA with a single component, couldn't it? From the PCA I could save the scores as

Re: [R] save scores from sem

I should have specified: lavaan is not familiar to me. I'm also not familiar enough with the sem package to tell you how to obtain the scores, but all information regarding the fit is in the object. With that, it shouldn't be too difficult to get the scores you want. This paper might give you some

Re: [R] save scores from sem

Dear Joris, thanks again for these very useful insights. Your point about PCA, sem, FA is clear to me now. And I understand what you say about it not being the point of sem to make corrections in a model and then try to save the scores for another analysis. But still I am wondering what would

[R] save scores from sem

Dear expeRts, sorry for such a newbie question - in PCA/factor analysis e.g. in SPSS it is possible to save scores from the factors. Is it analogously possible to save the implied scores from the latent variables in a measurement model or structural model e.g. using the sem or lavaan

Re: [R] save scores from sem

PCA and factor analysis is implemented in the core R distribution, no extra packages needed. When using princomp, they're in the object. pr.c - princomp(USArrests,scale=T) pr.c$scores # gives you the scores see ?princomp When using factanal, you can ask for regression scores or Bartlett