diff options
Diffstat (limited to 'man/one_box.Rd')
| -rw-r--r-- | man/one_box.Rd | 21 |
1 files changed, 11 insertions, 10 deletions
diff --git a/man/one_box.Rd b/man/one_box.Rd index cf8a91a..c675186 100644 --- a/man/one_box.Rd +++ b/man/one_box.Rd @@ -62,23 +62,24 @@ Create a time series of decline data 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") plot(pred_2) -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) } |