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")


Calvagone/campsis documentation built on April 17, 2024, 5:33 a.m.