diff options
Diffstat (limited to 'man/plot.mixed.mmkin.Rd')
| -rw-r--r-- | man/plot.mixed.mmkin.Rd | 13 |
1 files changed, 12 insertions, 1 deletions
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{ |