library(dplyr)
# Use default enrollment and calendar cut date
# for 50 events in the "Positive" stratum
x <- sim_pw_surv(
n = 200,
stratum = data.frame(
stratum = c("Positive", "Negative"),
p = c(.5, .5)
),
fail_rate = data.frame(
stratum = rep(c("Positive", "Negative"), 2),
period = rep(1, 4),
treatment = c(rep("control", 2), rep("experimental", 2)),
duration = rep(1, 4),
rate = log(2) / c(6, 9, 9, 12)
),
dropout_rate = data.frame(
stratum = rep(c("Positive", "Negative"), 2),
period = rep(1, 4),
treatment = c(rep("control", 2), rep("experimental", 2)),
duration = rep(1, 4),
rate = rep(.001, 4)
)
)
d <- get_cut_date_by_event(x |> filter(stratum == "Positive"), event = 50)
y <- cut_data_by_date(x, cut_date = d)
table(y$stratum, y$event)
Run the code above in your browser using DataLab