tests/testthat/test-npde.R

In nlmixr: Nonlinear Mixed Effects Models in Population PK/PD

nlmixrTest(
{
  context("npde runs successfully")

  one.compartment <- function() {
    ini({
      tka <- 0.45 # Log Ka
      tcl <- 1 # Log Cl
      tv <- 3.45 # Log V
      eta.ka ~ 0.6
      eta.cl ~ 0.3
      eta.v ~ 0.1
      add.sd <- 0.7
    })
    model({
      ka <- exp(tka + eta.ka)
      cl <- exp(tcl + eta.cl)
      v <- exp(tv + eta.v)
      d/dt(depot) <- -ka * depot
      d/dt(center) <- ka * depot - cl/v * center
      cp <- center/v
      cp ~ add(add.sd)
    })
  }

  expect_error(suppressWarnings(nlmixr(one.compartment, theo_sd,
                                       est = "focei",
                                       control = list(print = 0),
                                       table=tableControl(npde=TRUE)
                                       )), NA)

  expect_error(suppressWarnings(nlmixr(one.compartment, theo_sd,
                                       est="focei",
                                       control = list(print = 0)
                                       )) %>% addNpde(), NA)
},
test = "lvl2"
)