Type=RANDOM and standardization PreviousNext
Mplus Discussion > Growth Modeling of Longitudinal Data >
 Daniel E Bontempo posted on Wednesday, November 28, 2007 - 10:37 am
Hi -

I am seeking further insight into why standardized estimates and confidance intervals are not available when type=RANDOM. Is it a computational problem or is it in some way thoretically invalid?

We have a tricotomus item checklist with 5 factor structure and data at 4 waves (intervals vary enough to require inidivdual time basis). Because of the large number of indicators for each factor and other data conditions, we are not able to run a 5-factor multivariate growth model. What we have done is run pairwise growth models, getting the covariance of the slopes in each model. Then, for tabulation purposes,we standardize the covariances using the estimated slope variances. I'm not sure these are true correlations, but they are in a common metric we can compare.

My concern is that the reason standardized values are not available might speak to some gross invalidity to what we are doing. Any further discussion/comment would be appreciated.
 Linda K. Muthen posted on Wednesday, November 28, 2007 - 2:45 pm
The basic idea is that for TYPE=RANDOM, there is not a single variance/covariance matrix but rather one for each value of x. This makes in unclear how one would standardize.
 Daniel E Bontempo posted on Wednesday, November 28, 2007 - 3:32 pm
I see the problem.

Is there some better way to compare associations of slopes if they are not all in the same model?

Setting aside the curves of factors (as describe above). Say I calculate the manifest score of factors A & B and a covariate C at each occasion. I run two multivariate growth models with TYPE=RANDOM. One has trajectory of A with trajectory of C, while the 2nd has trajectory of B with trajectory of C. A and B are not on the same scale because each has a different number of items. IF I want to compare the relative associations A with C and B with C, does our claculation offer any value?
 Bengt O. Muthen posted on Thursday, November 29, 2007 - 8:51 am
Sounds like you are asking about two different things. You can certainly use a standardization of a covariance between two slope growth factors - that correlation is fine. In contrast your last paragraph seems to talk about regressions. But it confuses me. It first talks about factors (I assume DVs) A and B and a covariate C, so it sounds like regressing A or B on C - but then you say "trajectory A with trajectory of C", which I don't understand.
 Daniel E Bontempo posted on Thursday, November 29, 2007 - 9:59 am
Sorry for any confusion. We have done this calculation for a few different growth models and I had them collectively in mind. Let me try to clarify to get us all on the same age. The same underlying data is used and we have 4 waves.

In one case we had factors at each occasion with tricotomous indicators, linear growth model, obtained the covariance of the growth parameters. We did this pairwise because we could not get the 5-factor model to run. I called this the curves of factors model, or 2nd order growth. The growth parameters are not regressed on any covariates. We obtain the covariances of the growth parameters from MPlus and standardize them for tabulation using the estimated variances of the respective growth parameters.

In the 2nd application we used the manifest scores of the factors as well as the manifest score of another instrument. We did five multivariate (two linear growth curves) models using each factor score respectively with the 2nd instrument. (I distinguish this as 1st order growth model). In this model we do regress the growth parameters on covariates such as gender and IQ deficit. We obtain the growth parameter covariances (conditioned on the covariates) and standardize them with the estimated residual variances of the growth parameters.
 Daniel E Bontempo posted on Thursday, November 29, 2007 - 10:00 am
When I refer to "trajectory with" I am only talking about the growth parameter covariances (level-level, level-slope, slope-slope) obtained from our MPlus Model. I guess I, somewhat confusingly, distinguished A & B as factors and C as a covariate because A & B are two dimensions of the instrument with 5 factor structure, and we have some models that include the item level and estimate the factors. C on the other hand is a global health assessment and is always in our models as a summary score. Sorry for any confusion.

If I understand you, you see no problem with our standardization of the covariance.

Is this true in both applications, including the one where the growth parameters are conditioned on a few covariates?

How does this relate to Type=Random and Linda's initial response about individual covariance matrices?
 Bengt O. Muthen posted on Thursday, November 29, 2007 - 8:20 pm
Didn't this thread start with you using Type = Random for which standardizations are not produced? You work with growth modeling so perhaps you used Type=Random to handle individually-varying times of observation. Linda's answer concerned having a random slope in a regression of y on x in which case the y variance is a function of the x values and standardization is therefore not well defined. But it sounds like you have a growth model with regular regression on covariates and standardization would be straightforward and fine. As long as you don't use the standardized coefficients for comparisons across groups - but it sounds like you work with the same group and just different DVs.

At one point you mention standardizing with respect to the residual variance of a DV, but you want to use the full variance.
 Daniel E Bontempo posted on Friday, November 30, 2007 - 1:33 pm
Bengt, thanks for the followup. I am doing growth modeling, and we do have individual times of observation. We are using Type=Random, and if we ask for STD in the OUTPUT, we get the message that it is not available for Type=Random. Or, at least we do inthe growth model where the growth parameters are regressed on gender and IQ.

Yes. it is the same group and just different DV in the growth model. Glad this is not a problem as I did not think of another way.

Concerning your point about full variances, do you mean that I should run an unconditional model, note the vairances, and then use those to standardize the covariances reported in the conditional model? Residual variances are reported in the conditional model. Or, are the full variances in one of the tech outputs?
 Bengt O. Muthen posted on Friday, November 30, 2007 - 3:13 pm
Yes, Type = Random does not produce standardized values. When the only reason for Type = Random is due to "AT" (individually-varying times of observation) as in your case, it would in principle be possible to provide standardized results (current Mplus just shuts it off whenever RANDOM is specified) - in that case we don't have a random slope for the covariate (only for the time variable), which is the sticking point. To make a long story short, in your case standardization of coefficients for growth factors regressed on covariates is appropriate and you could compute these by hand given the output in Mplus (see standardization formulas in the new Mplus Version 5 User's Guide, chapter 17).

Regarding your last question, you should use the conditional run and just get the full variances from Tech4.
 Daniel E Bontempo posted on Friday, November 30, 2007 - 5:35 pm
Thanks again. This is so helpful.

Unfortunately TECH4 is also not available when Type=RANDOM. But, I do not mund doing the uncnditional run if that is what I need to do.

If I can ask, since the covariance is of the residualized levels and slopes, why is it necessary to use the full variances instead of the residual variances? Is this not conceptually a covariance of residuals which generally would be standardized with residual variances?

That aside, the formulas you mention refer to regression coefficients. No formula is explicitly given for covariances. I am dividing the covariance by the square root of the product of the variances. Correct?
 Bengt O. Muthen posted on Friday, November 30, 2007 - 6:01 pm
Talking about growth factors regressed on covariates - no need to approximate by using the unconditional run - unless you have many covariates - it can be hand calculated.

Talking about WITH parameters - yes, if you are interested in residual correlations you should certainly use the residual variances.

That's right - square root of the product of the variances.
Back to top
Add Your Message Here
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Options: Enable HTML code in message
Automatically activate URLs in message