set.seed(12345)
age <- rnorm(100, 50, 10)
size <- rexp(100,1)
stime <- rexp(100)
cens <- runif(100,.5,2)
sevent <- as.numeric(stime <= cens)
stime <- pmin(stime, cens)
c1 <- concordance.index(x=age, surv.time=stime, surv.event=sevent,
method="noether")
c2 <- concordance.index(x=size, surv.time=stime, surv.event=sevent,
method="noether")
cindex.comp(c1, c2)
Run the code above in your browser using DataLab