Simulated sample with exponentially or Weibull distributed time-to-event; log-hazard depends non-linearly on risk factors, and includes cross-terms.
simulate_crossterms(
N = 300,
observe_time = 10,
percentcensored = 0.75,
randomseed = NULL,
lambda = 0.1,
distr = "Exp",
rho_w = 1,
drop_out = 0.3
)
data frame; "time" and "event" columns describe survival outcome; predictors are "age", "sex", "hyp", "bmi"
sample size, 300 by default
study's observation time, 10 by default
expected number of non-events by observe_time, 0.75 by default (i.e. event rate is 0.25)
random seed for replication
baseline hazard rate, 0.1 by default
time-to-event distribution, "Exp" for exponential (default), "W" for Weibull
shape parameter for Weibull distribution, 0.3 by default
expected rate of drop out before observe_time, 0.3 by default
mydata <- simulate_crossterms()
head(mydata)
Run the code above in your browser using DataLab