Lois Downey posted on Thursday, June 18, 2009 - 11:03 am
Before the release of Mplus 5.21, I ran a survival analysis that produced a loglikelihood of -292.574 and a parameter estimate of -.638 (p=.007) for my predictor of interest. With Mplus 5.21, the same model generates a loglikelihood of +6346.678 and a parameter estimate of -.102 (p=.173) for my predictor of interest.
Most other models that I have compared from before and after the Mplus 5.21 release look identical. The substantial change in this model may reflect the fact that I have a relatively small sample (n=488) and a reasonably large number of free parameters (17). Do you find the difference in results believable, and would you advocate accepting the Mplus 5.21 result?
Lois Downey posted on Thursday, June 18, 2009 - 2:45 pm
As an amendment to my earlier posting, I would note that at least part of the problem may be the result of the fact that I have observations with survival time of 0. Do the two versions of 5.21 handle survival times of 0 differently?
Michael posted on Wednesday, March 07, 2012 - 12:52 pm
I am attempting to use a discrete-time survival model to examine initiation of substance use.
I have a single, dichotomous outcome variable (substance use initiation), assessed at 5 time points. I would like to include both time-invariant (e.g., gender) and time-varying (e.g., stressful experiences) covariates, also assessed at the 5 time points. The central question I am interested in is whether recent stressful experiences represent a risk factor for substance use initiation.
Following are the VARIABLE and MODEL statements I am working with.
VARIABLES ARE SUW1-SUW5 !Substance use initiation, coded 0 (no), 1 (yes), and 999 (missing or yes at a previous wave) Gender !coded 0 = female, 1 = male STW1-STW4 !Past year stressful experiences (continuous)
Hello, I may have a situation that is similar to Michael's. I have 5 time points of a discrete outcome, and 2 time varying covariate that corresponds to each of the 5 outcome time points. Here are my questions. 1) Is it possible to have interaction terms between the time varying covariates? e.g. OutcomeT2 on CovariateA_T1 CovariateB_T1 CovariateA_T1*CovariateB_T1;
2) There are theoretical reasons to believe that each time varying covariate should be correlated with itself at the next time point (CovariateA_T1 CovariateA_T2 etc). How would that be written into the survival analysis syntax?
What I would really like to try to do is have a path analysis that ends in a survival function, with the caveat that the path model varies over time. Can that be accomplished in mplus? Or would such a scenario be better modeled in mplus but using a different framework like SEM or a multilevel model with a within subjects time factor for the first level?