### simulating survival data ###
set.seed(123456)
n = 2000
p = 2
x = matrix(0, n, p)
colnames(x) = c("cov1", "cov2")
x[,1] = rnorm(n, 0, 2)
x[,2] = sample(c(0, 1), n, replace = TRUE)
# - true value of parameter
beta.true = c(0.5, 0.5)
alpha.true = 1.5
kappa.true = 0.02
cens <- c(30, 40)
simData <- simSurv(x, beta.true, alpha.true, kappa.true, cens)
y <- simData$y
delta <- simData$delta
survData <- as.data.frame(cbind(y, delta, x))
colnames(survData)[1:2] <- c("time", "event")
Run the code above in your browser using DataLab