library('CompAREdesign') # Load library knitr::opts_chunk$set(echo = TRUE) # Show code
This example is based on the data from the TUXEDO trial [1].
The Taxus Element versus Xience Prime in a Diabetic Population (TUXEDO)–India study is an investigator-initiated, multicenter, randomized clinical trial. The trial protocol is available with the full text of this article at NEJM.org. The trial was funded by Boston Scientific, the manufacturer of the paclitaxeleluting stent (Taxus Element).
To use the functions in CompAREdesign, we first of all need to include the information for the components of the composite endpoint.
## Probabilities of observing the event in control arm at the end of follow-up p0_e1 <- 0.059 # Ischemia-driven target-lesion revascularization p0_e2 <- 0.032 # Cardiac death or target-vessel MI ## Effect size (absolute reduction) for each endpoint AR_e1 <- -0.0196 # Ischemia-driven target-lesion revascularization AR_e2 <- -0.0098 # Cardiac death or target-vessel MI ## Correlation rho <- 0.4
Aiming to compare the gain in power of using the composite endpoint over the most relevant of its components. One can use the function ARE. As the ARE is greater than 1, we can state that the design using the composite endpoint is more efficient.
ARE_cbe(p0_e1 = p0_e1 , p0_e2 = p0_e2, eff_e1 = AR_e1 , eff_e2 = AR_e2, effm_e1 = "diff" , effm_e2 = "diff", effm_ce = "or", rho = rho)
If we can anticipate the baseline information of the composite components and expected effect sizes together with the correlation between the components, we can obtain the expected treatment effect on the composite endpoint (odds ratio, OR).
effectsize_cbe(p0_e1 = p0_e1 , p0_e2 = p0_e2, eff_e1 = AR_e1 , eff_e2 = AR_e2, effm_e1 = "diff" , effm_e2 = "diff", effm_ce = "or", rho = rho)
Lastly, we can compute the required sample size for the design using the composite endpoint, in this case obtaining 2644.
samplesize_cbe(p0_e1 = p0_e1 , p0_e2 = p0_e2, eff_e1 = AR_e1 , eff_e2 = AR_e2, effm_e1 = "diff" , effm_e2 = "diff", effm_ce = "or", rho = rho, alpha = 0.05, beta = 0.2)
rm(list = ls())
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.