From 48c463680b51fa767b4cd7bd62865f192d0354ac Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sat, 6 Feb 2021 18:30:32 +0100 Subject: Reintroduce interface to saemix Also after the upgrade from buster to bullseye of my local system, some test results for saemix have changed. --- man/plot.mixed.mmkin.Rd | 11 ++++++++++- 1 file changed, 10 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 87a82286..b1200729 100644 --- a/man/plot.mixed.mmkin.Rd +++ b/man/plot.mixed.mmkin.Rd @@ -27,7 +27,7 @@ ) } \arguments{ -\item{x}{An object of class \link{mixed.mmkin}, \link{nlme.mmkin}} +\item{x}{An object of class \link{mixed.mmkin}, \link{saem.mmkin} or \link{nlme.mmkin}} \item{i}{A numeric index to select datasets for which to plot the individual predictions, in case plots get too large} @@ -94,6 +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, 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