Variability levels"

In campsis: Generic PK/PD Simulation Platform CAMPSIS

if (campsis::onCran()) {
  cat("This vignette was not built on CRAN. Please check out the online version [here](https://calvagone.github.io/campsis.doc/articles/v02_uncertainties.html).")
  knitr::knit_exit()
}

library(campsis)

Simulate with all variabilities

Let's use a simple 1-compartment model with absorption compartment to illustrate the different levels of variabilities.

refModel <- model_suite$nonmem$advan2_trans2
refModel

We're going to use a very basic dataset. 1000 mg QD shall be administered for a week.

ds <- Dataset(25) %>%
  add(Bolus(time=0, amount=1000, ii=24, addl=6)) %>%
  add(Observations(times=seq(0,24*7,by=4)))

All subjects are different due to IIV and RUV.

results <- refModel %>% simulate(dataset=ds, seed=1)
spaghettiPlot(results, "CONC_ERR")

Simulate without RUV

Disabling RUV is done as follows:

model <- refModel %>% disable(c("RUV"))
model@parameters

In that case, CONC_ERR (the observed concentration) is identical as CONC (the model-simulated plasma concentration).

results <- model %>% simulate(dataset=ds, seed=1)
spaghettiPlot(results, "CONC_ERR")

spaghettiPlot(results, "CONC")

Simulate without RUV and IIV

Disabling RUV and IIV is done as follows:

model <- refModel %>% disable(c("IIV", "RUV"))
model@parameters

Now, the typical profile is shown for all subjects.

results <- model %>% simulate(dataset=ds, seed=1)
spaghettiPlot(results, "CONC_ERR")

campsis documentation built on April 12, 2025, 1:52 a.m.