Multiple sample path analysis with co... PreviousNext
Mplus Discussion > Multilevel Data/Complex Sample >
Message/Author
 Miha Moldovan posted on Tuesday, February 03, 2009 - 1:21 pm
Hello,

I was interested to run a path analysis comparing males versus females. I am also using complex weights and I wasn't sure how to write the syntax.
I have a PSU, Stratum and the sampling weight...is the syntax below correct? I wasn't sure what to consider PSU and stratum.

DATA:
FILE IS "E:\asci.dat";
NGROUPS = 2;
NOBSERVATION = 181 200;

VARIABLE:
NAMES ARE v1 v2 v3 v4 v5;
USEVARIABLES ARE v1 v2 v3 v4 v5;
WEIGHT IS v2; (is this the sampling weight from my data set?)
STRATIFICATION IS v4 (is this the PSU?)
CLUSTER IS v5 (is this the stratum?)


ANALYSIS:
TYPE IS COMPLEX;
ESTIMATOR IS MLR;
ITERATIONS = 1000;
CONVERGENCE = 0.00005;

I know the model part but just wanted to make sure I had the weighting part correct. Thank you.
 Linda K. Muthen posted on Tuesday, February 03, 2009 - 2:59 pm
CLUSTER is PSU and STRATIFICATION is stratum. See pages 455-460 of the Mplus User's Guide for a discussion of complex survey data features in Mplus.
 Miha Moldovan posted on Monday, February 09, 2009 - 5:40 pm
Thank you. I have another question if you could help me out. I have this syntax and it gives me lots of the same error message

Variable: USEVARIABLES
*** ERROR in VARIABLE command
Duplicate variable on NAMES list: WORRY
*** ERROR in VARIABLE command
Duplicate variable on NAMES list: RISK
*** ERROR in VARIABLE command
etc......

my syntax is as follows for a multi-sample path analysis with type complex.



TITLE: Lung Cancer

DATA:
FILE IS "E:\Lung cancer.dat";

VARIABLE:
NAMES ARE caseid fwgt wght2 know zknow hi zhi age hist1 hist2 risk worry incomes incomec
USEVARIABLES ARE caseid fwgt know hi age hist1 hist2 risk worry cigs race2 strat educ7;
MISSING ARE age(999) hist1(9) hist2(9) risk(9) worry(9) cigs(99999) race2(99) educ7(9999
GROUPING IS gender (0 = female 1 = male);
STRATIFICATION IS strat;
CLUSTER IS caseid;
WEIGHT IS fwgt;

ANALYSIS:
TYPE IS COMPLEX;
ESTIMATOR IS MLR;
ITERATIONS = 1000;
CONVERGENCE = 0.00005;

.....


Is it possible for you to tell me if I am doing something wrong in my syntax? thank you.
 Linda K. Muthen posted on Monday, February 09, 2009 - 6:38 pm
You need a semicolon at the end of the NAMES statement.
 Miha Moldovan posted on Monday, February 09, 2009 - 6:48 pm
In my syntax there is a semicolon at the end of that statement but because the statement is longer than 90 characters the software deleted parts of it. The same happened with MISSING ARE

any suggestions how to go about it?

thank you.
 Linda K. Muthen posted on Tuesday, February 10, 2009 - 5:48 am
You need to continue the record on a second line so you do not exceed 90 characters. If you cannot see how to do this and continue to have problems, send your input, data, output, and license number to support@statmodel.com.
 Miha Moldovan posted on Monday, February 23, 2009 - 2:20 pm
Hello,

I am using MLR and SB factor to compare models. I read somewhere that the fit indices also need to be adjusted as well. I looked on your website but I only found the calculation for the adjusted chi-square. Do you know if there is and where I could find the calculations for the other fit indices (CFI, RMSEA, etc)? Thank you.
 Linda K. Muthen posted on Tuesday, February 24, 2009 - 9:08 am
The chi-square statistic is adjusted so all statistics that use chi-square like RMSEA,CFI, etc. are already adjusted.
 Miha Moldovan posted on Wednesday, March 04, 2009 - 7:30 am
Dear Prof. Muthen,

I am conducting a path analysis with multiple group comparison. I know I have to report the unstandardized coefficients as opposed to the standardized ones but could you explain to me why that is? Also could you maybe provide me with a reference or two where I can learn more about this issue? Thank you.

Mihaela
 Linda K. Muthen posted on Wednesday, March 04, 2009 - 10:46 am
The basic issue with standardized coefficients in multiple group analysis is that for raw coefficients that are constrained to be equal across groups, the standardized coefficients are not equal because the standardization uses group-specific information. For more information, Google Sander Greenland who has written on this topic.
 Terry Ng-Knight posted on Wednesday, August 03, 2016 - 8:16 am
Hi,
I am running a similar model to those described above – a path model with stratification, cluster and weight variables. The main difference is that I have a binary dependent variable. I can run the model for the whole sample with MLR (with montecarlo integration), but have had to use WLSMV when I use the GROUPING command to assess gender differences. Unfortunately some key coefficients differ between the MLR and WLSMV models (I suspect due to the difference between FIML and pairwise present methods for each estimator). Is there any way I can run my model using a FIML based estimator and test gender differences?
Many thanks in advance!

Variable:
WEIGHT is Dweight;
CLUSTER is SampPSU;
STRATIFICATION is SampStra;
Names are
W7_wt ks4 zks4 Dweight SampPSU SampStra W3finwt female uni1dk uni2dk uni3dk uni1 uni2 uni3 eff1
eff2 eff3 atuni zeff1 zeff2 zeff3 zuni1 zuni2 zuni3 zks3 znssec3 zhiqual zincome zks2 russ prest zprest;
Missing are all (-9999) ;
CATEGORICAL=atuni;
USEVARIABLES ARE atuni zuni1 zuni2 zuni3 zeff1 zeff2 zeff3 znssec3 zhiqual zincome zks2 zks4;

Analysis:
TYPE IS COMPLEX;
ESTIMATOR IS MLR;
INTEGRATION=MONTECARLO;
 Bengt O. Muthen posted on Wednesday, August 03, 2016 - 2:27 pm
Stay in the MLR framework and use Type = Mixture with Knownclass to capture the groups.
 Terry Ng-Knight posted on Thursday, August 04, 2016 - 6:44 am
Thanks Bengt,
When I run the model with type=mixture and knownclass the models won't run - first I get a warning message saying it does not support INTEGRATION=MONTECARLO and when I try ALGORITHM=INTEGRATION as the warning messages suggest I get the following message:
*** FATAL ERROR
THIS MODEL CAN BE DONE ONLY WITH MONTECARLO INTEGRATION.

I can run the models if I do not specify my dependent variable as categorical - is this a viable option for multiple group testing?
 Bengt O. Muthen posted on Thursday, August 04, 2016 - 6:44 pm
Send the data, input, output and license number to Support.
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: