# NOT RUN {
## Generate simulated data from transformation model
datgen <- function(n) {
a <- -0.3
X <- rweibull(n, 2, 4) ## failure times
U <- rweibull(n, 2, 1) ## latent truncation time
T <- (1 + a) * U - a * X ## apply transformation
C <- rlnorm(n, .8, 1) ## censoring
dat <- data.frame(trun = T, obs = pmin(X, C), delta = 1 * (X <= C))
return(subset(dat, trun <= obs))
}
set.seed(123)
dat <- datgen(300)
with(dat, pmcc(trun, obs))
# }
Run the code above in your browser using DataLab