tests/testthat/test-rxode-issue-272.R

In nlmixr2/rxode2: Facilities for Simulating from ODE-Based Models

rxTest({
  test_that("issue RxODE#272 and RxODE#273", {
    m1 <- rxode2({
      KA <- 2.94E-01
      CL <- 1.86E+01
      V2 <- 4.02E+01
      Q <- 1.05E+01
      V3 <- 2.97E+02
      Kin <- 1
      Kout <- 1
      EC50 <- 200
      ## Added modeled bioavaiblity, duration and rate
      fdepot <- 1
      durDepot <- 8
      rateDepot <- 1250
      C2 <- centr / V2
      C3 <- peri / V3
      d / dt(depot) <- -KA * depot
      f(depot) <- fdepot
      dur(depot) <- durDepot
      rate(depot) <- rateDepot
      d / dt(centr) <- KA * depot - CL * C2 - Q * C2 + Q * C3
      d / dt(peri) <- Q * C2 - Q * C3
      d / dt(eff) <- Kin - Kout * (1 - C2 / (EC50 + C2)) * eff
      eff(0) <- 1
    })
    
    ev <- et(timeUnits = "hr") %>%
      et(amt = 10000, ii = 12, addl = 3) %>%
      et(time = 6, cmt = "-depot", evid = 2, ii = 12, addl = 3) %>%
      et(seq(0, 24, length.out = 100))
    
    expect_error(s <- rxSolve(m1, ev), NA)
    expect_false(s$C2[3] == 0) # 273
  })
})