Assignment 3

download data

library(dplyr)
dat <- readxl::read_excel("assignment3.xlsx") %>% 
  mutate(srv_time_mon = srv_time_mon+1)

Kaplan-Meier fit

library(survival)
fit <- survfit(Surv(srv_time_mon, STAT_REC) ~ Surgery, data = dat)

log-minus-log plot

log-minus-log plot for Leukemia dataset stratified by surgery

library(survminer)
ggsurvplot(
  fit,
  fun = "cloglog",
  risk.table = FALSE,
  linetype = 1:2
)

fit output

summary(fit)