Message/Author 


Hi, Is there a way I can make indicators on separate factors orthogonal to the other factors and their indicators using the monte carlo command? I'd like to produce a data set of an ideal CFA with items loading on intended factors at .6 and no modification indices from each factor to items that aren't supposed to load onto them (correlations between indicators on the same factor would be fine though). I thought adding f1 BY y1y6*.6 y8y36@0; f2 BY y7y12*.6 y1y6@0 y13y36@0; f3 BY y13y18*.6 y1y12@0 y19y36@0; f4 BY y19y24*.6 y1y18@0 y25y36@0; f5 BY y25y30*.6 y1y24@0 y31y36@0; f6 BY y31y36*.6 y1y30@0; and y1y6@.6; y7y12@.6; y13y18@.6; y19y24@.6; y25y30@.6; y31y36@.6; would do the trick, but I'm still getting small modification indices from every indicator to every factor. Thanks, Mike 


In any given sample you can get MIs that are not zero even when the population parameter value is zero. It is not clear if your residual variance statements like y1y6@.6; refer to Model Population or Model. If the former, that is not correct and would be another reason for nonzero MIs. To see how to set up Monte Carlo runs, see the Monte Carlo counterparts to the UG examples which are on the Mplus CD and also on the web site with the UG. 


Thank you for your support! I've been working at the models using the examples from ch 5 & 11. I've run into an issue that makes me unsure if I'm setting it up right. This is my syntax for a 6factor (6 indicators each) model with 6 bad crossloadings each (the .5 loadings). montecarlo: names = y1y36; nobs = 200; nreps = 1; save = highhighhigh.dat; model population: [y1y36*0] y1y36*1; f1 BY y1@1 y2y6*.6 y7y12@.5 y13y36*0; f2 BY y7@1 y8y12*.6 y13y18@.5 y1y6*0 y19y36*0; f3 BY y13@1 y14y18*.6 y19y24@.5 y1y12*0 y25y36*0; f4 BY y19@1 y20y24*.6 y25y30@.5 y1y18*0 y31y36*0; f5 BY y25@1 y26y30*.6 y31y36@.5 y1y24*0; f6 BY y31@1 y32y36*.6 y1y6@.5 y7y30*0; f1f6@1; f1 WITH f2f6*.3; f2 WITH f3f6*.3; f3 WITH f4f6*.3; f4 WITH f5f6*.3; f5 WITH f6*.3; y1y6*.36; model: f1 BY y1@1 y2y6*.6 y7y12@.5 y13y36*0; [repeats, same as above] f1 WITH f2f6*.3; [repeats, same as above] y1y6*.36; Analysis is done with another syntax file with simple syntax for the intended simple structure of the model (no crossloadings, only the 6 intended indicators loading onto each factor). 


(sorry, the syntax pushed me over the message length limit) Fit indices come out too good (CFI = .99, RMSEA = .026) for this model for me to think that the data are coming out as misspecified as I intend it. Is there an example in the manual/online of setting up models with crossloadings in the data that are not in the model? Thank you. 


In your model statement you should not fix the crossloadings as you do in for instance y7y12@.5. But then you say that analysis is done with another syntax file and it sounds like that's the one that think gets too good fit indices. We need to see that file. This sounds like it is best handled via support so please send it there. 

Back to top 