Add convertKinR0

nlmixr2 · Aug 22, 2024 · b6fdd58 · b6fdd58
1 parent 1bdeae5
commit b6fdd58
Show file tree

Hide file tree

Showing 5 changed files with 58 additions and 13 deletions.
diff --git a/NAMESPACE b/NAMESPACE
@@ -44,6 +44,7 @@ export(camelCombine)
 export(combinePaste2)
 export(convertEmax)
 export(convertEmaxHill)
+export(convertKinR0)
 export(convertLogLin)
 export(convertMM)
 export(convertQuad)

diff --git a/R/convertKinR0.R b/R/convertKinR0.R
@@ -34,7 +34,7 @@ convertKinR0 <- function(ui,
   }
   .modelLines <- .ui$lstExpr
   .w <- .whichDdt(.modelLines, R, start="", end="(0)")
-  if (!length(.w)) {
+  if (length(.w) != 1L) {
     stop(paste0("the model does not have the expected ",
                 R, "(0) expression"),
          call.=FALSE)

diff --git a/man/convertKinR0.Rd b/man/convertKinR0.Rd
diff --git a/man/modeldb.Rd b/man/modeldb.Rd
diff --git a/tests/testthat/test-convertKinR0.R b/tests/testthat/test-convertKinR0.R
@@ -0,0 +1,23 @@
+test_that("test the convertKinR0 function", {
+
+  m <- addIndirect(stim="in") |> convertKinR0()
+
+  expect_true("uR0" %in% names(m$theta))
+
+  m2 <- m |> rxode2::model(a=kin, append=TRUE) |>
+    rxode2::model(R(0) <- kout) |>
+    rxode2::model(d/dt(R) <- -k*R) |>
+    rxode2::model(-R0)
+
+  expect_error(convertKinR0(m2))
+
+  m <- addIndirect(stim="in") |> rxode2::model(-R(0))
+
+  expect_error(m |> convertKinR0())
+
+  m <-  addIndirect(stim="in") |>
+    rxode2::model(-d/dt(R))
+
+  expect_error(m |> convertKinR0())
+
+})