Improve examle code for `one_box`HEADmain

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)
 }