quiet |
@@ -267,36 +261,33 @@ using mmkin.
state.ini = c(parent = 100), fixed_initials = "parent", quiet = TRUE)
f_saem_p0_fixed <- saem(f_mmkin_parent_p0_fixed)
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:01 2021"
+#> [1] "Mon Jan 25 14:41:42 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:03 2021"
+#> [1] "Mon Jan 25 14:41:43 2021"
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:04 2021"
+#> [1] "Mon Jan 25 14:41:45 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:06 2021"
f_saem_fomc <- saem(f_mmkin_parent["FOMC", ])
+#> [1] "Mon Jan 25 14:41:46 2021"
f_saem_fomc <- saem(f_mmkin_parent["FOMC", ])
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:06 2021"
+#> [1] "Mon Jan 25 14:41:47 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:08 2021"
f_saem_dfop <- saem(f_mmkin_parent["DFOP", ])
+#> [1] "Mon Jan 25 14:41:49 2021"
f_saem_dfop <- saem(f_mmkin_parent["DFOP", ])
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:09 2021"
+#> [1] "Mon Jan 25 14:41:49 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:12 2021"
+#> [1] "Mon Jan 25 14:41:52 2021"
#> Package saemix, version 3.1.9000
#> please direct bugs, questions and feedback to emmanuelle.comets@inserm.fr
#> Likelihoods computed by importance sampling
#> AIC BIC
-#> 1 624.2484 622.2956
-#> 2 467.7096 464.9757
-#> 3 495.4373 491.9222
#> Error in compare.saemix(list(f_saem_sfo$so, f_saem_fomc$so, f_saem_dfop$so)): 'compare.saemix' requires at least two models.
#> Plotting convergence plots
#> Plotting individual fits
#> Simulating data using nsim = 1000 simulated datasets
@@ -333,13 +324,11 @@ using
mmkin.
f_mmkin_parent_tc <- update(f_mmkin_parent, error_model
= "tc")
f_saem_fomc_tc <- saem(f_mmkin_parent_tc["FOMC",
])
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:14 2021"
+#> [1] "Mon Jan 25 14:41:55 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:19 2021"
#> Likelihoods computed by importance sampling
#> AIC BIC
-#> 1 467.7096 464.9757
-#> 2 469.5208 466.3963
+#> [1] "Mon Jan 25 14:42:00 2021"
#> Error in compare.saemix(list(f_saem_fomc$so, f_saem_fomc_tc$so)): 'compare.saemix' requires at least two models.
#> Temporary DLL for differentials generated and loaded
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:22 2021"
+#> [1] "Mon Jan 25 14:42:02 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:27 2021"
f_saem_dfop_sfo <- saem(f_mmkin["DFOP-SFO", ])
+#> [1] "Mon Jan 25 14:42:07 2021"
f_saem_dfop_sfo <- saem(f_mmkin["DFOP-SFO", ])
#> Running main SAEM algorithm
-#> [1] "Mon Jan 11 12:42:28 2021"
+#> [1] "Mon Jan 25 14:42:08 2021"
#> ....
#> Minimisation finished
-#> [1] "Mon Jan 11 12:42:37 2021"
# We can use print, plot and summary methods to check the results
+#> [1] "Mon Jan 25 14:42:17 2021"
#> Kinetic nonlinear mixed-effects model fit by SAEM
#> Structural model:
@@ -408,8 +397,8 @@ using
mmkin.
#> saemix version used for fitting: 3.1.9000
#> mkin version used for pre-fitting: 0.9.50.4
#> R version used for fitting: 4.0.3
-#> Date of fit: Mon Jan 11 12:42:38 2021
-#> Date of summary: Mon Jan 11 12:42:38 2021
+#> Date of fit: Mon Jan 25 14:42:18 2021
+#> Date of summary: Mon Jan 25 14:42:18 2021
#>
#> Equations:
#> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -424,7 +413,7 @@ using
mmkin.
#>
#> Model predictions using solution type analytical
#>
-#> Fitted in 10.092 s using 300, 100 iterations
+#> Fitted in 9.954 s using 300, 100 iterations
#>
#> Variance model: Constant variance
#>
diff --git a/man/plot.mkinfit.Rd b/man/plot.mkinfit.Rd
index bf38c5a5..9fd4c038 100644
--- a/man/plot.mkinfit.Rd
+++ b/man/plot.mkinfit.Rd
@@ -12,7 +12,7 @@
fit = x,
obs_vars = names(fit$mkinmod$map),
xlab = "Time",
- ylab = "Observed",
+ ylab = "Residue",
xlim = range(fit$data$time),
ylim = "default",
col_obs = 1:length(obs_vars),
diff --git a/man/plot.mmkin.Rd b/man/plot.mmkin.Rd
index 998e25e1..c127be8c 100644
--- a/man/plot.mmkin.Rd
+++ b/man/plot.mmkin.Rd
@@ -10,6 +10,7 @@ of an mmkin object}
main = "auto",
legends = 1,
resplot = c("time", "errmod"),
+ ylab = "Residue",
standardized = FALSE,
show_errmin = TRUE,
errmin_var = "All data",
@@ -32,6 +33,8 @@ column.}
\code{\link{mkinresplot}}, or as squared residuals against predicted
values, with the error model, using \code{\link{mkinerrplot}}.}
+\item{ylab}{Label for the y axis.}
+
\item{standardized}{Should the residuals be standardized? This option
is passed to \code{\link{mkinresplot}}, it only takes effect if
\code{resplot = "time"}.}
diff --git a/tests/figs/evaluations-according-to-2015-nafta-guidance/nafta-sop-appendix-b.svg b/tests/figs/evaluations-according-to-2015-nafta-guidance/nafta-sop-appendix-b.svg
index 617e39b3..d63d8da2 100644
--- a/tests/figs/evaluations-according-to-2015-nafta-guidance/nafta-sop-appendix-b.svg
+++ b/tests/figs/evaluations-according-to-2015-nafta-guidance/nafta-sop-appendix-b.svg
@@ -47,7 +47,7 @@
Time
-
Observed
+
Residue
@@ -195,7 +195,7 @@
Time
-
Observed
+
Residue
@@ -343,7 +343,7 @@
Time
-
Observed
+
Residue
diff --git a/tests/figs/evaluations-according-to-2015-nafta-guidance/plot-nafta-analysis.svg b/tests/figs/evaluations-according-to-2015-nafta-guidance/plot-nafta-analysis.svg
index 5e98487b..283c5c9e 100644
--- a/tests/figs/evaluations-according-to-2015-nafta-guidance/plot-nafta-analysis.svg
+++ b/tests/figs/evaluations-according-to-2015-nafta-guidance/plot-nafta-analysis.svg
@@ -41,7 +41,7 @@
Time
-
Observed
+
Residue
@@ -189,7 +189,7 @@
Time
-
Observed
+
Residue
@@ -337,7 +337,7 @@
Time
-
Observed
+
Residue
diff --git a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-defaults.svg b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-defaults.svg
index 82577502..ef238a51 100644
--- a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-defaults.svg
+++ b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-defaults.svg
@@ -40,7 +40,7 @@
80
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-residuals-like-in-gmkin.svg b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-residuals-like-in-gmkin.svg
index 8d16bfd3..a84e170d 100644
--- a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-residuals-like-in-gmkin.svg
+++ b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-residuals-like-in-gmkin.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-sep-true.svg b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-sep-true.svg
index 327161da..9b16c583 100644
--- a/tests/figs/plotting/mkinfit-plot-for-focus-c-with-sep-true.svg
+++ b/tests/figs/plotting/mkinfit-plot-for-focus-c-with-sep-true.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/mmkin-plot-for-focus-c.svg b/tests/figs/plotting/mmkin-plot-for-focus-c.svg
index 1c697328..0f824e6c 100644
--- a/tests/figs/plotting/mmkin-plot-for-focus-c.svg
+++ b/tests/figs/plotting/mmkin-plot-for-focus-c.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
@@ -187,7 +187,7 @@
Time
-Observed
+Residue
@@ -327,7 +327,7 @@
Time
-Observed
+Residue
@@ -466,7 +466,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/mmkin-plot-for-sfo-focus-c-and-d.svg b/tests/figs/plotting/mmkin-plot-for-sfo-focus-c-and-d.svg
index 72496e86..40aa6bef 100644
--- a/tests/figs/plotting/mmkin-plot-for-sfo-focus-c-and-d.svg
+++ b/tests/figs/plotting/mmkin-plot-for-sfo-focus-c-and-d.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
@@ -183,7 +183,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/plot-err-for-focus-d.svg b/tests/figs/plotting/plot-err-for-focus-d.svg
index e7f05168..53049f6e 100644
--- a/tests/figs/plotting/plot-err-for-focus-d.svg
+++ b/tests/figs/plotting/plot-err-for-focus-d.svg
@@ -47,7 +47,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/plot-errmod-with-focus-c-tc.svg b/tests/figs/plotting/plot-errmod-with-focus-c-tc.svg
index 8adefc27..b0d84808 100644
--- a/tests/figs/plotting/plot-errmod-with-focus-c-tc.svg
+++ b/tests/figs/plotting/plot-errmod-with-focus-c-tc.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/plot-errmod-with-focus-d-obs-eigen.svg b/tests/figs/plotting/plot-errmod-with-focus-d-obs-eigen.svg
index afd3e064..d5e1b6b2 100644
--- a/tests/figs/plotting/plot-errmod-with-focus-d-obs-eigen.svg
+++ b/tests/figs/plotting/plot-errmod-with-focus-d-obs-eigen.svg
@@ -47,7 +47,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/plot-res-for-focus-c.svg b/tests/figs/plotting/plot-res-for-focus-c.svg
index 43b94055..226d9014 100644
--- a/tests/figs/plotting/plot-res-for-focus-c.svg
+++ b/tests/figs/plotting/plot-res-for-focus-c.svg
@@ -45,7 +45,7 @@
Time
-Observed
+Residue
diff --git a/tests/figs/plotting/plot-res-for-focus-d.svg b/tests/figs/plotting/plot-res-for-focus-d.svg
index ea3388ed..6504365c 100644
--- a/tests/figs/plotting/plot-res-for-focus-d.svg
+++ b/tests/figs/plotting/plot-res-for-focus-d.svg
@@ -47,7 +47,7 @@
Time
-Observed
+Residue
--
cgit v1.2.1