Dear Drs. Muthen, I have never used Mplus. Before I buy and learn how to use it I want to make sure that it can run the SEM analyses I need. I have a sample of 143 sites sampled at different environmental conditions and localities and in each site I have measured the abundance of three species. I wish to explore the direct and indirect effects on the relative abundance of the three species by testing the relations between these three species and ~14 different variables, and within the 14 different variables. Some of these variables are dichotomic, some ordinal and some continuous. Here are the questions: 1. If I consider the presence/absence of the three species some of my endogenous variables become dichotomic; is there a way to create a model including both endogenous and exogenous dicotomic variables together with continuous ones? 2. If I consider the relative abundance of the species, most of the sites do not have the species and only few have many, so the distribution of the endogenous variables is negative binomial. Is there a way to account for the discrepancy from normality in this case? 3. My main aim is to find the best model out of few candidates; can I calculate for the models that are fit to my previous questions AIC values? Thanks in advance, Hadas
Dear Dr. Muthen Thanks for your super-fast and helpful answer. I have read through the article, however I didn’t manage to fully understand how one can use nominal independent variable in the model. In particular, I didn’t understand how latent variables are used. I managed to run the model when I defined my independent variable as categorical but not as nominal. For example if x->y and x is nominal with three categories and y is binary categorical. How do I formulate it into the Mplus script? Thanks in advance, Hadas
Do you mean that after I create the two dummy variables I define these two new variables as categorical and run the analyses as I would with categorical variables? So, if I would have just 2 categories (for example female-0 and male-1) than I just treat the variable as categorical? Thanks, Hadas
The CATEGORICAL option is for dependent variables only. The dummy variables are independent variables. You simply use them on the right-hand side of an ON statement. In regression, independent variables can be binary or continuous. In both cases, they are treated as continuous.
When I use the dummy variables as independents I get two odd ratios for the relation of each new variable on the dependent variable. Should I average them?
Also, what should I do if the mediator is nominal. For example if x->m->y and x is categorical, m is nominal with three categories and y is binary categorical. It was hard for me to follow the relevant example in the article. Should it be the same?
Tracy Witte posted on Wednesday, March 20, 2013 - 11:12 am
I have noticed that with version 7, Mplus lets me specify binary predictor variables as categorical. In previous versions, I used to get an error message when I tried to do this, saying that the categorical descriptor only worked for dependent variables. Is this a new feature of version 7? If so, what is different about the modeling that allows this to be done?
No, the CATEGORICAL option should be used only for dependent variables. Please send the output where the error message does not show up to email@example.com.
Stephanie posted on Wednesday, January 15, 2014 - 8:36 am
I have a question on the one hand regarding observed exogenous variables and, on the other, the imputation of missing values in this variable.
I ran a sem (WLSMV because of a binary dependent variable) with three observed exogenous variables from which one is binary. This variable also includes missing values which I would like to impute.
To run the imputation and to avoid the delition of cases I additionaly included a line with "x1 x2 x3;" in the model command to mention the variance of all observed variables. All missings are defined correctly. Additionally, I use theta parameterization.
My question is, whether this is correct with a binary exogenous variable as my unstandardized coefficients for these variables are comparably rather high (e.g. 12.123) which seems to be implausible to me?
I thank you very much for your kind help!
Stephanie posted on Wednesday, January 15, 2014 - 8:55 am
I am sorry! I have forgotten to mention that I only have included "x1 x2 x3;" for all variables in the model that are observed exogenous and include missing values (NOT for all observed exogenous variables). One of them is the above mentioned binary exogenous variable. The other two variables are observed and continious and also include missing values.
Maybe I should also mention that my model includes five observed exogenous variables, in sum. But I have only mentioned those with missings in the "x1 x2 x3;" command.
You should not include observed exogenous variable variances in the MODEL command when using weighted least squares estimation. This is appropriate only for maximum likelihood estimation. When you do this, you should include variances for all observed exogenous variables not a subset of them.
Stephanie posted on Wednesday, January 15, 2014 - 1:09 pm
Thank you very much for your prompt answer.
Did I understand you right, that when using WLSMV I should include the variance for all five exogenous variables in the model command (-> "x1 x2 x3 x4 x5;") to get the missings imputed?
I have just tried this but unfortunately got three times the warning:
WARNING: VARIABLE x1 MAY BE DICHOTOMOUS BUT DECLARED AS CONTINUOUS.
No, for WLSMV you should never include the variances of the observed exogenous variables.
For maximum likelihood, if you include them, you should include all of them not a subset.
Stephanie posted on Wednesday, January 15, 2014 - 5:36 pm
Thank you once more for your answere. But may I ask you if there is a possibility to I avoid the deletion of cases when I would like to impute missing values in observed exogenous variables in an wlsmv?