Unidentified Model - Increase H1 iter... PreviousNext
Mplus Discussion > Multilevel Data/Complex Sample >
Message/Author
 Patchara Popaitoon posted on Wednesday, November 09, 2011 - 8:24 am
Dear Linda,

I got an error message from the analysis of multilevel data. There are several errors reported.

THE STANDARD ERRORS OF THE MODEL PARAMETER ESTIMATES MAY NOT BE
TRUSTWORTHY FOR SOME PARAMETERS DUE TO A NON-POSITIVE DEFINITE
FIRST-ORDER 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.143D-16. PROBLEM INVOLVING PARAMETER 215.

THE NONIDENTIFICATION IS MOST LIKELY DUE TO HAVING MORE PARAMETERS THAN THE
NUMBER OF CLUSTERS. REDUCE THE NUMBER OF PARAMETERS.


THE MODEL ESTIMATION TERMINATED NORMALLY

THE H1 MODEL ESTIMATION DID NOT CONVERGE. SAMPLE STATISTICS COULD NOT
BE COMPUTED. INCREASE THE NUMBER OF H1ITERATIONS.

***

I would like to know if the non-positive definite first-order derivative product matrix is the one identified as parameter 215.

Another problem is the model is unidentified due to the number of parameters over that of clusters. I can't reduce number of parameters. I remember you once suggested an analysis using Monte Carlo. Could you please suggest how to start.

Lastly, could you please advise how to increase the number H1 iteration as suggested by the system.

Thanks.
pat
 Linda K. Muthen posted on Wednesday, November 09, 2011 - 5:49 pm
The entire matrix is non-positive definite. It is because of parameter 215.

You can set up a Monte Carlo study using the sample size and parameter estimates from your study including the number of clusters. Then you can see if you can recover the parameter estimates and standard errors given that you have more parameters than clusters.

See the H1ITERATIONS option in the user's guide.
 Patchara Popaitoon posted on Thursday, November 10, 2011 - 10:22 pm
Dear Linda,

Thanks for your advise.

I got another message from the system that THE MLR STANDARD ERRORS COULD NOT BE COMPUTED. THE MLF STANDARD ERRORS WERE COMPUTED INSTEAD.

I would like to know if the results obtained given this error message are reliable.

Thanks.
Pat
 Linda K. Muthen posted on Friday, November 11, 2011 - 5:40 am
Yes, if you obtain standard errors in the results section, they are fine.
 Patchara Popaitoon posted on Friday, November 11, 2011 - 6:42 am
Dear Linda,

I tried to set up the Monte Carlo simulation for my study as suggested. I think I can follow the guideline provided in Example 12.11. I am however stuck by the term NCSIZES, CSIZES and NREP. I have the following details for my data analysis:

Type = TWOLEVEL
Number of Cluster = 213
Number of free parameters = 240
Observations = 1481

I think my NCSIZES is 213 but I don't understand how to identify the values for CSIZES and NREP.

Please advise. Thanks.

Pat
 Linda K. Muthen posted on Friday, November 11, 2011 - 11:36 am
NCSIZES specifies the number of unique cluster sizes. So if all clusters had ten observations, NCSIZES = 1, that is all clusters have the same size. CSIZES specifies the number of clusters and their sizes. CSIZES = 100 (100; says there are a hundred clusters and they are all size 10. See these options in the user's guide for further information.

The NREP option tells the number of replications you want in your simulation.
 Nadim Khatib posted on Tuesday, January 28, 2020 - 8:41 am
Hi, I'm trying to run a model with a number of data on the individual level and one indicator (IV) on the group level. I'm basically trying to predict individual outcomes using both individual and group-level data. but I keep receive this error message, even after increasing iterations. Can you please help?

here's the model:
MISSING = all (99.00, 999.00);
CLUSTER = site_id ;
Usevar = bel2 bel5 CSRVSC_2 CSRVCL_2 lead int awar YIA ;

Define:
Lead=(Lead2+Lead3+Lead4+Lead5)/4;
int = (Others1+Commun1+Commun3)/3;
awar = (Others2+Others5+Commun4)/3;
YIA = (YIA3m+YIA4m+YIA1m+YIA2m+YIA5m)/5;
ANALYSIS: TYPE = TWOLEVEL RANDOM ;
ESTIMATOR = MLR ;
Model:
%within%
Empower by lead int awar;
Heard by bel2 bel5;
Heard with YIA;
Empower on heard YIA;
CSRVCL_2 on empower ;
CSRVSC_2 on empower;
CSRVCL_2 with CSRVSC_2;
%Between%
CSRVCL_2 on YIA;
CSRVSC_2 on YIA;

And here's the error message:

THE H1 MODEL ESTIMATION DID NOT CONVERGE. CHI-SQUARE TEST AND SAMPLE STATISTICS COULD NOT BE COMPUTED. INCREASE THE NUMBER OF H1ITERATIONS.
 Bengt O. Muthen posted on Wednesday, January 29, 2020 - 5:23 pm
We need to see your full output - send to Support along with your license number.
Back to top
Add Your Message Here
Post:
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action: