Conditional Model PreviousNext
Mplus Discussion > Growth Modeling of Longitudinal Data >
Message/Author
 Carolin posted on Thursday, January 03, 2013 - 3:25 am
Hi,

I'm running a GMM wit covariates. I know that in TECH4 I can find the means of the growth factors - but how can I know if they are significant?

A second question: how can I interpret the influence of covariates on s and q? For example, depressive symptoms have a significant negative effect on the (negative) linear slope, but a positive significant effect on the quadratic term in youth with externalizing symptoms?

Thanks a lot for your help
 Bengt O. Muthen posted on Thursday, January 03, 2013 - 9:49 am
To see if the growth factor means are significant you can run an unconditional model (no covariates).

Perhaps your linear slope s is negative and the quadratic q positive so that the development first goes down and then up. The negative slope on a covariate for s means that as the covariate value increases s gets larger negative values. The positive slope on a covariate for q says that as the cov value increases q increases. So all in all a deeper dip and a stronger upswing as the covariate value increases.
 Carolin posted on Friday, January 04, 2013 - 7:43 am
Thanks a lot for your quick response and the very helpful interpretation!!!

Two related questions concerning the unconditional model: due to missings in my covariates I have a smaller sample size in the conditional model (in comparison with the unconditional model). Is this a problem when I refer to the significance of means in the unconditional model?
Furthermore, the sizes of classes change a little bit when I add the covariates. Is that a problem?

Thanks a lot again!
 Bengt O. Muthen posted on Monday, January 07, 2013 - 9:43 am
I would not worry about those minor issues. But if this bothers you, you may want to think about a new 3-step mixture approach described in Web Note 15 - see our home page.
 Carolin posted on Thursday, January 10, 2013 - 4:40 am
Thank you very much, you are very helpful.

I have another question concerning the log-likelihood difference testing. First, I analyzed different models (GMM) for girls and boys separately. Now I want to compare classes with multiple group analysis (two times log-likelihood difference). For example, I am interested if the intercept of class 1 in girls is equal to that of class 1 in boys.
My question: Is the H0 model the more restrictive one (constrain two intercepts to be equal) and the H1 model the less rectrictive (no equalities)? Or differently asked: which value do I have to subtract from which one?
 Linda K. Muthen posted on Thursday, January 10, 2013 - 10:11 am
Subtract the model with the most restrictions from the model with the least restrictions.
 Carolin posted on Thursday, January 10, 2013 - 11:49 pm
Thank you!!
 Carolin posted on Thursday, January 17, 2013 - 4:15 am
Hello,

it is me again. One more question: how can I constrain linear regressions of i on covariates to be equal across groups in multiple group analyses with knownclass? I am not sure about the syntax.

Or is that the default in Mplus
- if that is the case, how can I free them?

Thanks a lot for your help!
 Linda K. Muthen posted on Thursday, January 17, 2013 - 7:38 am
I believe this is the default. You can tell by checking your results or TECH1. You free them by mentioning the in the class-specific part of the MODEL command:

%c#1%

i ON x;

%c#2%

i ON x;
 Carolin posted on Thursday, January 24, 2013 - 5:56 am
Thank you very much!!!

Can I also free the multinomial logistic regression of latent class variable (c) on covariate x for group 1 (girls) and group 2 (boys) in multiple group analysis? At the moment, I can find only one OR for each class in comparison to the reference class in the results , but I expected to have different values for girls and boys...
How can I handle this (syntax)?

I would like to compare the odds ratios via chi-square difference testing... or isn't that possible (sorry for that silly question)


Thanks a lot again for your help!
 Linda K. Muthen posted on Thursday, January 24, 2013 - 11:58 am
You should be able to free c ON x in each known class.

You can use MODEL TEST to test if the multinomial regression coefficients are equal across known classes. You would need to define the odds ratios in MODEL CONSTRAINT to test if they are different.
 Carolin posted on Friday, January 25, 2013 - 4:35 am
Could you please help me with the syntax for
a) free c ON x in each known class
b) for defining the odds ratios in Model constraint?

Thank you very much
 Bengt O. Muthen posted on Friday, January 25, 2013 - 4:41 pm
Let's see if we can get a) going first.

Have you tried

classes = cg(2) c(2);

and then for

Model cg:
%cg#1%
c on x;
%cg#2%
c on x;
 Carolin posted on Monday, January 28, 2013 - 8:08 am
That worked!!! I tried that before but I did not include the line "Model cg" before the specific model commands.

Thank you very much!

So now we can continue with b)?
 Bengt O. Muthen posted on Monday, January 28, 2013 - 8:19 am
For b), see the UG chapt3er 14, Section Calculating Probabilities From Logistic...

It shows log odds and to get odds you just exponentiate, odds = exp(logodds).
 Carolin posted on Wednesday, April 17, 2013 - 6:43 am
Hello,

I had to focus on some other things in the meantime but now I am back on this.
Unfortunately, I still do not understand how to test if two odds ratios are equal in two known classes (even after I read the recommended section in chapter 14...).
Can you tell me the syntax I need?

I would be very happy about your help...
Thanks a lot
 Bengt O. Muthen posted on Wednesday, April 17, 2013 - 12:32 pm
In the Model command, for say class 1 you simply say

c ON x (b1);

and then in Model Constraint:

New(eb);
eb = exp(b1);

to get the odds ratio.
 Carolin posted on Thursday, April 18, 2013 - 7:51 am
Thank you very much!
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: