1 files changed, 84 insertions, 65 deletions
diff --git a/R/mkinresplot.R b/R/mkinresplot.R
index 974d0549..5377dbf2 100644
--- a/R/mkinresplot.R
+++ b/R/mkinresplot.R
@@ -1,65 +1,84 @@
-# Copyright (C) 2008-2014,2019 Johannes Ranke
-# Contact: jranke@uni-bremen.de
-
-# This file is part of the R package mkin
-
-# mkin is free software: you can redistribute it and/or modify it under the
-# terms of the GNU General Public License as published by the Free Software
-# Foundation, either version 3 of the License, or (at your option) any later
-# version.
-
-# This program is distributed in the hope that it will be useful, but WITHOUT
-# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
-# FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
-# details.
-
-# You should have received a copy of the GNU General Public License along with
-# this program. If not, see <http://www.gnu.org/licenses/>
-if(getRversion() >= '2.15.1') utils::globalVariables(c("variable", "residual"))
-
-mkinresplot <- function (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", 
-  col_obs = "auto", pch_obs = "auto",
-  frame = TRUE,
-  ...)
-{
-  obs_vars_all <- as.character(unique(object$data$variable))
-
-  if (length(obs_vars) > 0){
-      obs_vars <- intersect(obs_vars_all, obs_vars)
-  } else obs_vars <- obs_vars_all
-
-  residuals <- subset(object$data, variable %in% obs_vars, residual)
-
-  if (maxabs == "auto") maxabs = max(abs(residuals), na.rm = TRUE)
-
-  # Set colors and symbols
-  if (col_obs[1] == "auto") {
-    col_obs <- 1:length(obs_vars)
-  }
-
-  if (pch_obs[1] == "auto") {
-    pch_obs <- 1:length(obs_vars)
-  }
-  names(col_obs) <- names(pch_obs) <- obs_vars
-
-  plot(0, type = "n", frame = frame,
-       xlab = xlab, ylab = ylab,
-       xlim = xlim,
-       ylim = c(-1.2 * maxabs, 1.2 * maxabs), ...)
-
-  for(obs_var in obs_vars){
-    residuals_plot <- subset(object$data, variable == obs_var, c("time", "residual"))
-    points(residuals_plot, pch = pch_obs[obs_var], col = col_obs[obs_var])
-  }
-
-  abline(h = 0, lty = 2)
-
-  if (legend == TRUE) {
-    legend(lpos, inset = c(0.05, 0.05), legend = obs_vars,
-      col = col_obs[obs_vars], pch = pch_obs[obs_vars])
-  }
-}
+if(getRversion() >= '2.15.1') utils::globalVariables(c("variable", "residual"))
+
+#' Function to plot residuals stored in an mkin object
+#' 
+#' 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 \code{\link{plot.mkinfit}} using the
+#' argument \code{show_residuals = TRUE}.
+#' 
+#' @param object A fit represented in an \code{\link{mkinfit}} object.
+#' @param 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
+#' @param xlim plot range in x direction.
+#' @param xlab Label for the x axis. Defaults to "Time [days]".
+#' @param ylab Label for the y axis. Defaults to "Residual [\% of applied
+#'   radioactivity]".
+#' @param maxabs Maximum absolute value of the residuals. This is used for the
+#'   scaling of the y axis and defaults to "auto".
+#' @param legend Should a legend be plotted? Defaults to "TRUE".
+#' @param lpos Where should the legend be placed? Default is "topright". Will
+#'   be passed on to \code{\link{legend}}.
+#' @param col_obs Colors for the observed variables.
+#' @param pch_obs Symbols to be used for the observed variables.
+#' @param frame Should a frame be drawn around the plots?
+#' @param \dots further arguments passed to \code{\link{plot}}.
+#' @return Nothing is returned by this function, as it is called for its side
+#'   effect, namely to produce a plot.
+#' @author Johannes Ranke
+#' @seealso \code{\link{mkinplot}}, for a way to plot the data and the fitted
+#'   lines of the mkinfit object.
+#' @examples
+#' 
+#' model <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))
+#' fit <- mkinfit(model, FOCUS_2006_D, quiet = TRUE)
+#' mkinresplot(fit, "m1")
+#' 
+#' @export
+mkinresplot <- function (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", 
+  col_obs = "auto", pch_obs = "auto",
+  frame = TRUE,
+  ...)
+{
+  obs_vars_all <- as.character(unique(object$data$variable))
+
+  if (length(obs_vars) > 0){
+      obs_vars <- intersect(obs_vars_all, obs_vars)
+  } else obs_vars <- obs_vars_all
+
+  residuals <- subset(object$data, variable %in% obs_vars, residual)
+
+  if (maxabs == "auto") maxabs = max(abs(residuals), na.rm = TRUE)
+
+  # Set colors and symbols
+  if (col_obs[1] == "auto") {
+    col_obs <- 1:length(obs_vars)
+  }
+
+  if (pch_obs[1] == "auto") {
+    pch_obs <- 1:length(obs_vars)
+  }
+  names(col_obs) <- names(pch_obs) <- obs_vars
+
+  plot(0, type = "n", frame = frame,
+       xlab = xlab, ylab = ylab,
+       xlim = xlim,
+       ylim = c(-1.2 * maxabs, 1.2 * maxabs), ...)
+
+  for(obs_var in obs_vars){
+    residuals_plot <- subset(object$data, variable == obs_var, c("time", "residual"))
+    points(residuals_plot, pch = pch_obs[obs_var], col = col_obs[obs_var])
+  }
+
+  abline(h = 0, lty = 2)
+
+  if (legend == TRUE) {
+    legend(lpos, inset = c(0.05, 0.05), legend = obs_vars,
+      col = col_obs[obs_vars], pch = pch_obs[obs_vars])
+  }
+}