Function to plot residuals stored in an mkin object

Usage

mkinresplot(object,
    obs_vars = names(object$mkinmod$map),
    xlim = c(0, 1.1 * max(object$data$time)),
    xlab = "Time", ylab = "Residual",
    maxabs = "auto", legend = TRUE, lpos = "topright", ...)

Arguments

object
A fit represented in an mkinfit object.
obs_vars
A character vector of names of the observed variables for which residuals should be plotted. Defaults to all observed variables in the model
xlim
plot range in x direction.
xlab
Label for the x axis. Defaults to "Time [days]".
ylab
Label for the y axis. Defaults to "Residual [% of applied radioactivity]".
maxabs
Maximum absolute value of the residuals. This is used for the scaling of the y axis and defaults to "auto".
legend
Should a legend be plotted? Defaults to "TRUE".
lpos
Where should the legend be placed? Default is "topright". Will be passed on to legend.
...
further arguments passed to plot.

Description

This function plots the residuals for the specified subset of the observed variables from an mkinfit object. A combined plot of the fitted model and the residuals can be obtained using plot.mkinfit using the argument show_residuals = TRUE.

Value

Nothing is returned by this function, as it is called for its side effect, namely to produce a plot.

Examples

model <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))
Successfully compiled differential equation model from auto-generated C code.
fit <- mkinfit(model, FOCUS_2006_D, quiet = TRUE) mkinresplot(fit, "m1")