tests/testthat/test-derive_param_tte.R
In Roche-GSK/admiral: ADaM in R Asset Library
# derive_param_tte ----
## Test 1: new observations with analysis date are derived correctly ----
test_that("derive_param_tte Test 1: new observations with analysis date are derived correctly", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~LSTALVDT, ~TRTSDT, ~TRTSDTF,
"03", "Y", ymd("2021-08-21"), ymd("2021-08-21"), ymd("2021-08-10"), NA,
"04", "N", NA, ymd("2021-05-24"), ymd("2021-02-03"), NA
) %>%
mutate(STUDYID = "AB42")
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lstalv <- censor_source(
dataset_name = "adsl",
date = LSTALVDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST KNOWN ALIVE DATE",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR,
"03", ymd("2021-08-21"), 0L, "DEATH", "ADSL", "DTHDT",
"04", ymd("2021-05-24"), 1L, "LAST KNOWN ALIVE DATE", "ADSL", "LSTALVDT"
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "OS",
PARAM = "Overall Survival"
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT, STARTDTF = TRTSDTF), by = "USUBJID")
actual_output <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(death),
censor_conditions = list(lstalv),
source_datasets = list(adsl = adsl),
set_values_to = exprs(
PARAMCD = "OS",
PARAM = "Overall Survival"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 2: new parameter with analysis datetime is derived correctly ----
test_that("derive_param_tte Test 2: new parameter with analysis datetime is derived correctly", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~TRTSDTM, ~TRTSDTF, ~TRTSTMF,
"01", "Y", ymd("2021-06-12"), ymd_hms("2021-01-01 00:00:00"), "M", "H",
"02", "N", NA, ymd_hms("2021-02-03 10:24:00"), NA, NA,
"03", "Y", ymd("2021-08-21"), ymd_hms("2021-08-10 00:00:00"), NA, "H",
"04", "N", NA, ymd_hms("2021-02-03 10:24:00"), NA, NA,
"05", "N", NA, ymd_hms("2021-04-05 11:22:33"), NA, NA
) %>%
mutate(STUDYID = "AB42")
adrs <- tibble::tribble(
~USUBJID, ~AVALC, ~ADTM, ~ASEQ,
"01", "SD", ymd_hms("2021-01-03 10:56:00"), 1,
"01", "PR", ymd_hms("2021-03-04 11:13:00"), 2,
"01", "PD", ymd_hms("2021-05-05 12:02:00"), 3,
"02", "PD", ymd_hms("2021-02-03 10:56:00"), 1,
"04", "SD", ymd_hms("2021-02-13 10:56:00"), 1,
"04", "PR", ymd_hms("2021-04-14 11:13:00"), 2,
"04", "CR", ymd_hms("2021-05-15 12:02:00"), 3
) %>%
mutate(STUDYID = "AB42", PARAMCD = "OVR")
pd <- event_source(
dataset_name = "adrs",
filter = AVALC == "PD",
date = ADTM,
set_values_to = exprs(
EVENTDESC = "PD",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lastvisit <- censor_source(
dataset_name = "adrs",
date = ADTM,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST TUMOR ASSESSMENT",
SRCDOM = "ADRS",
SRCVAR = "ADTM"
)
)
start <- censor_source(
dataset_name = "adsl",
date = TRTSDTM,
censor = 1,
set_values_to = exprs(
EVENTDESC = "TREATMENT START",
SRCDOM = "ADSL",
SRCVAR = "TRTSDTM"
)
)
# nolint start
expected_output <- tibble::tribble(
~USUBJID, ~ADTM, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd_hms("2021-05-05 12:02:00"), 0L, "PD", "ADRS", "ADTM", 3,
"02", ymd_hms("2021-02-03 10:56:00"), 0L, "PD", "ADRS", "ADTM", 1,
"03", as_datetime(ymd("2021-08-21")), 0L, "DEATH", "ADSL", "DTHDT", NA,
"04", ymd_hms("2021-05-15 12:02:00"), 1L, "LAST TUMOR ASSESSMENT", "ADRS", "ADTM", NA,
"05", ymd_hms("2021-04-05 11:22:33"), 1L, "TREATMENT START", "ADSL", "TRTSDTM", NA
) %>%
# nolint end
mutate(
STUDYID = "AB42",
PARAMCD = "PFS",
PARAM = "Progression Free Survival"
) %>%
left_join(
select(adsl, USUBJID, STARTDTM = TRTSDTM, STARTDTF = TRTSDTF, STARTTMF = TRTSTMF),
by = "USUBJID"
)
actual_output <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDTM,
event_conditions = list(pd, death),
censor_conditions = list(lastvisit, start),
source_datasets = list(adsl = adsl, adrs = adrs),
create_datetime = TRUE,
set_values_to = exprs(
PARAMCD = "PFS",
PARAM = "Progression Free Survival"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 3: no records for subjects not in ADSL, multiple events/cens ----
test_that("derive_param_tte Test 3: no records for subjects not in ADSL, multiple events/cens", {
# For subject 01 both events are at the same date, for subject 04 both censorings
# are at the same date
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~RSPDT,
"01", "Y", ymd("2021-05-05"), ymd("2021-03-04"),
"02", "N", NA, NA,
"03", "Y", ymd("2021-08-21"), NA,
"04", "N", NA, ymd("2021-05-15"),
"05", "N", NA, NA
) %>%
mutate(STUDYID = "AB42")
adrs <- tibble::tribble(
~USUBJID, ~AVALC, ~ADT, ~ASEQ,
"01", "SD", ymd("2021-01-03"), 1,
"01", "PR", ymd("2021-03-04"), 2,
"01", "PD", ymd("2021-05-05"), 3,
"02", "PD", ymd("2021-02-03"), 1,
"04", "SD", ymd("2021-02-13"), 1,
"04", "SD", ymd("2021-04-14"), 2,
"04", "CR", ymd("2021-05-15"), 3
) %>%
mutate(STUDYID = "AB42", PARAMCD = "OVR")
pd <- event_source(
dataset_name = "adrs",
filter = AVALC == "PD",
date = ADT,
set_values_to = exprs(
EVENTDESC = "PD",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lastvisit <- censor_source(
dataset_name = "adrs",
date = ADT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST TUMOR ASSESSMENT",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
first_response <- censor_source(
dataset_name = "adsl",
date = RSPDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "FIRST RESPONSE",
SRCDOM = "ADSL",
SRCVAR = "RSPDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd("2021-05-05"), 0L, "PD", "ADRS", "ADTM", 3,
"04", ymd("2021-05-15"), 1L, "LAST TUMOR ASSESSMENT", "ADRS", "ADTM", 3,
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "DURRSP",
PARAM = "Duration of Response"
) %>%
left_join(
select(adsl, USUBJID, STARTDT = RSPDT),
by = "USUBJID"
)
actual_output <- derive_param_tte(
dataset_adsl = filter(adsl, !is.na(RSPDT)),
start_date = RSPDT,
event_conditions = list(pd, death),
censor_conditions = list(first_response, lastvisit),
source_datasets = list(adsl = adsl, adrs = adrs),
set_values_to = exprs(
PARAMCD = "DURRSP",
PARAM = "Duration of Response"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 4: error is issued if DTC variables specified for date ----
test_that("derive_param_tte Test 4: error is issued if DTC variables specified for date", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ,
"01", "2021-01-03T10:56", 1,
"01", "2021-03-04", 2,
"01", "2021", 3
) %>%
mutate(STUDYID = "AB42")
ttae <- event_source(
dataset_name = "ae",
date = AESTDTC,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_error(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
)
),
class = "assert_date_var"
)
})
## Test 5: by_vars parameter works correctly ----
test_that("derive_param_tte Test 5: by_vars parameter works correctly", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# nolint start
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~PARCAT2, ~PARAMCD,
"01", ymd("2021-01-01"), 0L, "AE", "AE", "AESTDTC", 3, "Flu", "TTAE2",
"02", ymd("2021-01-30"), 1L, "END OF TRT", "ADSL", "TRTEDT", NA, "Flu", "TTAE2",
"01", ymd("2021-03-04"), 0L, "AE", "AE", "AESTDTC", 2, "Cough", "TTAE1",
"02", ymd("2021-01-30"), 1L, "END OF TRT", "ADSL", "TRTEDT", NA, "Cough", "TTAE1"
) %>%
# nolint end
mutate(
STUDYID = "AB42",
PARCAT1 = "TTAE",
PARAM = paste("Time to First", PARCAT2, "Adverse Event")
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT), by = "USUBJID")
expect_dfs_equal(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 6: an error is issued if some of the by variables are missing ----
test_that("derive_param_tte Test 6: an error is issued if some of the by variables are missing", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEBODSYS, AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 7: errors if all by vars are missing in all source datasets ----
test_that("derive_param_tte Test 7: errors if all by vars are missing in all source datasets", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEBODSYS),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 8: errors if PARAMCD and by_vars are not one to one ----
test_that("derive_param_tte Test 8: errors if PARAMCD and by_vars are not one to one", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 9: errors if set_values_to contains invalid expressions ----
test_that("derive_param_tte Test 9: errors if set_values_to contains invalid expressions", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = past("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 10: error is issued if parameter code already exists ----
test_that("derive_param_tte Test 10: error is issued if parameter code already exists", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd("2021-01-01"), 0L, "AE", "AE", "AESTDTC", 3,
"02", ymd("2021-02-03"), 1L, "END OF STUDY", "ADSL", "EOSDT", NA
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT), by = "USUBJID")
expect_error(
derive_param_tte(
expected_output,
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
)
),
class = "assert_param_does_not_exist"
)
})
## Test 11: ensuring ADT is not NA because of missing start_date ----
test_that("derive_param_tte Test 11: ensuring ADT is not NA because of missing start_date", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~LSTALVDT,
"01", NA, ymd("2022-08-10"),
"02", NA, ymd("2022-09-12"),
"03", ymd("2020-10-13"), ymd("2022-07-21")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESEQ, ~ASTDT,
"01", 1, ymd("2020-08-10"),
"02", 2, ymd("2020-08-15"),
"03", 3, ymd("2020-12-10"),
) %>%
mutate(STUDYID = "AB42")
eos <- censor_source(
"adsl",
date = LSTALVDT,
set_values_to = exprs(
EVNTDESC = "Last Known Alive Date",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
ttae <- event_source(
dataset_name = "adae",
date = ASTDT,
set_values_to = exprs(
EVNTDESC = "Any Adverse Event",
SRCDOM = "ADAE",
SRCVAR = "AEDECOD",
SRCSEQ = AESEQ
)
)
actual_output <- derive_param_tte(
dataset_adsl = adsl,
source_datasets = list(adae = ae, adsl = adsl),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
set_values_to = exprs(
PARAMCD = "ANYAETTE",
PARAM = "Time to any first adverse event"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT, ~STARTDT,
"01", "Any Adverse Event", "ADAE", "AEDECOD", 1, 0L, ymd("2020-08-10"), NA,
"02", "Any Adverse Event", "ADAE", "AEDECOD", 2, 0L, ymd("2020-08-15"), NA,
"03", "Any Adverse Event", "ADAE", "AEDECOD", 3, 0L, ymd("2020-12-10"), ymd("2020-10-13") # nolint
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "ANYAETTE",
PARAM = "Time to any first adverse event"
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 12: test dataset and dynamic byvars populated ----
test_that("derive_param_tte Test 12: test dataset and dynamic byvars populated", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVNTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVNTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
actual_output <- derive_param_tte(
dataset = adsl %>% select(STUDYID, USUBJID) %>% mutate(PARAMCD = "XYZ"),
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
)
expected_output <- bind_rows(
adsl %>% select(STUDYID, USUBJID) %>% mutate(PARAMCD = "XYZ"),
tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT,
"01", "AE", "AE", "AESTDTC", 2, 0L, ymd("2021-03-04"),
"02", "END OF TRT", "ADSL", "TRTEDT", NA_real_, 1L, ymd("2021-01-30"),
) %>%
mutate(
STUDYID = "AB42",
STARTDT = if_else(USUBJID == "01", ymd("2020-12-06"), ymd("2021-01-16")),
PARAMCD = "TTAE1",
PARAM = "Time to First Cough Adverse Event",
PARCAT1 = "TTAE",
PARCAT2 = "Cough"
),
tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT,
"01", "AE", "AE", "AESTDTC", 3, 0L, ymd("2021-01-01"),
"02", "END OF TRT", "ADSL", "TRTEDT", NA_real_, 1L, ymd("2021-01-30"),
) %>%
mutate(
STUDYID = "AB42",
STARTDT = if_else(USUBJID == "01", ymd("2020-12-06"), ymd("2021-01-16")),
PARAMCD = "TTAE2",
PARAM = "Time to First Flu Adverse Event",
PARCAT1 = "TTAE",
PARCAT2 = "Flu"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 13: error if dataset_name not in source_datsets ----
test_that("derive_param_tte Test 13: error if dataset_name not in source_datsets", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~LSTALVDT, ~TRTSDT, ~TRTSDTF,
"03", "Y", ymd("2021-08-21"), ymd("2021-08-21"), ymd("2021-08-10"), NA,
"04", "N", NA, ymd("2021-05-24"), ymd("2021-02-03"), NA
) %>%
mutate(STUDYID = "AB42")
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lstalv <- censor_source(
dataset_name = "adls",
date = LSTALVDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST KNOWN ALIVE DATE",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(death),
censor_conditions = list(lstalv),
source_datasets = list(adsl = adsl),
set_values_to = exprs(
PARAMCD = "OS",
PARAM = "Overall Survival"
)
),
error = TRUE
)
})
## Test 14: detects duplicates in input datasets via pipeline functions ----
test_that("derive_param_tte Test 14: detects duplicates in input datasets via pipeline functions", {
# Define ADSL dataset
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>% mutate(STUDYID = "AB42")
# Define AE dataset with duplicates
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-03", 3, "Flu"
) %>% mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
# Define event and censor sources
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# Run derive_param_tte and check for warning
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
)
})
## Test 15: using order for resolving ties ----
test_that("derive_param_tte Test 15: using order for resolving ties", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
# Sort the input AE dataset in descending order by AESEQ
# to confirm that the order argument re-sorts it correctly.
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AESER, ~AEDECOD,
"01", "2021-01-03", 2, "Y", "Cough",
"01", "2021-01-03", 1, "Y", "Flu",
"01", "2021-01-20", 3, "N", "Headache"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
) %>%
arrange(desc(AESEQ)) # Intentionally sort descending to test the order argument
result <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "Serious AE",
SRCSEQ = AESEQ
),
filter = AESER == "Y",
order = exprs(AESEQ) # Should re-sort so that AESEQ=1 (Flu) is chosen on tie
)),
censor_conditions = list(censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(EVENTDESC = "End of Study")
)),
set_values_to = exprs(
PARAMCD = "TTSAE",
PARAM = "Time to First Serious AE"
),
source_datasets = list(adsl = adsl, ae = ae)
)
# Check that for USUBJID = "01", the first serious AE selected is the one with AESEQ = 1 (Flu),
# despite the input AE data initially being arranged to show AESEQ=2 (Cough) first.
selected_seq <- result %>%
filter(USUBJID == "01", PARAMCD == "TTSAE") %>%
pull(SRCSEQ)
expect_equal(selected_seq, 1, info = "The order argument should ensure AE with AESEQ=1
is chosen on tie.")
})
## Test 16: produces consistent results regardless of input sort order ----
test_that("derive_param_tte Test 16: produces consistent results regardless of input sort order", {
# Define ADSL dataset
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>% mutate(STUDYID = "AB42")
# Define AE dataset with duplicates
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-03", 3, "Flu"
) %>% mutate(STUDYID = "AB42", AESTDT = ymd(AESTDTC))
# Define event and censor sources
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
order = exprs(AESTDT, AESEQ),
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ # Ensure AESEQ is included here
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# Run derive_param_tte with sorted AE dataset
result_sorted <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = arrange(ae, AESEQ)),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
# Run derive_param_tte with reverse-sorted AE dataset
result_unsorted <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = arrange(ae, desc(AESEQ))),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
expect_equal(result_sorted, result_unsorted)
})
# list_tte_source_objects ----
## Test 17: error is issued if package does not exist ----
test_that("list_tte_source_objects Test 17: error is issued if package does not exist", {
expect_snapshot(
list_tte_source_objects(package = "tte"),
error = TRUE
)
})
## Test 18: expected objects produced ----
test_that("list_tte_source_objects Test 18: expected objects produced", {
expected_output <- tibble::tribble(
~object, ~dataset_name, ~filter,
"ae_ser_event", "adae", quote(TRTEMFL == "Y" & AESER == "Y"),
"ae_gr2_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "2"),
"ae_sev_event", "adae", quote(TRTEMFL == "Y" & AESEV == "SEVERE"),
"ae_gr4_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "4"),
"ae_gr3_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "3"),
"lastalive_censor", "adsl", NULL,
"ae_event", "adae", quote(TRTEMFL == "Y"),
"death_event", "adsl", quote(DTHFL == "Y"),
"ae_gr35_event", "adae", quote(TRTEMFL == "Y" & ATOXGR %in% c("3", "4", "5")),
"ae_wd_event", "adae", quote(TRTEMFL == "Y" & AEACN == "DRUG WITHDRAWN"),
"ae_gr1_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "1"),
"ae_gr5_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "5")
) %>%
mutate(
date = case_when(
object == "lastalive_censor" ~ "LSTALVDT",
object == "death_event" ~ "DTHDT",
TRUE ~ "ASTDT"
),
censor = if_else(object == "lastalive_censor", 1, 0),
filter = as.character(filter),
censor = as.integer(censor)
)
observed_output <- list_tte_source_objects(package = "admiral") %>%
select(object, dataset_name, filter, date, censor)
expect_dfs_equal(expected_output, observed_output, keys = c("object"))
})
Roche-GSK/admiral documentation built on April 14, 2025, 12:36 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
# derive_param_tte ----
## Test 1: new observations with analysis date are derived correctly ----
test_that("derive_param_tte Test 1: new observations with analysis date are derived correctly", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~LSTALVDT, ~TRTSDT, ~TRTSDTF,
"03", "Y", ymd("2021-08-21"), ymd("2021-08-21"), ymd("2021-08-10"), NA,
"04", "N", NA, ymd("2021-05-24"), ymd("2021-02-03"), NA
) %>%
mutate(STUDYID = "AB42")
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lstalv <- censor_source(
dataset_name = "adsl",
date = LSTALVDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST KNOWN ALIVE DATE",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR,
"03", ymd("2021-08-21"), 0L, "DEATH", "ADSL", "DTHDT",
"04", ymd("2021-05-24"), 1L, "LAST KNOWN ALIVE DATE", "ADSL", "LSTALVDT"
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "OS",
PARAM = "Overall Survival"
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT, STARTDTF = TRTSDTF), by = "USUBJID")
actual_output <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(death),
censor_conditions = list(lstalv),
source_datasets = list(adsl = adsl),
set_values_to = exprs(
PARAMCD = "OS",
PARAM = "Overall Survival"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 2: new parameter with analysis datetime is derived correctly ----
test_that("derive_param_tte Test 2: new parameter with analysis datetime is derived correctly", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~TRTSDTM, ~TRTSDTF, ~TRTSTMF,
"01", "Y", ymd("2021-06-12"), ymd_hms("2021-01-01 00:00:00"), "M", "H",
"02", "N", NA, ymd_hms("2021-02-03 10:24:00"), NA, NA,
"03", "Y", ymd("2021-08-21"), ymd_hms("2021-08-10 00:00:00"), NA, "H",
"04", "N", NA, ymd_hms("2021-02-03 10:24:00"), NA, NA,
"05", "N", NA, ymd_hms("2021-04-05 11:22:33"), NA, NA
) %>%
mutate(STUDYID = "AB42")
adrs <- tibble::tribble(
~USUBJID, ~AVALC, ~ADTM, ~ASEQ,
"01", "SD", ymd_hms("2021-01-03 10:56:00"), 1,
"01", "PR", ymd_hms("2021-03-04 11:13:00"), 2,
"01", "PD", ymd_hms("2021-05-05 12:02:00"), 3,
"02", "PD", ymd_hms("2021-02-03 10:56:00"), 1,
"04", "SD", ymd_hms("2021-02-13 10:56:00"), 1,
"04", "PR", ymd_hms("2021-04-14 11:13:00"), 2,
"04", "CR", ymd_hms("2021-05-15 12:02:00"), 3
) %>%
mutate(STUDYID = "AB42", PARAMCD = "OVR")
pd <- event_source(
dataset_name = "adrs",
filter = AVALC == "PD",
date = ADTM,
set_values_to = exprs(
EVENTDESC = "PD",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lastvisit <- censor_source(
dataset_name = "adrs",
date = ADTM,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST TUMOR ASSESSMENT",
SRCDOM = "ADRS",
SRCVAR = "ADTM"
)
)
start <- censor_source(
dataset_name = "adsl",
date = TRTSDTM,
censor = 1,
set_values_to = exprs(
EVENTDESC = "TREATMENT START",
SRCDOM = "ADSL",
SRCVAR = "TRTSDTM"
)
)
# nolint start
expected_output <- tibble::tribble(
~USUBJID, ~ADTM, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd_hms("2021-05-05 12:02:00"), 0L, "PD", "ADRS", "ADTM", 3,
"02", ymd_hms("2021-02-03 10:56:00"), 0L, "PD", "ADRS", "ADTM", 1,
"03", as_datetime(ymd("2021-08-21")), 0L, "DEATH", "ADSL", "DTHDT", NA,
"04", ymd_hms("2021-05-15 12:02:00"), 1L, "LAST TUMOR ASSESSMENT", "ADRS", "ADTM", NA,
"05", ymd_hms("2021-04-05 11:22:33"), 1L, "TREATMENT START", "ADSL", "TRTSDTM", NA
) %>%
# nolint end
mutate(
STUDYID = "AB42",
PARAMCD = "PFS",
PARAM = "Progression Free Survival"
) %>%
left_join(
select(adsl, USUBJID, STARTDTM = TRTSDTM, STARTDTF = TRTSDTF, STARTTMF = TRTSTMF),
by = "USUBJID"
)
actual_output <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDTM,
event_conditions = list(pd, death),
censor_conditions = list(lastvisit, start),
source_datasets = list(adsl = adsl, adrs = adrs),
create_datetime = TRUE,
set_values_to = exprs(
PARAMCD = "PFS",
PARAM = "Progression Free Survival"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 3: no records for subjects not in ADSL, multiple events/cens ----
test_that("derive_param_tte Test 3: no records for subjects not in ADSL, multiple events/cens", {
# For subject 01 both events are at the same date, for subject 04 both censorings
# are at the same date
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~RSPDT,
"01", "Y", ymd("2021-05-05"), ymd("2021-03-04"),
"02", "N", NA, NA,
"03", "Y", ymd("2021-08-21"), NA,
"04", "N", NA, ymd("2021-05-15"),
"05", "N", NA, NA
) %>%
mutate(STUDYID = "AB42")
adrs <- tibble::tribble(
~USUBJID, ~AVALC, ~ADT, ~ASEQ,
"01", "SD", ymd("2021-01-03"), 1,
"01", "PR", ymd("2021-03-04"), 2,
"01", "PD", ymd("2021-05-05"), 3,
"02", "PD", ymd("2021-02-03"), 1,
"04", "SD", ymd("2021-02-13"), 1,
"04", "SD", ymd("2021-04-14"), 2,
"04", "CR", ymd("2021-05-15"), 3
) %>%
mutate(STUDYID = "AB42", PARAMCD = "OVR")
pd <- event_source(
dataset_name = "adrs",
filter = AVALC == "PD",
date = ADT,
set_values_to = exprs(
EVENTDESC = "PD",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lastvisit <- censor_source(
dataset_name = "adrs",
date = ADT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST TUMOR ASSESSMENT",
SRCDOM = "ADRS",
SRCVAR = "ADTM",
SRCSEQ = ASEQ
)
)
first_response <- censor_source(
dataset_name = "adsl",
date = RSPDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "FIRST RESPONSE",
SRCDOM = "ADSL",
SRCVAR = "RSPDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd("2021-05-05"), 0L, "PD", "ADRS", "ADTM", 3,
"04", ymd("2021-05-15"), 1L, "LAST TUMOR ASSESSMENT", "ADRS", "ADTM", 3,
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "DURRSP",
PARAM = "Duration of Response"
) %>%
left_join(
select(adsl, USUBJID, STARTDT = RSPDT),
by = "USUBJID"
)
actual_output <- derive_param_tte(
dataset_adsl = filter(adsl, !is.na(RSPDT)),
start_date = RSPDT,
event_conditions = list(pd, death),
censor_conditions = list(first_response, lastvisit),
source_datasets = list(adsl = adsl, adrs = adrs),
set_values_to = exprs(
PARAMCD = "DURRSP",
PARAM = "Duration of Response"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 4: error is issued if DTC variables specified for date ----
test_that("derive_param_tte Test 4: error is issued if DTC variables specified for date", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ,
"01", "2021-01-03T10:56", 1,
"01", "2021-03-04", 2,
"01", "2021", 3
) %>%
mutate(STUDYID = "AB42")
ttae <- event_source(
dataset_name = "ae",
date = AESTDTC,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_error(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
)
),
class = "assert_date_var"
)
})
## Test 5: by_vars parameter works correctly ----
test_that("derive_param_tte Test 5: by_vars parameter works correctly", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# nolint start
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~PARCAT2, ~PARAMCD,
"01", ymd("2021-01-01"), 0L, "AE", "AE", "AESTDTC", 3, "Flu", "TTAE2",
"02", ymd("2021-01-30"), 1L, "END OF TRT", "ADSL", "TRTEDT", NA, "Flu", "TTAE2",
"01", ymd("2021-03-04"), 0L, "AE", "AE", "AESTDTC", 2, "Cough", "TTAE1",
"02", ymd("2021-01-30"), 1L, "END OF TRT", "ADSL", "TRTEDT", NA, "Cough", "TTAE1"
) %>%
# nolint end
mutate(
STUDYID = "AB42",
PARCAT1 = "TTAE",
PARAM = paste("Time to First", PARCAT2, "Adverse Event")
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT), by = "USUBJID")
expect_dfs_equal(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 6: an error is issued if some of the by variables are missing ----
test_that("derive_param_tte Test 6: an error is issued if some of the by variables are missing", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEBODSYS, AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 7: errors if all by vars are missing in all source datasets ----
test_that("derive_param_tte Test 7: errors if all by vars are missing in all source datasets", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEBODSYS),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 8: errors if PARAMCD and by_vars are not one to one ----
test_that("derive_param_tte Test 8: errors if PARAMCD and by_vars are not one to one", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 9: errors if set_values_to contains invalid expressions ----
test_that("derive_param_tte Test 9: errors if set_values_to contains invalid expressions", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = past("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
),
error = TRUE
)
})
## Test 10: error is issued if parameter code already exists ----
test_that("derive_param_tte Test 10: error is issued if parameter code already exists", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eos <- censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF STUDY",
SRCDOM = "ADSL",
SRCVAR = "EOSDT"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~ADT, ~CNSR, ~EVENTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ,
"01", ymd("2021-01-01"), 0L, "AE", "AE", "AESTDTC", 3,
"02", ymd("2021-02-03"), 1L, "END OF STUDY", "ADSL", "EOSDT", NA
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
) %>%
left_join(select(adsl, USUBJID, STARTDT = TRTSDT), by = "USUBJID")
expect_error(
derive_param_tte(
expected_output,
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = "TTAE",
PARAM = "Time to First Adverse Event"
)
),
class = "assert_param_does_not_exist"
)
})
## Test 11: ensuring ADT is not NA because of missing start_date ----
test_that("derive_param_tte Test 11: ensuring ADT is not NA because of missing start_date", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~LSTALVDT,
"01", NA, ymd("2022-08-10"),
"02", NA, ymd("2022-09-12"),
"03", ymd("2020-10-13"), ymd("2022-07-21")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESEQ, ~ASTDT,
"01", 1, ymd("2020-08-10"),
"02", 2, ymd("2020-08-15"),
"03", 3, ymd("2020-12-10"),
) %>%
mutate(STUDYID = "AB42")
eos <- censor_source(
"adsl",
date = LSTALVDT,
set_values_to = exprs(
EVNTDESC = "Last Known Alive Date",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
ttae <- event_source(
dataset_name = "adae",
date = ASTDT,
set_values_to = exprs(
EVNTDESC = "Any Adverse Event",
SRCDOM = "ADAE",
SRCVAR = "AEDECOD",
SRCSEQ = AESEQ
)
)
actual_output <- derive_param_tte(
dataset_adsl = adsl,
source_datasets = list(adae = ae, adsl = adsl),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eos),
set_values_to = exprs(
PARAMCD = "ANYAETTE",
PARAM = "Time to any first adverse event"
)
)
expected_output <- tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT, ~STARTDT,
"01", "Any Adverse Event", "ADAE", "AEDECOD", 1, 0L, ymd("2020-08-10"), NA,
"02", "Any Adverse Event", "ADAE", "AEDECOD", 2, 0L, ymd("2020-08-15"), NA,
"03", "Any Adverse Event", "ADAE", "AEDECOD", 3, 0L, ymd("2020-12-10"), ymd("2020-10-13") # nolint
) %>%
mutate(
STUDYID = "AB42",
PARAMCD = "ANYAETTE",
PARAM = "Time to any first adverse event"
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 12: test dataset and dynamic byvars populated ----
test_that("derive_param_tte Test 12: test dataset and dynamic byvars populated", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-01", 3, "Flu"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVNTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVNTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
actual_output <- derive_param_tte(
dataset = adsl %>% select(STUDYID, USUBJID) %>% mutate(PARAMCD = "XYZ"),
dataset_adsl = adsl,
by_vars = exprs(AEDECOD),
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(
PARAMCD = paste0("TTAE", as.numeric(as.factor(AEDECOD))),
PARAM = paste("Time to First", AEDECOD, "Adverse Event"),
PARCAT1 = "TTAE",
PARCAT2 = AEDECOD
)
)
expected_output <- bind_rows(
adsl %>% select(STUDYID, USUBJID) %>% mutate(PARAMCD = "XYZ"),
tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT,
"01", "AE", "AE", "AESTDTC", 2, 0L, ymd("2021-03-04"),
"02", "END OF TRT", "ADSL", "TRTEDT", NA_real_, 1L, ymd("2021-01-30"),
) %>%
mutate(
STUDYID = "AB42",
STARTDT = if_else(USUBJID == "01", ymd("2020-12-06"), ymd("2021-01-16")),
PARAMCD = "TTAE1",
PARAM = "Time to First Cough Adverse Event",
PARCAT1 = "TTAE",
PARCAT2 = "Cough"
),
tibble::tribble(
~USUBJID, ~EVNTDESC, ~SRCDOM, ~SRCVAR, ~SRCSEQ, ~CNSR, ~ADT,
"01", "AE", "AE", "AESTDTC", 3, 0L, ymd("2021-01-01"),
"02", "END OF TRT", "ADSL", "TRTEDT", NA_real_, 1L, ymd("2021-01-30"),
) %>%
mutate(
STUDYID = "AB42",
STARTDT = if_else(USUBJID == "01", ymd("2020-12-06"), ymd("2021-01-16")),
PARAMCD = "TTAE2",
PARAM = "Time to First Flu Adverse Event",
PARCAT1 = "TTAE",
PARCAT2 = "Flu"
)
)
expect_dfs_equal(
actual_output,
expected_output,
keys = c("USUBJID", "PARAMCD")
)
})
## Test 13: error if dataset_name not in source_datsets ----
test_that("derive_param_tte Test 13: error if dataset_name not in source_datsets", {
adsl <- tibble::tribble(
~USUBJID, ~DTHFL, ~DTHDT, ~LSTALVDT, ~TRTSDT, ~TRTSDTF,
"03", "Y", ymd("2021-08-21"), ymd("2021-08-21"), ymd("2021-08-10"), NA,
"04", "N", NA, ymd("2021-05-24"), ymd("2021-02-03"), NA
) %>%
mutate(STUDYID = "AB42")
death <- event_source(
dataset_name = "adsl",
filter = DTHFL == "Y",
date = DTHDT,
set_values_to = exprs(
EVENTDESC = "DEATH",
SRCDOM = "ADSL",
SRCVAR = "DTHDT"
)
)
lstalv <- censor_source(
dataset_name = "adls",
date = LSTALVDT,
censor = 1,
set_values_to = exprs(
EVENTDESC = "LAST KNOWN ALIVE DATE",
SRCDOM = "ADSL",
SRCVAR = "LSTALVDT"
)
)
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(death),
censor_conditions = list(lstalv),
source_datasets = list(adsl = adsl),
set_values_to = exprs(
PARAMCD = "OS",
PARAM = "Overall Survival"
)
),
error = TRUE
)
})
## Test 14: detects duplicates in input datasets via pipeline functions ----
test_that("derive_param_tte Test 14: detects duplicates in input datasets via pipeline functions", {
# Define ADSL dataset
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>% mutate(STUDYID = "AB42")
# Define AE dataset with duplicates
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-03", 3, "Flu"
) %>% mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
)
# Define event and censor sources
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# Run derive_param_tte and check for warning
expect_snapshot(
derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = ae),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
)
})
## Test 15: using order for resolving ties ----
test_that("derive_param_tte Test 15: using order for resolving ties", {
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-02-03")
) %>%
mutate(STUDYID = "AB42")
# Sort the input AE dataset in descending order by AESEQ
# to confirm that the order argument re-sorts it correctly.
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AESER, ~AEDECOD,
"01", "2021-01-03", 2, "Y", "Cough",
"01", "2021-01-03", 1, "Y", "Flu",
"01", "2021-01-20", 3, "N", "Headache"
) %>%
mutate(
STUDYID = "AB42",
AESTDT = ymd(AESTDTC)
) %>%
arrange(desc(AESEQ)) # Intentionally sort descending to test the order argument
result <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(event_source(
dataset_name = "ae",
date = AESTDT,
set_values_to = exprs(
EVENTDESC = "Serious AE",
SRCSEQ = AESEQ
),
filter = AESER == "Y",
order = exprs(AESEQ) # Should re-sort so that AESEQ=1 (Flu) is chosen on tie
)),
censor_conditions = list(censor_source(
dataset_name = "adsl",
date = EOSDT,
censor = 1,
set_values_to = exprs(EVENTDESC = "End of Study")
)),
set_values_to = exprs(
PARAMCD = "TTSAE",
PARAM = "Time to First Serious AE"
),
source_datasets = list(adsl = adsl, ae = ae)
)
# Check that for USUBJID = "01", the first serious AE selected is the one with AESEQ = 1 (Flu),
# despite the input AE data initially being arranged to show AESEQ=2 (Cough) first.
selected_seq <- result %>%
filter(USUBJID == "01", PARAMCD == "TTSAE") %>%
pull(SRCSEQ)
expect_equal(selected_seq, 1, info = "The order argument should ensure AE with AESEQ=1
is chosen on tie.")
})
## Test 16: produces consistent results regardless of input sort order ----
test_that("derive_param_tte Test 16: produces consistent results regardless of input sort order", {
# Define ADSL dataset
adsl <- tibble::tribble(
~USUBJID, ~TRTSDT, ~TRTEDT, ~EOSDT,
"01", ymd("2020-12-06"), ymd("2021-03-02"), ymd("2021-03-06"),
"02", ymd("2021-01-16"), ymd("2021-01-20"), ymd("2021-02-03")
) %>% mutate(STUDYID = "AB42")
# Define AE dataset with duplicates
ae <- tibble::tribble(
~USUBJID, ~AESTDTC, ~AESEQ, ~AEDECOD,
"01", "2021-01-03", 1, "Flu",
"01", "2021-03-04", 2, "Cough",
"01", "2021-01-03", 3, "Flu"
) %>% mutate(STUDYID = "AB42", AESTDT = ymd(AESTDTC))
# Define event and censor sources
ttae <- event_source(
dataset_name = "ae",
date = AESTDT,
order = exprs(AESTDT, AESEQ),
set_values_to = exprs(
EVENTDESC = "AE",
SRCDOM = "AE",
SRCVAR = "AESTDTC",
SRCSEQ = AESEQ # Ensure AESEQ is included here
)
)
eot <- censor_source(
dataset_name = "adsl",
date = pmin(TRTEDT + days(10), EOSDT),
censor = 1,
set_values_to = exprs(
EVENTDESC = "END OF TRT",
SRCDOM = "ADSL",
SRCVAR = "TRTEDT"
)
)
# Run derive_param_tte with sorted AE dataset
result_sorted <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = arrange(ae, AESEQ)),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
# Run derive_param_tte with reverse-sorted AE dataset
result_unsorted <- derive_param_tte(
dataset_adsl = adsl,
start_date = TRTSDT,
event_conditions = list(ttae),
censor_conditions = list(eot),
source_datasets = list(adsl = adsl, ae = arrange(ae, desc(AESEQ))),
set_values_to = exprs(PARAMCD = "TTAE"),
check_type = "warning"
)
expect_equal(result_sorted, result_unsorted)
})
# list_tte_source_objects ----
## Test 17: error is issued if package does not exist ----
test_that("list_tte_source_objects Test 17: error is issued if package does not exist", {
expect_snapshot(
list_tte_source_objects(package = "tte"),
error = TRUE
)
})
## Test 18: expected objects produced ----
test_that("list_tte_source_objects Test 18: expected objects produced", {
expected_output <- tibble::tribble(
~object, ~dataset_name, ~filter,
"ae_ser_event", "adae", quote(TRTEMFL == "Y" & AESER == "Y"),
"ae_gr2_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "2"),
"ae_sev_event", "adae", quote(TRTEMFL == "Y" & AESEV == "SEVERE"),
"ae_gr4_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "4"),
"ae_gr3_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "3"),
"lastalive_censor", "adsl", NULL,
"ae_event", "adae", quote(TRTEMFL == "Y"),
"death_event", "adsl", quote(DTHFL == "Y"),
"ae_gr35_event", "adae", quote(TRTEMFL == "Y" & ATOXGR %in% c("3", "4", "5")),
"ae_wd_event", "adae", quote(TRTEMFL == "Y" & AEACN == "DRUG WITHDRAWN"),
"ae_gr1_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "1"),
"ae_gr5_event", "adae", quote(TRTEMFL == "Y" & ATOXGR == "5")
) %>%
mutate(
date = case_when(
object == "lastalive_censor" ~ "LSTALVDT",
object == "death_event" ~ "DTHDT",
TRUE ~ "ASTDT"
),
censor = if_else(object == "lastalive_censor", 1, 0),
filter = as.character(filter),
censor = as.integer(censor)
)
observed_output <- list_tte_source_objects(package = "admiral") %>%
select(object, dataset_name, filter, date, censor)
expect_dfs_equal(expected_output, observed_output, keys = c("object"))
})
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.