scholarly journals Impact of a non-constant baseline hazard on detection of time-dependent treatment effects: a simulation study

2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Kim Jachno ◽  
Stephane Heritier ◽  
Rory Wolfe

Abstract Background Non-proportional hazards are common with time-to-event data but the majority of randomised clinical trials (RCTs) are designed and analysed using approaches which assume the treatment effect follows proportional hazards (PH). Recent advances in oncology treatments have identified two forms of non-PH of particular importance - a time lag until treatment becomes effective, and an early effect of treatment that ceases after a period of time. In sample size calculations for treatment effects on time-to-event outcomes where information is based on the number of events rather than the number of participants, there is crucial importance in correct specification of the baseline hazard rate amongst other considerations. Under PH, the shape of the baseline hazard has no effect on the resultant power and magnitude of treatment effects using standard analytical approaches. However, in a non-PH context the appropriateness of analytical approaches can depend on the shape of the underlying hazard. Methods A simulation study was undertaken to assess the impact of clinically plausible non-constant baseline hazard rates on the power, magnitude and coverage of commonly utilized regression-based measures of treatment effect and tests of survival curve difference for these two forms of non-PH used in RCTs with time-to-event outcomes. Results In the presence of even mild departures from PH, the power, average treatment effect size and coverage were adversely affected. Depending on the nature of the non-proportionality, non-constant event rates could further exacerbate or somewhat ameliorate the losses in power, treatment effect magnitude and coverage observed. No single summary measure of treatment effect was able to adequately describe the full extent of a potentially time-limited treatment benefit whilst maintaining power at nominal levels. Conclusions Our results show the increased importance of considering plausible potentially non-constant event rates when non-proportionality of treatment effects could be anticipated. In planning clinical trials with the potential for non-PH, even modest departures from an assumed constant baseline hazard could appreciably impact the power to detect treatment effects depending on the nature of the non-PH. Comprehensive analysis plans may be required to accommodate the description of time-dependent treatment effects.

2020 ◽  
Vol 29 (12) ◽  
pp. 3525-3532
Author(s):  
Thomas J Prior

Clinical trials in oncology often involve the statistical analysis of time-to-event data such as progression-free survival or overall survival to determine the benefit of a treatment or therapy. The log-rank test is commonly used to compare time-to-event data from two groups. The log-rank test is especially powerful when the two groups have proportional hazards. However, survival curves encountered in oncology studies that differ from one another do not always differ by having proportional hazards; in such instances, the log-rank test loses power, and the survival curves are said to have “non-proportional hazards”. This non-proportional hazards situation occurs for immunotherapies in oncology; immunotherapies often have a delayed treatment effect when compared to chemotherapy or radiation therapy. To correctly identify and deliver efficacious treatments to patients, it is important in oncology studies to have available a statistical test that can detect the difference in survival curves even in a non-proportional hazards situation such as one caused by delayed treatment effect. An attempt to address this need was the “max-combo” test, which was originally described only for a single analysis timepoint; this article generalizes that test to preserve type I error when there are one or more interim analyses, enabling efficacious treatments to be identified and made available to patients more rapidly.


2020 ◽  
Vol 41 (Supplement_2) ◽  
Author(s):  
C Perego ◽  
M Sbolli ◽  
C Specchia ◽  
C Oriecuia ◽  
G Peveri ◽  
...  

Abstract Background The hazard ratio (HR) is the most common measure used to quantify treatment effects in heart failure (HF) clinical trials. However, the HR is only valid when the proportional hazards assumption is plausible, and the HR may be difficult to interpret for clinicians and laypeople. Restricted mean survival time (RMST), defined as the average time-to-event before a specific timepoint, is an intuitive summary of group-wise survival. The difference between two RMSTs measures treatment effects without model assumptions and may communicate more clinically interpretable results. Purpose To evaluate statistical and clinical properties of RMST-based statistics applied to clinical trial data for treatments of HF with reduced ejection fraction. Methods Patient time-to-event data was reconstructed from the published primary and secondary outcome Kaplan-Meier curves from landmark HF clinical trials. We estimated the RMST-differences between treatment groups as a measure of treatment effect with published data, and compared statistical testing results and effect size values to HR analysis results. Results We analyzed 7 HF clinical trials, including data from a total of 27,845 patients (Table 1). RMST should be interpreted as the average number of months that the outcome is avoided over the study period. As examples: On average, treatment with enalapril for 12 months extended each patient's life by 2.2 months compared to placebo, and treatment with spironolactone for 34 months extended each patient's life by 2.2 months compared to placebo. Conclusions RMST-difference test statistic has identical statistical conclusions as HRs but provided an intuitive estimate of each treatment effect. RMST-based data can potentially be used to better communicate treatment effects to patients, to assist in patient-preference discussions and shared decision-making Funding Acknowledgement Type of funding source: None


Author(s):  
Sean Wharton ◽  
Arne Astrup ◽  
Lars Endahl ◽  
Michael E. J. Lean ◽  
Altynai Satylganova ◽  
...  

AbstractIn the approval process for new weight management therapies, regulators typically require estimates of effect size. Usually, as with other drug evaluations, the placebo-adjusted treatment effect (i.e., the difference between weight losses with pharmacotherapy and placebo, when given as an adjunct to lifestyle intervention) is provided from data in randomized clinical trials (RCTs). At first glance, this may seem appropriate and straightforward. However, weight loss is not a simple direct drug effect, but is also mediated by other factors such as changes in diet and physical activity. Interpreting observed differences between treatment arms in weight management RCTs can be challenging; intercurrent events that occur after treatment initiation may affect the interpretation of results at the end of treatment. Utilizing estimands helps to address these uncertainties and improve transparency in clinical trial reporting by better matching the treatment-effect estimates to the scientific and/or clinical questions of interest. Estimands aim to provide an indication of trial outcomes that might be expected in the same patients under different conditions. This article reviews how intercurrent events during weight management trials can influence placebo-adjusted treatment effects, depending on how they are accounted for and how missing data are handled. The most appropriate method for statistical analysis is also discussed, including assessment of the last observation carried forward approach, and more recent methods, such as multiple imputation and mixed models for repeated measures. The use of each of these approaches, and that of estimands, is discussed in the context of the SCALE phase 3a and 3b RCTs evaluating the effect of liraglutide 3.0 mg for the treatment of obesity.


2018 ◽  
Vol 15 (5) ◽  
pp. 489-498 ◽  
Author(s):  
Jennifer G Le-Rademacher ◽  
Ryan A Peterson ◽  
Terry M Therneau ◽  
Ben L Sanford ◽  
Richard M Stone ◽  
...  

Background/aims The goal of this article is to illustrate the utility of multi-state models in cancer clinical trials. Our specific aims are to describe multi-state models and how they differ from standard survival methods, to illustrate how multi-state models can facilitate deeper understanding of the treatment effect on multiple paths along the disease process that patients could experience in cancer clinical trials, to explain the differences between multi-state models and time-dependent Cox models, and to briefly describe available software to conduct such analyses. Methods Data from 717 newly diagnosed acute myeloid leukemia patients who enrolled in the CALGB 10603 trial were used as an illustrative example. The current probability-in-state was estimated using the Aalen–Johansen estimator. The restricted mean time in state was calculated as the area under the probability-in-state curves. Cox-type regression was used to evaluate the effect of midostaurin on the various clinical paths. Simulation was conducted using a newly constructed shiny application. All analyses were performed using the R software. Results Multi-state model analyses of CALGB 10603 suggested that the overall improvement in survival with midostaurin seen in the primary analysis possibly resulted from a higher complete remission rate in combination with a lower risk of relapse and of death after complete remission in patients treated with midostaurin. Simulation results, in a three-state illness-death without recovery model, demonstrate that multi-state models and time-dependent Cox models evaluate treatment effects from different frameworks. Conclusion Multi-state models allow detailed evaluation of treatment effects in complex clinical trial settings where patients can experience multiple paths between study enrollment and the final outcome. Multi-state models can be used as a complementary tool to standard survival analyses to provide deeper insights to the effects of treatment in trial settings with complex disease process.


Author(s):  
Suzanne Freeman ◽  
Nicola Cooper ◽  
Alex Sutton ◽  
Michael Crowther ◽  
James Carpenter ◽  
...  

IntroductionSynthesis of clinical effectiveness is a well-established component of health technology assessment (HTA) combining data from multiple trials to obtain an overall pooled estimate of clinical effectiveness, which may inform an associated economic evaluation. Time-to-event outcomes are often synthesized using effect measures from Cox proportional hazards models assuming a constant hazard ratio over time. However, where treatment effects vary over time an assumption of proportional hazards is not always valid. Several methods have been proposed for synthesizing time-to-event outcomes in the presence of non-proportional hazards. However, guidance on choosing between these methods and the implications for HTA is lacking.MethodsWe applied five methods for estimating treatment effects from time-to-event outcomes, which relax the proportional hazards assumption to a network of melanoma trials, reporting overall survival: restricted mean survival time, an accelerated failure time generalized gamma model, piecewise exponential, fractional polynomial and Royston-Parmar models. We conducted a simulation study to compare these five methods. Simulated individual patient data was generated from a mixture Weibull distribution assuming a treatment-time interaction. Each simulated meta-analysis consisted of five trials with varying numbers of patients and length of follow-up across trials. For each model fitted to each dataset, we calculated the restricted mean survival time at the end of observed follow-up and following extrapolation to a 20-year time horizon.ResultsAll models fitted the melanoma data reasonably well with some variation in the treatment rankings and differences in the survival curves. The simulation study demonstrated the potential for different conclusions from different modelling approaches.ConclusionsThe restricted mean survival time, generalized gamma, piecewise exponential, fractional polynomial and Royston-Parmar models can all accommodate non-proportional hazards and differing lengths of trial follow-up within an evidence synthesis of time-to-event outcomes. Further work is needed in this area to extend the simulation study to the network meta-analysis setting and provide guidance on the key considerations for informing model choice for the purposes of HTA.


Sign in / Sign up

Export Citation Format

Share Document