reg <- new_regimen(
amt = 100,
n = 3,
interval = 12,
type = "infusion",
t_inf = 2
)
par <- list(CL = 5, V = 50)
omega <- c(0.1, 0.0, 0.1)
t_obs <- c(2, 48)
test_that("return_event_table=TRUE returns an appropriate event table", {
evtab1 <- sim_ode(
mod_1cmt_iv,
parameters = par,
regimen = reg,
t_obs = t_obs,
return_event_table = TRUE
)
expect_equal(
evtab1,
data.frame(
t = c(0, 2, 12, 14, 24, 26, 48, 48),
dose = c(100, 0, 100, 0, 100, 0, 0, 0),
type = c(1, 1, 1, 1, 1, 1, 0, 0),
dum = c(0, 1, 0, 1, 0, 1, 0, 0),
dose_cmt = c(1, 1, 1, 1, 1, 1, 0, 0),
t_inf = c(2, 0, 2, 0, 2, 0, 0, 0),
evid = c(1, 2, 1, 2, 1, 2, 0, 0),
bioav = c(1, 0, 1, 0, 1, 0, 0, 0),
rate = c(50,-50, 50,-50, 50,-50, 0, 0),
obs_type = c(0, 1, 0, 0, 0, 0, 1, 1)
)
)
})
test_that("return_event_table=TRUE returns an appropriate event table with covariate", {
covs <- list(CRCL = new_covariate(value = c(70, 80), t = c(0, 24)), WT = new_covariate(70))
evtab2 <- sim_ode(
mod_1cmt_iv,
parameters = par,
regimen = reg,
t_obs = t_obs,
covariates = covs,
return_event_table = TRUE
)
expect_equal(
evtab2,
structure(
list(
t = c(0, 0, 2, 12, 14, 24, 24, 26, 48),
dose = c(0, 100, 0, 100, 0, 0, 100, 0, 0),
type = c(0, 1, 1, 1, 1, 0, 1, 1, 0),
dum = c(0, 0, 1, 0, 1, 0, 0, 1, 0),
dose_cmt = c(0, 1, 1, 1, 1, 0, 1, 1, 0),
t_inf = c(0, 2, 0, 2, 0, 0, 2, 0, 0),
evid = c(2, 1, 2, 1, 2, 2, 1, 2, 0),
bioav = c(1, 1, 0, 1, 0, 1, 1, 0, 0),
rate = c(0, 50,-50, 50,-50, 0, 50,-50, 0),
cov_CRCL = c(70, 70, 70, 70, 70, 80, 80, 80, 80),
cov_t_CRCL = c(0, 0, 0, 0, 0, 24, 24, 24, 24),
gradients_CRCL = c(0.416666666666667, 0.416666666666667, 0.416666666666667, 0.416666666666667, 0.416666666666667, 0, 0, 0, 0),
cov_WT = c(70, 70, 70, 70, 70, 70, 70, 70, 70),
cov_t_WT = c(0, 0, 0, 0, 0, 0, 0, 0, 0),
gradients_WT = c(0, 0, 0, 0, 0, 0, 0, 0, 0),
obs_type = c(0, 0, 1, 0, 0, 0, 0, 0, 1)
),
row.names = c(2L, 1L, 3L, 4L, 5L,
6L, 7L, 8L, 9L),
class = "data.frame"
)
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.