From d28ce9f8ad6f9573e403ebd8eb637ecd5e5b0e02 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Mon, 21 Dec 2020 06:02:10 +0100
Subject: plot.mixed: Possibility to overlay predictions

---
 man/plot.mixed.mmkin.Rd | 13 ++++++++++++-
 1 file changed, 12 insertions(+), 1 deletion(-)

(limited to 'man/plot.mixed.mmkin.Rd')

diff --git a/man/plot.mixed.mmkin.Rd b/man/plot.mixed.mmkin.Rd
index b90a4b3a..c7b2344f 100644
--- a/man/plot.mixed.mmkin.Rd
+++ b/man/plot.mixed.mmkin.Rd
@@ -12,6 +12,7 @@
   xlab = "Time",
   xlim = range(x$data$time),
   resplot = c("predicted", "time"),
+  pred_over = NULL,
   ymax = "auto",
   maxabs = "auto",
   ncol.legend = ifelse(length(i) <= 3, length(i) + 1, ifelse(length(i) <= 8, 3, 4)),
@@ -45,6 +46,9 @@ variables in the model.}
 \item{resplot}{Should the residuals plotted against time or against
 predicted values?}
 
+\item{pred_over}{Named list of alternative predictions as obtained
+from \link{mkinpredict} with a compatible \link{mkinmod}.}
+
 \item{ymax}{Vector of maximum y axis values}
 
 \item{maxabs}{Maximum absolute value of the residuals. This is used for the
@@ -90,8 +94,15 @@ plot(f[, 3:4], standardized = TRUE)
 f_nlme <- nlme(f, control = list(pnlsMaxIter = 120, tolerance = 1e-3))
 plot(f_nlme)
 
-f_saem <- saem(f)
+f_saem <- saem(f, transformations = "saemix")
 plot(f_saem)
+
+# We can overlay the two variants if we generate predictions
+pred_nlme <- mkinpredict(dfop_sfo,
+  f_nlme$bparms.optim[-1],
+  c(parent = f_nlme$bparms.optim[[1]], A1 = 0),
+  seq(0, 180, by = 0.2))
+plot(f_saem, pred_over = list(nlme = pred_nlme))
 }
 }
 \author{
-- 
cgit v1.2.1