certara/tidyvpc

ALQ is not supported for binless()

Closed this issue · 1 comments

When using binless(), if specifying alq data (either alone or in addition to blq data), the vpc$pctalq data frame is not created for the binless() method e.g.,

library(tidyvpc)
obs_data <- obs_data[MDV == 0]
sim_data <- sim_data[MDV == 0]

# Add LLOQ/ULOQ
obs_data$LLOQ <- obs_data[, ifelse(GENDER == "F", 50, 25)]
obs_data$ULOQ <- obs_data[, ifelse(GENDER == "F", 125, 80)]

vpc  <- 
  observed(obs_data, x = TIME, y = DV ) |> 
  simulated(sim_data, y = DV) |>
  censoring(blq = DV < LLOQ, lloq = LLOQ,
             alq = DV > ULOQ, uloq = ULOQ) |> 
  stratify(~ GENDER)  |>
  binless() |>
  vpcstats()

plot(vpc, censoring.type = "both")

Closed by #50