External Monte Carlo Sets PreviousNext
Mplus Discussion > Exploratory Factor Analysis >
Message/Author
 Eric Turkheimer posted on Monday, January 15, 2007 - 7:14 am
I would like to generate monte carlo data sets using an external program and then model them in mplus. Is there a way to accomplish this?

Thanks,
Eric
 Linda K. Muthen posted on Monday, January 15, 2007 - 7:57 am
See Step 2 of Example 11.6.
 Eric Turkheimer posted on Monday, January 15, 2007 - 11:35 am
Thanks, that helps. I now have...

TITLE: Montecarlo of snap;
DATA: FILE = "z:\projects\snap archive\boot.txt";
TYPE = MONTECARLO;
VARIABLE: NAMES = y1-y285;
categorical are y1-y285;
MISSING = .;

ANALYSIS: TYPE = efa 1 2;

OUTPUT: TECH9;

The job runs through the (currently 3) datasets that I specify, but the output doesn't contain anything about the EFA, it ends after computing the correlation matrix. TECH9 lists each iteration, but no errors. I also note that each dataset runs fine on its own.

Here is the tail end of the output:


SAMPLE TETRACHORIC CORRELATIONS
Y281 Y282 Y283 Y284 Y285
________ ________ ________ ________ ________
Y282 0.012
Y283 -0.032 -0.012
Y284 -0.002 -0.013 0.007
Y285 0.024 -0.049 0.023 0.003


Errors for replication with data file z:\projects\snap archive\boot1.dat:


Errors for replication with data file z:\projects\snap archive\boot2.dat:


Errors for replication with data file z:\projects\snap archive\boot3.dat:


If that isn't enough to figure anything out, I can attach data.

Thanks,
Eric
 Linda K. Muthen posted on Monday, January 15, 2007 - 11:58 am
Please send your input, data, output, and license number to support@statmodel.com so I can see what is happening.
 Eric Turkheimer posted on Wednesday, January 17, 2007 - 8:30 am
Linda, by way of follow-up, I thought I might try the same problem using the Mplus monte carlo data generation rather than doing it internally. That would entail running a step one model to simply estimate all 285 thresholds, which could then be used as the basis for a monte carlo run using the population statement.

I guess my question is, how can I run a model to do nothing but estimate and save the thresholds for a whole lot of categorical variables?
 Linda K. Muthen posted on Wednesday, January 17, 2007 - 9:05 am
I based this on the Monte Carlo counterpart Example 3.12:

MODEL POPULATION:
[u1$1*-1 u1$2*1];
[u2$1*0];
[u3$1*-.5 u3$2*.5 u3$3*1.5];
u1-u3*1;
 Hsien-Yuan Hsu posted on Tuesday, July 15, 2008 - 3:24 pm
Dr. Muthen

I use MONTECARLO to generate 100 dataset.
Then, I follow Example 11.6.

I want to run a EFA model.
I write:
ANALYSIS: TYPE = EFA 1 2;

Mplus doesn't show any results.

Then I try to run a CFA model.
I just write:
model:
f1 by y1-y10;

I get the results.

Q: Does Mplus allow us to get EFA results based on 100 dataset ?

Thanks.

Mark
 Linda K. Muthen posted on Tuesday, July 15, 2008 - 3:34 pm
EFA and Monte Carlo are not allowed because of problems with the order of the factors changing.
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: