Introduction

Setup

Load the package hce and check the version:

library(hce)
packageVersion("hce")
#> [1] '0.6.5'

For citing the package run citation("hce") (Samvel B. Gasparyan 2024).

Definitions

Prioritized outcome composite endpoints (POCE) are a general class of endpoints combining different clinical outcomes of patients into a composite so as to preserve their different natures. A particular case of these endpoints is the hierarchical composite endpoint (HCE). It is evaluated in a fixed follow-up period and accounts for the patient’s clinically most important outcome for the analysis. HCEs are analyzed using win odds and other win statistics (Samvel B. Gasparyan et al. 2023).

Examples

Here we provide examples of HCE using in clinical trials from different therapeutic areas. General considerations for creating HCEs can be found in Samvel B. Gasparyan et al. (2022).

COVID-19

The DARE-19 (M. Kosiborod et al. 2021; M. N. Kosiborod et al. 2021) trial used an HCE to assess outcomes in patients hospitalized for COVID-19 and treated for 30 days. The COVID-19 HCE is presented below. It combines death, in hospital organ dysfunction events with clinical status at Day 30 for patients alive, still hospitalized but without previous organ dysfunction events, and hospital discharge as the most favorable outcome for patients discharging without organ dysfunction events and being alive at Day 30.

Below a higher category signifies a better outcome. Patients are ranked into one and only one category based on their clinically most severe event. For example, patients experiencing an in-hospital new or worsening organ dysfunction event then dying will be included in the category I.

#>   Order                                               Category
#> 1     I                                                  Death
#> 2    II More than one new or worsened organ dysfunction events
#> 3   III            One new or worsened organ dysfunction event
#> 4    IV          Hospitalized at the end of follow-up (Day 30)
#> 5     V                 Discharged from hospital before Day 30

Patients in the category I are compared using the timing of the event, with an earlier event being a worse outcome (are assigned a lower rank). Similarly, in the category III the timing of the event is used for ranking patients within this category. In the category II patients are compared using the number of events with a higher number signifying a worse outcome. Patients in the category IV - hospitalized at the end of follow-up without previous worsening events - are further ranked according to oxygen support requirements at the hospital (IV.1 on high flow oxygen devices, IV.2 requiring supplemental oxygen, IV.3 not requiring supplemental oxygen, with a higher rank being a better outcome). Patients in the category V are compared using the timing of the event, but, the hospital discharge being a favorable outcome, here the earlier event signifies a better outcome than the late event (reverse of the ranking in categories I and III).

The simplest case of a COVID-19 HCE is an endpoint with ordinal scale outcomes assessed at a given timepoint. The endpoint uses 1-8 categories for assessing the physical limitations of hospitalized patients with COVID-19 after 15 or 30 days of treatment. But unlike the DARE-19 HCE, within each category it does not use the timing of events to reduce the ties in a paiwrise comparison of patients in the active group with patients in the control group. See, for example, COVID-19 and COVID-19b datasest ordinal scale outcomes (Beigel et al. 2020).

table(COVID19)
#>          GROUP
#> TRTP        1   2   3   4   5   6   7   8
#>   Active   34  95  28  58  38  14 117 157
#>   Placebo  58 121  24  60  33   8 102 115

The function hce::summaryWO() provides the number of wins, losses, and ties by categories. We can calculate the probability of ties from the provided numbers.

COVID19HCE <- hce(GROUP = COVID19$GROUP, TRTP = COVID19$TRTP)
SUM <- summaryWO(COVID19HCE, ref = "Placebo")$summary
SUM$Ptie <- round(SUM$TIE/SUM$TOTAL, 2)
SUM
#>   TRTP    WIN   LOSS   TIE  TOTAL        WR        WO Ptie
#> 1    A 135744  97143 48974 281861 1.3973627 1.3173641 0.17
#> 2    P  97143 135744 48974 281861 0.7156338 0.7590916 0.17

Kidney HCE

The kidney HCE defined in Heerspink et al. (2023) has the following ordinal outcomes (for the review of the topic see Little et al. (2023)).

#>   Order                                      Category
#> 1     I                                         Death
#> 2    II            Dialysis or kidney transplantation
#> 3   III         Sustained GFR < 15 ml/min per 1.73 m2
#> 4    IV Sustained GFR decline from baseline of >= 57%
#> 5     V Sustained GFR decline from baseline of >= 50%
#> 6    VI Sustained GFR decline from baseline of >= 40%
#> 7   VII                          Individual GFR slope

The dataset KHCE contains data on a kidney HCE outcomes

dat <- KHCE
Order <- c("Death (adj)", "Chronic dialysis (adj) >=90 days", 
           "Sustained eGFR<15 (mL/min/1.73 m2)", "Sustained >=57% decline in eGFR", 
           "Sustained >=50% decline in eGFR", "Sustained >=40% decline in eGFR", "eGFR slope")   
dat$GROUP <- factor(dat$GROUP, levels = Order)
table(dat$GROUP, dat$TRTP)
#>                                     
#>                                        A   P
#>   Death (adj)                         40  50
#>   Chronic dialysis (adj) >=90 days    17  29
#>   Sustained eGFR<15 (mL/min/1.73 m2)  16  28
#>   Sustained >=57% decline in eGFR      2   9
#>   Sustained >=50% decline in eGFR      7  22
#>   Sustained >=40% decline in eGFR     36  34
#>   eGFR slope                         632 578

This dataset is derived from ADSL which contains baseline characteristics,

ADSL dataset structure
ID TRTPN EGFRBL STRATAN
1 1 44 1
2 1 31 3
3 1 28 4
4 1 61 4
5 1 37 4
6 2 35 1

ADLB laboratory measurements of kidney function,

ADLB dataset structure
ID TRTPN AVAL ADAY AVISITN PARAM PARAMCD PARAMN
1 1 40 357 7 eGFR measurements eGFR 7
1 1 49 119 5 eGFR measurements eGFR 7
1 1 41 238 6 eGFR measurements eGFR 7
1 1 36 469 8 eGFR measurements eGFR 7
1 1 40 721 10 eGFR measurements eGFR 7
1 1 34 602 9 eGFR measurements eGFR 7

and ADET for the time-to-event outcomes with their timing.

ADET dataset structure
ID AVAL PARAM PARAMCD PARAMN TRTPN
11 359 Sustained >=50% decline in eGFR EGFR50 5 2
11 467 Sustained >=57% decline in eGFR EGFR57 4 2
11 359 Sustained >=40% decline in eGFR EGFR40 6 2
11 841 Sustained eGFR<15 (mL/min/1.73 m2) EGFR15 3 2
15 962 Death (adj) DTHADJ 1 1
21 737 Sustained >=40% decline in eGFR EGFR40 6 1

For the detailed derivation see the Technical Appendix in Heerspink et al. (2023).

Heart Failure

In the Heart Failure population (see Kondo et al. (2023)) the following HCE was considered

#>   Order                                        Category
#> 1     I                            Cardiovascular death
#> 2    II Total (first and recurrent) HF hospitalizations
#> 3   III                          Total urgent HF visits
#> 4    IV           Improvement/deterioration in KCCQ-TSS

References

Beigel, John H, Kay M Tomashek, Lori E Dodd, Aneesh K Mehta, Barry S Zingman, Andre C Kalil, Elizabeth Hohmann, et al. 2020. “Remdesivir for the Treatment of Covid-19.” New England Journal of Medicine 383 (19): 1813–26. https://doi.org/10.1056/NEJMoa2007764.
Gasparyan, Samvel B. 2024. hce: Design and Analysis of Hierarchical Composite Endpoints. CRAN: The Comprehensive R Archive Network, R Package, Version 0.6.5. https://CRAN.R-project.org/package=hce.
Gasparyan, Samvel B, Joan Buenconsejo, Elaine K Kowalewski, Jan Oscarsson, Olof F Bengtsson, Russell Esterline, Gary G Koch, Otavio Berwanger, and Mikhail N Kosiborod. 2022. “Design and Analysis of Studies Based on Hierarchical Composite Endpoints: Insights from the DARE-19 Trial.” Therapeutic Innovation & Regulatory Science 56 (5): 785–94. https://doi.org/10.1007/s43441-022-00420-1.
Gasparyan, Samvel B, Elaine K Kowalewski, Joan Buenconsejo, and Gary G Koch. 2023. “Hierarchical Composite Endpoints in COVID-19: The DARE-19 Trial.” In Case Studies in Innovative Clinical Trials, 95–148. Chapman; Hall/CRC. https://doi.org/10.1201/9781003288640-7.
Heerspink, Hiddo L, Niels Jongs, Patrick Schloemer, Dustin J Little, Meike Brinker, Christoph Taso, Martin Karpefors, et al. 2023. “Development and Validation of a New HCE for Clinical Trials of Kidney Disease Progression.” Journal of the American Society of Nephrology 34 (12): 2025–38. https://doi.org/10.1681/ASN.0000000000000243.
Kondo, Toru, Samvel B Gasparyan, Pardeep S Jhund, Olof Bengtsson, Brian L Claggett, Rudolf A de Boer, Adrian F Hernandez, et al. 2023. “Use of Win Statistics to Analyze Outcomes in the DAPA-HF and DELIVER Trials.” NEJM Evidence 2 (11): EVIDoa2300042. https://doi.org/10.1056/EVIDoa2300042.
Kosiborod, Mikhail N, Russell Esterline, Remo HM Furtado, Jan Oscarsson, Samvel B Gasparyan, Gary G Koch, Felipe Martinez, et al. 2021. “Dapagliflozin in Patients with Cardiometabolic Risk Factors Hospitalised with COVID-19 (DARE-19): A Randomised, Double-Blind, Placebo-Controlled, Phase 3 Trial.” The Lancet Diabetes & Endocrinology 9 (9): 586–94. https://doi.org/10.1016/S2213-8587(21)00180-7.
Kosiborod, Mikhail, Otavio Berwanger, Gary G Koch, Felipe Martinez, Omar Mukhtar, Subodh Verma, Vijay Chopra, et al. 2021. “Effects of Dapagliflozin on Prevention of Major Clinical Events and Recovery in Patients with Respiratory Failure Because of COVID-19: Design and Rationale for the DARE-19 Study.” Diabetes, Obesity and Metabolism 23 (4): 886–96. https://doi.org/10.1111/dom.14296.
Little, Dustin J, Samvel B Gasparyan, Patrick Schloemer, Niels Jongs, Meike Brinker, Martin Karpefors, Christoph Taso, et al. 2023. “Validity and Utility of a Hierarchical Composite Endpoint for Clinical Trials of Kidney Disease Progression: A Review.” Journal of the American Society of Nephrology 34 (12): 1928–35. https://doi.org/10.1681/ASN.0000000000000244.