

Monte Carlo ordinal CFA y* total vari... 

Message/Author 

Doug Snee posted on Tuesday, June 16, 2009  9:46 am



Hello, I just want to make sure that I haven't lost my mind. I'm doing some Monte Carlo analyses on CFAs with 2 correlated factors (4 items each), where each indicator is ordinal. One indicator for each factor crossloads on the other factor. Each factor variance is fixed to 1.0. For the population model, I use the following in the relevant part of the syntax: MODEL population: f1@1 f2@1; f1 WITH f2@.3; f1 BY u1u4@.7 u5@.35; f2 BY u5u8@.7 u4@.35; u1u3@.51 u6u8@.51 u4u5@.2405; (I specify thresholds next in the syntax because I'm simulating ordinal indicators with latent y*s. I also tell it to make them categorical earlier in the syntax.) I just want to verify that this indeed should cause each y* variable to have a variance of 1.0 (in the population model). For u1 for instance, .7^2 + .51 = 1.0. For u4, .7^2 + .35^2 + 2*(.7*.3*.35) + .2405 = 1.0. I am not crazy, right? 


This seems correct. 

Doug Snee posted on Tuesday, June 16, 2009  11:26 pm



Thank you, Linda! I was afraid I was losing it. One related thing: If a particular estimated solution results in a variance >1 for a y*, as indicated by the above formulas, does this correspond to a situation in which Mplus flags the solution as inadmissable? In other words, are the y* automatically scaled to have variance = 1 during estimation, meaning that a variance greater than 1 as indicated by the above formulas would indicate a negative error variance? Doug 


A variance greater than one indicates a negative residual variance. 

Doug Snee posted on Wednesday, June 17, 2009  6:04 pm



Thank you, Linda. One other thing that I realized is that somehow Mplus isn't producing the nonconvergence indicator variable. (On either p. 480 or 488 of the version 3 manual, at the bottom of the page, it mentions that a nonconvergence indicator variable is produced as part of the Monte Carlo output.) The variable somehow just doesn't appear in the results that are output from each MC run. (I'm using version 3). Perhaps I need to enter some command to get it? Or it doesn't come out if you request the TECH9 output? (I did) 


Please send the output, saved data, and your license number to support@statmodel.com. 

Back to top 

