inst/doc/v05_parameters.R

In campsismod: Generic Implementation of a PK/PD Model

## ----echo=FALSE, warning=FALSE, message=FALSE, results='hide'-----------------
EXPORT_PNG <- FALSE

## ----echo=TRUE, warning=FALSE, message=FALSE----------------------------------
library(campsismod)

## -----------------------------------------------------------------------------
model <- model_suite$pk$`2cpt_fo`
model

## -----------------------------------------------------------------------------
model %>% find(Theta("CL"))
model %>% find(Omega("KA"))
model %>% find(Sigma("RUV_FIX"))

## -----------------------------------------------------------------------------
model@parameters %>% get_by_index(Theta(index = 5))
model@parameters %>% get_by_index(Omega(index = 1, index2 = 1))
model@parameters %>% get_by_index(Sigma(index = 1, index2 = 1))

## -----------------------------------------------------------------------------
thetaCL <- model %>% find(Theta("CL"))
thetaCL@value

## -----------------------------------------------------------------------------
theta <- Omega(name = "TEST", index = 1, index2 = 1, value = 15, type = "cv%")
theta@value # 15 is returned

theta_standardised <- theta %>% standardise() # Conversion to variance
theta_standardised
theta_standardised@value

## -----------------------------------------------------------------------------
model <- model %>% replace(Theta("KA", value = 2)) # Previous value for KA was 1
model <- model %>% replace(Omega("CL", value = 20, type = "cv%")) # Previous value was a 25% CV
model

## -----------------------------------------------------------------------------
model <- model %>% delete(Theta("KA"))
model <- model %>% delete(Omega("CL"))
model

## -----------------------------------------------------------------------------
tryCatch(
  {
    validObject(model, complete = TRUE)
  },
  error = function(msg) {
    print(msg)
  }
)