Workpackage 2

Published

October 20, 2022

Package 1B

Code
library(PKPDsim)
library(PKPDplot)
library (ggplot2)
pk1 <- new_ode_model(model = "pk_2cmt_iv") # ceftazidime
pk2 <- new_ode_model(model = "pk_2cmt_iv") # avibactam
p <- list(
  CL = 6.95,
  V  = 10.5,
  Q  = 31.5,
  V2 = 7.57)

p2 <- list(
  CL = 10.2,
  V  = 11.1,
  Q  = 5.44,
  V2 = 6.91)
n <- 100

omega <- cv_to_omega(list(CL=0.423, V=1.05, Q=0.315, V2=1.10))

omega2 <- cv_to_omega(list(CL=0.591, V=1.07, Q=1.22, V2=6.91))
r1 <- new_regimen(
  amt = 2000,
  interval = 8,
  n = 5,
  type = "infusion",
  t_inf = c(2,2,2,2,2)
)

r2 <- new_regimen(
  amt = 500,
  interval = 8,
  n = 5,
  type = "infusion",
  t_inf = c(2,2,2,2,2)
)
dat <- sim(
  ode = pk1,
  parameters = p,
  omega = omega,
  n_ind = 20,
  res_var = list(prop=0.1, add=7.5),
  regimen = r1, 
  covariates = list("SCR" = new_covariate(value = 80)), output_include = list(parameters=TRUE, variable=TRUE, covariates=TRUE)
)

dat2 <- sim(
  ode = pk2,
  parameters = p2,
  omega = omega2,
  n_ind = 20,
  res_var = list(prop=0.17, add=0.5),
  regimen = r2, 
  covariates = list("SCR" = new_covariate(value = 80)), output_include = list(parameters=TRUE, variable=TRUE, covariates=TRUE)
)

plot(dat, show_population = list(obs = TRUE, spaghetti = FALSE, ci = TRUE, 
                                 median =TRUE, regimen = TRUE), ci = c(0.05, 0.95))