diff options
Diffstat (limited to 'R/twa.R')
| -rw-r--r-- | R/twa.R | 22 |
1 files changed, 11 insertions, 11 deletions
@@ -21,13 +21,22 @@ #' pred_0 <- one_box(10) #' plot(pred_0) #' -#' # Use a fitted mkinfit model +#' # Use a custom mkin model #' require(mkin) +#' SFO_SFO <- mkinmod( +#' parent = mkinsub("SFO", to = "m1"), +#' m1 = mkinsub("SFO")) +#' c_0 = c(parent = 100, m1 = 0) +#' deg_parms = c(k_parent = 0.15, f_parent_to_m1 = 0.5, k_m1 = 0.01) +#' pred_sfo_sfo <- one_box(SFO_SFO, odeparms = deg_parms, ini = c_0) +#' plot(pred_sfo_sfo) +#' +#' # Use a fitted mkin model #' fit <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE) #' pred_1 <- one_box(fit) #' plot(pred_1) #' -#' # Use a model with more than one observed variable +#' # Use a fitted model with a metabolite #' m_2 <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO")) #' fit_2 <- mkinfit(m_2, FOCUS_2006_D, quiet = TRUE) #' pred_2 <- one_box(fit_2, ini = "model") @@ -86,15 +95,6 @@ one_box.character <- function(x, ini = 1, parms, ..., #' @param solution_type Will be passed to [mkinpredict] #' @importFrom mkin mkinpredict #' @export -#' @examples -#' library(mkin) -#' SFO_SFO <- mkinmod( -#' parent = mkinsub("SFO", to = "m1"), -#' m1 = mkinsub("SFO")) -#' c_0 = c(parent = 100, m1 = 0) -#' deg_parms = c(k_parent = 0.15, f_parent_to_m1 = 0.5, k_m1 = 0.01) -#' sfo_sfo_box <- one_box(SFO_SFO, odeparms = deg_parms, ini = c_0) -#' head(sfo_sfo_box) one_box.mkinmod <- function(x, ini = c(1, rep(0, length(x$diffs) - 1)), odeparms, |