diff options
| author | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2013-02-16 21:15:52 +0000 |
|---|---|---|
| committer | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2013-02-16 21:15:52 +0000 |
| commit | 329dadc5fd557cbeec0d1f3a6db57b6b4d6f41d4 (patch) | |
| tree | 8159a9ba673fc6c150a92748e142697c845ece1d /R | |
| parent | e6ba7af29db678ef437e7f8301299ec18f4c66a4 (diff) | |
- 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
Diffstat (limited to 'R')
| -rw-r--r-- | R/mkinfit.R | 2 | ||||
| -rw-r--r-- | R/mkinresplot.R | 2 | ||||
| -rw-r--r-- | R/plot.mkinfit.R | 4 |
3 files changed, 4 insertions, 4 deletions
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 <http://www.gnu.org/licenses/>
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 |