# create a trial_sequence-class object
trial_itt_dir <- file.path(tempdir(), "trial_itt")
dir.create(trial_itt_dir)
trial_itt <- trial_sequence(estimand = "ITT") |>
set_data(data = data_censored) |>
set_outcome_model(adjustment_terms = ~ x1 + x2)
trial_itt_csv <- set_expansion_options(
trial_itt,
output = save_to_csv(file.path(trial_itt_dir, "trial_csvs")),
chunk_size = 500
) |>
expand_trials()
# load_expanded_data default behaviour returns all trial_periods and doesn't sample
load_expanded_data(trial_itt_csv)
# load_expanded_data can subset the data before sampling
load_expanded_data(
trial_itt_csv,
p_control = 0.2,
period = 1:20,
subset_condition = "followup_time %in% 1:20 & x2 < 1",
)
# delete after use
unlink(trial_itt_dir, recursive = TRUE)
Run the code above in your browser using DataLab