tests/multiple_regimens.R

In ronkeizer/PKPDsim: Simulate ODE models

library(testit)
library(PKPDsim)
Sys.setenv("R_TESTS" = "")

cov_table <- data.frame(WT = rnorm(10, 70, 5))
multi_regs <- list()
for(i in seq(cov_table$WT)) {
  multi_regs[[i]] <- new_regimen(amt = 10 * cov_table$WT[i], interval = 12, type = "infusion")
}
class(multi_regs) <- "regimen_multiple"

mod <- new_ode_model(code = "
  dAdt[1] = -(CL/V) * A[1];
", obs = list(cmt = 1, scale = "V"),
   covariates = list(WT = new_covariate(70)),
   cpp_show_code=F
)
par <- list(CL = 5, V = 50)
reg <- new_regimen(amt = 2000, interval = 24, type = "infusion")
covariates = list(WT = new_covariate(1))
res <- sim_ode(ode = mod,
               parameters = par,
               covariates = covariates,
               regimen = multi_regs,
               only_obs = TRUE
               #regimen = reg
               )

assert("all 10 subjects", length(unique(res$id)) == 10)
assert("no NAs", any(!is.na(res$y)))