Message/Author 

Anonymous posted on Friday, September 03, 2004  11:42 am



I try to use your article Muthen & Currant (1997) "General Longitudinal Model.." as dipicted in Figure 5. I am using the following set up but it doesn' t work. S2 it is supposed to be the added growth. What changes need to be done in this program. Thanks VARIABLE: NAMES ARE id y1 y2 y3 grp; USEVARIABLES y1y3 grp; MISSING ARE .; GROUPING is grp(1=LGM 0=ADD); ANALYSIS: TYPE=MEANSTRUCTURE; MODEL: i s  y1@0 y2@1 y2@2; MODEL ADD: i s y1@0 y2@1 y3@2; s2  y1@0 y2@1 y3@2; i s WITH s2@0; 

bmuthen posted on Saturday, September 04, 2004  5:05 pm



See the handout for "Day 2" of our short courses  see bottom of the web page http://www.statmodel.com/trainhandouts.html 


I am also interested in the model depicted in Muthen & Curran (1997) and not quite sure how to implement it in Mplus. In the above reply did you mean the Day 2 handout from the twoday workshop on multilevel modeling in Mplus? I can't find the needed input there. Thanks! Rick 


TITLE: Aggressive behavior intervention growth model n = 111 for control group n = 75 for tx group MODEL: i s q  y1@0 y2@1 y3@2 y4@3 y5@5 y6@7 y7@9 y8@11; i t  y1@0 y2@1 y3@2 y4@3 y5@5 y6@7 y7@9 y8@11; [y1y8] (1); !alternative growth !model [i@0]; !parameterization i (2); s (3); i WITH s (4); [s] (5); [q] (6); t@0 q@0; q WITH i@0 s@0 t@0; y1y7 PWITH y2y8; t ON i; MODEL control: [s] (5); [q] (6); t ON i@0; [t@0]; 


thanks Linda! As always, you guys rock! What confuses me is that it looks like the mean of t is set to 0 in the control group [t@0] but the variance is not set to zero in the control group so doesn't this factor still exist in the control group? I am also confused as to why the t@0 and q@0 statements are included in the overall model statement, perhaps I am misunderstanding these statements as I thought they set the variance of these two growth factors equal to 0 (and I can understand why one would do this in the control group for which the Muthen & Curran Fig 5 shows that t does not exist but am not sure why one would state this in the overall model section so that it would apply for the treatment group too) 


The overall part of the model sets the *residual* variance of t at zero. This setting therefore holds in both the treatment and control group and says that in the treatment group t regressed on i has Rsquare 1 (so t varies as a function of only i)  it just turned out that this residual variance was not needed. Since the residual variance isn't mentioned in the control group, it is still fixed at zero in the control group. There was no indication of a need for a q growth factor variance so q@0 in the overall part of the model makes q a fixed effect. 


thanks I think that helps Bengt. Just to make sure, it is the "t on i@0" statement in the control model combined with the residual variance of t being set to zero in both groups that sets the variance of t to zero in the control group, right? and combined with setting its mean to equal zero in the control group makes this factor specific to the treatment group? And the t on i statement models the treatmentinitial staus interaction discussed on pp. 379  381 in Muthen & Curran? What if one didn't model the treatmentinitial status interaction  how would one make t specific to the treatment group in this case? 


Answering your questions in turn: Yes. Yes  this makes t vanish in the control group. Yes. Yes. 


great! thanks again! my last question was not yesno though. It was, if one didn't model the treatmentinitial status interaction and so didn't have the t on i@0 constraint in the control model, how would one make t specific to the treatment group? That is, how does one make t vanish in the control group without the t oni@0 contraint (which wouldn't be relevant if one didn't model the treatmentinitial status interaction)? 


Having t@0; in the overall and [t@0]; for controls, eliminates t for controls (a variable with zero variance and mean zero is zero all the time). 


when I have t@0 in the overall model (for a similar model I am currently conducting power analyses on)I get 2 warnings that the PSI matrix is not positive definite in either group with the problem involving variable t. If I put the t@0 only in the model for controls I get the same warning for the control group. Are these warnings ignorable? 


You need to send your input, data, output and license number to support@statmodel.com so we can see what you are talking about. 

Back to top 