simulate.pkmod: Simulate method for pkmod objects

In tci: Target Controlled Infusion (TCI)

Simulate method for pkmod objects

Description

Simulate observations from a pkmod object

Usage

## S3 method for class 'pkmod'
simulate(
  object,
  nsim = 1,
  seed = NULL,
  ...,
  inf,
  tms,
  obs_cmpt = 1,
  resp_bounds = NULL
)

Arguments

object	An object with class "pkmod" generated by 'pkmod'
nsim	Number of observations to simulate at each time point. Defaults to 1.
seed	An integer used to initialize the random number generator.
...	Arguments passed to 'update.pkmod'.
inf	A matrix of infusion rates with columns 'begin', 'end', and 'inf_rate'. This can be created manually, by 'inf_manual', or by 'inf_tci'.
tms	Times at which to simulate observations.
obs_cmpt	Integer value indicating compartment in which observations are taken. Overridden if a PD model is included.
resp_bounds	Optional vector of two values indicating minimum and maximum values possible for the response.

Examples

# simulate data from a 2 compartment model with multiplicative error
dose <- inf_manual(inf_tms = c(0,0.5,4,4.5,10), inf_rate = c(100,0,80,0,0))
my_mod <- pkmod(pars_pk = c(CL = 10, V1 = 10, Q2 = 4, V2 = 30))
inf <- inf_tci(my_mod, target_vals = c(2,3,4,4), target_tms = c(0,2,3,10), "plasma")
simulate(my_mod, nsim = 3, inf = inf, tms = c(1,2,4,6,10), sigma_mult = 0.2)
# simulate with PD model
my_mod_pd <- pkmod(pars_pk = c(v1 = 8.995, v2 = 17.297, v3 = 120.963, cl = 1.382,
q2 = 0.919, q3 = 0.609, ke0 = 1.289),
pars_pd = c(c50 = 2.8, gamma = 1.47, gamma2 = 1.89, e0 = 93, emx = 93),
pdfn = emax, pdinv = emax_inv, ecmpt = 4)
simulate(my_mod_pd, inf = dose, tms = c(1.5,2.5,3))

tci documentation built on Aug. 15, 2022, 9:09 a.m.