Hi. I'm a doctoral student in nutrition epidemiology. I am relatively new to SEM, but am hoping to use it for my dissertation to test mediation of the relationship between body mass index (BMI) and Coronary Heart Disease (CHD) hazard. I have continuous time-to-event (censored) CHD data, and BMI and all potential mediators are continuous, though I have some dichotomous exogenous covariates, such as gender. I understand this is doable in Mplus, but I've seen no papers testing mediation of a continuous time to event outcome and only a few papers at all on survival analysis in SEM. But I'm more familiar with the medical literature, so I'm probably not doing a very good job of searching other literature. So my questions are: a) can someone confirm that testing mediation of censored continuous time to event is doable (in Mplus), and b) does anyone know of any papers on similar analyses? Any advice would be very much appreciated. Thanks.
I'm afraid I confused the issue by talking about mediation of my time to event variable. What I meant was that my time to event variable (CHD) is my outcome variable, with body mass index as my exposure and mediators such as continuous systolic blood pressure and glucose. I don't need to model any latent variables, as all of my variables are observed, but I do want to use path analysis to assess mediation with one exposure, one time-to-event outcome and multiple mediators in the model at once. Is that not possible?
Jill McClain posted on Tuesday, February 05, 2008 - 10:49 am
I have another question on this same topic. If I run the analysis I've suggested above, will I get fit indices such as one normally gets for SEM, such as RMSEA, CFI, etc.? I'm specifically interested in comparing the fit of the same mediation model in race and gender groups (I have Whites and African Americans of both genders). I believe the CAIC and BIC would be appropriate for comparing across different populations, but will they be output with a Cox regression?
P.S. I tried to find the answer to this in the user guide rather than bothering you, but I can't find it. Sorry if I missed it. Assuming I decide to go forward with this analysis, I will acquire Mplus and learn how to use it, which will then prevent my having to ask such basic questions.
No, you will not obtain chi-square and related fit statistics. You will obtain a loglikelihood value, AIC, BIC, and sample-size adjusted BIC. For future reference, the way you could have known this is by looking at the output from the continuous-time survival examples.
We propose to test a similar type of mediation model (time-to-event outcome and single/multiple mediators).
A few follow up questions regarding the test of the mediated (indirect) effect:
1) Can Mplus compute the estimate of the indirect effect and its standard error when the outcome is a censored, event time? 2) While the test of mediation is well understood for normally distributed variables, does it make sense to use the product of coefficients (path a x path b) method to test the indirect effect, when these estimates correspond to OLS and Cox regression, respectively?
Any insight or a sense for the appropriateness of the mediation test under this setup would be greatly appreciated!
Kesinee posted on Tuesday, August 23, 2011 - 2:21 pm
Hi, I have a model with continuous survival time (age at diagnosis), m1, m2, m3, m4 are continuous (observed). While case is the event classified as a binary (0,1) or ordinary (0,1,2). My code: USEV are case x1 x2 x3 m1 m2 m3 m4 t; survival = t ; timecensored = case (0=not, 1=yes); missing values all (-9); Model: m1 m2 m3 on x1 x2 x3; M4 on x1 x2 x3 m1 m2; T on x1 x2 x3 m1 m2 m3 m4;
Questions: 1) Is it correct for the code above? What is the estimator? 2) Can we use this model when the event is ordinary? Thank you very much.
Kesinee posted on Wednesday, August 24, 2011 - 7:30 am
One more question Can Mplus 6.11 handle with late entry?
1) I would say that it all looks correct. The estimator is ML (actually MLR) for Cox proportional hazard model.
2) If the event is ordinary 0=not diagnosed, 1= diagnosed, 2=?. You might have to use competing risk for that or simply treat 2=diagnosed as well.
3) Late entry should generally be treated the same way unless at the time of entry it is already diagnosed which would mean that it is left censored and that is currently not available in Mplus 6.11
Kesinee posted on Wednesday, August 24, 2011 - 4:57 pm
1)For the ordinal outcome; 0=no disease, 1= susceptible, 2= definite, so how to use competing risk. I’m sorry I just first use Mplus. Could you please give me some advice where I can find the example for this? 2)I use age at diagnosis as a time scale. However, we recruited participants with different age and all of them free from disease at the time of entry, in SAS or R, we can use some command to account for this different such as use (exit*entry)*disease (0). I do not know that how I can do with Mplus. Since It makes more precise than use a study time and control for age.
This SAS program does implement left censoring but this is not appropriate in your case because ... the "entry" variable specifies when the subject first got exposed to the "risk" (and not the age you enrolled him/her in the study), your subjects were exposed to the risk even before you enrolled them. Here is an example of how to use the SAS option correctly. This is from
"Another useful application of the counting process formulation is delayed entry of subjects into the risk set. For example, in studying the mortality of workers exposed to a carcinogen, the survival time is chosen to be the worker’s age at death by malignant neoplasm. Any worker joining the workplace at a later age than a given event failure time is not included in the corresponding risk set. The variables of a worker consist of Entry (age at which the worker entered the workplace), Age (age at death or age censored), Status (an indicator of whether the observation time is censored, with the value 0 identifying a censored time), and X1 and X2 (explanatory variables thought to be related to survival). The specification for such an application is as follows:
I have a multiple mediation model with continuous survival as the outcome and dichotomous mediators. I know Mplus uses MLR estimation for Cox proportional hazards models, but is there a way I can use WLS estimation for the mediators in the path analysis? I was hoping to avoid the problem with using logistic regression for the dichotomous mediators. If not, do you have any suggestions for running this analysis? Thanks for your help!
Dear professors, I am a doctoral student and I´ve started to test the demo version of Mplus to see if it is the software I could purchase and use for the study I need to conduct. I would like to do a path analysis, with continuous survival as the outcome and categorical variables (with 2 or more categories) as mediators. I have three questions:
1- Is this implementation correct (for a the test with only one mediator)?
VARIABLE: NAMES = t wc cs tc; SURVIVAL = t(ALL); TIMECENSORED = tc (0=NOT 1=RIGHT); CATEGORICAL = cs; ANALYSIS: ALGORITHM = INTEGRATION; BASEHAZARD = OFF; MODEL: t ON wc cs; cs ON wc;
2- Can Mplus 7 handle left censoring in Cox? 3- Can direct and indirect effects be estimated in Mplus 7 for time-to-event outcomes? Thank you very much in advance for your help! Luciana