Variability levels"

In campsis: Generic PK/PD Simulation Platform Campsis

if (campsis::on_cran()) {
  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.

ref_model <- model_suite$nonmem$advan2_trans2
ref_model

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 <- simulate(model = ref_model, dataset = ds, seed = 1)
spaghetti_plot(results, "CONC_ERR")

Simulate without RUV

Disabling RUV is done as follows:

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

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

results <- simulate(model = model, dataset = ds, seed = 1)
spaghetti_plot(results, "CONC_ERR")

spaghetti_plot(results, "CONC")

Simulate without RUV and IIV

Disabling RUV and IIV is done as follows:

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

Now, the typical profile is shown for all subjects.

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

campsis documentation built on Aug. 5, 2026, 9:07 a.m.