Message/Author 

Anonymous posted on Monday, March 22, 2004  6:20 am



The univariate theory of LoMendellRubin test for normal mixture is in article 'Testing the numer of components in a normal mixture(Biometrika 2001). MPLUS calculate the test for multivariate cases. How it's configured and where I could find related theory. 

bmuthen posted on Monday, March 22, 2004  1:37 pm



The multivariate case is a direct generalization of the univariate case. There is no writing on this. See also Vuong's original article that LMR refer to. 


Are there any Mplus examples which show how to test for best fit between a latent class model (LCM) and a latent trait model (LTM). In other words, could one use a set of four indicators and test whether a two class LCM provided better fit than a simple CFA where the indicators simply loaded onto a continuous latent variable (LTM). My understand thus far is that a fair test would involve a LTM with two graded levels versus a LCM with two classes 


See the following paper: Muthén, B. & Asparouhov, T. (2006). Item response mixture modeling: Application to tobacco dependence criteria. Addictive Behaviors, 31, 10501066. 

Paul Widdop posted on Tuesday, March 04, 2008  1:21 pm



Good Afternoon, I am looking at participation in different sporting activities. I have participation habits (binary yes no) for 15 activities, but believe there to be a latent element that explains the types of individuals that partake. I am using the following script...... Title: Stata2Mplus conversion for p:\methodology\sport.dta Data: File is p:\methodology\sport.dat ; Variable: Names are swimming snooker darts football fishing outdoor wintersp water tennis badmin squash cycling fitness cricket golf horserid yoga tenpin jog; Missing are all (9999) ; USEVARIABLES swimming snooker darts football outdoor wintersp water tennis badmin cycling fitness cricket golf tenpin jog ethnic; CATEGORICAL swimming snooker darts football newoutdo wintersp water tennis badmin cycling fitness cricket golf tenpin jog ethnic; CLASSES = C (4); Analysis: Type = mixture; Model: %Overall% Output: Tech11; ....... However, in previous work of this nature scholars have placed into the model a local dependency between the swimming and fitness activities, but I am unsure of how do this, I was wondering if anybody here could point me in the right direction. Cheers Paul 


See Example 7.16. 

Paul Widdop posted on Tuesday, March 04, 2008  2:53 pm



Thanks Linda, Paul 


Good Evening Linda, Sorry to bother you again. I have ran my model using the example 7.16, and does improve my model fit. Unfortunately I am a novice with Mplus, when I run the model with a local dependancy, the model results in the output identify thresholds for my latent classes but does not give me results in the probability scale. I guess my question is can I get these results in the probability scale from the model in 7.16. Many thanks, Paul 


You will not get the results in probability scale when numerical integration is used. You will need to compute the probabilities yourself for the items not included in the local dependency. For the items included in the local dependency, computing the probabilities would require numerical integration. 

Paul Widdop posted on Wednesday, March 19, 2008  2:43 pm



Thanks Linda, I have computed the probabilities for the items not included in the local dependency, but unsure of the method for calculating the probabilities of the items included in the local dependency. Especially using numerical integration, can I do this in MPLUS? Many thanks Paul 


These cannot be computed by hand and Mplus does not compute them. 

Rachel posted on Thursday, January 29, 2009  10:36 am



Hi! I have a question. How can I compare fitnesses between two nonnested models? One has 3 classes and another has no class(=1class). I think these models are not nested. Can I use Vuong test for these models by Mplus? thanks. 


First of all, I would use a stepwise approach, namely, testing your single class model against two classes. Then three classes against two and so on. You can use BIC or aBIC to compare nonnested class models regardless of the stepwise approach I suggested. However, Vuong test (as far as I know) requires a stepwise approach as described above, because it tests a kclass model against a (k1)class model. You can request Vuongtest via "tech11". 


Is it possible to test differences between, say, a 3class model and 1class model? I understand that the stepwise approach has been recommended, but I have a marginally significant 3 class model (using the Vuong test) and a nonsignificant 2 class model. I would like to show that the 3class model fits better than the 1class. Any help or references regarding this would be appreciated. 


Regarding Vuong test: it can happen that it's pvalue bounces across various class models. However, once Vuong test has rejected let's say two classes you can irgnore significant pvalues for 3 or 4 classes. That's why a stepwise approach is recommended when using tech11. In your case vuong test clearly points to a single class model. If you don't want to rely on tech11 alone you can use BIC to compare a single class model vs. 3 classes. You should also consider using BLRT (tech14). But here I would also test in a stepwise fashion. 


I agree. And there is not a way in Mplus to test a 3class model against a 1class model using Tech11 or Tech14. 


Hi, I am using Mplus to run some MC simulations with misspecified mixtures and I have a quick question related to the VuoungLoMendelRubin test. The test requires evaluation of a weighted sum of chisquare(1) variables, the distribution of which is stated by Lo Mendel and Rubin (2001). This same test distribution can ostensibly be used to for other hypotheses about strictly nested and possibly misspecified models. I am wondering if there is a way to ask Mplus to evaluate the weighted sum chisquare for an arbitrary pair of models. Alternatively, given either the partial derivatives of the log likelihoods or the weights themselves, can I use Mplus to compute the weighted chisquare probabilities? If possible I would like to do this within Mplus since I am running a large number of simulations. Thanks for your time. 


On a related note, does Mplus output the firstorder derivatives of the loglikelihood evaluated at the ML estimates for each observation? And the Hessian of the loglikelihood at the ML estimates? I am having a hard time finding these options if they exist. Thanks again 


No, their is no way currently to use weighted sum chisquare for an arbitrary pair of models. The only thing that comes to mind is the LOGLIKELIHOOD setting of the OUTLIER option. Mplus does not give the Hessian of the loglikelihood at the ML estimates. 


Thanks for your reply. Guess I'll have to do it the hard way : ) 

Back to top 