Message/Author 

Chris Stride posted on Wednesday, October 19, 2005  7:45 am



Hi I am currently trying to fit a SEM which contains interaction terms between continuous latent variables. I have upgraded from v2 to v3.13 and have hence successfully used the next XWITH command to compute my interaction term. However, when I run the analyses, the only fit indices in the output are the AIC and BIC. No chisquared stat, CFI, TLI, RMSEA. From this I have 2 questions... i) Is it possible to get these indices for this type of analyses? ii) And a sort of related question  can a SEM which contains such an interaction effect between factors be considered as nested w.r.t.o. the model that contains just the main effects? And if so, is there a way of testing for whether the model is significantly improved by fitting the interaction? cheers Chris 


i) No, this is not possible. These fit statistics are available only for models where means, variances, and covariances are sufficient statistics for model estimation. They are printed whenever they are available. ii) I think these models would be considered nested. It is the case that 2 times the loglikelihood differenece is distributed as chisquare and can be used to test the difference between two nested models. 


Dear All, I am currently doing a similar analysis and have two broad questions: 1. Is it possible to get the effect size of interaction between two continuous latent variables? It seems for RANDOM model the Rsq is not available. How can I solve this? What's the problem if I use the Rsq change from means interaction? 2. May I have references on interaction between continuous latent variables? 


Many thanks. Xiaojing 


1. To me, effect size concerns a standardized mean difference between two conditions (treatmentcontrol, malefemale, etc). But it sounds like you use this term for Rsquare. Rsquare when predictors involve products of latent variables is less straightforward because the variance and covariance related to a product of two variables is involved. Rsquare for the means would not be the same. I would simply see if the influence of the product term is significant. 2. Our User's Guide refers to KleinMoosbrugger in Psychometrika 2000 


Many thanks. 


Dear Bengt, is it just as difficult to calculate Rsquare if the predictors involve a quadratic or cubic latent variable? Thanks in advance. 


Yes, this is more involved also. 


Thank you for your reply. Nevertheless, is it possible to calculate Rsquare by hand in these cases with the help of the MplusOutput? Maybe do you have any references? 


I don't know. You might want to check the Aiken and West book on regression analysis with observed variable interactions or contact Andreas Klein. 

Steve posted on Friday, July 12, 2013  4:14 am



Dear Linda, I have read the above discussion (and recognizing that it is a bit dated now), I am wondering if it is now possible somehow to estimate model fit statistics (chisquared stat, CFI, TLI, RMSEA) when including an interaction term in the model. I have run my analysis and this information is not included in the output. I would like to do a chisquare difference test after adding in the interaction variable. Mplus is specifically referenced as being able to do this in a recent paper (Graves, Sarkis & Zhu, 2013, p. 86, Journal of Environmental Psychology), and I would like to replicate this process if you could please explain how to get this output. Many thanks. 


They probably did difference testing using 2 times the loglikelihood difference which is distributed as chisquare. This is the same as the ztest of the interaction. 

Steve posted on Friday, July 12, 2013  1:07 pm



Dear Linda, Thank you very much for this information. Have a nice weekend. 


Hello, I am a bit confused about how to compare an SEM model with an interaction to a "nested" model without an interaction using the log likelihood value. Specifically, I only get an Ho loglikelihood value in my output for the model with the interaction and do not get a scaling factor for my model without the interactions. Could you please clarify which of these I should compare? 


You should get an H0 loglikelihood and scaling correction factor in both outputs. This is what you compare. It is not really necessary, however, because the ztest of the interaction gives the same information. 


Thank you Linda. 


Hi Linda, Please can you elaborate on your post on October 09, 2013? "You should get an H0 loglikelihood and scaling correction factor in both outputs. This is what you compare. It is not really necessary, however, because the ztest of the interaction gives the same information." Does that mean that the significant level of the interaction term (Est./S.E.) would provide enough information for reporting the interaction effect? If not, can you let me know where I can find "the ztest of the interaction" you referred to? I am comparing two nested models using H0 log likelihood. I specified a model with only direct effects and the other with the direct effects and the interaction term: Null model H0 value 4103.707 free parameters 70 Interaction model (the interaction term is significant, p < 0.05) H0 value 4103.948 free parameters 71 2*(4103.948(4103.707))=0.482, which is not significant at the 0.05 level for 1 df. Does that mean that the model is not significantly different from the null model, i.e. the fit is not significantly worsen when adding the interaction term, and therefore the interaction model can be accepted? 


The ztest is in the third column of the results. It is the parameter estimate divided by its standard error. 


Hi Linda, Thanks. So does it mean that the ztest for the interaction term would be sufficient when reporting the interaction effect? Does the log likelihood value comparison between two nested models (without vs. with interaction term) necessary when the ztest is significant? if so, can you advise as to whether my reasoning in the previous post is appropriate? Thanks. 


The loglikelihood test with one degree of freedom is the same as the ztest. You don't need both. 


Dear Drs. Muthen, For a latent interaction model, I have learned that it is possible to compute a chisquare statistic by comparing the loglikelihood of interaction model versus independence model. May I know if it is possible to hand compute other fit indices (e.g., TLI, CFI) based on this particular chisquare? 


I don't think it is known how well or how poorly CFI would work in this context. 


Hi Dr. Muthen, Thank you for your reply. I am currently trying to fit my a model, which contains two interaction terms between continuous latent variables using the XWITH command to compute my interaction term. When I run the analyses, the only fit indices in the output are the AIC and BIC. No chisquared stat, CFI, TLI, RMSEA. Do I only report the AIC and BIC in my result section? Thanks. 


Right. But see also Section 1.2 Model Testing in the FAQ on our website: Latent variable interactions 

Back to top 