Anonymous posted on Wednesday, June 09, 2004 - 7:10 pm
Hello. I have 2 variables (y1 & y2) measured at 4 time points. At each time point, multiple regression analyses show that the relationship between y1 & y2 is an inverted U shape, i.e., y1 is highest when y2 is at a moderate level, y1 is low when y2 is either low or high.
Growth models show that y1 did not change over time, i.e., there is significant variation among our subjects in the latent intercept factor (i1) but no variation in the latent linear slope factor. On the other hand, both the intercept (i2) and linear slope factors (s2) of y2 have significant variation among subjects. There is no evidence of a quadratic growth factor for either y1 or y2.
Growth modelling results also show that there is no correlation between i1 and i2, and between i1 and s2. This does not contradict my cross-sectional analyses. However, is there any way that I can model the data longitudinally to reflect the cross-sectional relationship?
I think this will capture the nonlinear relationship between i1 and i2.
Anonymous posted on Thursday, June 10, 2004 - 8:52 am
My growth models has random times of observation. I got the following error message when I did the above - *** ERROR in Model command The model specified with the following set of MODEL statements is not supported for TYPE = RANDOM: I2 | RESILI1 AT AGE1 I1SQ | I1 XWITH I1 I2 ON I1SQ
My program is listed below - i1 s1| contol1-contol4 at age1-age4; i1sq | i1 XWITH i1; i2| resili1-resili4 at age1-age4; resili1 resili2 resili3 resili4 (1); [contol1-contol4@0]; [i1* s1* i2*]; i1 with s1; i2 ON i1 i1sq;
I ran a parallel growth curve model that used all categorical indicators. The latent factors were not specified to be categorical factors, so the analysis continued using numerical integration on the 4 dimensions.
After the solution was found, I noticed an input error in my model. After I corrected the error, I tried to re-run the analysis. However, the analysis would not and I received an error message: "THERE IS NOT ENOUGH MEMORY SPACE TO RUN THE PROGRAM ON THE CURRENT INPUT FILE...etc."
First, what would the suggested pc-hardware specifications be for the analysis? Second, any idea how the analysis would run once, but not again? Third, what would the pros and cons be (or reference) for running the analysis using Monte Carlo as the integration method?
bmuthen posted on Wednesday, March 23, 2005 - 7:48 am
The combination of many dimensions of integration and a large sample can cause the MEMORY SPACE message - although it sounds like you haven't installed the latest update, version 3.12, which I think gives a more explicit message.
I can't say why it worked the first time and not thereafter - would have to have exact details on the two runs.
To bring down the computational demands, you can try working with fewer integration points to give an approximate ML solution, e.g.
integration = 7;
- or, as you say, use MONTECARLO integration, which might be quicker still. They should give similar results.
Kersi Antia posted on Wednesday, August 17, 2005 - 11:39 am
Hi Bengt and Linda:
I was wondering whether it is possible to specify and estimate a Growth Mixture Model (GMM) for a parallel process model. I have data on 166 salespeoples' monthly sales performance and corresponding usage of a web-based training system and have been able to specify a parallel process growth model thus far (sales performance and website usaage as the growth processes modeled). I was wondering whether it is mathematically possible to have finite classes identified wrt performance and usage with distinct though overlapping covariates, and whether MPLUS can handle this.
fati posted on Friday, November 25, 2005 - 1:06 pm
HI Bengt and Linda; I ran a growth model for two parallel processes for continuous outcome, I have doing the same thing as example 6.13 in M-plus user's guide. 1-if I want to detect the effect of first processes to the second process (cross-lagged analysis), how can I add to my program the regression between the first and the second processe? 2-it is sufficient to make only the regression between the intercepts and slope in order to mesure this effect? 3- it is logic to make the regression between intercept and slope for the same processes? 4-can I do the same analysis for the dichotomuos outcome? thank you for your help in advance
bmuthen posted on Saturday, November 26, 2005 - 6:32 am
1. Do you want the second process to be influenced by the first, or each influencing the other? In the second case you say s2 on i1; s1 on i2;
2. That depends on your application.
3. Yes, this can make sense.
FATI posted on Monday, November 28, 2005 - 6:25 am
thank you for your response;
I want just to be sure that I have understand. in basic analysis (FOR EXAMPLE, LINEAR REGRESSION ANALYSIS) if we want to see the effect the processe to the another, I make the regression between two variables (PROCESSES) (for example, detresse and satisfaction), then if I make the regression between slopes and intercepts in Grwowth model analysis for two parallel processes, that is replace the regression between variables in linear regression , if yes, can you tell me the reference to this in order to cite this in my article?
Growth factors are continuous latent variables. If you regress one on the other, it is a simple linear regression. I don't know of any special citation.
fati posted on Tuesday, November 29, 2005 - 5:43 am
Excuse me I think that my message wasn't clear enough, when we do two processes growth model, we do a regression between a slope and intercept and not between factors directly, my question is that is sufficient to mesure the effect between the factors (latent variables) or I must add a regression between factors in my model, for example in user's guide example 6.13 (page 89) how can I mesure the effects between the first processe (score of the Y1j) et the second ( score of the Y2j), do I add a regression between the Yij or it is sufficient to use the regression between intercept and slope, I think that is a basic question, I don't have enough information to the growth modeling, thanks you
I am conducting a mediation analysis in which both the mediator and outcome are repeatedly measured (5 waves). I suppose the parallel growth model is what I'm trying to do. The outcome variable is a categorical measure with linear growth over time. I used a multiple indicator growth model for the mediator measure, because that is a multivariate measure with missing values. To assess the mediation effect, we evaluate the path from intervention (0/1 dummy) to the linear trend of the mediator, to the linear trend of the outcome. From analysis results, the mediation path is significant. My question is, I'm not able to obtain fit statistics, such as chi-square, CFI, TLI, and RMSEA. How do I get these values in Mplus?
I'm having trouble with estimation, probably due to a small sample (as usual!), but perhaps due to my own ignorance. I have 2 symptoms each measured 7x (e.g., sleep disturbance & morning fatigue), and I want to estimate a parallel process growth model for them with MLR along the lines of ex6.13. The simple growth model for one indicates a significant linear and quadratic slope, the other just a linear slope. When I try to include the quadratic slope in the PP model, I get "NO CONVERGENCE. SERIOUS PROBLEMS IN ITERATIONS" and "THE MODEL MAY NOT BE IDENTIFIED" sorts of warnings. I've tried setting covariances between the uninteresting or nonsignificant latent variables to zero (e.g., i1 with s1, i1 with q1, s1 with q1) with no success. So, here I am "hat in hand" looking for suggestions!
Note: really small sample: 60 paired observations of patients and caregivers; I count 34 or so free parameters without fixing any. Recognizing that this is REALLY a small sample for so many estimates, I'm only looking for some fixed effects to explore/describe these two processes, so I'd be happy fixing some random effects to zero (like random slopes) if it would still allow estimation of the crossed regressions (e.g., caregiver morning fatigue on patient sleep disturbance).
Is it possible to do a twolevel parallel process model? I tried it with a community level variable (n=60) and an individual level variable (n=3636) where individuals have been sampled from those communities. Here is my code: Analysis: TYPE = MISSING h1 TWOLEVEL; ALGORITHM=INTEGRATION; MITERATIONS=10000;
I ran a parallel growth curve model for two latent variables of authoritative teaching based on four items each measured at three occasions. N= 900. I have missing data. I used a step by step procedure to ensure that the CFA for each variable at each occasion and for the growth models separately yielded good fit. Also the parallel model has good fit. When adding covariates this works well for observed variables (e.g. seniority - 4 categories and type of intervention -2 categories). However, when adding latent variables as covariates there is an error message saying that the model may not be identified:
The standard errors of the model parameter estimates could not be computed. The model may not be identified. Check your model. Problem involving parameter X.
When I include one latent covariate (inovat) based on 7 item, parameter x is the regression of s(lope)2 on s(slope)1. However, if I include the same latent covariate as well as one or two observed variables (intervention and gender), the problem appears in the regression of inovat on s1.
Do you have any suggestions on what the problem might be and how I can solve it?
Hello, I have two questions; 1) I want to test mediating effect with parallel process LGM, but i have a clustered data. I have 3-levels, namely 1)time, 2)children, 3)school.
I could not find any previous example, article or any message about multilevel parallel process LGM for mediation analysis (in the Mplus discussion topics also). Did anybody apply this kind of analysis before? I am not sure if it is possible because even without parallel process i had some problems with multilevel LGM part. Mplus did not run without giving any error message for some models.
Is it possible to do multilevel parallel process LGM for mediation and moderation analysis with Mplus?
I appreciate it a lot if you could give some references on that topic?
This should be possible in Mplus. 3-level analysis (2 levels in Mplus) has no extra problem with parallel processes. You should explore each process separately to see if you have enough school variation (and enough schools) to capture the parameters of the school level. I don't know any references. We have 3-level growth model examples in our growth modeling teaching in Topic 8 (coming up this March).
Hello, I am running a parallel-process model on cognitive and physical health outcomes, and have a question regarding the set-up of interactive effects between latent and observed variables. The observed variable is an interactive term between two categorical variables. When setting up the interaction between the latent variable and the observed interaction, do I need to include the main effects for the observed interaction?
Thank you for the input. I have two questions related to the proposed parallel process model: 1) Can you provide a source that explains how to probe interactions between a continuous latent variable and a continuos observed variable? 2) I would like to extend this parallel process model to account for NMAR data on the outcome variables. I have tried to implement Diggle-Kenward and Wu-Carroll selection adjustment and pattern-mixture adjustment to the I and S of both outcomes, but the models would not converge. I am able to successfully implement these NMAR models when only adjusting one of the outcome variables.
What potential biases would be introduced to a parallel-process model when only controlling for data NMAR on one of the outcome trajectories?
2) If the missingness is the same for the 2 processes, then NMAR should be done for both. But perhaps NMAR for one is better than MAR - but hard to argue for in a pub. There shouldn't be a problem doing NMAR for 2 processes.
Regarding 1) and in reference to the relevant slides for Topic 3, in parallel-process output, I am not given the unit variances for the initial status of the outcome variables, only the mean and standard errors for the estimated intercepts. Can these be used like the unit variances to standardize and aid in interpreting the interaction?
Also, in the handout and video you center the intercept for the single growth process at zero using the statement [i@0 s];. In the parallel-process framework, would I want to do this for the intercept of each process? In addition, what does including the slope (s) after i@0 invoke?
Regarding your second paragraph, yes, you would do this for the processes that you want to have an interaction with [i] for. You ask what [s] means inside [i@0 s] - it means that the mean of s is freely estimated.
Regarding your first paragraph, I don't follow you. The unit variance for the initial status can be obtained from the standardized solution. You cannot standardize by using standard errors.
Hello again, I am hoping to clarify the restrictions I am using on the parameters across dropout patterns in my pattern mixture parallel process model. I am using example 4 from Section 5.1.1 from the 2010 manuscript on non-ignorable dropout as a reference.
When specifying classes = c(1); and type = mixture; and %overall% in the model statement, are the means for the un-identified parameters held equal to those of the pattern corresponding to dropping out one time point later for all un-identified parameters (for both trajectories in the parallel process model)?
I would be happy to post more of my syntax if that would help. Thanks.
Is the default assumption for the pattern mixture model provided in example 5.1.1 that the means for the un-identified parameters are held equal to the pattern corresponding to dropping out one time point later? In the Enders piece that accompanies your work on NMAR in Psychological Methods, he places a number of identifying restrictions on the pattern mixture model to produce three latent sub-groups based on patterns of missing data.
I want to be clear that without these identifying restrictions, the default in mplus is that the means of the un-identified parameters are held equal to the pattern corresponding to dropout at t + 1.
No, it is not the default. The analyst has to apply the restrictions he/she deems appropriate since there are many possible choices as seen in the Demitras et al article I refer to. I just happen to like the approach I used in my Psych Meth article.
Using the Psychological Methods article and corresponding syntax (5.1.1), please tell me if this is correct:
The identifying restrictions constrain the effect of the slope on d1-d2 to be equal, and the effect of d1-d3 on q to be equal. For the slope constraint, the mean of patterns corresponding to d1-d2 are held equal to the mean of the pattern corresponding to dropout at d3. For the quadratic constraint, the mean of patterns corresponding to d1-d3 are held equal to the mean of the pattern corresponding to dropout at d4.
i-q on d4-d5; i on d1 d2 d3; s on d3; s on d1 (1); s on d2 (1); q on d1 (2); q on d2 (2); q on d3 (2);
Hello, I am running a parallel process latent growth model with two dichotomous sets of repeated measures and wlsmv estimation. The univariate latent growth models show that each of these curves have intercepts, linear slopes, and quadratic slopes (one curve has no variance in the quadratic). I was hoping to be able to regress the linear slope of one curve on the other to say something about how the two processes are related over time (I expect that growth in one of the processes is driving growth in the other), but I'm not quite sure how the quadratic slopes should be modeled to fit into this interpretation, especially since one of the quadratic slopes has zero variance and is thus not correlated with the other intercepts and slopes. Any thoughts or advice would be most appreciated!
It is hard to work with a model that's predicting from a linear and a quadratic slope because those 2 growth factors are not separately interpretable. In contrast, I can imagine that the intercept of one process influences the 2 slope factors of the other.
It is ok to have predictors of a growth factor with zero residual variance - you may still find a signficant predictor. This happens quite frequently.
I am currently struggling with the fit of a parallel proces growth model. Both seperate growth models have a good fit, but when we put them together, the fit of the "global" model is not good. Do you have suggestions to solve this problem, or the reasons for the occurence of this problem? Is it possible that the two outcomes are too highly correlated? (correlation of the two slopes was .70).
Thank you for your response. We have tried to improve the model fit by following the modification indices. However, the correlation between slopes remains extremely high (.90). Can we report this, or is this comparable to multicollinearity?
I'm working on a growth model of moderated mediation with a 3-wave dataset and I have some questions:
1. My 2 interaction terms are between the intercept of one variable and the slope of another one. Should I specify this as: Interaction | s1 XWITH i6 ?
2. Many of my slopes should not follow necessarily an 1-2-3 fashion, but most likely 1-2-1 or 1-2-2. If for 1-2-3 I use 0-1-2 as factor loadings, what should the factor loadings be for 1-2-1 and 1-2-2?
3. Some of my variables should not have slope but only intercept. Does it simply suffice to enter s@0 in the model?
Any help would be very much appreciated, thanx in advance! Paris
1. Using 0 1 * as factor loadings results in the following error: Error in parsing line: "X_T1@0X_T2@1 X_T3@*"
2. If next to i and s I also specify a "q" (quadratic term), don't I actually test the above question?
3. One term of my interactions is a variable measured three times which I assume to be stable over time. There are two ways to model that: i) Specify its i and s and constrain s to be equal to 0 or ii) Make a latent factor with the three measurments loading onto the factor. What method is the most preferable?
It is a good idea to test each model separately before estimating a model with both processes. If the mean of the slope growth factor is not signficant, this means that there is no development over time for the process.
I have a follow-up question on the parallel processes LGCM. One of the model does not show a sig. mean growth factor, and the other one does. When analyzing them in the parallel process, both slopes are sig. correlated as well as their intercepts. I am just wondering whether the results are valid for interpretation.
I ran a parallel process LGM with both processes having positive linear slope, a negative quadratic slope, and a positive cubic term. I found that the linear terms covary significantly and positively with one another, which means that individuals who show faster gains over time in x show faster gains over time in y as well. The quadratic terms covary significantly and negatively with one another, but I am having some difficulty trying to interpret this. Would you be able to advice me on this? Thank you.
It is hard to give a separate interpretation of linear and quadratic growth factors since they interact.
Neha Goyal posted on Monday, July 14, 2014 - 4:19 pm
I have two outcomes which I seperately modeled. The best unconditional model for one variable was a piecewise model with one linear piece from 0-8 months and second linear piece from 8-60 months. Furthermore, I had to set the slope variance of the first linear piece to zero since I encountered a non-positive definite covariance matrix. The best unconditional model for the second variable was a piecewise model with one linear piece from 0-18 months and the second linear piece from 18-60 months. In other words, they have different growth trajectories. My questions are
1) If the variance of one piece for a variable is zero, does it even make sense to do a parallel process model. 2) How do you do a parallel process model with two variables with different piecewise growth trajectories, if at all?
1) You still have variance in the intercept and therefore you want to correlate them over processes.
2) Just correlate all the growth factor that have variances.
Myung Lee posted on Monday, December 22, 2014 - 3:44 pm
I estimated a parallel processes latent growth curve model with two continuous outcomes measured at three time points. I would like to study the correlation between the two slope factors. In the model, slope loadings on the first factor were specified as linear (0, 1, 2) and those on the second factor were specified as nonlinear using a free time score at T2 (0, *, 1).
Given that the two slope factors have different shapes of the growth trajectories, would I still be able to interpret the slope factor correlation (e.g., r = .70) as: the rate of change in the first outcome from T1 to T3 traveled together with the rate of change in the second outcome from T1 to T3?
Hi, when I tried to run a growth model with an ordinal variable (scaled 1-6) using the categorical option I got the following error message: The categorical variables in the growth model do not have the same number of categories. Use the CATEGORICAL option to allow the number of categories to differ for maximum likelihood estimation. Problem with: I S Do you know what went wrong?
Another question is if it is possible to test 4 potential mediators (2 continous, 2 ordinal) between an intervention (binary, randomized) and a continous outcome variable within one model using parallel process latent growth curve modeling (do you know any example)?
I am working with a parallel process latent growth curve model, with two continuous outcomes measured at four time points. Reviewers have requested that I provide an analysis of how the common variance in the two outcomes is changing across time and how variables predict this common variance. Are you aware of an analysis in Mplus that could answer this question?
Thank you for your reply. From the Estimated Model part of the output I can see the covariances of outcomes 1 and 2 at each time point. I notice that the covariances are decreasing over time. Is there a way to model the predictors of the covariance at each time point? Also, is it possible to model the predictors of the non-shared or error variance in order to understand why the trajectories become increasingly differentiated over time?
Liu BAI posted on Friday, November 03, 2017 - 8:02 am
Dear Linda and Bengt, I'm running a parallel growth curve model to test if two people change together. My model statements are as follows and I have 2 questions. MODEL: ii is| x1976@0x1977@firstname.lastname@example.org@5.63;
First Question: When I test if one's intercept could impact the other one's slope using: is on mi; ms on ii;
I get a warning saying that: "THE LATENT VARIABLE COVARIANCE MATRIX (PSI) IS NOT POSITIVE DEFINITE." But if I change it to WITH statement as above, the model fits good with no warning. I wonder why this happened. And, do I have to use ON statement if I'm not interested in the effects from one to the other?
Second Question: I wonder if it is possible to get the covariances between the two slope for each dyad, because I'm interested how this co-change cause influence or influenced by other variables? I tried SAVE = FSCORES, but it only gave the variance for each variable no the covariances.
Q1: Send the output with the warning to Support along with your license number.
Q2: You don't have to use ON.
Q3: Use TECH4.
Liu BAI posted on Friday, November 03, 2017 - 6:16 pm
Thank you for your reply! About my third question. I think I didn't make it clear enough.
I have two people in the same family and their cortisol variables. I want to see if these two family members have similar cortisol patterns across the day, therefore I used the parallel growth curve model.
Now I wonder if I can get a score measuring the synchrony between the two family members ( (i.e., each family has their only score)? I'm think the the slope covariances for each family. Maybe something like Empirical Bayes covariance. I wonder if it is possible?
In your model, the covariance between the slopes is not a random parameter that varies across units (families). It is assumed to be a regular, fixed parameter. Therefore, you cannot get an estimate of a covariance score for each family.