Improve examle code for `one_box`HEADmain

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
 ```