Admiral Bds

library(admiral)
library(dplyr)
library(lubridate)
library(pharmaversesdtm)

# Load source domain — replace with vs/lb per dataset being derived
vs  <- pharmaversesdtm::vs
adsl <- <loaded ADSL dataset>

# Remove DOMAIN to avoid conflicts in derive_vars_merged() calls
vs <- select(vs, -DOMAIN)

# REVIEW: Confirm which population flags and ADSL variables are required by
#   the ADaM spec for this dataset. Add or remove from new_vars accordingly.
advs <- vs |>
  derive_vars_merged(
    dataset_add = adsl,
    by_vars     = exprs(STUDYID, USUBJID),
    new_vars    = exprs(TRTSDT, TRTEDT, TRT01P, TRT01PN, TRT01A, TRT01AN,
                        SAFFL, ITTFL)
  )

# REVIEW: PARAMCD mapping must match the ADaM spec parameter list exactly.
#   Confirm VSTESTCD values in VS and align with ADaM PARAMCD conventions.
#   Remove parameters not in scope for this study.
param_lookup <- tibble::tribble(
  ~VSTESTCD, ~PARAMCD, ~PARAM,                    ~PARAMN,
  "SYSBP",   "SYSBP",  "Systolic Blood Pressure",  1L,
  "DIABP",   "DIABP",  "Diastolic Blood Pressure",  2L,
  "PULSE",   "PULSE",  "Pulse Rate",                3L,
  "WEIGHT",  "WEIGHT", "Weight",                    4L,
  "HEIGHT",  "HEIGHT", "Height",                    5L,
  "TEMP",    "TEMP",   "Temperature",               6L
)

advs <- advs |>
  derive_vars_merged_lookup(
    dataset_add = param_lookup,
    by_vars     = exprs(VSTESTCD),
    new_vars    = exprs(PARAMCD, PARAM, PARAMN)
  ) |>
  filter(!is.na(PARAMCD))   # drop records for out-of-scope tests

advs <- advs |>
  mutate(
    AVAL  = VSSTRESN,    # numeric result in standard units
    AVALC = VSSTRESC,    # character result (for non-numeric or verbatim)
    AVALU = VSSTRESU     # analysis value units
  )

advs <- advs |>
  # ADT: analysis date from VSDTC
  derive_vars_dt(
    dtc             = VSDTC,
    new_vars_prefix = "A",
    date_imputation = "first",
    flag_imputation = "auto"
  ) |>
  # ADY: study day relative to TRTSDT
  derive_vars_dy(
    reference_date = TRTSDT,
    source_vars    = exprs(ADT)
  )

# REVIEW: Visit map must come from the ADaM spec. The example below is
#   illustrative. Confirm VISIT names and AVISITN codes against the study CRF
#   and ADaM spec before use.
visit_map <- tibble::tribble(
  ~VISIT,        ~AVISIT,       ~AVISITN,
  "SCREENING 1", "Screening",   -1L,
  "BASELINE",    "Baseline",     0L,
  "WEEK 2",      "Week 2",       2L,
  "WEEK 4",      "Week 4",       4L,
  "WEEK 8",      "Week 8",       8L,
  "WEEK 16",     "Week 16",     16L,
  "WEEK 26",     "Week 26",     26L
)

advs <- advs |>
  derive_vars_merged_lookup(
    dataset_add = visit_map,
    by_vars     = exprs(VISIT),
    new_vars    = exprs(AVISIT, AVISITN)
  )

# REVIEW: Baseline window definition is protocol-specific. Confirm whether
#   the baseline is the last pre-dose record (ADT <= TRTSDT), last on-or-before
#   treatment start, or a specific visit (e.g. DAY 1 only). Adjust filter below.
advs <- advs |>
  restrict_derivation(
    derivation = derive_var_extreme_flag,
    args = params(
      by_vars  = exprs(STUDYID, USUBJID, PARAMCD, BASETYPE),
      order    = exprs(ADT, AVISITN),
      new_var  = ABLFL,
      mode     = "last"
    ),
    filter = ADT <= TRTSDT & !is.na(AVAL)
  )

advs <- advs |>
  derive_var_base(
    by_vars  = exprs(STUDYID, USUBJID, PARAMCD, BASETYPE),
    source_var = AVAL,
    new_var    = BASE
  ) |>
  derive_var_base(
    by_vars    = exprs(STUDYID, USUBJID, PARAMCD, BASETYPE),
    source_var = AVALC,
    new_var    = BASEC
  )

advs <- advs |>
  derive_var_chg() |>    # CHG = AVAL - BASE
  derive_var_pchg()      # PCHG = CHG / BASE * 100; NA if BASE = 0 or NA

# REVIEW: ANL01FL definition is protocol-specific. The condition below
#   (on-treatment, non-baseline) is a common starting point. Confirm against
#   the SAP before use.
advs <- advs |>
  restrict_derivation(
    derivation = derive_var_extreme_flag,
    args = params(
      by_vars  = exprs(STUDYID, USUBJID, PARAMCD),
      order    = exprs(ADT, AVISITN),
      new_var  = ANL01FL,
      mode     = "last"
    ),
    filter = ADT >= TRTSDT & !is.na(AVAL) & is.na(DTYPE)
  )

# Uniqueness assertion — adjust key variables per spec
stopifnot(
  advs |>
    filter(is.na(DTYPE)) |>
    count(STUDYID, USUBJID, PARAMCD, AVISITN, VSPOS) |>
    filter(n > 1) |>
    nrow() == 0
)