From 329dadc5fd557cbeec0d1f3a6db57b6b4d6f41d4 Mon Sep 17 00:00:00 2001 From: jranke Date: Sat, 16 Feb 2013 21:15:52 +0000 Subject: - Some maintenance fixes - Start of the examples vignette git-svn-id: svn+ssh://svn.r-forge.r-project.org/svnroot/kinfit/pkg/mkin@61 edb9625f-4e0d-4859-8d74-9fd3b1da38cb --- R/mkinfit.R | 2 +- R/mkinresplot.R | 2 +- R/plot.mkinfit.R | 4 ++-- 3 files changed, 4 insertions(+), 4 deletions(-) (limited to 'R') diff --git a/R/mkinfit.R b/R/mkinfit.R index f574b6e..411cb40 100644 --- a/R/mkinfit.R +++ b/R/mkinfit.R @@ -158,7 +158,7 @@ mkinfit <- function(mkinmod, observed, points(subset(observed, name == obs_var, c(time, value)), pch = pch_obs[obs_var], col = col_obs[obs_var]) } - matlines(out_plot$time, out_plot[-1], co = col_obs, lty = lty_obs) + matlines(out_plot$time, out_plot[-1], col = col_obs, lty = lty_obs) legend("topright", inset=c(0.05, 0.05), legend=obs_vars, col=col_obs, pch=pch_obs, lty=1:length(pch_obs)) } diff --git a/R/mkinresplot.R b/R/mkinresplot.R index 6392794..b348278 100644 --- a/R/mkinresplot.R +++ b/R/mkinresplot.R @@ -19,7 +19,7 @@ # this program. If not, see mkinresplot <- function (object, obs_vars = vector(), - xlab = "Time [days]", ylab = "Residual [% of applied radioactivity]", + xlab = "Time", ylab = "Residual", maxabs = "auto", legend= TRUE, lpos = "topright", ...) { obs_vars_all <- as.character(unique(object$data$variable)) diff --git a/R/plot.mkinfit.R b/R/plot.mkinfit.R index 1858aa8..b47678b 100644 --- a/R/plot.mkinfit.R +++ b/R/plot.mkinfit.R @@ -1,8 +1,8 @@ plot.mkinfit <- function(x, fit = x, xlab = "Time", ylab = "Observed", - xlim = range(fit$data$time), ylim = range(fit$data$observed, na.rm = TRUE), + xlim = range(fit$data$time), ylim = c(0, max(fit$data$observed, na.rm = TRUE)), col_obs = 1:length(fit$mkinmod$map), - pch_obs = col_obs, lty_obs = rep(1, length(fit$mkinmod$map)), + pch_obs = col_obs, lty_obs = 1, add = FALSE, legend = !add, ...) { solution_type = fit$solution_type -- cgit v1.2.1