inst/models/example_iov.R

In tdmore-dev/tdmore: In silico evaluation of precision dosing

library(tdmore)
nlmixr::nlmixrUI(function(){
  ini({
    TVV1 <- 24.4;
    TVV2 <- 7.01;
    TVQ <- 4.97;
    TVCL <- 9.87;
    ECL ~ 0.194 # This value corresponds to OMEGA_CL (44% SD)
    EV1 ~ 0.287 # This value corresponds to OMEGA_V1 (54% SD)
    ECL_IOV ~ 0.1 #31%
    EV1_IOV ~ 0.1 #31%
    EPS_PROP <- 0.371 # Proportional error (37% SD)
  })
  model({
    CL <- TVCL * exp(ECL + ECL_IOV)
    V1 <- TVV1 * exp(EV1 + EV1_IOV)
    V2 <- TVV2
    Q <- TVQ
    K12 <- Q/V1
    K21 <- Q/V2

    d/dt(center) = - CL/V1 * center - K12*center + K21 * periph
    d/dt(periph) = K12*center - K21 * periph

    CONC = center / V1
    CONC ~ prop(EPS_PROP) # Proportional error linked to the PK model
  })
}) %>% tdmore(iov=c("ECL_IOV", "EV1_IOV"))