Plot a fitted nonlinear mixed model obtained via an mmkin row object

# S3 method for nlme.mmkin
plot(
  x,
  i = 1:ncol(x$mmkin_orig),
  main = "auto",
  legends = 1,
  resplot = c("time", "errmod"),
  standardized = FALSE,
  show_errmin = TRUE,
  errmin_var = "All data",
  errmin_digits = 3,
  cex = 0.7,
  rel.height.middle = 0.9,
  ymax = "auto",
  ...
)

Arguments

x	An object of class `nlme.mmkin`
i	A numeric index to select datasets for which to plot the nlme fit, in case plots get too large
main	The main title placed on the outer margin of the plot.
legends	An index for the fits for which legends should be shown.
resplot	Should the residuals plotted against time, using `mkinresplot`, or as squared residuals against predicted values, with the error model, using `mkinerrplot`.
standardized	Should the residuals be standardized? This option is passed to `mkinresplot`, it only takes effect if `resplot = "time"`.
show_errmin	Should the chi2 error level be shown on top of the plots to the left?
errmin_var	The variable for which the FOCUS chi2 error value should be shown.
errmin_digits	The number of significant digits for rounding the FOCUS chi2 error percentage.
cex	Passed to the plot functions and `mtext`.
rel.height.middle	The relative height of the middle plot, if more than two rows of plots are shown.
ymax	Maximum y axis value for `plot.mkinfit`.
...	Further arguments passed to `plot.mkinfit` and `mkinresplot`.

Value

The function is called for its side effect.

Author

Johannes Ranke

Examples

ds <- lapply(experimental_data_for_UBA_2019[6:10],
 function(x) subset(x$data[c("name", "time", "value")], name == "parent"))
f <- mmkin("SFO", ds, quiet = TRUE, cores = 1)
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0195
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.011
#plot(f) # too many panels for pkgdown
plot(f[, 3:4])
library(nlme)
f_nlme <- nlme(f)

#plot(f_nlme) # too many panels for pkgdown
plot(f_nlme, 3:4)