Message/Author 

jsandber posted on Thursday, May 03, 2001  6:34 pm



Hi; I'm contemplating running a model as in the subject line, with either a latent class or other latent variable as the random effect(s) to capture heterogeneity. I was wondering if this was possible with Mplus, using a mixture of categorical and continuously measured covariates (this won't work with Jeroen Vermunt's LEM). I don't know how it would work, just hoping. Also, if estimating such a model is possible, is the user manual available for download so I can research how in more detail? 


Yes, you can use Mplus Version 2 to do discretetime survival analysis with a combination of categorical and continuous covariates. You can use a latent class variable to capture heterogeneity, while random effects (continuous latent variables) are not available in survival analysis with the current Mplus. The User's Guide has an example that shows how to do this. The User's Guide is not available for download. There is a forthcoming paper by Muthen & Masyn on this topic that I'd be happy to send as soon as I put the finishing touches on it. 

jsandber posted on Friday, May 04, 2001  8:16 pm



Thanks so much for your response. The latent class approach is what I prefered anyway so this would be ideal. I would definitely like a copy of the paper when it's ready. Thanks again Jack Sandberg 

Anonymous posted on Saturday, October 27, 2001  3:58 am



do you mean the discretetime survival analysis is not available in version 1? 


Yes. It requires TYPE=MIXTURE MISSING and MISSING was not available for mixture models in Version 1. 

David Bard posted on Sunday, June 18, 2006  2:22 am



I have a question related to identification of singleclass discretetime survival analyses. I was playing around with the examples in the manual using trial and error proof of principle to determine identification for various frailty specifications. I was a little confused and amazed that when I estimated an unstructured DTSA (i.e., unique thresholds for each DT variable) without covariates, estimation of a latent variance term appeared to be identified (provided I used threshold starting values that were correctly ordered and of reasonable magnitude). I should mention that while a replicable error variance did converge during estimation of this particular model, it was small and nonsignificant. However, my initial expectation was that the model would not converge at all. I've run some quick simulations and discovered this result is contingent on the presence of left censoring. Is this estimation of latent heterogeneity meaningful in this case or simply an artifact of the left censoring? In a related question/concern, when left censoring occurs, the Pearson and Likelihood Chisquare values do not equal 0 when estimating the unstructured DTSA without covariates nor latent heterogeneity specified. Are these chisquare values interpretable/meaningful in this situation or should they be adjusted for this type of censoring? Thank you. 


I don't know the answer to this question. I will try to find the answer and post it if I do. 


I would like to do something very similar to example 3.9 random coefficient regression. Only difference is that Y is binary. But the s with y part doesn't work because they have to be continuous. Is a similar model with binary outcome possible? Thanks! Kanako Ishida 


Yes, this is possible. See the partial input below that uses Example 3.5: CATEGORICAL IS u1; ANALYSIS: ESTIMATOR = ML; TYPE=RANDOM; MODEL: s  u1 ON x1 ; f BY s@1 u1*; f@1; [f@0]; u1 s ON x3; The residual covariance between s and u1 is specified using the BY option. 


Thank you, Linda. I ran the model, but the integration points are 200, it gave me a fatal error message due to memory shortage. On p.386, it says that integration points of 15 are already very heavy, so unfortunately, it seems that this model is not possible to run. I tried 1. "integration=montecarlo," 2. ignore the complex design of my survey data, 3. considerably reduce the sample size (I have 300,000 observations), and 4. reduce the number of variables to 3 as in the example, just to see if the model works, but none jointly or individually worked. Is there anything else I can do? 


Please send your input, data, output, and license number to support@statmodel.com. I would have to see exactly what you are doing to comment further. 


I know that 1 and 2 work individually, but 3, which is the combination of 1 and 2 doesn't work. I am not sure if I did wrong, or it is not possible to do it in mplus. 1. introduce random slope for x1 and residual covariance. CATEGORICAL IS u1; ANALYSIS: TYPE=RANDOM; MODEL: %within% s1  u1 ON x1; f BY s@1 u1*; f@1;[f@0]; u1 ON x2 x3 z1 z2; s1 ON x3 z1; 2. introduce random slope for X3 and twolevel structure WITHIN are x2 x3 x1; BETWEEN ARE z1 z2; CATEGORICAL IS u1; ANALYSIS: TYPE=TWOLEVEL RANDOM; MODEL: %within% u1 ON x1 x2; s2  u1 ON x3; %between% u1 ON z1; s2 on z2; 3. 1&2 %within% s1  u1 on x1; f by s1@1 u1*;f@1;[f@0]; u1 ON x2 x3; s1 ON x3; s2  u1 on x3; %between% u1 ON z1; s1 ON z1; s2 ON z2; Except X1 and Z2, all variables are dichotomous. The error messages that I get from 3 are; 1.u1 is a withinlevel variable and cannot be used on the between level variables. 2.s1 is a betweenlevel variable and cannot be used on the withinlevel variables. 3.The regression of u1 on X3 cannot be estimated because it defines random effect s2. Thanks! 


1. This should not happen unless u1 is on the within = list. 2. A random slope can vary on both within and between if you use the * approach described in the V5 UG in the random slope section of chapter 16 on pages 548549. 3. You cannot define a regression coefficient both as random and as fixed. 

Back to top 