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
mkinfitobject. - 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")
See also
mkinplot, for a way to plot the data and the fitted lines of the
mkinfit object.