Simulate event times for each individual from a piecewise constant hazard model
sim_pw_const_haz(n, hazard_periods = NULL, hazard_values)
Arguments
- n
Number of individuals
- hazard_periods
Vector of break points between time periods with separate constant
hazards, e.g., c(6,8) defines [0,6), [6,8), [8, infinity). Leave as NULL if
defining only one hazard period.
- hazard_values
Vector of constant hazard values associated with the time intervals
Value
Vector of simulated times from the time-to-event distribution
Examples
tte_dat <- sim_pw_const_haz(n = 100000, hazard_periods = c(6, 8), hazard_values = c(0.1, 0.1, 0.1))
hist(tte_dat, breaks = 100, main = "Event Time Distribution", xlab = "Event Time")
