require(survival)
set.seed(8)
n <- 300
Z <- rnorm(n)
X <- rnorm(n, mean = Z)
time <- rexp(n, rate = exp(X + Z + X * Z)) # survival time
C <- rexp(n, rate = exp(X + Z + X * Z)) # censoring time
U <- pmin(time, C) # time at risk
D <- as.numeric(time < C) # event indicator
dd <- data.frame(Z, X, U, D)
x <- standardize_coxph(
formula = Surv(U, D) ~ X + Z + X * Z,
data = dd, values = list(X = seq(-1, 1, 0.5)), times = c(2,3,4)
)
tidy(x)
Run the code above in your browser using DataLab