# Sample Code for integrated risk difference (IRD) source("FUN-IRD-ver001.R") library(survival) D=cbind(pbc$time, as.numeric(pbc$status==2), as.numeric(pbc$trt==1)) D=D[1:312,] FUN.IRD(data=D,points=c(1000,2000,3000),t0=0,t1=3000)