A Univariable and Multivariable Analysis of Right Censored Time-to-event Data Based on Restricted Mean Survival Time: A Combination with Traditional Survival Methods.

Qiao Huang,Jun Lyv,Bing-hui Li,Lin-lu Ma,Tong Deng,Xian-Tao Zeng

DOI: https://doi.org/10.21203/rs.2.18724/v1

2020-01-01

Abstract:Background: Many traditional survival analyses and restricted mean survival time (RMST) based analyses have been proposed for dealing with right censored time-to-event data. It is necessary to sort out the conditions and relationship among these methods for instruction. Comparison between hazard ratio (HR) and RMST in a study may promote appropriate understanding and application of the two effect measures.Method : Traditional survival analyses and RMST-based analyses were integrated into a flowchart and applied for univariable and multivariable analyses, RMST-based analyses included group comparison for difference in RMST, pseudo-value (PV) regressions, inverse probability of censoring probability (IPCW) regressions with group-specific weights and individual weights with homogeneity of censoring mechanism assumption. 4405 osteosarcomas were captured from Surveillance, Epidemiology and End Results Program Database. HR and difference in RMST (survival time lost or gain, STL or STG) were considered as effect measures and the effect of all included covariates on overall survival were reported for comparison. The relationship between difference in RMST and HR were explored when proportional hazard (PH) assumption was and was not met, respectively. Results: In group comparison and univariable regressions, using difference in RMST calculated by Kaplan-Meier methods as reference, PV regressions (R2=0.99) and 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. Conclusion: The flowchart may be intuitive and helpful to instruct appropriate use of RMST based and traditional methods. PH assumption and homogeneity of censoring mechanism assumption will determine the appropriate selection of these method. HR and difference in RMST can be reported comprehensively.