data(bmt)
bmt$cluster <- sample(1:100,408,replace=TRUE)
out1 <- cif(Event(time,cause)~+1,data=bmt,cause=1)
out2 <- cif(Event(time,cause)~+1+cluster(cluster),data=bmt,cause=1)
par(mfrow=c(1,2))
plot(out1,se=TRUE)
plot(out2,se=TRUE)
Run the code above in your browser using DataLab