weibull_rmst: RMST using Weibull fit

In reconstructKM: Reconstruct Individual-Level Data from Published KM Plots

Description

RMST for time-to-event data under parametric Weibull fit for data in each arm separately. Also can provide CI for RMST estimate and difference in RMST.

Usage

weibull_rmst(num_boots = 1000, dat, tau, alpha, find_pval = FALSE, seed = NULL)

Arguments

num_boots	Number of bootstrap iterations
dat	Data frame of time-to-event data which MUST have the columns 'time', 'arm', and 'status
tau	How long of a follow-up to consider, i.e. we integrate the survival functions from 0 to tau
alpha	Confidence interval is given for (alpha/2, 1-alpha/2) percentiles
find_pval	Boolean, if TRUE then does bootstrap under the null to find p-value of mean difference and RMST difference
seed	For reproducibility

Value

A list including out_tab (estimate and CI in both arms), trt_rmst, pbo_rmst, diff_rmst, trt_CI, pbo_CI, diff_CI. Assumes trt coded as arm 1 and placebo coded as arm 0.

Examples

time <- rexp(100)
status <- rbinom(n=100, prob=0.5, size=1)
arm <- c( rep(1, 50), rep(0, 50))
dat <- data.frame(time=time, status=status, arm=arm)
weibull_rmst(dat=dat, tau=1, alpha=0.05, num_boots=200)

reconstructKM documentation built on Nov. 25, 2020, 5:08 p.m.