Function to plot residuals stored in an mkin object

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.

Usage

1
2
3
4
5
  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.

Value

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

Author(s)

Johannes Ranke

See Also

mkinplot, for a way to plot the data and the fitted lines of the mkinfit object.

Examples

1
2
3
model <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))
fit <- mkinfit(model, FOCUS_2006_D, quiet = TRUE)
mkinresplot(fit, "m1")

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.