A monte carlo simulation study with c... PreviousNext
Mplus Discussion > Categorical Data Modeling >
Message/Author
 WJCAO posted on Tuesday, September 20, 2011 - 7:16 am
Hi,
I'm trying to use Montecarlo in Mplus to generate 40 categorical variables with 6 factors.Then MLR were used to analysis these data.But I got error:
*** ERROR in MODEL command
EFA factors are not allowed with ALGORITHM=INTEGRATION.
EFA factors are declared with (*label).
Montecarlo:
names are y1-y40;
nobservations = 500;
nreps = 500;
seed = 12345;
generate = y1-y40(3 p);
categorical = y1-y40;

Model population:
f1 by y1-y5* .4;
f2 by y6-y13* .5;
f3 by y14-y21* .6;
f4 by y22-y27* .7;
f5 by y28-y33* .8;
f6 by y34-y40* .9;
y1-y40* 1;
f1-f6*1;

ANALYSIS: ESTIMATOR = MLR;
Model:
f1 by y1-y5* .4 y6-y40*0 (*1);
f2 by y1-y5*0 y6-y13* .5 y13-y40*0(*1);
f3 by y1-y13*0 y14-y21* .6 y21-y40*0(*1);
f4 by y1-y21*0 y22-y27* .7 y28-y40*0(*1);
f5 by y1-y27*0 y28-y33* .8 y34-y40*0(*1);
f6 by y1-y33*0 y34-y40* .9(*1);
y1-y40* 1;
f1-f6*1;
OUTPUT: TECH9;
 Linda K. Muthen posted on Tuesday, September 20, 2011 - 10:17 am
I think this can be done only with the WLSMV estimator. Try that.
 WJCAO posted on Tuesday, September 20, 2011 - 6:27 pm
Thank you.But I know we can use MLR with categorical outcomes for regular EFA or for CFA.And in this simmulation, the aim of our research is to recover the difference between MLR and WLSMV. So, under this scenario, what can I do?
 Linda K. Muthen posted on Wednesday, September 21, 2011 - 7:37 am
You would need to generate the data and save the data sets in one step and analyze each data set separately and put the results together yourself. Mplusautomation may be helpful. See the left-hand column fo the website for further information.
 WJCAO posted on Wednesday, September 21, 2011 - 8:18 am
Thank you. I have do as you said.But I cannot generate data sucessfully.And the
EXAMPLE 11.5 in Mplus User's Guide cannot do too. The data I got is not a four category indicators.

Montecarlo:
names are y1-y40;
nobservations = 500;
nreps = 1;
seed = 12345;
generate = y1-y40(3 p);
categorical = y1-y40;
save=M1.dat;

Model population:
f1 by y1-y5* .4;
f2 by y6-y13* .5;
f3 by y14-y21* .6;
f4 by y22-y27* .7;
f5 by y28-y33* .8;
f6 by y34-y40* .9;
y1-y40* 1;
f1-f6@1;
 Linda K. Muthen posted on Wednesday, September 21, 2011 - 2:40 pm
You need to give threshold population parameter values in MODEL POPULATION. See mcex7.6 where three-category variables are generated.
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: