1
|
Zhang Y, Schaubel DE. Semiparametric Additive Modeling of the Restricted Mean Survival Time. Biom J 2024; 66:e202200371. [PMID: 39149839 DOI: 10.1002/bimj.202200371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2022] [Revised: 04/06/2024] [Accepted: 04/15/2024] [Indexed: 08/17/2024]
Abstract
Analysis of the restricted mean survival time (RMST) has become increasingly common in biomedical studies during the last decade as a means of estimating treatment or covariate effects on survival. Advantages of RMST over the hazard ratio (HR) include increased interpretability and lack of reliance on the often tenuous proportional hazards assumption. Some authors have argued that RMST regression should generally be the frontline analysis as opposed to methods based on counting process increments. However, in order for the use of the RMST to be more mainstream, it is necessary to broaden the range of data structures to which pertinent methods can be applied. In this report, we address this issue from two angles. First, most of existing methodological development for directly modeling RMST has focused on multiplicative models. An additive model may be preferred due to goodness of fit and/or parameter interpretation. Second, many settings encountered nowadays feature high-dimensional categorical (nuisance) covariates, for which parameter estimation is best avoided. Motivated by these considerations, we propose stratified additive models for direct RMST analysis. The proposed methods feature additive covariate effects. Moreover, nuisance factors can be factored out of the estimation, akin to stratification in Cox regression, such that focus can be appropriately awarded to the parameters of chief interest. Large-sample properties of the proposed estimators are derived, and a simulation study is performed to assess finite-sample performance. In addition, we provide techniques for evaluating a fitted model with respect to risk discrimination and predictive accuracy. The proposed methods are then applied to liver transplant data to estimate the effects of donor characteristics on posttransplant survival time.
Collapse
Affiliation(s)
- Yuan Zhang
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Douglas E Schaubel
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| |
Collapse
|
2
|
Yu Z, Geng X, Li Z, Zhang C, Hou Y, Zhou D, Chen Z. Time-varying effect in older patients with early-stage breast cancer: a model considering the competing risks based on a time scale. Front Oncol 2024; 14:1352111. [PMID: 39015489 PMCID: PMC11249566 DOI: 10.3389/fonc.2024.1352111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 06/10/2024] [Indexed: 07/18/2024] Open
Abstract
Background Patients with early-stage breast cancer may have a higher risk of dying from other diseases, making a competing risks model more appropriate. Considering subdistribution hazard ratio, which is used often, limited to model assumptions and clinical interpretation, we aimed to quantify the effects of prognostic factors by an absolute indicator, the difference in restricted mean time lost (RMTL), which is more intuitive. Additionally, prognostic factors of breast cancer may have dynamic effects (time-varying effects) in long-term follow-up. However, existing competing risks regression models only provide a static view of covariate effects, leading to a distorted assessment of the prognostic factor. Methods To address this issue, we proposed a dynamic effect RMTL regression that can explore the between-group cumulative difference in mean life lost over a period of time and obtain the real-time effect by the speed of accumulation, as well as personalized predictions on a time scale. Results A simulation validated the accuracy of the coefficient estimates in the proposed regression. Applying this model to an older early-stage breast cancer cohort, it was found that 1) the protective effects of positive estrogen receptor and chemotherapy decreased over time; 2) the protective effect of breast-conserving surgery increased over time; and 3) the deleterious effects of stage T2, stage N2, and histologic grade II cancer increased over time. Moreover, from the view of prediction, the mean C-index in external validation reached 0.78. Conclusion Dynamic effect RMTL regression can analyze both dynamic cumulative effects and real-time effects of covariates, providing a more comprehensive prognosis and better prediction when competing risks exist.
Collapse
Affiliation(s)
- Zhiyin Yu
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Xiang Geng
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Zhaojin Li
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Chengfeng Zhang
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Yawen Hou
- Department of Statistics and Data Science, School of Economics, Jinan University, Guangzhou, China
| | - Derun Zhou
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Zheng Chen
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| |
Collapse
|
3
|
Bai F, Yang X. Semiparametric estimation of restricted mean survival time as a function of restriction time. Stat Med 2023; 42:5389-5404. [PMID: 37737510 DOI: 10.1002/sim.9918] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 08/28/2023] [Accepted: 09/13/2023] [Indexed: 09/23/2023]
Abstract
The restricted mean survival time (RMST) is an appealing measurement in clinical or epidemiological studies with censored survival outcome and receives a lot of attention in the past decades. It provides a useful alternative to the Cox model for evaluating the covariate effect on survival time. The covariate effect on RMST usually varies with the restriction time. However, existing methods cannot address this problem properly. In this article, we propose a semiparametric framework that directly models RMST as a function of the restriction time. Our proposed model adopts a widely-used proportional form, enabling the estimation of RMST predictions across an interval using a unified model. Furthermore, the covariate effect for multiple restriction time points can be derived simultaneously. We develop estimators based on estimating equations theories and establish the asymptotic properties of the proposed estimators. The finite sample properties of the estimators are evaluated through extensive simulation studies. We further illustrate the application of our proposed method through the analysis of two real data examples. Supplementary Material are available online.
Collapse
Affiliation(s)
- Fangfang Bai
- School of Statistics, University of International Business and Economics, Beijing, China
| | - Xiaoran Yang
- School of Statistics, University of International Business and Economics, Beijing, China
| |
Collapse
|
4
|
Yang Z, Zhang C, Hou Y, Chen Z. Analysis of dynamic restricted mean survival time based on pseudo-observations. Biometrics 2023; 79:3690-3700. [PMID: 37337620 DOI: 10.1111/biom.13891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 05/07/2023] [Accepted: 06/01/2023] [Indexed: 06/21/2023]
Abstract
In clinical follow-up studies with a time-to-event end point, the difference in the restricted mean survival time (RMST) is a suitable substitute for the hazard ratio (HR). However, the RMST only measures the survival of patients over a period of time from the baseline and cannot reflect changes in life expectancy over time. Based on the RMST, we study the conditional restricted mean survival time (cRMST) by estimating life expectancy in the future according to the time that patients have survived, reflecting the dynamic survival status of patients during follow-up. In this paper, we introduce the estimation method of cRMST based on pseudo-observations, the statistical inference concerning the difference between two cRMSTs (cRMSTd), and the establishment of the robust dynamic prediction model using the landmark method. Simulation studies are conducted to evaluate the statistical properties of these methods. The results indicate that the estimation of the cRMST is accurate, and the dynamic RMST model has high accuracy in coefficient estimation and good predictive performance. In addition, an example of patients with chronic kidney disease who received renal transplantations is employed to illustrate that the dynamic RMST model can predict patients' expected survival times from any prediction time, considering the time-dependent covariates and time-varying effects of covariates.
Collapse
Affiliation(s)
- Zijing Yang
- Stomatological Hospital, School of Stomatology, Southern Medical University, Guangzhou, China
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Chengfeng Zhang
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| | - Yawen Hou
- Department of Statistics and Data Science, School of Economics, Jinan University, Guangzhou, China
| | - Zheng Chen
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, China
| |
Collapse
|
5
|
Vilain-Abraham FL, Tavernier E, Dantan E, Desmée S, Caille A. Restricted mean survival time to estimate an intervention effect in a cluster randomized trial. Stat Methods Med Res 2023; 32:2016-2032. [PMID: 37559486 DOI: 10.1177/09622802231192960] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/11/2023]
Abstract
For time-to-event outcomes, the difference in restricted mean survival time is a measure of the intervention effect, an alternative to the hazard ratio, corresponding to the expected survival duration gain due to the intervention up to a predefined time t*. We extended two existing approaches of restricted mean survival time estimation for independent data to clustered data in the framework of cluster randomized trials: one based on the direct integration of Kaplan-Meier curves and the other based on pseudo-values regression. Then, we conducted a simulation study to assess and compare the statistical performance of the proposed methods, varying the number and size of clusters, the degree of clustering, and the magnitude of the intervention effect under proportional and non-proportional hazards assumption. We found that the extended methods well estimated the variance and controlled the type I error if there was a sufficient number of clusters (≥ 50) under both proportional and non-proportional hazards assumption. For cluster randomized trials with a limited number of clusters (< 50), a permutation test for pseudo-values regression was implemented and corrected the type I error. We also provided a procedure to estimate permutation-based confidence intervals which produced adequate coverage. All the extended methods performed similarly, but the pseudo-values regression offered the possibility to adjust for covariates. Finally, we illustrated each considered method with a cluster randomized trial evaluating the effectiveness of an asthma-control education program.
Collapse
Affiliation(s)
| | - Elsa Tavernier
- INSERM, SPHERE, U1246, Tours University, Nantes University, Tours, France
| | - Etienne Dantan
- INSERM, SPHERE, U1246, Nantes University, Tours University, Nantes, France
| | - Solène Desmée
- INSERM, SPHERE, U1246, Tours University, Nantes University, Tours, France
| | - Agnès Caille
- INSERM, SPHERE, U1246, Tours University, Nantes University, Tours, France
| |
Collapse
|
6
|
Zhang C, Li Z, Yang Z, Huang B, Hou Y, Chen Z. A Dynamic Prediction Model Supporting Individual Life Expectancy Prediction Based on Longitudinal Time-Dependent Covariates. IEEE J Biomed Health Inform 2023; 27:4623-4632. [PMID: 37471185 DOI: 10.1109/jbhi.2023.3292475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
In the field of clinical chronic diseases, common prediction results (such as survival rate) and effect size hazard ratio (HR) are relative indicators, resulting in more abstract information. However, clinicians and patients are more interested in simple and intuitive concepts of (survival) time, such as how long a patient may live or how much longer a patient in a treatment group will live. In addition, due to the long follow-up time, resulting in generation of longitudinal time-dependent covariate information, patients are interested in how long they will survive at each follow-up visit. In this study, based on a time scale indicator-restricted mean survival time (RMST)-we proposed a dynamic RMST prediction model by considering longitudinal time-dependent covariates and utilizing joint model techniques. The model can describe the change trajectory of longitudinal time-dependent covariates and predict the average survival times of patients at different time points (such as follow-up visits). Simulation studies through Monte Carlo cross-validation showed that the dynamic RMST prediction model was superior to the static RMST model. In addition, the dynamic RMST prediction model was applied to a primary biliary cirrhosis (PBC) population to dynamically predict the average survival times of the patients, and the average C-index of the internal validation of the model reached 0.81, which was better than that of the static RMST regression. Therefore, the proposed dynamic RMST prediction model has better performance in prediction and can provide a scientific basis for clinicians and patients to make clinical decisions.
Collapse
|
7
|
Mao L. Nonparametric inference of general while-alive estimands for recurrent events. Biometrics 2023; 79:1749-1760. [PMID: 35731993 PMCID: PMC9772359 DOI: 10.1111/biom.13709] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 06/16/2022] [Indexed: 12/24/2022]
Abstract
Measuring the treatment effect on recurrent events like hospitalization in the presence of death has long challenged statisticians and clinicians alike. Traditional inference on the cumulative frequency unjustly penalizes survivorship as longer survivors also tend to experience more adverse events. Expanding a recently suggested idea of the "while-alive" event rate, we consider a general class of such estimands that adjust for the length of survival without losing causal interpretation. Given a user-specified loss function that allows for arbitrary weighting, we define as estimand the average loss experienced per unit time alive within a target period and use the ratio of this loss rate to measure the effect size. Scaling the loss rate by the width of the corresponding time window gives us an alternative, and sometimes more photogenic, way of showing the data. To make inferences, we construct a nonparametric estimator for the loss rate through the cumulative loss and the restricted mean survival time and derive its influence function in closed form for variance estimation and testing. As simulations and analysis of real data from a heart failure trial both show, the while-alive approach corrects for the false attenuation of treatment effect due to patients living longer under treatment, with increased statistical power as a result. The proposed methods are implemented in the R-package WA, which is publicly available from the Comprehensive R Archive Network (CRAN).
Collapse
Affiliation(s)
- Lu Mao
- Department of Biostatistics and Medical Informatics, University of Wisconsin-Madison, Madison, WI, 53792, USA
| |
Collapse
|
8
|
Zhang C, Huang B, Wu H, Yuan H, Hou Y, Chen Z. Restricted mean survival time regression model with time-dependent covariates. Stat Med 2022; 41:4081-4090. [PMID: 35746886 PMCID: PMC9545070 DOI: 10.1002/sim.9495] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2021] [Revised: 05/09/2022] [Accepted: 05/16/2022] [Indexed: 11/06/2022]
Abstract
In clinical or epidemiological follow‐up studies, methods based on time scale indicators such as the restricted mean survival time (RMST) have been developed to some extent. Compared with traditional hazard rate indicator system methods, the RMST is easier to interpret and does not require the proportional hazard assumption. To date, regression models based on the RMST are indirect or direct models of the RMST and baseline covariates. However, time‐dependent covariates are becoming increasingly common in follow‐up studies. Based on the inverse probability of censoring weighting (IPCW) method, we developed a regression model of the RMST and time‐dependent covariates. Through Monte Carlo simulation, we verified the estimation performance of the regression parameters of the proposed model. Compared with the time‐dependent Cox model and the fixed (baseline) covariate RMST model, the time‐dependent RMST model has a better prediction ability. Finally, an example of heart transplantation was used to verify the above conclusions.
Collapse
Affiliation(s)
- Chengfeng Zhang
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, People's Republic of China
| | - Baoyi Huang
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, People's Republic of China
| | - Hongji Wu
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, People's Republic of China
| | - Hao Yuan
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, People's Republic of China
| | - Yawen Hou
- Department of Statistics, School of Economics, Jinan University, Guangzhou, People's Republic of China
| | - Zheng Chen
- Department of Biostatistics, School of Public Health (Guangdong Provincial Key Laboratory of Tropical Disease Research), Southern Medical University, Guangzhou, People's Republic of China
| |
Collapse
|
9
|
Chen X, Harhay MO, Li F. Clustered restricted mean survival time regression. Biom J 2022. [PMID: 35593026 DOI: 10.1002/bimj.202200002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Revised: 03/23/2022] [Accepted: 04/18/2022] [Indexed: 11/05/2022]
Abstract
For multicenter randomized trials or multilevel observational studies, the Cox regression model has long been the primary approach to study the effects of covariates on time-to-event outcomes. A critical assumption of the Cox model is the proportionality of the hazard functions for modeled covariates, violations of which can result in ambiguous interpretations of the hazard ratio estimates. To address this issue, the restricted mean survival time (RMST), defined as the mean survival time up to a fixed time in a target population, has been recommended as a model-free target parameter. In this article, we generalize the RMST regression model to clustered data by directly modeling the RMST as a continuous function of restriction times with covariates while properly accounting for within-cluster correlations to achieve valid inference. The proposed method estimates regression coefficients via weighted generalized estimating equations, coupled with a cluster-robust sandwich variance estimator to achieve asymptotically valid inference with a sufficient number of clusters. In small-sample scenarios where a limited number of clusters are available, however, the proposed sandwich variance estimator can exhibit negative bias in capturing the variability of regression coefficient estimates. To overcome this limitation, we further propose and examine bias-corrected sandwich variance estimators to reduce the negative bias of the cluster-robust sandwich variance estimator. We study the finite-sample operating characteristics of proposed methods through simulations and reanalyze two multicenter randomized trials.
Collapse
Affiliation(s)
- Xinyuan Chen
- Department of Mathematics and Statistics, Mississippi State University, Mississippi State, MS, USA
| | - Michael O Harhay
- Clinical Trials Methods and Outcomes Lab, PAIR (Palliative and Advanced Illness Research) Center, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA.,Department of Biostatistics, Epidemiology and Informatics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Fan Li
- Department of Biostatistics, Yale University School of Public Health, New Haven, CT, USA.,Center for Methods in Implementation and Prevention Science, Yale University, New Haven, CT, USA
| |
Collapse
|
10
|
Zhao L. Deep Neural Networks For Predicting Restricted Mean Survival Times. Bioinformatics 2021; 36:5672-5677. [PMID: 33399818 PMCID: PMC8023687 DOI: 10.1093/bioinformatics/btaa1082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 11/30/2020] [Accepted: 12/16/2020] [Indexed: 11/14/2022] Open
Abstract
Restricted mean survival time (RMST) is a useful summary measurement of the time-to-event data, and it has attracted great attention for its straightforward clinical interpretation. In this article, I propose a deep neural network model that directly relates the RMST to its baseline covariates for simultaneous prediction of RSMT at multiple times. Each subject's survival time is transformed into a series of jackknife pseudo observations and then used as quantitative response variables in a deep neural network model. By using the pseudo values, a complex survival analysis is reduced to a standard regression problem, which greatly simplifies the neural network construction. By jointly modelling RMST at multiple times, the neural network model gains prediction accuracy by information sharing across times. The proposed network model was evaluated by extensive simulation studies and was further illustrated on three real datasets. In real data analyses, I also used methods to open the blackbox by identifying subject-specific predictors and their importance in contributing to the risk prediction. AVAILABILITY AND IMPLEMENTATION The source code is freely available at http://github.com/lilizhaoUM/DnnRMST. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Lili Zhao
- Department of Biostatistics, School of Public Health, University of Michigan, Ann Arbor, MI 48105, USA
- To whom correspondence should be addressed.
| |
Collapse
|