mkin (version 0.9-31)

mkinresplot: 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

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.

See Also

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

Examples

data <- mkin_wide_to_long(schaefer07_complex_case, time = "time")
model <- mkinmod(
  parent = list(type = "SFO", to = c("A1", "B1", "C1"), sink = FALSE),
  A1 = list(type = "SFO", to = "A2"),
  B1 = list(type = "SFO"),
  C1 = list(type = "SFO"),
  A2 = list(type = "SFO"))
fit <- mkinfit(model, data, plot=TRUE)
mkinresplot(fit, "A1")