I have a question regarding continuous-time survival analysis in MPlus. I am looking at survival in a sample of elderly over several years, based on gender. I have 3 variables...
FEMALE is coded 0=male, 1=female. DIED indicates whether the person died (0=living at end of study, 1=not living). LASTTIME is the years post-baseline at the time of death (DIED=1) or the final assessment (DIED=0).
The MPLUS input is...
USEVARIABLES ARE FEMALE DIED LASTTIME; SURVIVAL = LASTTIME (ALL); TIMECENSORED = DIED (1 = NOT 0 = RIGHT); ANALYSIS: TYPE=MISSING; BASEHAZARD = OFF; MODEL: LASTTIME ON FEMALE;
I obtain an estimate for LASTTIME ON FEMALE of -0.787
I know that women were less likely to die.
Given that DIED is not coded as in the Mplus users guide (where 0=not right censored, 1=right censored), I want to make certain that I interpret the output correctly.
Is it the log of the odds-ratio for survival for females relative to males? Is it the log of the odds of death for females to males? Is it a predicting years post-baseline at death/final assessment? Or something else?
In essence, (1) should I recode my indicator for time censoring and (2) how should I interpret the estimate for the FEMALE effect.
Thanks for your swift reply. I have a following question about the interpretation of the results. Let's consider my own data in which the continuous covariate is externalizing problems at age 3 and in which the critical event is parental divorce between age 3 and age 12. When I regress divorce (age of child in months at divorce) on ext, then I get an estimate of 0.021 (p=0.002). So, the odds of the probability of divorce is 1,021 times higher when children have a one unit increase in ext. Does this mean that children with 10 units increase have a 1.021^10=1.23=23% more chance that there parents will divorce? Or does it mean that ext is related to EARLIER divorce? Please correct me if my interpretation is wrong. Is the estimate similar to the hazard ratio? (sorry if this is a stupid question, but I am new to survival analysis). Thanks alot in advance.
The log odds increase for a 10 unit increase in x (externalizing) is 10*0.021. If you are new to survival analysis, I recommend warmly the pedagogical book by Singer & Willet, Applied Longitudinal Data Analysis. They discuss things not only in terms of formulas but also in carefully worded words, including how to interpret results.
where "p1" is the label for the parameter in the model. This gives you the hr estimate and SE from which you can compute the CI you asked for.
Or, you can get CI limits for p1 by +-1.96*SE(p1) and then exponentiate these limits to be used for the hr CI (the latter is better - see the logistic regression literature).
Reza ELAIDI posted on Wednesday, April 02, 2014 - 3:31 am
Dear Pr Muthen, I'd like to obtain the confidence interval for the hazard ratio (HR) of a covariate X1 involved in an interaction with another covariate X2. If beta1 and beta2 Cox regression estimates, then the CI for HR is exp(beta1 + beta2 +/-1.96(variance(beta1)+variance(beta2)+2(covariance(beta1*beta2)). I tried to get the covariance using TECH3 output but results are very different from what I obtain using SAS CovB option in PROC PHREG despite the same estimates are obtained with Mplus and SAS. Is TECH3 the right option to get the variance/covariance of estimates ? What could be the reason for difference with SAS ? Thank you in advance for your assistance
TECH3 is the right option. Perhaps SAS is using ML SEs and Mplus MLR SEs? If not, please send both outputs, data, and license number to Support.
Reza ELAIDI posted on Thursday, April 03, 2014 - 5:25 am
Dear Pr Muthen, You are right, SAS uses ML. The problem I face now is that I'm unable to match the estimates covariances in the TECH3 file with those of the output matrix, despite I used the variable numbering of TECH1. How is TECH3 file organized ? Than k you again for your help.
TECH3 contains the variances and covariances of the parameters. This will not match the results in the output. The first entry in TECH3 is the variance of parameter 1 from TECH1. The second parameter is the covariance between parameters 1 and 2. The third parameter is the variance of parameter 2, etc.
Laura Brown posted on Thursday, February 15, 2018 - 3:40 am
Hi, I'm following the advice from Bengt's post (Thurs Nov 13, 2008, 10:47) to calculate CIs around a hazard ratio but the values I get are non-sensical. The model gives an estimate of 0.485 for the effect of my continuous exposure on duration (SE 0.118, p-value <0.001) and an HR of 1.624 (SE 0.192, p<0.001). Following the two methods suggested in Bengt's post, I get the following CI estimates:
Method 1 can't be right because HRs and their CIs can't be negative values. Method 2 can't be right because the low p-value dictates that the relationship is significant and so the CI shouldn't cross 1. I sense-checked this by running the same cox model in Stata, which gives the CIs automatically. I got an estimate of 0.485 (SE 0.113, p-value <0.001, LCI 0.263, UCI 0.707) and an HR of 1.624 (SE 0.184, p<0.001, LCI 1.300, UCI 2.028). Neither of the methods provided gives me the correct confidence intervals here either: