# NOT RUN {
T <- rexp(30) ## event times
Z <- rbinom(30, 1, 0.5) ## treatment assignment
C <- rexp(30) ## drop-out times
data <- data.frame(time=pmin(T,C), status=T<=C, Z=Z)
x <- permIPT(Surv(time, status) ~ Z, data)
summary(x)
# }
Run the code above in your browser using DataLab