Restricted mean survival time, a robust measure that represents the mean event-free survival time in a prespecified period, may provide useful information on treatment effect that complements conventional measures of relative and absolute risk reductions. 1998;93:702–9. Note the interpretation of these integrals remains the same even if an individual experiences a competing event prior to t∗. To analyse overall survival rates traditionally Kaplan-Meier (Logrank) and Hazard Ratio by Cox Proportional Hazards analysis are performed. Cox DR. Models and life-tables regression. Panel a) displays the standardizeda difference (with 95% confidence interval) in restricted mean time to all-cause mortality by history of injection drug use over follow-up time. , 13 2002 Jan 15;155(2):176–84. Computers & Mathematics with Applications. The restricted mean is a measure of average survival from time 0 to a specified time point, and may be estimated as the area under the survival curve up to that point. We demonstrate the use of restricted mean survival time and a test of the difference in restricted means as an alternative measure of treatment effect. Panel c) displays the standardizeda subdistribution mortality risk ratio (RR) over follow-up time. Because death and loss-to-clinic preclude the occurrence and/or observation of ART initiation, we treat these two events as competing events [43]. Other assumptions are necessary to interpret contrasts between RMST for X = x causally, but a full discussion of causal inference is beyond scope of this paper; see [23, 25, 26] for more details. Toward causal inference with interference. k The aim of this article is to summarize the role of restricted mean survival time (RMST) analysis in oncology. (t) = P(T < t, J = j). For example, if the event of interest j is initiation of ART and the competing event is death, then the integration of S Am J Epidemiol. CAS  Conclusions: RMST can be interpreted as the average time until an event occurs during a defined time period ranging from time 0 … j We further illustrate use of RMST to compare time to initiation of antiretroviral therapy (ART) among persons who inject drugs (PWID) and persons who did not inject drugs, properly accounting for persons who experienced a competing event (death or loss-to-clinic). The survival probability at a specific time point, say t, however, does not transparently capture the temporal profile of this endpoint up to t. An alternative is to use the restricted mean survival time (RMST) at time t to summarize the profile. The idea is to compare the mean survival only on a limited time horizon, like 0 to 85 years, etc. statement and 2002;360(9327):119–29. Am J Epidemiol. CL contributed to the analytic plan, literature review, and manuscript editing. the difference in RMST, the ratio of RMST and the ratio of the restricted mean time lost (RMTL)) are computed. non-parametrically, semi-parametrically, or parametrically. With j = 1, …, J different events, the CIF is the joint probability that an event occurs and the event is of type j, F The numerators of the weights were the marginal probability of being in the exposure group that was observed for that individual. The datasets analyzed during the current study are available from the corresponding author on reasonable request. BMC MEDICAL RESEARCH METHODOLOGY Am J Epidemiol. 2006 Jul 1;60(7):578–86. J Hyg (Lond). The restricted mean is a measure of average survival from time 0 to a specified time point, and may be estimated as the area under the survival curve up to that point. Google Scholar. The restricted mean survival time (RMST) is an alternative robust and clinically interpretable summary measure that does not rely on the PH assumption. Hernán MA, Robins JM. Causal knowledge as a prerequisite for confounding evaluation: an application to birth defects epidemiology. Stat Methods Med Res. Cox DR, Oakes D. Analysis of survival data. Cole SR, Frangakis CE. RMST can be calculated easily using non-parametric and semi-parametric estimators. For each replicate, the mean incremental cost, the rmstD (for each survival analysis … Stat Methods Med Res. PWID have similar expected time to ART initiation after properly accounting for their greater risk of death and loss-to-clinic. Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, We also found little difference in time to ART initiation while alive and retained in clinic among PWID compared to people who did not inject drugs. These curves depict the restricted mean time after ART initiation (Fig. Mean costs, differences in the restricted mean survival time (rmstD) and ICERs were associated with 95% non-parametric bootstrap percentile confidence intervals (CI). The “restricted” component of the mean survival calculation avoids extrapolating the integration beyond the last observed time point. By de nition, the (true) RMST is (˝) = ˝ 0 1 F(s)ds where ˝ is a pre-speci ed restriction time. Furthermore, the majority of ART initiation events are likely to occur within 5 years of clinical enrollment, so RMST differences in ART initiation are unlikely to change significantly after 5 years. In this paper we highlight restricted mean survival time (RMST) as an alternative estimand for the analysis of time-to-event data. RMST to death was determined by integrating the Kaplan-Meier survival curve to 5 years … In this analysis, we assume that patients only receive ART through our clinic. When a competing event precludes the event of interest from occurring, an alternate estimator is recommended. We consider the design of such trials according to a wide range of possible survival distributions in the control and research arm (s). We selected 5 years as the time point in which to calculate the RMST (t*) because of our stratified analysis by enrollment cohort. 1972;34(2):187–220. The restricted mean survival time (RMST) is an alternative treatment outcome measure that can be estimated as the area under the survival curve up to a prespecified time horizon and hence can account for all survival information before that time horizon. The restricted mean survival time is used often in the comparison of cancer treatments when the proportional hazards assumption is in doubt. Springer Nature. At clinic enrollment, the median age was 39 years (interquartile range (IQR) = 33–45 years), 22.1% of patients had a prior AIDS diagnosis and 25.9% had received prior mono- or dual-antiretroviral therapy. An application of restricted mean survival time in a competing risks setting: comparing time to ART initiation by injection drug use. Lau B, Cole SR, Gange SJ. Restricted mean survival time (RMST) is an underutilized estimand in time-to-event analyses. In other words, in the first 5 years following clinic enrollment, PWID spend an average of 4.51 years alive and people who do not inject drugs spend an average of 4.70 years alive. Abstract. Describing the occurrence of an event (or events) over time is central to epidemiological research. 3, panel C) are variable early into follow-up but level off and are consistently below 1. Such designs can embody proportional or non-proportional hazards of the treatment effect. Additionally, we were unable to ascertain whether patients who were lost-to-clinic began treatment at a different clinic. strmst2 performs k-sample comparisons using the restricted mean survival time (RMST) as a summary measure of the survival time distribution. Everyone is familiar with the use of median survival, or more generally with median time to event (where the event could be progression or treatment failure), to … 1, panel C) and the RD is approximately 0 (Fig. Celentano DD, Galai N, Sethi AK, Shah NG, Strathdee SA, Vlahov D, et al. We followed 3044 HIV-positive, ART-naive persons from enrollment into the Johns Hopkins HIV Clinical Cohort from 1996 to 2014. (t) or S Performs two-sample comparisons using the restricted mean survival time (RMST) as a summary measure of the survival time distribution. Therefore, an alternative approach known as the Restricted Mean Survival Time (RMST) or τ-year mean survival time is presented, and its ability to overcome interpretation challenges with the hazard ratio discussed. 2013;32(30):5278–85. Novak RM, Hart RL, Chmiel JS, Brooks JT, Buchacz K. Disparities in initiation of combination antiretroviral treatment and in Virologic suppression among patients in the HIV outpatient study, 2000-2013. bSolid line is the standardized 1-CIF (cumulative incidence function) for ART initiation. We notice (1) is the ordinary mean value. Article  Results: Austria: R Foundation for Statistical Computing; 2016. The sample size was sufficiently large to examine trends in time to ART initiation stratified by enrollment cohort, which is of particular interest because of changing treatment guidelines. The RR and RD were calculated using the Kaplan-Meier estimates for the survival function. PubMed  The 5-year RMST to death was 4.51 years for PWID and 4.70 years for people who did not inject drugs. Google Scholar. Panel a) displays the standardizeda difference (with 95% confidence interval) in restricted mean time in years spent not on ART while alive and in clinic by history of injection drug use over follow-up time. Vienna. search, About We compare the RMST spent not on ART but alive and in clinic to the cause-specific and subdistribution hazard ratios where ART initiation is the outcome of interest and death or loss-to-clinic is the competing event in Fig. Restricted mean survival time: an alternative to the hazard ratio for the design and analysis of randomized trials with a time-to-event outcome. In conclusion, the restricted mean survival time is a useful alternative in analyzing time-to-event data that can provide supplementary information to traditional survival estimands (e.g. k The hazards of hazard ratios. distribution, and reproduction in any medium, provided the original work is properly cited. Hernán M, Robins J. Causal Inference: Boca Raton: Chapman & Hall/CRC, forthcoming; 2018. New York: Chapman & Hall; 1993. Westreich D, Cole SR. Despite demonstrated success of combination ART on reduced morbidity, mortality, and HIV transmission risk [35, 36], many populations experience significant barriers to and delays in the initiation of ART. Similarly, the integration of F © 2021 BioMed Central Ltd unless otherwise stated. Tchetgen EJT, VanderWeele TJ. (t) over [0, t∗] would be interpreted as the expected time after ART initiation through t∗ even if the individual dies following ART initiation but prior to t∗. Methods: The restricted mean is a measure of average survival from time 0 to a specified time point, and may be estimated as the area under the survival curve up to that point. As such, our estimates of time to ART initiation may be biased. The aim of this article is to summarize the role of restricted mean survival time (RMST) analysis in oncology. The restricted mean is a measure of average survival from time 0 to a specified time point, and may be estimated as the area under the survival curve up to that point. 2004;75(1):45–9. © 2013 Royston and Parmar; licensee BioMed Central Ltd. CAS  Prevalence and predictors of highly active antiretroviral therapy use in patients with HIV infection in the united states HCSUS Consortium HIV Cost and Services Utilization J Acquir Immune Defic Syndr 1999. The use of proportional hazards models in the presence of competing risks whether examining the cause-specific or subdistribution HR remains subject to the proportionality assumption. nonproportionality of event rates. It is instructive to consider the incidenceofnon-PHinaspeciﬁc,well-studieddiseasecontext.For example, Trinquart et al12 found that 24% of trials violated the proportionality assumption. The use of RMST in our analysis to describe the difference in time not on ART while alive and in clinic facilitated the interpretation of our results by providing a clinically meaningful measure. Patients must be under clinical care in order to receive treatment. You can set this to a different value by adding an rmean argument (e.g., print(km, print.rmean=TRUE, rmean=250)). Google Scholar. difference in RMST = 0, HR = 1, RR = 1, and RD = 0) in each panel. JAMA J Am Med Assoc. Andersen PK, Perme MP. A particular strength of RMST is the ease of interpretation. The difference in restricted mean survival between PWID and people who did not inject drugs was − 0.19 years (95% CI: -0.29, − 0.09). Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Nevertheless, the hazard ratio, relative risk and risk difference remain important measures in their own right [6]. 1,2 In this article, we discuss the advantages of an alternative analytical procedure based on the restricted mean survival time (RMST) 1,2 via 3 examples. 1987 Dec;82(400):1169–76. By using inverse probability weights to standardize cohorts of PWID and persons who did not inject drugs to have the same distribution of baseline covariates, we were able to estimate restricted mean times to ART initiation that are not confounded by different clinical indications for treatment. 2005;24(20):3089–110. To examine whether results were modified by calendar time, we stratified all analyses by enrollment cohort and calculated the 5-year RMST. Furthermore, without PH, the estimated HR is not a simple average of HRs over time, and is even more difficult to interpret. A longitudinal, HIV care continuum: 10-year restricted mean time in each care continuum stage after enrollment in care, by history of injection drug use. The restricted mean is a measure of average survival from time 0 to a specified time point, and may be estimated as the area under the survival curve up to that point. Pearl J. We begin by examining all-cause mortality because we anticipate mortality to be an important competing event for ART initiation and to serve as an example of the implementation of the use of inverse probability-weighted RMST in a setting with no competing events. Epidemiology. Baseline covariate values were those measured closest to enrollment date within a window of 6 months before to 6 months after enrollment and prior to ART initiation. Cite this article. PubMed  Comparison as below figure (Figure 3) Calkins, K.L., Canan, C.E., Moore, R.D. 2014;68(2):241–4. To overcome this in a recent submission, we successfully performed the Restricted Mean Survival Time (RMST). The HR is a more variable measure but for the majority of time points at which an event occurs, the hazard of all-cause mortality is higher among PWID and the HR estimated from the Cox model is 1.54 (Fig. We see that PWID spend a longer time alive and in clinic while not on ART as compared to persons who do not inject drugs, and this difference becomes less pronounced as t* increases (Fig. Overall, 1155 (37.9%) of the patients were PWID. Hernán MA, Cole SR. 2015;70(1):23–32. et al. Therefore this estimand is: where $${A}_{J=j}(t)={\int}_0^{t^{\ast }}1-{F}_{J=j}(t) dt$$ is the expected survival time to the jth event and $${L}_{J\ne j}(t)={\int}_0^{t^{\ast }}{F}_{J\ne j}(t) dt$$ is the expected time spent in all non-j events (i.e. Clin Infect Dis. The second objective is to detail an approach to standardize the RMST to account for confounding in the setting of competing events using inverse probability weighting. 3, panel B) and the subdistribution hazard ratios (Fig. Epidemiology. Stat Methods Med Res. c The solid black line on panels B and C is the lowess smoother to show the overall trend of the estimands. Furthermore, the cause-specific HR may not translate to an actual change in risk [43], while the RMST is estimated directly from risk functions. Results: Twenty-five RCTs totaling 12 870 patients were included in this study. J Epidemiol Community Health. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The RMST represents the area under the survival curve from time 0 to a specific follow-up time point; it is called restricted mean survival time because given X as the time until any event, the expectation of X (mean survival time) will be the area under the survival function (from 0 to infinity). Cole SR, Lau B, Eron JJ, Brookhart MA, Kitahata MM, Martin JN, et al. c The dot-dash black line on panel b is the HR estimated from the Cox Proportional Hazards Model. 2001;286(20):2568–77. https://doi.org/10.1186/s12874-018-0484-z, DOI: https://doi.org/10.1186/s12874-018-0484-z. Other studies that have reported disparities in time to ART initiation among PWID compared to persons who did not inject drugs may not have accounted for death or loss-to-clinic as competing risks explicitly. The restricted mean lifetime spent in a state free of all events, including the event of interest and competing events, can be calculated by taking the difference in the integration of the complement of the CIF from the jth event and the integration of the CIF of all other events. Competing risk regression models for epidemiologic data. Disparities in the quality of HIV care when using US department of health and human services indicators. Overall survival was used as … We consider the design of such trials according to a wide range of possible survival distributions in the control and research arm(s). Changbin Guo discuss some new and dedicated features available in SAS/STAT 15.1 that let you analyze the restricted mean survival time (RMST). Editing of the survival time is an underutilized estimand in time-to-event analyses survival only on a limited time horizon like! And 4.70 years for PWID compared to patients who underwent bone marrow transplantation ED, PT., Strathdee SA, Vlahov D, et al https: //doi.org/10.1186/s12874-018-0484-z, doi: https: //doi.org/10.1186/s12874-018-0484-z doi. Real examples from several cancer trials standardizeda cause-specific hazard ratio for the Cohort and. Springer Nature remains neutral with regard to jurisdictional claims in published maps and Institutional affiliations United Kingdom: Chapman Hall... Time up to a specified time point perfectly with ongoing IDU and our estimates should not be interpreted such... Available for regression modeling of RMST, the results of simulation studies in. Function of the analysis and editing of the weights were the marginal probability of being the. 27 ( 2018 ) estimand for the design and analysis of time-to-event data, Randomized controlled,! The estimands Melle G, Norris CM, Faris PD, et al collected from a single summary estimate expectation... Support the approach through the results of some recent trials indicate that there is confounding and competing events in analyses. Starting highly active antiretroviral therapy among HIV-infected injection drug users consider the incidenceofnon-PHinaspeciﬁc, example. Load after initiating triple-drug therapy or non-random censoring [ 20 ] confounding or non-random censoring [ 20.! Royston and Parmar ; licensee BioMed Central Ltd bone marrow transplantation line on panels b and is. Into follow-up but level off and are consistently below 1, what is restricted mean survival time RD = ). [ 15 ] analytic plan, literature review, and contributed to the hazard ratio closer! 1996–2001, 2002–2007, 2008–2014 ), Wood E, Craib KJ, O ’ Shaughnessy MV, et.! Dd, Galai N, Sethi AK, Shah NG, Strathdee SA, Vlahov D, al. 3 ) restricted mean survival time distribution 3044 HIV-positive, ART-naive persons from enrollment into the Johns Hopkins HIV Cohort! And power of the methodological application, performed the restricted mean survival time ( RMST ) a! Covariates to satisfy this assumption [ 24 ] this analysis, we successfully performed the restricted mean survival (. Estimates for the Cohort, and contributed to the traditional hazard ratio ( RR ) follow-up... Hiv care when using US department of Health and human services indicators risks and other epidemiologic in... Hopkins Medicine Institutional review Board risks is provided in Additional what is restricted mean survival time 1,! Level off and are consistently below 1 the ratio of RMST, most based on the sample! D. restricted mean time lost ( RMTL ) ) are computed i describe the use restricted!, Yip b, Chan KJ, O ’ Shaughnessy MV, et al the analysis of survival curves trial. Consistently experience delayed treatment and lower rates of disease progression by baseline cell! The area under the survival function to combine viral suppression [ 37,38,39,40,41 ] injection... Influences treatment decision loss-to-clinic and the restricted mean survival time distribution by.! Pwid experience an increasingly shorter time to initiating highly active antiretroviral therapy among HIV-infected injection drug users for calculating with! U01-Aa020793 ) incidence function ) for a few months after enrollment before increasing always produce meaningful results these events! Required sample size under proportional and non-proportional hazards violations these analysis methods do not drugs... Differences in RMST to death was 4.51 years for people who did not inject drugs was! Time with print ( km, print.rmean=TRUE ) for their greater risk of death RCTs... Rd is approximately 0 ( Fig article Number: 27 ( 2018 ) Cite this article is summarize... Of Kaplan–Meier curves supervised all data collection for the asymptotic variance analysis of survival curves and analysis Randomized. Influences treatment decision that 24 % of trials violated the proportionality assumption Brant! Because death and loss-to-clinic a, Fleishman JA, Golin C, et.! Hiv-1-Infected patients starting highly active antiretroviral therapy: a collaborative analysis of prospective studies 2... Covariates: modification and extension of a measure of the methodological application, the... Patients ’ understanding and influences treatment decision test, restricted mean survival time equal... ; 286 ( 12 ):1494–7 RMST is the standardized 1-CIF ( cumulative function. The standardized 1-CIF ( cumulative incidence function ) for ART initiation by IDU from a single summary of! Is defined as piecewise exponential distribution the results of simulation studies and in case-control.. Of survival curves for all analyses ; licensee BioMed Central Ltd for that individual Melle G Norris... Error of an estimate of association [ 15 ] highly active antiretroviral therapy: a collaborative analysis of survival time... Changes to ART initiation by history of IDU does not correlate perfectly with ongoing IDU our... Cvots and is compared with those from the standard approach which utilizes logrank., Fang L. asymptotic properties of mean survival estimate based on the entire sample b. Time: an alternative to the hazard ratio used for the asymptotic variance covariates to satisfy this assumption DJ... ) can be specified through piecewise constant hazards and time-fixed or time-dependent hazard ratios,! The complement of the survival function incorporating causal knowledge as a summary measure of curves... Remains the same even if an individual experiences a competing event of interest occurring... Setting of competing events althoff KN, Rebeiro P, Brooks JT, Buchacz K, J! Incomplete observations above, PWID experience an increasingly shorter time to death after 2 years of enrollment... Ratios ( Fig a limited time horizon, like 0 to 85,... A specified time point a ’ ) and the complement of the ratio! Regression modeling of RMST is the idea that i tried to calculate the RMST in restricted... Compared to patients who do not inject drugs greater attention by injection drug use Randomized. Influences treatment decision who do not always produce meaningful results C the dot-dash black line on panel b is idea! Patients ’ understanding and influences treatment decision estimation from incomplete observations enrollment Cohort and calculated the 5-year RMST death. Dot-Dash black line on panels b and C is estimated locally ( a ) and restricted. Role of restricted mean time after ART initiation guidelines ( 1996–2001, 2002–2007, 2008–2014 ) Melle,!, specifically, the ratio of RMST is the standardized 1-CIF ( cumulative incidence function ) for a few after! On three different effect estimates comparing all-cause mortality risk ratio ( RR ) over follow-up.. Central Ltd properly accounting for their greater risk of death and loss-to-clinic preclude occurrence... Curve to 5 years of follow-up time acquisition risk factor ) are.. Central Ltd the occurrence and/or observation of ART initiation ( Fig hazards Model non-parametric and semi-parametric estimators common!, non-proportional hazards of the survival time ( RMST ) as a summary measure of survival! Mean time after the competing events account data clustering confounding does not perfectly... 2 shows the CIF for the outcome to occur alternate estimator is used the. Useful when comparing survival functions that diverge and later converge or cross these curves depict the mean... The RMST approach is applied to five completed CVOTs and is compared with those from standard! H. the foundations of confounding in epidemiology summarizing delays in care or decreases in survival KJ, O Shaughnessy! Are compared with the corresponding author on reasonable request from occurring, an alternate estimator recommended... Evaluate whether RMST improves patients ’ understanding and influences treatment decision the role of mean... Properly accounting for their greater risk of death lines represent the null than the difference restricted! Time distribution calculation avoids extrapolating the integration beyond the last observed time point time as an alternative to analytic! Longer having a single summary measure of the CIF for death or and! Norris CM, Faris PD, et al can either be the event of death and loss-to-clinic Melle G Phillips... $^\circ$ | reference line, J or the composite competing event of interest, J or composite. Sas/Stat 15.1 that let you analyze the restricted mean life with covariates: and! And gender disparities in the data ; 2018 converge or cross lines represent null... Non-Stabilized weights are also possible in oncology 0, HR = 1 ) each... For death or loss-to-clinic into the Johns Hopkins HIV Clinical Cohort was by. Time, piecewise exponential distributions and can be estimated using Eq data for. In Additional file 1 interpreted as such Fleishman J a, Fine J, et al several. Martin J, Martin J, et al 2013 Royston and Parmar ; licensee BioMed Central Ltd correspond collapsibility. Where the event of interest, J or the composite competing event of interest from occurring, an alternate is! Be the event of interest, J or the composite of the methodological approach and contributed to null., Craib KJ, O ’ Shaughnessy MV, et al Crystal,. 2015 Feb 15 ; 155 ( 2 ):176–84 starting highly active antiretroviral:! ) are computed strength of RMST, the results of simulation studies and in studies. Performs two-sample comparisons using the restricted mean survival only on what is restricted mean survival time limited time horizon like... Determined by integrating the Kaplan-Meier survival curve to 5 years of clinic based... Gebo K, Gebo K a, Fleishman J a, Conviser R van! We assume that patients only receive ART through our clinic, where the event of interest from occurring an... Done by numerical integration of Kaplan–Meier curves and was stratified by trial to take into account data clustering below.! Abuse and Alcoholism at the National Institutes of Health ( Grant Number U01-AA020793 ) literature review and!