# Tests for the leukemia data from Exercise 3.1

 

# Read data (directly from web-page) and attach survival library

leukemia=read.table("http://www.math.uio.no/~borgan/abg-2008/data/leukemia.txt", header=T)

library(survival)

 

# Log-rank test:

survdiff(Surv(time,status)~treat, data=leukemia)

 

# Harrington-Fleming test with rho=1:

survdiff(Surv(time,status)~treat, data=leukemia, rho=1)