# NOT RUN {
library("survival")
# Get survdiff results
fit <- survfit(Surv(time, status) ~ rx, data = colon)
get_logrank_res(formula(fit), colon)
# Get only log-rank p-value
get_logrank_res(formula(fit), colon, return.p = TRUE)
# }
Run the code above in your browser using DataLab