A univariable and multivariable analysis of time-to-event data based on restricted mean survival time: A combination with traditional survival methods.

Abstract

Background Hazard ratio is considered as an appropriate effect measure of time-to-event data. However, hazard ratio is only valid when proportional hazards (PH) assumption is met. The use of the restricted mean survival time (RMST) is proposed and recommended without limitation of PH assumption. Method 4405 osteosarcomas were captured from Surveillance, Epidemiology and End Results Program Database. Traditional survival analyses and RMST-based analyses were integrated into a flowchart and applied for univariable and multivariable analyses, using hazard ratio (HR) and difference in RMST (survival time lost or gain, STL or STG) as effect measures. The relationship between difference in RMST and HR were explored when PH assumption was and was not met, respectively. Results In univariable analyses, using difference in RMST calculated by Kaplan-Meier methods as reference, pseudo-value regressions (R2=0.99) and inverse probability of censoring probability (IPCW) regressions with group-specific weights (R2=1.00) provided more consistent estimation on difference in RMST than IPCW with individual weights (R2=0.09). In multivariable analysis, age (HR:1.03, STL: 3.86 months), diagnosis in 1970~1980s (HR:1.39 STL:27.49 months), metastasis (HR:4.47, STL: 202 months), surgery (HR:0.58, SLG:35.55 months) and radiation (HR:1.46, SLT:44.65 months), met PH assumption and were main independent factors for overall survival. In both univariable and multivariable variables, a robust negative logarithmic linear relationship between HRs estimated by Cox regression and differences in RMST by pseudo-value regressions was only observed when PH assumption was hold (Difference in RMST = -109.3✕ln (HR) - 0.83, R² = 0.97, and Difference = -127.7✕ln (HR) – 9.49, R² = 0.93, respectively.) Conclusion The flowchart will be intuitive and helpful to instruct appropriate use of RMST based and traditional methods. RMST-based methods provided an absolute effect measure to inspect effects of covariates on survival time and promote evidence communication with HR. Difference in RMST should be reported with hazard ratio routinely.

Comparison of survival analysis approaches to modelling age at first sex among youth in Kisesa Tanzania

“…The pseudo values (PV) method uses jack-knife leave -one-out estimation to generate pseudo values, which are then used to model the effects of covariates on the outcome of interest through a generalized estimating equation. The PV method is not constrained by the assumption of homogeneity of the censorship mechanism [31,32].…”