From e959fde98f95f3595e01490b67892678bbcd1b27 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 7 May 2014 14:47:28 +0200 Subject: Fork the gmkin GUI from mkin. See ChangeLog for details --- R/mkinresplot.R | 53 ----------------------------------------------------- 1 file changed, 53 deletions(-) delete mode 100644 R/mkinresplot.R (limited to 'R/mkinresplot.R') diff --git a/R/mkinresplot.R b/R/mkinresplot.R deleted file mode 100644 index 07bd7df..0000000 --- a/R/mkinresplot.R +++ /dev/null @@ -1,53 +0,0 @@ -# Copyright (C) 2008-2014 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 -if(getRversion() >= '2.15.1') utils::globalVariables(c("variable", "residual")) - -mkinresplot <- function (object, - obs_vars = names(object$mkinmod$map), - xlab = "Time", ylab = "Residual", - maxabs = "auto", legend= TRUE, lpos = "topright", ...) -{ - 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) - - col_obs <- pch_obs <- 1:length(obs_vars) - names(col_obs) <- names(pch_obs) <- obs_vars - - plot(0, xlab = xlab, ylab = ylab, - xlim = c(0, 1.1 * max(object$data$time)), - 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]) - } -} -- cgit v1.2.1