I ran an unconditional growth model with both a linear and a quadratic growth factor. My primary outcome is categorical, and treating as such in the input file using the code below. Using the code below the intercept was not estimated. I was curious why this would happen since it is estimated when I treat the categorical observed variable as continuous?
The growth model paramaterization with categorical outcomes is to hold thresholds equal and not estimate the mean of the intercept growth factor. The two equivalent but different growth model parameterizations are shown in Chapter 16 where growth models are discussed. This one is used with categorical outcomes because it offers more flexibility with ordinal variables.
I have a small sample (N=86) and 3 time points with binary outcomes. I plan to run a LCGM to find trajectory classes. Time intervals between T1 and T2 are varying across individuals between 10 and 50 years, between T2 and T3 it is the same for each person (14 years).
My idea is to center time at T2 (slope loading=0), using a loading of 1 for T3 and allow the loading of T1 to be free.
Example for a 3-class solution:
classes = c(3); categorical = t1 t2 t3; missing are t1 t2 t3 (99);
thank you for your response. When comparing with model 8.9, is it because of only 3 time points? Is there another way to test for or to confirm trajectories of these binary data (diagnosis vs. none) across three time points?
Jon Heron posted on Thursday, July 29, 2010 - 8:59 am
the method I obscurely referred to as "AT" is shown in example 6.12. Ignore the stuff with "ST" in it.
This allows for varying ages at each observation, e.g. you might have an intended data collection when the subjects are age 12 but due to a number of reasons the actual age is scattered across 11.5 - 12.5.
You can get an equivalent model by stacking the dataset into long format and fitting an MLwiN-style growth model.
This would seem a reasonable approach if your ages were spread in the way I describe, however I am a little bit wary in your case as you say the spread of ages is 10-50 years. can anything be modelled in a linear fashion across such a range?
Thank you...okay, I reshaped the dataset to long format and now it contains ID, "age" (0 at baseline) and individually varying ages for t1 and t2 as well as the associated diagnosis status (0/1). There is no way to run such a model in MPlus? What about STATA? Or did I misunderstand you?
Jon Heron posted on Monday, August 02, 2010 - 2:58 am
a3 is the time-varying covariate. This might be useful. I've fitted a growth model to 3 measures of weight (wt1-wt3) with variable times of measurement. Both these bits of syntax give the same result (wide and then long format). Should be simple to extend to multiple classes.
VARIABLE: names= wt1 wt2 wt3 agewks1 agewks2 agewks3; missing are all (-9999) ; tscores = agewks1 agewks2 agewks3; ANALYSIS: proc = 2 ; type=random ; MODEL: i s | wt1 wt2 wt3 at agewks1 agewks2 agewks3; i s ; i with s ; wt1 wt2 wt3 (equalvar) ;
VARIABLE: names= id wt agewks ; missing are all (-9999) ; cluster=id ; within=agewks ; ANALYSIS: type= twolevel random ; algorithm = integration ; integration=100 ; miterations=1000 ; MODEL: %within% s | wt on agewks ; %between% wt with s ;
model: i s | dia0 dia1 dia2 at zeit0 zeit1 zeit2; i s ; i with s ; dia0 dia1 dia2 (equalvar) ;
...I get the error message, that "proc" is unknown. Furthermore, how can I extend to classes (type=mixture?)? Is (equalvar) part of the model? didn't find any explanation about it.
Jon Heron posted on Tuesday, August 03, 2010 - 5:32 am
1] Proc problem. Have you got a really old version of Mplus - v4 perhaps??
2] Equalvar. Anything in round brackets implies a constraint. This is simply saying that the residual variances are constrained equal across time. Not necessary for the model, but needed to obtain the same answer as long-format
3] Add "type = mixture;" to the analysis, "classes = c(2);" to the variable section and "%overall%" immediately after "model:"
Thank you, Jon! I hope that's okay that I ask all these things but I've just begun to work with MPlus.
I followed all your suggestions but two errors occured:
1. tscores: are these relevant for the model? tscores are not compatible with mixture models. 2. the known error message: Categorical variable DIA1 contains less than 2 categories. ...although I doublechecked it!
Jon Heron posted on Thursday, August 05, 2010 - 4:54 am
1] Yes, you need them for allowing varying times of observation if you're modelling in wide-format. Sounds like long-format the only option then.
2] Has your dataset been read in correctly? If you fit a very simple model e.g.