Message/Author 

Stat_love posted on Monday, November 14, 2005  6:28 am



After fitting the model, my result have this kinds of message in output. (1) ONE OR MORE PARAMETERS WERE FIXED TO AVOID SINGULARITY OF THE INFORMATION MATRIX. THE SINGULARITY IS MOST LIKELY BECAUSE THE MODEL IS NOT IDENTIFIED, OR BECAUSE OF EMPTY CELLS IN THE JOINT DISTRIBUTION OF THE CATEGORICAL VARIABLES IN THE MODEL. THE FOLLOWING PARAMETERS WERE FIXED:4 > this mean that some parameteres are not meanningless isn't it? so I didn't use the anyresult after showing this message. But not sure, could you give me some comment? (2) In output, Tests of model fit part didn't give the pvalue but, Chisquare test of model fit...part give two pvalue like pearson chisqure and Likelihood ratio chisqure. Some output give different result. one is significant the other is nosignificant. How to decise test of model fit result? (3)In model results, they give estimates and S.E. est./S.E., how to get the Pvalue? just calcuate ourself? or some option? Thank you. p.s. I have a Mplus 'User's Guide' But, if I am a new user in Mplus, I have lots of Questions about output, error message, and interpretation...etc. Could you give me some comment? good book or something.? 

bmuthen posted on Monday, November 14, 2005  5:58 pm



(1) you probably have a large threshold that is fixed. Such a large thresholds means that the probability of item endorsement is essentially 1. This helps interpretation and is not a problem. (2) The Pearson and LR chisquares are for the frequency table. They are only trustworthy if the table has no small or zero cell frequencies and they typically disagree when neither can be trusted. With many items, many cells are zero. In such cases you should instead go by other means of testing described in the literature (see e.g. Muthen 2004 in the Kaplan handbook on the Mplus web site) and in our courses (BIC, loglikelihood differences, etc). (3) yes calculate yourself. I simply note if I have significance or not. We have not yet written the Mplus books that we have in mind. Until then you can watch the web videos or come to our courses. 

John posted on Sunday, December 11, 2005  6:07 am



Hello, I am new to LCA multilevel analysis. My question is: how can one reduce the number of parameters in a multilevel LCA? Having read some papers and the mplus manual, I tried to estimate a LCA with 5 dichotomous indicators for 3000 individuals in 20 Groups. For 2, 3 and four classes, I get a normal output with the AIC /BIC needed for model comparison. However, when trying to estimate a LCA with 5 classes, I get the following message: THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NONPOSITIVE DEFINITE FIRSTORDER DERIVATIVE PRODUCT MATRIX. THIS MAY BE DUE TO THE STARTING VALUES BUT MAY ALSO BE AN INDICATION OF MODEL NONIDENTIFICATION. THE CONDITION NUMBER IS 0.459D16. PROBLEM INVOLVING PARAMETER 16. THE NONIDENTIFICATION IS MOST LIKELY DUE TO HAVING MORE PARAMETERS THAN THE NUMBER OF CLUSTERS. REDUCE THE NUMBER OF PARAMETERS. many thanks for your efforts 


With 5 binary latent class indicators and 5 classes, the degrees of freedom are 2. This is a necessary but not sufficient condition for model identification. I am not sure if this model is identified. I think the major problem is that you have 29 parameters and only 20 clusters. I think trying to get 5 classes out of this data may be asking too much of it. 

Tim Victor posted on Thursday, July 19, 2007  8:04 pm



We have been modeling a binary outcome as a function of a continuous predictor, age, and higher order polynomials age, up to the 4th order. Given the shape of the distribution of outcomes across age, the distribution appears to be mixture of 2 or maybe 3 classes. How might such a model be fit using Mplus? Thanks so much 


You could use a mixture regression model as shown in Examples 7.1 and 7.2 except your dependent variable would be specified as categorical rather than continuous or count. 

Min Jung Kim posted on Thursday, February 07, 2008  1:44 pm



I am running a mixture model using continuous outcome variables repeatedly measured for 4 times. I don't know whether I need to fix intercepts (or means) of these outcome variables to zero. It seems that the default is to freely estimate intercepts of the outcome variables but to constrain those parameters to be equal across classes. Do I understand correctly? Could you clarify for me when to fix and when to release these intercepts. Also, what does mean fixing intercepts to zero? FYI, I copied my model command below. Thank you very much for your help in advance. MODEL: %OVERALL% I BY aggt2aggt5@1 ; S BY aggt2@0 aggt3@1 aggt4@2 aggt5@3 ; Q BY aggt2@0 aggt3@1 aggt4@4 aggt5@9; [aggt2aggt5@0 I* S* Q*] ; 


There are two ways a growth model can be parameterized. One is to fix the intercepts to zero and estimate the means of the growth factors. The other is to hold the intercepts equal, fix the intercept growth factor mean to zero, and estimate the other growth factor means. There is a table in Chapter 16 that shows the Mplus defaults for the special  growth model language. 

Back to top 