diff options
Diffstat (limited to 'docs/reference/one_box.md')
| -rw-r--r-- | docs/reference/one_box.md | 36 |
1 files changed, 13 insertions, 23 deletions
diff --git a/docs/reference/one_box.md b/docs/reference/one_box.md index f956b8f..89ccecf 100644 --- a/docs/reference/one_box.md +++ b/docs/reference/one_box.md @@ -87,39 +87,29 @@ 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")) +#> Temporary DLL for differentials generated and loaded +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")) #> Temporary DLL for differentials generated and loaded fit_2 <- mkinfit(m_2, FOCUS_2006_D, quiet = TRUE) #> Warning: Observations with value of zero were removed from the data pred_2 <- one_box(fit_2, ini = "model") plot(pred_2) - -library(mkin) -SFO_SFO <- mkinmod( - parent = mkinsub("SFO", to = "m1"), - m1 = mkinsub("SFO")) -#> Temporary DLL for differentials generated and loaded -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) -#> Time Series: -#> Start = c(0, 1) -#> End = c(0, 6) -#> Frequency = 100 -#> parent m1 -#> 0.00 100.00000 0.00000000 -#> 0.01 99.85011 0.07494003 -#> 0.02 99.70045 0.14976024 -#> 0.03 99.55101 0.22446081 -#> 0.04 99.40180 0.29904192 -#> 0.05 99.25281 0.37350376 ``` |