Message/Author 


Hi, When assessing the fit of a GMM, is it possible to test whether a slope or intercept variance should be free, or held equal to zero, using a chisquare difference test on the difference in the loglikelihoods? I.e. if the less restricted model M1 is (with likelihood L1), k=3; VARIABLE: NAMES are t1 t2 t3 t4 t5; CATEGORICAL = t1 t2 t3 t4 t5; CLASSES = c(3); ANALYSIS: TYPE = MIXTURE; ALGORITHM = integration; STARTS 200; MODEL: %OVERALL% i s  t1@0 t2@1 t3@2 t4@3 t5@4; S20; And the more restricted model M2 (likelihood M2)is the same except; MODEL: %OVERALL% i s  t1@0 t2@1 t3@2 t4@3 t5@4; IS@0 Then can we use a chisquare difference test on L1L2 with 2df? 


I would not do a difference test in this way with a mixture model. I would use MODEL TEST to test whether coefficients are equal. I would not test if a variance is zero. This is on the border of the permissible parameter space. 


Could you please clarify what exactly you recommend testing using MODEL TEST? 


Holding the parameters equal may change the classes from where they are not equal. 


I see. So to determine whether or not I should set the slope variance to zero (as in an LCGA), or let it be free, would the following be the correct use of MODEL TEST? VARIABLE: NAMES are t1 t2 t3 t4 t5; CATEGORICAL = t1 t2 t3 t4 t5; CLASSES = c(3); ANALYSIS: TYPE = MIXTURE; ALGORITHM = integration; STARTS 200; MODEL: %OVERALL% i s  t1@0 t2@1 t3@2 t4@3 t5@4; MODEL TEST: s = 0; ? 


What you specify in MODEL TEST is no different from what is given in the results section in the third column. Note that testing a variance against zero may not be correct as zero in on the boundary of the admissible parameter space. You can use BIC. Look at the BIC with variance of zero versus BIC when the variance is free. 


Ok thankyou, I'm very grateful for the advice. 

Back to top 