library(survival)

source('../../Data/get_bile_data.R')
DF = get_bile_data()

plot(survfit(Surv(Time, Censored) ~ Treatment, data=DF),
     col=c('blue', 'red'), xlim=c(0, 1500),
     lwd=2,
     mark.time=TRUE,
     cex.lab=1.5,
     xlab='Time (days)',
     ylab='Survival fraction')
legend(500, 0.8, c('Radiation + 5-FU (Treatment=0)', 'Placebo (Treatment=1)'), lwd=2, col=c('red', 'blue'))
grid()


## Log-rank survival difference:
##
res = survdiff(Surv(Time, Censored) ~ Treatment, data=DF)
print(res)

res = survdiff(Surv(Time, Censored) ~ Treatment, data=DF, rho=1)
print(res)