# NOT RUN {
library(survival)
dat <- pbc[!is.na(pbc$trt),]
time <- dat$time/365.25
status <- (dat$status == 2) + 0
arm <- (dat$trt == 2) + 0
KM2NNT(time, status, arm, tau = NULL, confint = 0.95, digits = 3)
# }
Run the code above in your browser using DataLab