inst/models/two_error_models.R

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

nlmixr::nlmixrUI(function(){
  ini({
    TVKA <- 3.7
    TVV1 <- 61
    TVQ <- 10
    TVCL <- 3.7
    ECL ~ 0.0784 #ETA1, 28%
    EV1 ~ 0.0361 #ETA2, 19%
    EPS_PROP <- 0.001 # Proportional error, 10% SD
    EPS_PROP2 <- 0.001 # Proportional error, 10% SD
  })
  model({
    KA <- TVKA
    CL <- TVCL * (WT/70)^0.75 * exp(ECL)
    V1 <- TVV1 * exp(EV1)
    V2 <- V1
    Q <- TVQ
    K12 <- Q/V1
    K21 <- Q/V2

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

    CONC_C = center / V1
    CONC_C ~ prop(EPS_PROP) | center

    CONC_P = periph / V2
    CONC_P ~ prop(EPS_PROP2) | periph
  })
})