Feature request: New logslope & c1 impute_start methods #350

NAMESPACE

@@ -90,8 +90,10 @@ export(PKNCA.choose.option)
 export(PKNCA.options)
 export(PKNCA.options.describe)
 export(PKNCA.set.summary)
+export(PKNCA_impute_method_start_c1)
 export(PKNCA_impute_method_start_cmin)
 export(PKNCA_impute_method_start_conc0)
+export(PKNCA_impute_method_start_logslope)
 export(PKNCA_impute_method_start_predose)
 export(PKNCAconc)
 export(PKNCAdata)

NEWS.md

@@ -6,6 +6,7 @@ the dosing including dose amount and route.
 
 # PKNCA 0.11.0.9000
 
+* New default PKNCA impute start methods for C1 and IV bolus logslope to add for intervals without a start concentration
 * PKNCA options can now use `tmax` as a reference for BLQ handling by using new names in the `conc.blq` argument (`before.tmax`,`after.tmax`)
 * PKNCA will now give an error when there are unexpected interval columns.
   The `keep_interval_cols` option can be used to mitigate this error.

R/impute.R

@@ -120,6 +120,46 @@ PKNCA_impute_method_start_predose <- function(conc, time, start, end, conc.group
   ret
 }
 
+#' @describeIn PKNCA_impute_method Imputes based on a logarithmic slope on the
+#'   first points a time zero concentration (usually in IV bolus at dose time)
+#' @export
+PKNCA_impute_method_start_logslope <- function(conc, time, start, end, ..., options = list()) { # nolint
+
+  ret <- data.frame(conc = conc, time = time)
+  mask_start <- time %in% start
+  if (!any(mask_start)) {
+    all_concs <- conc[time >= start  &  time <= end]
+    all_times <- time[time >= start  &  time <= end]
+    if (!all(is.na(all_concs))) {
+      c0 <- PKNCA::pk.calc.c0(all_concs, all_times, method = "logslope")
+      if (!is.na(c0)) {
+        ret <- rbind(ret, data.frame(time = start, conc = c0))
+        ret <- ret[order(ret$time), ]
+      }
+    }
+  }
+  ret
+}
+
+#' @describeIn PKNCA_impute_method Shift the following concentration to a start
+#'   to become the time zero concentration (rarely used; non-monodecay IV bolus)
+#' @return A data frame with imputed start concentration.
+#' @export
+PKNCA_impute_method_start_c1 <- function(conc, time, start, end, ..., options = list()) { # nolint
+  ret <- data.frame(conc = conc, time = time)
+  mask_start <- time %in% start
+  if (!any(mask_start)) {
+    all_concs <- conc[time >= start  &  time <= end]
+    all_times <- time[time >= start  &  time <= end]
+    if (!all(is.na(all_concs))) {
+      c1 <- all_concs[which.min(all_times)]
+      ret <- rbind(ret, data.frame(time = start, conc = c1))
+      ret <- ret[order(ret$time), ]
+    }
+  }
+  ret
+}
+
 #' Separate out a vector of PKNCA imputation methods into a list of functions
 #'
 #' An error will be raised if the functions are not found.

tests/testthat/test-impute.R

@@ -92,6 +92,58 @@ test_that("PKNCA_impute_method_start_cmin", {
 
 })
 
+test_that("PKNCA_impute_method_start_logslope", {
+  # No imputation when start is in the data
+  expect_equal(
+    PKNCA_impute_method_start_logslope(conc = 3:1, time = 0:2, start = 0, end = 2),
+    data.frame(conc = 3:1, time = 0:2)
+  )
+  # Impute when start is not in the data
+  expect_equal(
+    PKNCA_impute_method_start_logslope(conc = 3:1, time = 1:3, start = 0, end = 3),
+    data.frame(conc = c(4.5, 3:1), time = 0:3),
+    ignore_attr = TRUE
+  )
+  # Data outside the interval are ignored (before interval)
+  expect_equal(
+    PKNCA_impute_method_start_logslope(conc = c(0, 2:1), time = c(-1, 1:2), start = 0, end = 2),
+    data.frame(conc = c(0, 4, 2:1), time = c(-1, 0, 1:2)),
+    ignore_attr = TRUE
+  )
+  # No modification if no C1 -> C2 decline in samples
+  expect_equal(
+    PKNCA_impute_method_start_logslope(conc = c(1, 1, 1), time = 1:3, start = 0, end = 3),
+    data.frame(conc = c(1, 1, 1), time = 1:3),
+    ignore_attr = TRUE
+  )
+  # No modification if C1 = C2 in samples
+  expect_equal(
+    PKNCA_impute_method_start_logslope(conc = c(3, 3, 1), time = 1:3, start = 0, end = 3),
+    data.frame(conc = c(3, 3, 1), time = 1:3),
+    ignore_attr = TRUE
+  )
+})
+
+test_that("PKNCA_impute_method_start_c1", {
+  # No imputation when start is in the data
+  expect_equal(
+    PKNCA_impute_method_start_c1(conc = 1:3, time = 0:2, start = 0, end = 2),
+    data.frame(conc = 1:3, time = 0:2)
+  )
+  # Impute when start is not in the data
+  expect_equal(
+    PKNCA_impute_method_start_c1(conc = 1:3, time = 1:3, start = 0, end = 3),
+    data.frame(conc = c(1, 1:3), time = 0:3),
+    ignore_attr = TRUE
+  )
+  # Data outside the interval are ignored (before interval)
+  expect_equal(
+    PKNCA_impute_method_start_c1(conc = 1:3, time = c(-1, 1:2), start = 0, end = 2),
+    data.frame(conc = c(1, 2, 2:3), time = c(-1, 0, 1:2)),
+    ignore_attr = TRUE
+  )
+})
+
 test_that("PKNCA_impute_fun_list", {
   expect_equal(
     PKNCA_impute_fun_list(NA_character_),