A Novel Stratified Analysis Method for Testing and Estimating Overall Treatment Effects on Time-to-Event Outcomes Using Average Hazard with Survival Weight

Zihan Qian,Lu Tian,Miki Horiguchi,Hajime Uno
DOI: https://doi.org/10.48550/arXiv.2404.00788
2024-04-01
Abstract:Given the limitations of using the Cox hazard ratio to summarize the magnitude of the treatment effect, alternative measures that do not have these limitations are gaining attention. One of the recently proposed alternative methods uses the average hazard with survival weight (AH). This population quantity can be interpreted as the average intensity of the event occurrence in a given time window that does not involve study-specific censoring. Inference procedures for the ratio of AH and difference in AH have already been proposed in simple randomized controlled trial settings to compare two groups. However, methods with stratification factors have not been well discussed, although stratified analysis is often used in practice to adjust for confounding factors and increase the power to detect a between-group difference. The conventional stratified analysis or meta-analysis approach, which integrates stratum-specific treatment effects using an optimal weight, directly applies to the ratio of AH and difference in AH. However, this conventional approach has significant limitations similar to the Cochran-Mantel-Haenszel method for a binary outcome and the stratified Cox procedure for a time-to-event outcome. To address this, we propose a new stratified analysis method for AH using standardization. With the proposed method, one can summarize the between-group treatment effect in both absolute difference and relative terms, adjusting for stratification factors. This can be a valuable alternative to the traditional stratified Cox procedure to estimate and report the magnitude of the treatment effect on time-to-event outcomes using hazard.
Methodology
What problem does this paper attempt to address?
The problems that this paper attempts to solve mainly focus on several key limitations encountered in the existing stratified analysis using the Cox proportional - hazards model: 1. **Strong model assumptions**: The traditional stratified Cox regression model requires that the proportional - hazards assumption be maintained in each stratum, that is, the hazard ratio of the two treatment groups must remain constant over the follow - up time. If this assumption does not hold, the estimated Cox hazard ratio will be affected by the study - specific censoring time distribution, making the results difficult to interpret or generalize to future populations. 2. **Lack of within - group absolute risk**: Although this approach is advantageous from a statistical perspective, it poses challenges in practical applications. Because the effect of clinical interventions cannot be evaluated solely by the hazard ratio, but also needs to consider the absolute risk of the control group. Providing the absolute risk of each group is crucial for assisting clinical evaluation. However, it is usually difficult to convert the hazard ratio in the Cox proportional - hazards model into an absolute risk difference, except in specific cases. 3. **Homogeneity assumption of cross - stratum effects**: Many stratified analysis methods, including the well - known Cochran - Mantel - Haenszel method, are based on the assumption that each stratum estimates the same underlying quantity or effect. This means that the treatment effects (such as hazard ratio and odds ratio) remain consistent across all strata. However, this assumption does not always hold, which constitutes a significant limitation in interpreting the magnitude of the estimated treatment effect. To solve these problems, the authors propose a new stratified analysis method based on Average Hazard (AH). This method adjusts the stratification factors through standardization techniques and can simultaneously provide summaries of the absolute and relative differences between groups, thus providing a valuable alternative to the traditional stratified Cox process for estimating and reporting the magnitude of treatment effects for risk - based time - to - event outcomes. This method not only overcomes the limitations of traditional methods but also provides more robust and richer information on treatment effects, which is helpful for making risk - based treatment decisions.