Message/Author 

RDU posted on Monday, December 01, 2008  5:54 pm



I was looking at the online technical appendix and happened upon the formula used to compute factor scores for CFA/SEM models. After looking at this, I came up with two questions. 1.) It appears that Mplus takes the posterior moment approach, which uses the posterior disttribution or conditional distribution of the latent factors (Y) given the observed data (X): f(YX) Furthermore, the means of the posterior distribution are used as estimates of the factor scores. Am I correct in my assumption here? 2.) I know that SPSS uses Bartlett's method and the AndersonRubin method to compute factor scores. Are these methods different from what is used in Mplus because they use the alternative solution as opposed to the posterior moment or Bayesian solution? Thanks, and I hope all of this is clear. 


1) For continuous outcomes and for categorical outcomes with weighted least squares estimation, Mplus uses the maximum a posteriori method. For continuous outcomes, this is also called the regression method  probably the most commonly used method for factor score estimation. For categorical outcomes this is also called MAP (e.g. in IRT). For ML with categorical and other nonnormal outcomes, Mplus uses EAP, the expected a posterio method. 2) Those methods are different. 

RDU posted on Tuesday, December 02, 2008  9:31 am



Thank you for the quick response. As a followup question, for ML with continuous outcomes, does Mplus use EAP, MAP, or something different? Best, RDU 


MAP  the regression method. 


Dear Drs. Muthén, is it possible to estimate factor scores in Mplus 5.21 using weighted likelihood estimation (WLEs)? If so, what would one have to specify in the input file? Thanks for your help! 


Sampling weights can be used in the estimation of the model from which factor scores are derived. 


Dear Drs. Muthén, I have estimated factor scores for a measurement model with 7 indicators. All load a common factor (i. e. F1 BY a b c d e f g;). Additionally, I have established a nested factor for four indicators (i. e. F2 BY d e f g;) with F1 and F2 not being correlated (F1 WITH F2@0;). When I correlate the factor scores for F1 and F2, the correlation amounts to .289. Why would that happen? Is there a way to specify that their correlation is 0? Thanks for your help! 


Factor scores and factors are not the same so this correlation is not surprising. You cannot specify the correlation is zero for the purposes of generating factor scores. 


Hello, I have been provided with factor scores for my categorical variable factors, whose items were all measured on 5point likert scales. While I appreciate the utility of these and will be using them for the correlational analyses I need to perform, I'm also expected to provide raw means and standard deviations for these factors. Can you please tell me how I rescale these factor scores to have the same Ms and SDs as my original scale? Thanks, Heather 


You won't get factor scores coming out scored as 5point Likert scales. Factors are specified as continuous normal variables. With continuous outcomes the mean of the factor is also different from the mean of the factor indicators due to the inclusion of measurement intercepts which pick up the indicator means. Also, means and variances of estimated factor scores are not the same as those of the estimated model parameters. Furthermore, the estimated factor score do not have the same correlations with other variables as the true factors do. These are some of the reasons behind the need to do a 1step SEM instead of the multistep procedure it appears you have in mind. Only if you have many good indicators would the estimated factor scores work like true factor scores. It doesn't help to rescale the estimated factor scores. It sounds like you have categorical factor indicators which adds another layer of complexity  you can't put a continuous factor score into a categorical scale. 


Hello, I am using multilevel data (students nested in classrooms nested in teachers) and have computed the factor score of a classroomlevel scale using TYPE=COMPLEX. I am hoping to use these factor scores as independent variables in an HLM model. I understand from reading Linda's response to others that the factor scores MPLUS generates are not standardized but do have a mean of zero. I am wondering how to interpret any significant results for this scale in my HLM models because I am not sure what the units of the factor score are and whether it would be correct to say that they are standardized around the grand mean. For example, can I say that a score that is 1 standard deviation above the grand mean is associated with an X standard deviation change in my dependent variable (this variable is standardized around the grand mean)? Thanks for your help! 


You can standardize your estimated factor scores before you use them in your model. But why do you take this 3step approach (FAfactorscoresHLM)? You can do this in 1 step in Mplus. 

deana desa posted on Wednesday, January 15, 2014  1:40 am



Does Mplus uses MAP or EAP for computing factor scores when MLR is used? 


We use MAP for MLR continuous and EAP for MLR categorical. 

Elina Dale posted on Wednesday, January 22, 2014  9:16 pm



Dear Dr. Muthen, In the writeup by Dr. Skrondal on MPlus website, it says the following: "With continuous variables,Mplus estimates factor scores as the maximum of the posterior distribution of the factor, which is the same as the Regression Method for factor score estimation.With this method, using factor scores as predictors gives unbiased regression slopes, but using factor scores as dependent variables gives biased slopes." Does the last part of the statement, i.e. using factor scores as predictors gives unbiased slopes, also apply to MAP method used in MPlus with categorical outcomes? Thank you! 


With MAP, slopes will be biased if they are used as independent or dependent variables. 


I have a question regarding factor scores from a CFA. I remember reading something on the discussion board that they are centered? So will have a mean of zero. In my analyses, the factors were based on original variables that were difference scores so zero was meaningful in interpreting the direction of the effect. Am I correct in thinking that the zero point for the factor scores cannot be interpreted the same way? 


Factor scores are not centered, that is, Mplus does not subtract their means. Their means may not be zero (although typically close to it) because even if the modelestimated factor mean parameter is zero the estimated factor scores don't behave exactly the same way. 


Hi Dr. Muthen, Thanks for the clarification. Just to confirm, am I correct in thinking that the factors scores cannot be interpreted the same way as the difference scores for the original variables? Specifically, for the difference scores zero indicated that there was no change between the two scores used to create the difference score? Or since the factors were created from multiple difference scores,can I assume that a factor score of zero indicates on average no change? 


If all your factor indicators have zero sample means (as it sounds like your difference score variables do), then you would expect the estimates of your factor indicator intercepts to be zero (your factor mean parameters are fixed at zero as the default) to give estimated indicator means of zero, and you would expect your estimated factor scores to have close to zero means in your sample. The factor measures what the indicators have in common. So I would say that factor scores of zero can be interpreted as an overall difference score. 


Dear all, I wondered if there is an easy way to calculate Bartlett Factor scores in Mplus for a hierarchical CFA model. I understand that the method that is described in the User Guide is to obtain Regression factor scores. Thanks in advance 


Not easy to do with Mplus but with some matrix algebra calculations outside Mplus it should not be hard. 


I have done CFA with categorical indicators, and I would also like to generate factor scores from this to use in future regression analyses. When generating the factor scores, does it make a difference whether the CFA was run as unstandardized vs. standardized? In other words, will the factor scores still be valid if I do CFA with the default (first path fixed at 1) as opposed to fixing the variances? Also, do factor scores assume that measurement is invariant throughout the sample/data used in the CFA such that it is based on the overall relationship between the items and the factors? 


Q1: The factor scores will come out in different metrics. Q2: Yes, they will be valid but come out in different metrics. Q3: Yes, everyone is assumed to follow the same model. 


I am also using complex survey data in my analysis (which has stratification, clustering, and weights). When calculating the factor scores to use in future analysis, should I run a regular CFA as opposed to a complex CFA if I will be accounting for the survey design in the future analysis? (i.e. would calculating the factor scores using Type=complex and then analyzing the scores later as complex survey data duplicate the designbased analysis?) 


No, you should use Type=Complex when getting the factor scores. Otherwise, the model estimates that the scores are influenced by won't be the right ones. 


Dear professors, I have two questions regarding the metric of factor scores and latent factor means. First, I have read in other posts that the metric of latent factor means will be identical to the metric of the indicator with fixed factor loading (though centered around zero). What if there is no fixed factor loading, but the variance of the factor is fixed instead? What determines the metric then? Second, is the metric of the factor scores similarly scaled? Thank you in advance. 


Q1Q2: Mean zero, var 1. Q3: No. 


Thank you for your quick reply. A followup to my last question: how are factor scores scaled? How to interpret a difference of e.g. 0.5 in factor scores? 


You see the factor score mean and variance  that gives you the scale for interpretations. 

Back to top 