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

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

rxTest({
  model <- rxode2({
    C2 <- centr / V2
    C3 <- peri / V3
    #
    d / dt(depot) <- -KA * depot
    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
    ##
    ## central
    KA <- 2.94E-01
    CL <- 1.86E+01
    V2 <- 4.02E+01
    ## peripheral
    Q <- 1.05E+01
    V3 <- 2.97E+02
    ##
    ## effects
    Kin <- 1
    Kout <- 1
    EC50 <- 200
  })

  event_table <- rxode2::eventTable() %>%
    add.dosing(dose = 10000, nbr.doses = 10, dosing.interval = 12) %>%
    add.sampling(0:240)

  test_that("Bad solve raises error, RxODE#5", {
    tf <- tempfile()
    sink(tf)
    expect_error(
      model %>% solve(c(KA = -10), event_table, maxsteps = 5000),
      rex::rex("could not solve the system")
    )
    sink()
    unlink(tf)
  })
})